材料加工中的数值模拟方法-微观组织数值模拟ppt课件

材料加工过程的数值模拟微观组织数值模拟(V),任课教师:王锦程,Office：公字楼216TelO)Email：jchwang,TheCahn-HilliardEquation,(forconservedquantities),TheGinzbug-Landau(Allen-Cahn)Equation,(fornon-conservedquantities),描述同构相变过程中，只涉及到成分场的变化，成分场变化采用Cahn-Hilliard方程描述：,M为化学迁移率，一般为系统成分和温度的函数。这类模型可描述同构相变过程,系统化学自由能密度,梯度能系数,弹性能,Ni-Al合金系统析出过程,(1,2,3):(1,1,1)0,(-1,-1,1)0,(-1,1,-1)0,(1,-1,-1)0,Landau多项式自由能形式,三个序参量场,一个序参量场,CALPHAD自由能形式,高斯散度定理,平衡方程（物质或能量守恒）,微分方程,随时间的变化率,通量的散度,内部的源或汇,化学势,变分是求泛函的极值,广义菲克定律,化学势梯度，扩散的驱动力,守恒序参量,非守恒序参量,场量变化速率与驱动力成正比,序参量可看成广义坐标，能量对坐标的变分导数，类似于能量对广义坐标的微分导数，可看成广义力，即驱动力。粘性弛豫：运动速度（广义坐标的变化率）与驱动力成正比（最简单的情形）。,相场模型建立步骤：,分析问题的物理背景，确定描述该体系组织演变的序参量根据所确定的序参量，建立描述体系各状态的统一自由能形式根据序参量特性，确定相应的动力学方程确定模型中的各种参数确定是否考虑系统中可能存在的各向异性确定是否需要考虑系统的随机涨落选择合适的数值算法求解偏微分方程组,相场法应用实例（II）,纯物质凝固相场模型的建立,序参量确定自由能密度函数确定相场动力学方程温度场方程相场参数确定（渐近分析法，平衡解法）各向异性噪声,对于一个封闭空间Rd，材料由界面(t)Rd-1分割为固/液两态。对于纯物质有,枝晶生长：典型的Stefan问题（自由边界问题）,远场条件,尖锐界面模型,需跟踪界面,kineticmobility,传热方程,界面守恒,Gibbs-Thomson关系（曲率效应）,Gibbs-Thompson系数,纯物质相场模型包含两个场变量：相场和温度场(无浓度场)。,自由能密度写为序参量的级数展开形式：,序参量平衡值的确定:,自由能密度函数的构造(1),在凝固点Tm进行泰勒展开：,其中：,当T=Tm时：,在T=TM时，固相与液相的自由能相等,选取：,Thistermtiltsthefreeenergywellsup/down,allowingfordiscontinuouschangeinasTgoesthrough1stordertransitionatTM。,最终得到：,极小值在：,由于液相的自由能与序参量无关，可以将其作为参考态，令其取值为0，可将f(,T)重新写为如下形式：,The“Double-Well”FreeEnergyDensity,是一个双阱曲线，其在处，取最小值，对应着固相和液相。,选择如此形式要满足：当时，其导数，从而保证曲线独立于温度，在处（即固相和液相）取极小值，而使体系处于稳态；同时满足时，时,Needathermodynamicstatefunctionintermsof,and,sothat,solid,liquid,degreeofordering,andapenaltyforisothermallytransformingfromsolidtoliquid,自由能函数的构造(2),AssemblingtheFreeEnergyFunction,“classic”Gibbsfreeenergyoftransformation,canbesettozero,ifliqidcomponentAistakenasthereferencestate,故,TrendsintheFreeEnergy,liquid,solid,系统自由能,liquid,Interfacepenalty,without,solid,with,相场动力学方程,Solids=1,liquidL=0,propagationofinterface,根据标准热传导方程:,考虑固/液界面熵的变化：,温度场方程,This(phenomenological)termaccountsforthereleaseoflatentheatastheinterfaceadvances,Thisisaflux-conservativeequation,i.e.itconservesheatFluxintheabsenceoflatentheat,W,Fourier定律,纯物质凝固相场模型,Karma纯物质凝固相场模型,A.KarmaandW.J.RappelPRE1998,57:4323,固相：,液相：,Thephase-field(orderparameter)andtemperaturefieldapproachsharp-interfaceprofileswhenw0,渐近性分析的目的是将相场模型与尖锐界面模型联系起来，从而证明相场法在解决自由边界问题的适用性。分析的原理是选取界面厚度远小于某一物理尺度，尺度上的巨大差异可以采用扰动展开法来解决。,相场模型参数获取方法：渐近分析,尖锐界面模型,相场模型,？,渐近分析,Interfacekineticcoefficient,Capillarylength,Notethedetailsofa1anda2dependontheformofh().SeeA.KarmaandW.R.Rappel,Phys.Rev.E(1996),相场模型的极限形式,一维稳态解在相场方程中引入的各个系数都不是一个可测量的物理参数。为了解这些设定参数与宏观可测量的物理参量间的关系，考虑一维相变系统，其在区域(-,0)是固