材料科学计算机模拟技术,计算机模拟技术在材料科学中的应用

292广西大学学报(自然科学版)第26卷　的成功是对相变的理解.例如,固体的结晶有序,合金的成份有序或铁磁体的磁化.这种模拟属于所谓“物质的平衡态”,也就是物质从头至尾已弛豫至与环境达到热平衡和化学平衡.但是,实际许多工艺上情况是远离平衡的,例如,在铸造、焊接、拉丝和施压等情况下,平衡统计力学是不合适的.最近十年期间,非平衡过程的理论和这些过程的数学建模技术已经取得很大进步.随巨型计算机的出现,用于规则的结晶固体的模拟计算,已经达到了定量预测的能力.最新的进展表明有可能以相似的精度描述诸如缺陷附近的晶体形变、表面和晶粒边界的非规则图像.这些新方法甚至有可能用以研究物质的亚稳态或严重无序状态.

表1　计算机的模拟层次、空间尺度及模拟对象

模拟层次

电子层次

原子分子层次

微观层次

宏观层次空间尺度011nm～1nm1nm～10nm～1Λm>1模拟对象电子结构结构、力学性能、、烧结、33.1.因此,研究材料的电子结构及相关物性,对从微,也能为探索改善材料力学性能的可能途径提供指导.基于量子力学第一原理的局部密度函数(LDF)理论上的各种算法[5](LMTO,

-SW,LKKR等)已能够计算材料的电子结构及一些基本物理性能,包括晶界—非FLAPW,SCF-XΑ

晶—自由表面与断纹面—杂质—缺陷等各类原子组态的电子结构、相结构稳定性、点和切变面缺陷能量、理想解能量、原子键强及热力学函数等,这使得在实验和理论之间的比较不再局限于依靠经验或半经验参量势函数的计算模式.

3.2　原子模拟技术

按照获得原子位形或微观状态的方法,对于完整和非完整晶体的结构、动力学和热力学性质,有几种可行的模拟方法,如分子动力学方法(MD),蒙特卡罗方法(MC),最小能量法[7](EM)等.分子动力学的目标是研究体系中与时间和温度有关的性质而不只是静力学模拟中研究的构型方面.分子动力学方法是求解运动方程(如牛顿方程、哈密顿方程或拉格朗日方程),通过分析系统中各粒子的受力情况,用经典或量子的方法求解系统中各粒子在某时刻的位置和速度,来确定粒子的运动状态.蒙特卡罗方法是根据待求问题的变化规律,人为地构造出一个合适的概率模型,依照该模型进行大量的统计试验,使它的某些统计参量正好是待求问题的解.最小能量法是利用计算机计算晶体的能量,通过调整原子的位置、调整原子间的化学键长和键角得到最可能的结构,使其系统能量下降,达到最小,所计算的能量值与实验结果相比较,可达到相当精确的程度.

3.3　连续介质模型的模拟方法

为处理宏观问题,常用的方法主要包括传统的有限差分法、有限元法、边界元法等.例如,对材料研究中的传热温度场、传质扩散等问题都可借助这些方法进行求解.此外,对于某些连续的材料微观物理演变过程,也可以在对空间和时间的离散化处理的基础上,采用一定的算法对其进行数据模拟,如对材料的显微组织转变过程、晶粒或第二相粒子长大过程等现象的数值模拟.3.4　综合化模拟方法

综合模拟技术是近年来兴起并蓬勃发展的一类新技术.综合化的含义主要体现在研究方法和研究对象的空间尺度两个方面,前者除发展全新技术外,还包括将原有的基于交互作用势函数的原子模拟技术、从第一原理出发的各种计算技术、连续介质模型、离散化数值计算这三类技术相结合的模拟技术;后者或是直接研究介于原子尺度和宏观尺度之间中间尺度(1～100Λm)的材料结构与性能,或是将不同尺度的材料行为联系起来作为统一体加以研究,特别是如何将不同层次的研究联系起来,已成为材料模拟领域最富挑战性的重点课题.