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

上传人(卖家):三亚风情 文档编号:2256907 上传时间:2022-03-27 格式:PPT 页数:27 大小:1.41MB
下载 相关 举报
相场法数值模拟选编课件.ppt_第1页
第1页 / 共27页
相场法数值模拟选编课件.ppt_第2页
第2页 / 共27页
相场法数值模拟选编课件.ppt_第3页
第3页 / 共27页
相场法数值模拟选编课件.ppt_第4页
第4页 / 共27页
相场法数值模拟选编课件.ppt_第5页
第5页 / 共27页
点击查看更多>>
资源描述

1、相相 场场 法法 数数 值值 模模 拟拟phase-field modeling内容相场法数值模拟相场法数值模拟介绍介绍 (Introduction)相场变量相场变量(Phase-field variables)热力学势函数热力学势函数(thermodynamic energy functional)相场方程相场方程(Phase field equations)一、介绍一、介绍 相场模型是一种建立在热力学基础上,考相场模型是一种建立在热力学基础上,考虑有序化势与热力学驱动力的综合作用来虑有序化势与热力学驱动力的综合作用来建立相场方程描述系统演化动力学的模型。建立相场方程描述系统演化动力学的模型

2、。核心思想核心思想引入一个或多个连续变化的序参量,用弥散引入一个或多个连续变化的序参量,用弥散界面模型代替传统的尖锐界面来描述界面界面模型代替传统的尖锐界面来描述界面尖锐界面与弥散界面尖锐界面与弥散界面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 Co

