相场法数值模拟选编课件.ppt

相相 场场 法法 数数 值值 模模 拟拟phase-field modeling内容相相场场法法数数值值模模拟拟介绍介绍 (Introduction)相场变量相场变量(Phase-field variables)热力学势函数热力学势函数(thermodynamic energy functional)相场方程相场方程(Phase field equations)一、介绍一、介绍•相场模型是一种建立在热力学基础上，考相场模型是一种建立在热力学基础上，考虑有序化势与热力学驱动力的综合作用来虑有序化势与热力学驱动力的综合作用来建立相场方程描述系统演化动力学的模型建立相场方程描述系统演化动力学的模型核心思想核心思想引入一个或多个连续变化的序参量，用弥散引入一个或多个连续变化的序参量，用弥散界面模型代替传统的尖锐界面来描述界面界面模型代替传统的尖锐界面来描述界面尖锐界面与弥散界面尖锐界面与弥散界面sharp-interface versus diffuse-interface图图1（（a）性能不连续）性能不连续 （（b）性能连续）性能连续N. Moelans, B. Blanpain , P. Wollants, "An introduction to phase-field modeling of microstructure evolution", CALPHAD -- Computer Coupling of Phase Diagrams and Thermochemistry, 32, 268-294, 2008 尖锐界面尖锐界面α相扩散到相扩散到β相的溶质相的溶质扩散方程扩散方程 相场法是以相场法是以Ginzburg-Landau理论为基础理论为基础, , 用微用微分方程来体现扩散、有序分方程来体现扩散、有序化势和热力学驱动的综合化势和热力学驱动的综合作用作用, , 它是建立在它是建立在Ginzburg-Landau 唯象理论之上的唯象理论之上的一种近代方法。

一种近代方法基本类型基本类型1.1.连续相场法：连续相场法：扩散方程扩散方程 驰豫方程驰豫方程2.2.微观相场法微观相场法：：实际是实际是 Cahn-Hilliard 方程的微方程的微 观离散格点形式观离散格点形式Kha- chatuyran 引入微观场引入微观场, 用于描述由原子占据晶用于描述由原子占据晶 格位置的几率作为场变格位置的几率作为场变 量来描述微结构变化量来描述微结构变化于志生, 刘平, 龙永强.基于Ginzburg-Landau 理论的相场法研究进展[J].材料热处理技术,2008,37(16):94~98相场法原理相场法原理1.1.计算量巨大，可模拟的计算量巨大，可模拟的尺度较小（最大可达几十尺度较小（最大可达几十个微米）个微米）2.2.相场参数不容易确定相场参数不容易确定1.1.通过相场与温度场、溶质通过相场与温度场、溶质场及其它外部场的耦合，能场及其它外部场的耦合，能有效地将微观与宏观尺度结有效地将微观与宏观尺度结合起来2.2.由于不需要追踪晶界位置由于不需要追踪晶界位置能方便处理晶界上溶质聚集能方便处理晶界上溶质聚集和第二相析出问题，并能将和第二相析出问题，并能将晶界能和晶界迁移率的各向晶界能和晶界迁移率的各向异性方便地考虑进去，还能异性方便地考虑进去，还能够较大程度避免点阵的各向够较大程度避免点阵的各向异性。

