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

《材料加工过程的数值模拟》——微观组织数值模拟(II)任课教师:王锦程Office：公字楼216Tel：[phone](O)Email：[e-mail]微观组织数值模拟简介计算材料学简介相变理论基础热力学基础基础知识不同空间尺度下的组织结构MolecularDynamicsMonteCarloMethodCellularAutomataPhaseFieldModelsFrontTrackingMethodLevelSetmethod……微观组织的数值模拟方法分子动力学（MolecularDynamics）分子动力学是在原子、分子水平上求解多体问题的重要计算机模拟方法，可预测纳米尺度上的材料动力学特性。在分子动力学中，粒子的运动行为是通过经典的Newton运动方程所描述。通过求解所有粒子的运动方程，可以模拟与原子运动路径相关的基本过程。根据各个粒子运动的统计分析，即可推知体系的各种性质。如可能的构型、热力学性质、分子的动态性质、溶液中的行为，各种平衡态性质等。1.给定条件参数(温度、粒子数、时间等)2.体系初始化(初始位置和速度)3.计算作用于所有粒子上的力4.解牛顿运动方程，计算短时间内(TimeStep)粒子的新位置5.计算粒子新的速度和加速度6.重复3-5直至体系达到平衡，获得新的原子位置7.继续计算取得足够的数据，分析数据得到体系的统计性质。分子动力学模拟的一般步骤amoleculardynamicssimulationoflasermelting缺点:1.准确的原子间相互作用势获取困难2.空间、时间尺度均非常小（原子振动频率，ps）3.计算量巨大Monte-Carlo方法蒙特卡罗(MonteCarlo)方法，也称为计算机随机模拟方法，是以概率统计理论为基础的一种方法。二十世纪四十年代中期，由于科学技术的发展和电子计算机的发明，蒙特卡罗方法作为一种独立的方法被提出来，并首先在第二次世界大战美国进研制原子弹的曼哈顿计划中得到了应用。NicholasMetropolisEnricoFermi物理学家JohnvonNeumann计算机科学家StanislawMarcinUlam数学家Buffon投针问题1777年法国科学家蒲丰提出一种计算圆周π的方法：将长为l的一根针任意投到地面上，用针与一组相间距离为a(l＜a）的平行线相交的频率代替概率P，再利用准确的关系式P=(2l)/(pa)求出p值laayq2/0aypq0qsin2lyaladlPppqqp2)2/(2/sin(0）实验者年份投计次数π的实验值沃尔弗(Wolf)185050003.1596斯密思(Smith)185532043.1553福克斯(Fox)189411203.1419拉查里尼(Lazzarini)190134083.1415929基本思想：将求解问题转化为某随机事件A出现的概率，通过某种“实验”方法，得出A事件出现的概率，进而得到问题的解。1.建立构造一个与所求问题相关的概率过程;2.实现从已知概率分布抽样;3.建立各种估计量，从中得到问题的解。优点①能够逼真地描述具有随机性质事物的特点及物理实验过程。②受几何条件限制小。③收敛速度与问题的维数无关。④误差容易确定。⑤程序结构简单，易于实现。缺点①收敛速度慢。②误差具有概率性。③进行模拟的前提是各输入变量相互独立。蒙特卡罗方法在金融工程学，宏观经济学，计算物理学(如粒子输运计算、量子热力学计算、空气动力学计算)等领域应用广泛。在模拟晶粒的长大方面，目前最为广泛使用的是基于Potts模型的蒙特卡罗方法。MCgraingrowthmodel–polyhedralmicrostructuresMCsolidificationmodel–Two-phasemicrostructuresCellularAutomata（元胞自动机）元胞自动机法是一种用来描述复杂系统在离散空间-时间上演化规律的数学算法。元胞在某一时间步的状态转变由一定的演化规则来决定，并且这种转变是随时间推移对体系各元胞同步进行的。元胞的状态受其相邻元胞状态的影响，同时也影响着相邻元胞的状态。局部之间相互作用，相互影响，通过一定的规则变化而整合成一总体行为。PhaseFieldMethod相场法相场模型是一种建立在热力学基础上，考虑有序化势与热力学驱动力的综合作用来建立相场方程描述系统演化动力学的模型。其核心思想是引入一个或多个连续变化的序参量，用弥散界面模型代替传统的尖锐界面来描述界面。相场法的不足是计算量巨大，可模拟的尺度较小（最大可达几十个微米）。将固相和液相中的热、质传输方程与固/液界面上的边界条件耦合起来，建立一个统一的控制方程。原理简单，但需跟踪界面，计算步骤烦琐，难以处理枝晶生长时诸如枝晶分支合并等复杂的界面形貌演变及三维模拟的扩展，因此限制了它的发展。FrontTrackingMethod界面前沿跟踪法1988年，Osher和Sethian最先提出了水平函数的概念基本思想是将闭合轮廓表示为高维曲面等值点的集合，在一系列内力和外力的作用下，通过演化水平函数并跟踪它的零水平得到轮廓的演化过程.LevelSetMethod水平函数调整法(水平集)LevelSetMethod水平函数调整法界面运动方程||0tV(,)(,)0(,)dxtxxtxdxtxSigneddistance利用φ=0描述界面，定义φ0为液相且其值为至界面的距离，φ0为固相且其绝对值为至界面的距离。V为界面移动速度向量在水平函数调整法中，水平函数代替了相场变量，它表示距固/液界面的距离，液相为正，固相为负，而界面为零。根据水平函数求解速度场，确定新界面的位置；继之，根据新的界面位置调整原有水平函数。通过不断的迭代，求解出界面形貌的演化。水平函数调整法虽避免了相场法的渐近分析，但未解决相场法计算量巨大的根本问题。LevelSetMethod水平函数调整法..,.(2002)YTKimNGoldenfield&(1998)KarmaRapel0100200300400050100150200250300350400~3000hoursonDECAlpha:~20MeshelementsizeCPUT12~minuteonaGHzPCOurdiffusedinterfacemodelwithtrackingofinterfacePhasefieldmodelwithouttrackingofinterface:1Meshelementsize:~270nodeno:~160000nodenoTrackinginterfacemakesthedifference!1.S.Chen,B.Merriman,S.Osher,P.Smereka,J.Comp.Phys.135,8–29(1997)2.S.Osher,R.P.Fedkiwy,J.Comp.Phys.169,463–502(2001)3.LijianTan,NicholasZabaras,J.Comp.Phys.211(2006)36–63Computationtime:2dayswith8nodes(16CPUs).基础知识微观组织数值模拟简介计算材料学简介相变理论基础热力学基础L.D.Landau(1908-1968):1962NobelPrize“forhispioneeringtheoriesforcondensedmatter,especiallyliquidhelium二战后最杰出的物理学家，理论多面手(等离子物理，流体力学，核物理，量子场论和天体物理)凝聚态物理理论：二级相变理论，超导理论，超流理论和Fermi液体理论.概念:元激发，序参量和对称破缺朗道相变理论为了对连续相变进行理论分析，朗道提出了序参量的概念，认为连续相变的特征是物质有序程度的改变及与之相伴随的物质对称性质的变化。通常在临界温度以下的相，对称性较低，有序度较高，序参量非零；临界温度以上的相，对称性较高，有序性较低，序参量为零。随着温度的降低，序参量在临界点连续地从零变到非零。强调对称性的重要性，对称性的存在与否是不容模棱两可的，高对称性相中某一对称元素突然消失，就对应于相变的发生，导致低对称相的出现。是一种平均场理论。在临界点附近不一定适用（涨落很大且互相关联）朗道相变理论核心：对称破缺（SymmetryBroken）序参量:低温有序相的一个标志，描述偏离对称的性质和程度。为某个物理量的平均值，可以是标量、矢量、复数或更加复杂的量。随对称性的不同，它在高温时为零，而低温下取有限值，在Tc处转变，对称破缺意味着序参量不为零的有序相的出现。序参量可为标量、矢量、张量或复数。精髓：将对称破缺这一概念引入了相变理论，将序参量不为零的有序相的出现和母相的对称性下降联系在一起。自由能作为序参量的函数将自由能对序参量作级数展开：24024(,)T由于反演对称，不含奇次幂项高于相变温度时，＝0使系统自由能达到极小；低于相变温度时，0使系统自由能达到极小。Landau理论的具体表达TTcTTcT=Tc基础知识微观组织数值模拟简介计算材料学简介相变理论基础热力学基础相平衡0ixGiii21根据热力学原理，体系在恒温恒压下达到平衡的一般条件要求体系总的自由能达到最小值Gmin，即或组元在各相中的化学位相等，即)AAAdxdGxG1)rjjjijxGxG2i一般来说，Gibbs自由能表达式的建立主要有两种方法：一是直接提出描述体系Gibbs自由能的热力学模型；二是由实验数据经数学拟合得到数学表达式，并赋予物理意义。恒压下，Gibbs自由能是温度和成分的函数，任一多组元溶体相的Gibbs自由能，可表示为：refidEmixGGGG构成溶体相的纯组元对Gibbs自由能的贡献来自理想混合熵对Gibbs自由能的贡献过剩Gibbs自由能，表示溶液偏离理想溶液的程度。构成溶体相的纯组元对Gibbs自由能的贡献，是纯组元Gibbs自由能的线性叠加，相当于纯组元之间的简单机械混合，0refiiiGxG纯组元i的标准Gibbs自由能lnididmixmixiiiGTSRTxx因此，Gibbs自由能的求解关键是建立过剩Gibbs自由能的表达式而理想混合熵对Gibbs自由能的贡献可表示为EG0,,11lnCCEmagmjiiiiimjmiiGxGRTyxxGGRef.GibbsenergyIdealmixingenergyExcessfreeenergyMagneticcontributioniiimin1CiiGGijnnPTiinG,,EquilibriumcriteriaCALPHADapproachCALPHAD:Theonlymethodtodealwithmulti-componentandmulti-phasesystem!ExperimentaldataCrystalStructureThermodynamicdataphasediagramdata理想溶体模型规则溶体模型亚规则溶体模型缔合溶液模型亚点阵模型)GxxTRGxGexmiiigmiioi11ln机械混合理想溶液真实体系niiiexxxLxxG02121)(TLiii常数，则为规则溶液相互作用参数DisorderedphaseCuAuNiAlFCCBinaryAlloyCu3AuNi3Al亚点阵模型认为晶格是由几个亚点阵相互穿插构成，粒子在每个亚点阵中随机混合。))A,BC,DacA，B表示同处于一个亚点阵的两种组元；C，D表示