材料加工中的数值模拟方法-微观组织数值模拟(2)

1、材料加工过程的数值模拟微观组织数值模拟(II)任课教师: 王锦程Office：公字楼216Tel O)Email ： 微观组织数值模拟简介微观组织数值模拟简介 计算材料学简介计算材料学简介 相变理论基础相变理论基础 热力学基础热力学基础l 基础知识不不同同空空间间尺尺度度下下的的组组织织结结构构p Molecular Dynamicsp Monte Carlo Method p Cellular Automatap Phase Field Modelsp Front Tracking Methodp Level Set methodp 微观组织的数值模拟方法分子动力

2、学（Molecular Dynamics）v分子动力学是在原子、分子水平上求解多体问题的重要计算机模拟方法，可预测纳米尺度上的材料动力学特性。v在分子动力学中，粒子的运动行为是通过经典的Newton运动方程所描述。通过求解所有粒子的运动方程，可以模拟与原子运动路径相关的基本过程。v根据各个粒子运动的统计分析，即可推知体系的各种性质。如可能的构型、热力学性质、分子的动态性质、溶液中的行为，各种平衡态性质等。给定条件参数(温度、粒子数、时间等)体系初始化(初始位置和速度) 计算作用于所有粒子上的力解牛顿运动方程，计算短时间内(Time Step)粒子的新位置计算粒子新的速度和加速度重复3-5直至体

3、系达到平衡，获得新的原子位置继续计算取得足够的数据，分析数据得到体系的统计性质。分子动力学模拟的一般步骤a molecular dynamics simulation of laser melting 缺点: 准确的原子间相互作用势获取困难空间、时间尺度均非常小（原子振动频率，ps）1.计算量巨大Monte-Carlo方法p 蒙特卡罗(Monte Carlo)方法，也称为计算机随机模拟方法，是以概率统计理论为基础的一种方法。二十世纪四十年代中期，由于科学技术的发展和电子计算机的发明，蒙特卡罗方法作为一种独立的方法被提出来，并首先在第二次世界大战美国进研制原子弹的曼哈顿计划中得到了应用。 Nic

4、holas MetropolisEnrico Fermi物理学家John von Neumann计算机科学家Stanislaw Marcin Ulam数学家Buffon投针问题 1777年法国科学家蒲丰提出一种计算圆周的方法：将长为l的一根针任意投到地面上，用针与一组相间距离为a( la）的平行线相交的频率代替概率P，再利用准确的关系式P=(2l)/(pa)求出p值laayq2/0ay pq0qsin2ly aladlPppqqp2)2/(2/sin(0）实验者年份投计次数的实验值沃尔弗(Wolf)185050003.1596斯密思(Smith)185532043.1553福克斯(Fox)18

5、9411203.1419拉查里尼(Lazzarini) 190134083.1415929基本思想：将求解问题转化为某随机事件A出现的概率，通过某种“实验”方法，得出A事件出现的概率，进而得到问题的解。建立构造一个与所求问题相关的概率过程;实现从已知概率分布抽样;1.建立各种估计量，从中得到问题的解。优点能够逼真地描述具有随机性质事物的特点及物理实验过程。受几何条件限制小。收敛速度与问题的维数无关。误差容易确定。程序结构简单，易于实现。缺点收敛速度慢。误差具有概率性。进行模拟的前提是各输入变量相互独立。p 蒙特卡罗方法在金融工程学，宏观经济学，计算物理学(如粒子输运计算、量子热力学计算、空气动

6、力学计算)等领域应用广泛。p 在模拟晶粒的长大方面，目前最为广泛使用的是基于Potts模型的蒙特卡罗方法。MC grain growth model polyhedral microstructuresMC solidification model Two-phase microstructuresCellular Automata（元胞自动机） 元胞自动机法是一种用来描述复杂系统在离散空间-时间上演化规律的数学算法。 元胞在某一时间步的状态转变由一定的演化规则来决定，并且这种转变是随时间推移对体系各元胞同步进行的。 元胞的状态受其相邻元胞状态的影响，同时也影响着相邻元胞的状态。局部之间相互作

7、用，相互影响，通过一定的规则变化而整合成一总体行为。Phase Field Method 相场法p 相场模型是一种建立在热力学基础上，考虑有序化势与热力学驱动力的综合作用来建立相场方程描述系统演化动力学的模型。p 其核心思想是引入一个或多个连续变化的序参量，用弥散界面模型代替传统的尖锐界面来描述界面。p 相场法的不足是计算量巨大，可模拟的尺度较小（最大可达几十个微米）。 p 将固相和液相中的热、质传输方程与固/液界面上的边界条件耦合起来，建立一个统一的控制方程。p 原理简单，但需跟踪界面，计算步骤烦琐，难以处理枝晶生长时诸如枝晶分支合并等复杂的界面形貌演变及三维模拟的扩展，因此限制了它的发展。

8、 Front Tracking Method 界面前沿跟踪法p 1988年，Osher和Sethian最先提出了水平函数的概念p 基本思想是将闭合轮廓表示为高维曲面等值点的集合，在一系列内力和外力的作用下，通过演化水平函数并跟踪它的零水平得到轮廓的演化过程.Level Set Method 水平函数调整法(水平集)Level Set Method 水平函数调整法界面运动方程| 0tV( , ) ( , ) 0 ( , ) d x txx txd x txSigned distance利用 = 0 描述界面，定义 0 为液相且其值为至界面的距离， TcT TcT =Tc 基础知识 微观组织数值模

9、拟简介 计算材料学简介 相变理论基础 热力学基础相平衡相平衡0ixGiii 21根据热力学原理，体系在恒温恒压下达到平衡的一般条件要求体系总的自由能达到最小值Gmin，即 或组元在各相中的化学位相等，即)AAAdxdGxG1)rjjjijxGxG2ip 一般来说，Gibbs自由能表达式的建立主要有两种方法：一是直接提出描述体系Gibbs自由能的热力学模型；二是由实验数据经数学拟合得到数学表达式，并赋予物理意义。p 恒压下，Gibbs自由能是温度和成分的函数，任一多组元溶体相的Gibbs自由能，可表示为：refidEmixGGGG 构成溶体相的纯组元对Gibbs自由能的贡献来自理想混合熵对Gib