异性相场模型相场模型优缺点优缺点该方法自提出该方法自提出后，迅速成为后，迅速成为微观组织模拟微观组织模拟的热点的热点二、相场变量二、相场变量(phase-field variables)保守场保守场指那些满足局域守恒条件的场变量指那些满足局域守恒条件的场变量如人们最熟悉的浓度序参量如人们最熟悉的浓度序参量c c非保非保守守场场指那些不满足局域守恒条件的场变量指那些不满足局域守恒条件的场变量如长程序参量如长程序参量η保守场保守场(conserved variables)假设假设C组分体系组分体系成分变量成分变量图图2 2 两种不同组成区域两种不同组成区域非保守场非保守场(non-conserved variables)序参量序参量(order parameters)图图3 3 反相位结构反相位结构非保守场非保守场(non-conserved variables)图图4 4 立方结构转化成四方结构有三个等同取向立方结构转化成四方结构有三个等同取向非保守场非保守场(non-conserved variables)相场量相场量(phase-fields)两相两相多相多相三、热力学势函数三、热力学势函数(thermodynamic energy functional)经典热力学经典热力学 体积自由能体积自由能 界面能界面能 弹性应变能弹性应变能 电磁相互作用能电磁相互作用能(bulk free energy) (interfacial energy) (elastic strain energy) 相场法相场法-热力学热力学ϵϵ 和和 ——梯度能量系数梯度能量系数ϵ均质与非均质体系均质与非均质体系Homogeneous versus heterogeneous systemsNele Moelans.Phase field method to simulate microstructural evolution (June 2004) 图图5 5 自由能与浓度的关系自由能与浓度的关系固相固相转变转变 反相位结构反相位结构 (anti-phase domain structure) 立方转变为四方相立方转变为四方相 (cubic to tetragonal transformation)1 1、固态相变、固态相变- -对称性降低对称性降低各向异性各向异性　　　　界面能各向异性通过序参量的梯度项引入到自界面能各向异性通过序参量的梯度项引入到自由能表达式中，如：由能表达式中，如：ϵAnisotropyElastic misfit energy弹性失配能弹性失配能Nele Moelans.Phase field method: from fundamental theories to a phenomenological simulation method (June 2003) 图图6 多相结构应变理论多相结构应变理论bcd弹性失配能弹性失配能有了弹性应变，就可以求得体系总的弹性形变能：有了弹性应变，就可以求得体系总的弹性形变能：根据胡克定律可得：根据胡克定律可得：2 2、凝固、凝固- -单相场变量单相场变量ϵ均质自由能密度均质自由能密度等温凝固，假设摩尔体积不变，即组成梯度项不考虑则等温凝固，假设摩尔体积不变，即组成梯度项不考虑则=0=0ϵ图图7 7 组成自由能曲线组成自由能曲线图图8 8 双阱势函数与插值函数双阱势函数与插值函数图图9 9 两相体系均质自由能三维图两相体系均质自由能三维图非等温凝固非等温凝固非等温凝固用熵函数来表示，避免表达式中出非等温凝固用熵函数来表示，避免表达式中出现温度现温度ϵ四、相场方程四、相场方程(phase-field equations)Cahn-Hilliard 方程方程Ginzburg-Landau 方程方程数值解数值解(Numerical solution)图图10 10 相场变量在空间和时间离散化相场变量在空间和时间离散化1、有限差分方法、有限差分方法finite differences2、自适应有限元法、自适应有限元法finite elements based on adaptive3、谱方法、谱方法spectral methods计计算算方方法法代代入入总结总结将系统总能量写成所有取向将系统总能量写成所有取向场变量及其它们梯度的函数场变量及其它们梯度的函数相场方程相场方程构造局域自构造局域自由能密度由能密度数值模拟数值模拟结束语结束语　　相场模拟通过微积分放映扩散、有序化　　相场模拟通过微积分放映扩散、有序化以及热力学驱动力的综合作用。

相场方程的以及热力学驱动力的综合作用相场方程的解可以描述系统中固液界面的状态、曲率以解可以描述系统中固液界面的状态、曲率以及界面的移动把相场方程与外场（如温度及界面的移动把相场方程与外场（如温度场等）耦合，则可以模拟多晶的真实形貌场等）耦合，则可以模拟多晶的真实形貌　　相场法数值模拟从单相场发展到多相场，　　相场法数值模拟从单相场发展到多相场，从没有流场到包含流场，其数学模型越来越从没有流场到包含流场，其数学模型越来越接近真实模拟过程接近真实模拟过程人有了知识，就会具备各种分析能力，明辨是非的能力所以我们要勤恳读书，广泛阅读，古人说“书中自有黄金屋通过阅读科技书籍，我们能丰富知识，培养逻辑思维能力；通过阅读文学作品，我们能提高文学鉴赏水平，培养文学情趣；通过阅读报刊，我们能增长见识，扩大自己的知识面有许多书籍还能培养我们的道德情操，给我们巨大的精神力量，鼓舞我们前进。