3、upling of Phase Diagrams and Thermochemistry, 32, 268-294, 2008 尖锐界面尖锐界面cDtc2)()(int,int,cccDtc2相扩散到相扩散到相的溶质相的溶质扩散方程扩散方程cDtc2cDtc211int,int,)(rcDrcDcc)()(int,int,cc 相场法是以相场法是以Ginzburg-Landau理论为基础理论为基础, , 用微用微分方程来体现扩散、有序分方程来体现扩散、有序化势和热力学驱动的综合化势和热力学驱动的综合作用作用, , 它是建立在它是建立在Ginzburg-Landau 唯象理论之上的唯象理论之上的

4、一种近代方法。一种近代方法。基本类型基本类型1.1.连续相场法:连续相场法:扩散方程扩散方程 驰豫方程驰豫方程2.2.微观相场法微观相场法:实际是实际是 Cahn-Hilliard 方程的微方程的微 观离散格点形式。观离散格点形式。Kha- chatuyran 引入微观场引入微观场, 用于描述由原子占据晶用于描述由原子占据晶 格位置的几率作为场变格位置的几率作为场变 量来描述微结构变化量来描述微结构变化于志生, 刘平, 龙永强.基于Ginzburg-Landau 理论的相场法研究进展J.材料热处理技术,2008,37(16):9498相场法原理相场法原理1.1.计算量巨大,可模拟的计算量巨大,

5、可模拟的尺度较小(最大可达几十尺度较小(最大可达几十个微米)。个微米)。2.2.相场参数不容易确定。相场参数不容易确定。1.1.通过相场与温度场、溶质通过相场与温度场、溶质场及其它外部场的耦合,能场及其它外部场的耦合,能有效地将微观与宏观尺度结有效地将微观与宏观尺度结合起来。合起来。2.2.由于不需要追踪晶界位置由于不需要追踪晶界位置能方便处理晶界上溶质聚集能方便处理晶界上溶质聚集和第二相析出问题,并能将和第二相析出问题,并能将晶界能和晶界迁移率的各向晶界能和晶界迁移率的各向异性方便地考虑进去,还能异性方便地考虑进去,还能够较大程度避免点阵的各向够较大程度避免点阵的各向异性。异性。相场模型相场

6、模型优缺点优缺点该方法自提出该方法自提出后,迅速成为后,迅速成为微观组织模拟微观组织模拟的热点的热点二、相场变量二、相场变量(phase-field variables)指那些满足局域守恒条件的场变量指那些满足局域守恒条件的场变量如人们最熟悉的浓度序参量如人们最熟悉的浓度序参量c c指那些不满足局域守恒条件的场变量指那些不满足局域守恒条件的场变量如长程序参量如长程序参量保守场保守场(conserved variables)假设假设C组分体系组分体系摩尔浓度ic摩尔分数ixiVimVinrdxVrdc1mtotCiiCiiVVncx1, 111nntotiix miiiVxVcn成分变量成分变量

7、图图2 2 两种不同组成区域两种不同组成区域非保守场非保守场(non-conserved variables)表示无序域表示有序域或011序参量序参量(order parameters)图图3 3 反相位结构反相位结构非保守场非保守场(non-conserved variables)图图4 4 立方结构转化成四方结构有三个等同取向立方结构转化成四方结构有三个等同取向非保守场非保守场(non-conserved variables)在固液界面表示在液相中表示在固相中1001相场量相场量(phase-fields)两相两相多相多相kppk, 01rk1kk,在系统中任一点相,相应的变量三、热力学势

8、函数三、热力学势函数(thermodynamic energy functional)经典热力学经典热力学fyselbulkFFFFFint 体积自由能体积自由能 界面能界面能 弹性应变能弹性应变能 电磁相互作用能电磁相互作用能(bulk free energy) (interfacial energy) (elastic strain energy) 相场法相场法-热力学热力学rdxxfxxfxFVkkkBkBVkBkBkB220)(2)(2),(),(),(? 和和 梯度能量系数梯度能量系数krdxxfxFVkkkBkBkB220int)(2)(2),(),(?均质与非均质体系均质与非均质

9、体系Homogeneous versus heterogeneous systemsNele Moelans.Phase field method to simulate microstructural evolution (June 2004) 图图5 5 自由能与浓度的关系自由能与浓度的关系固相固相转变转变 反相位结构反相位结构 (anti-phase domain structure) 立方转变为四方相立方转变为四方相 (cubic to tetragonal transformation)42max0043204121)(4)()(ffDCBAffdis3312314312321614

10、121),(kikikiCBAf1 1、固态相变、固态相变- -对称性降低对称性降低各向异性各向异性界面能各向异性通过序参量的梯度项引入到自界面能各向异性通过序参量的梯度项引入到自由能表达式中,如:由能表达式中,如:rdrrrxrxxfFVjiplkjiilikijkljBiBijkB31,1,31,02121),(?AnisotropyElastic misfit energy弹性失配能弹性失配能Nele Moelans.Phase field method: from fundamental theories to a phenomenological simulation method

11、(June 2003) 000ijelijijelijijelijtotij图图6 多相结构应变理论多相结构应变理论bcd弹性失配能弹性失配能有了弹性应变,就可以求得体系总的弹性形变能:有了弹性应变,就可以求得体系总的弹性形变能:rdrrrCFrrrCrrCrelklelijvijklelkltotklijklelklijklij)()()(21)()()()()()(0ijBxijBBBkkijBxijdxdaarxtrxtrxrrxrBkB1)(),(),()()()(0,0ij根据胡克定律可得:根据胡克定律可得:3111311130000000000000000003212 2、凝固、凝

12、固- -单相场变量单相场变量rdxTxfFVBB22*0)(2)(2),(?均质自由能密度均质自由能密度等温凝固,假设摩尔体积不变,即组成梯度项不考虑等温凝固,假设摩尔体积不变,即组成梯度项不考虑则则=0=0?)10156()()(30)(),()(),()(1(),(11)()(),(),(23*22*0pgpTxfpTxfpTxfxxggTxfTxfBBBpABABBpB)()(双阱势函数图图7 7 组成自由能曲线组成自由能曲线图图8 8 双阱势函数与插值函数双阱势函数与插值函数图图9 9 两相体系均质自由能三维图两相体系均质自由能三维图非等温凝固非等温凝固非等温凝固用熵函数来表示,避免表

13、达式中出非等温凝固用熵函数来表示,避免表达式中出现温度现温度亥姆霍兹能内能000002200)(2)(2),(feTsferdxxesSVBB?四、相场方程四、相场方程(phase-field equations)kkkkkkkfLtrFLttr0),(),(),(),(),(),(),(10trxxxfMtrxxFMttrxVBBkBBkBBmCahn-Hilliard 方程方程Ginzburg-Landau 方程方程数值解数值解(Numerical solution)图图10 10 相场变量在空间和时间离散化相场变量在空间和时间离散化1、有限差分方法、有限差分方法finite differ

14、ences2、自适应有限元法、自适应有限元法finite elements based on adaptive3、谱方法、谱方法spectral methods计算方法计算方法代入代入总结总结将系统总能量写成所有取向将系统总能量写成所有取向场变量及其它们梯度的函数场变量及其它们梯度的函数相场方程相场方程构造局域自构造局域自由能密度由能密度数值模拟数值模拟结束语结束语相场模拟通过微积分放映扩散、有序化相场模拟通过微积分放映扩散、有序化以及热力学驱动力的综合作用。相场方程的以及热力学驱动力的综合作用。相场方程的解可以描述系统中固液界面的状态、曲率以解可以描述系统中固液界面的状态、曲率以及界面的移动

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

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 办公、行业 > 各类PPT课件(模板)
版权提示 | 免责声明

1,本文(相场法数值模拟选编课件.ppt)为本站会员(三亚风情)主动上传,163文库仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。
2,用户下载本文档,所消耗的文币(积分)将全额增加到上传者的账号。
3, 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(发送邮件至3464097650@qq.com或直接QQ联系客服),我们立即给予删除!


侵权处理QQ:3464097650--上传资料QQ:3464097650

【声明】本站为“文档C2C交易模式”,即用户上传的文档直接卖给(下载)用户,本站只是网络空间服务平台,本站所有原创文档下载所得归上传人所有,如您发现上传作品侵犯了您的版权,请立刻联系我们并提供证据,我们将在3个工作日内予以改正。


163文库-Www.163Wenku.Com |网站地图|