10、bs自由能的贡献过剩Gibbs自由能，表示溶液偏离理想溶液的程度。p 构成溶体相的纯组元对Gibbs自由能的贡献，是纯组元Gibbs自由能的线性叠加，相当于纯组元之间的简单机械混合，0refiiiGxG纯组元i 的标准Gibbs自由能lnididmixmixiiiGT SRTxx p 因此，Gibbs自由能的求解关键是建立过剩Gibbs自由能的表达式p 而理想混合熵对Gibbs自由能的贡献可表示为EG0,11lnCCEmagm jiiiiim jmiiGx GRT yxxGGRef. Gibbs energyIdeal mixing energyExcess free energyMagnet

11、ic contributioniii min1CiiGGijnnPTiinG,Equilibrium criteria CALPHAD approachCALPHAD: The only method to deal with multi-component and multi-phase system!Experimental dataCrystal Structure Thermodynamic dataphase diagram data 理想溶体模型理想溶体模型 规则溶体模型规则溶体模型 亚规则溶体模型亚规则溶体模型 缔合溶液模型缔合溶液模型 亚点阵模型亚点阵模型 )GxxTRGxGe

12、xmiiigmiioi11ln机械混合理想溶液真实体系niiiexxxLxxG02121)(TLiii常数，则为规则溶液相互作用参数Disordered phaseCuAuNiAlFCC Binary AlloyCu3Au Ni3Alp 亚点阵模型认为晶格是由几个亚点阵相互穿插构成，粒子在每个亚点阵中随机混合。) )A,BC,DacA，B表示同处于一个亚点阵的两种组元；C，D表示同处于另外一个亚点阵的两种组元；a和c表示两个亚点阵的结点数的比例。FCC_A1(disorder) L10(order) L12(order) F(order) (A,B)1(A,B)0.25(A,B)0.75(A,

13、B)0.5(A,B)0.5(A,B)0.5(A,B)0.25(A,B)0.251234iiiiyyyy1234iiiiyyyy1234iiiiyyyy1234iiiiyyyyp 各亚点阵之间的相互作用可以忽略不计，过剩Gibbs自由能是描述同一亚点阵内组元的相互作用，此相互作用与其它亚点阵内组元的种类无关。p 亚点阵模型中每一亚点阵的结点数目可以相同，也可以完全不同。p 亚点阵的性质完全不同，每一组元只能进入一种亚点阵，而不能进入另一种亚点阵。p 计算时，在每一个亚点阵内实行的是正规溶液近似。双亚点阵模型) AABBABABABCCDDCDCDCDnxnxyynna acnna acnxnxy

14、ynnc acnnc acsssiiissiinnyNn点阵分数组元i在亚点阵s中所占据的结点数亚点阵s中所有结点数) )A,BC,Dac(A,B )a ( A,B )ccacycaayxcacycaayxBBBAAA p 由于亚点阵模型假设各亚点阵之间的相互作用可以忽略不计，一个亚点阵内组元的相互作用与其它亚点阵内组元的种类无关。将每一个亚点阵中不同组元混合的理想混合熵相加，就可得到理想混合熵对Gibbs自由能的贡献。)lnlnlnlnidmixAABBCCDDGRT a yyyyc yyyyp Gibbs自由能的参考面是由每一个亚点阵中只有一种组元存在时的状态所定义的。是AaCc、 BaC

15、c、 AaDc、 BaDc等化合物Gibbs自由能的机械混合0000acacacacrefACA CBCB CADA DBDB DGy y Gy y Gy y Gy y G ) )A,BC,Dacp 亚点阵模型的过剩Gibbs自由能是描述同一亚点阵内组元的相互作用对理想溶体的偏差。如果一个亚点阵上只有一种组元，则这个亚点阵上的过剩Gibbs自由能为零。当亚点阵中有两种以上的组元时，可以按正规溶液模型计算其过剩Gibbs自由能。), :, : ,: ,EABCA B CDA B DCDAA C DBB C DGy yy Ly Ly yy Ly L )0000, :, : ,: , lnlnlnln acacacacrefidEACA CBCB CADA DBDB DAABBCCDDABCA B CDA B DCDAAC DBB C DGGGGy y Gy y Gy y Gy y GRT a