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

材料加工过程的数值模拟微观组织数值模拟(II),任课教师:王锦程,Office：公字楼216TelO)Email：jchwang,微观组织数值模拟简介计算材料学简介相变理论基础热力学基础,基础知识,不同空间尺度下的组织结构,MolecularDynamicsMonteCarloMethodCellularAutomataPhaseFieldModelsFrontTrackingMethodLevelSetmethod,微观组织的数值模拟方法,分子动力学（MolecularDynamics）,分子动力学是在原子、分子水平上求解多体问题的重要计算机模拟方法，可预测纳米尺度上的材料动力学特性。在分子动力学中，粒子的运动行为是通过经典的Newton运动方程所描述。通过求解所有粒子的运动方程，可以模拟与原子运动路径相关的基本过程。,根据各个粒子运动的统计分析，即可推知体系的各种性质。如可能的构型、热力学性质、分子的动态性质、溶液中的行为，各种平衡态性质等。,给定条件参数(温度、粒子数、时间等)体系初始化(初始位置和速度)计算作用于所有粒子上的力解牛顿运动方程，计算短时间内(TimeStep)粒子的新位置计算粒子新的速度和加速度重复3-5直至体系达到平衡，获得新的原子位置继续计算取得足够的数据，分析数据得到体系的统计性质。,分子动力学模拟的一般步骤,amoleculardynamicssimulationoflasermelting,缺点:准确的原子间相互作用势获取困难空间、时间尺度均非常小（原子振动频率，ps）计算量巨大,Monte-Carlo方法,蒙特卡罗(MonteCarlo)方法，也称为计算机随机模拟方法，是以概率统计理论为基础的一种方法。二十世纪四十年代中期，由于科学技术的发展和电子计算机的发明，蒙特卡罗方法作为一种独立的方法被提出来，并首先在第二次世界大战美国进研制原子弹的曼哈顿计划中得到了应用。,NicholasMetropolis,EnricoFermi物理学家,JohnvonNeumann计算机科学家,StanislawMarcinUlam数学家,Buffon投针问题1777年法国科学家蒲丰提出一种计算圆周的方法：将长为l的一根针任意投到地面上，用针与一组相间距离为a(la）的平行线相交的频率代替概率P，再利用准确的关系式P=(2l)/(pa)求出p值,基本思想：将求解问题转化为某随机事件A出现的概率，通过某种“实验”方法，得出A事件出现的概率，进而得到问题的解。建立构造一个与所求问题相关的概率过程;实现从已知概率分布抽样;建立各种估计量，从中得到问题的解。,优点能够逼真地描述具有随机性质事物的特点及物理实验过程。受几何条件限制小。收敛速度与问题的维数无关。误差容易确定。程序结构简单，易于实现。,缺点收敛速度慢。误差具有概率性。进行模拟的前提是各输入变量相互独立。,蒙特卡罗方法在金融工程学，宏观经济学，计算物理学(如粒子输运计算、量子热力学计算、空气动力学计算)等领域应用广泛。在模拟晶粒的长大方面，目前最为广泛使用的是基于Potts模型的蒙特卡罗方法。,MCgraingrowthmodelpolyhedralmicrostructures,MCsolidificationmodelTwo-phasemicrostructures,CellularAutomata（元胞自动机）,元胞自动机法是一种用来描述复杂系统在离散空间-时间上演化规律的数学算法。元胞在某一时间步的状态转变由一定的演化规则来决定，并且这种转变是随时间推移对体系各元胞同步进行的。元胞的状态受其相邻元胞状态的影响，同时也影响着相邻元胞的状态。局部之间相互作用，相互影响，通过一定的规则变化而整合成一总体行为。,PhaseFieldMethod相场法,相场模型是一种建立在热力学基础上，考虑有序化势与热力学驱动力的综合作用来建立相场方程描述系统演化动力学的模型。其核心思想是引入一个或多个连续变化的序参量，用弥散界面模型代替传统的尖锐界面来描述界面。相场法的不足是计算量巨大，可模拟的尺度较小（最大可达几十个微米）。,将固相和液相中的热、质传输方程与固/液界面上的边界条件耦合起来，建立一个统一的控制方程。原理简单，但需跟踪界面，计算步骤烦琐，难以处理枝晶生长时诸如枝晶分支合并等复杂的界面形貌演变及三维模拟的扩展，因此限制了它的发展。,FrontTrackingMethod界面前沿跟踪法,1988年，Osher和Sethian最先提出了水平函数的概念基本思想是将闭合轮廓表示为高维曲面等值点的集合，在一系列内力和外力的作用下，通过演化水平函数并跟踪它的零水平得到轮廓的演化过程.,LevelSetMethod水平函数调整法(水平集),LevelSetMethod水平函数调整法,界面运动方程,Signeddistance,利用=0描述界面，定义0为液相且其值为至界面的距离，Tc,TTc,T=Tc,基础知识微观组织数值模拟简介计算材料学简介相变理论基础热力学基础,相平衡,根据热力学原理，体系在恒温恒压下达到平衡的一般条件要求体系总的自由能达到最小值Gmin，即或组元在各相中的化学位相等，即,一般来说，Gibbs自由能表达式的建立主要有两种方法：一是直接提出描述体系Gibbs自由能的热力学模型；二是由实验数据经数学拟合得到数学表达式，并赋予物理意义。恒压下，Gibbs自由能是温度和成分的函数，任一多组元溶体相的Gibbs自由能，可表示为：,构成溶体相的纯组元对Gibbs自由能的贡献,来自理想混合熵对Gibbs自由能的贡献,过剩Gibbs自由能，表示溶液偏离理想溶液的程度。,构成溶体相的纯组元对Gibbs自由能的贡献，是纯组元Gibbs自由能的线性叠加，相当于纯组元之间的简单机械混合，,纯组元i的标准Gibbs自由能,因此，Gibbs自由能的求解关键是建立过剩Gibbs自由能的表达式,而理想混合熵对Gibbs自由能的贡献可表示为,Ref.Gibbsenergy,Idealmixingenergy,Excessfreeenergy,Magneticcontribution,Equilibriumcriteria,CALPHADapproach,CALPHAD:Theonlymethodtodealwithmulti-componentandmulti-phasesystem!,Experimentaldata,CrystalStructure,Thermodynamicdataphasediagramdata,理想溶体模型规则溶体模型亚规则溶体模型缔合溶液模型亚点阵模型,常数，则为规则溶液,相互作用参数,Disorderedphase,CuAuNiAl,FCCBinaryAlloy,Cu3AuNi3Al,亚点阵模型认为晶格是由几个亚点阵相互穿插构成，粒子在每个亚点阵中随机混合。,A，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,B)0.5(A,B)0.5,(A,B)0.5(A,B)0.25(A,B)0.25,各亚点阵之间的相互作用可以忽略不计，过剩Gibbs自由能是描述同一亚点阵内组元的相互作用，此相互作用与其它亚点阵内组元的种类无关。亚点阵模型中每一亚点阵的结点数目可以相同，也可以完全不同。亚点阵的性质完全不同，每一组元只能进入一种亚点阵，而不能进入另一种亚点阵。计算时，在每一个亚点阵内实行的是正规溶液近似。,双亚点阵模型,点阵分数,组元i在亚点阵s中所占据的结点数,亚点阵s中所有结点数,(A,B)a(A,B)c,由于亚点阵模型假设各亚点阵之间的相互作用可以忽略不计，一个亚点阵内组元的相互作用与其它亚点阵内组元的种类无关。将每一个亚点阵中不同组元混合的理想混合熵相加，就可得到理想混合熵对Gibbs自由能的贡献。,Gibbs自由能的参考面是由每一个亚点阵中只有一种组元存在时的状态所定义的。是AaCc、BaCc、AaDc、BaDc等化合物Gibbs自由能的机械混合,亚点阵模型的过剩Gibbs自由能是描述同一亚点阵内组元的相互作用对理想溶体的偏差。如果一个亚点阵上只有一种