混合物质原子结构的自洽场计算

基于含温有界Hartree-Fock-Slater相对论自洽场平均原子模型,应用体积相加原理,计算混合物质的原子结构.给出了求解混合物的Thomas-Fermi方程以及自洽场原子结构的方法,并就一些计算结果作了讨论. 关键词： 混合物 原子结构 Thomas-Fermi模型     

随机扰动下三维流体界面不稳定性的并行计算

对三维流体界面不稳定性的数值模拟引进了新的数值计算方法,并在MPI并行计算环境下进行了数值模拟.利用LevelSet方法确定界面位置,零水平集对应界面位置.对应离散LevelSet方程和界面两侧的两套Euler方程,借助于Ghost网格方法来完成离散.对最后网格点上的两套状态量的辨认依赖于该点的LevelSet值的符号.并进行了数值计算.     

三维流体界面不稳定性的Ghost方法

Euler方法对大多数气体流动处理得很好,但在临近材料界面的区域出现非物理振荡。Lagarange方法对多材料面适用性很好,但它在处理大变形和大多数气体流动的涡结构却在自身的困难。理想的稳定性性强的方法是结合两种方法的长处,利用Level Set方法来确定界面位置,利用Ghost方法为确定临近界面的网格点的状态量,采用单侧界面的熵插值技巧来捕捉适宜的界面边界条件,并与TVD／AC格式进行了数值比较。     

二维分子动力学程序(MDP)的并行与优化

基于消息传递并行程序设计平台ＭＰＩ，结合当前微处理器的高性能特征，探讨了二维分子动力学数值模拟程序的优化和并行，且具体应用到高速碰撞物理模型问题。其中，串行优化的性能提高了两倍，并行效率在由８台高性能微机构成的微机机群上，均大于９０％。     

冲击作用下金属表面微喷射的分子动力学模拟

利用二维分子动力学程序,结合类紧束缚杂化多体势,研究冲击载荷下金属表面包含沟槽型缺陷的微喷射动力学过程.在产生微喷射以后,结果表明材料内部传播着两种波系:反射稀疏波和二次加载压缩波,其中,反射稀疏波波面与沟槽形状相似,而二次压缩波波面随沟槽夹角的变化而变化;并形成了两个压强区:负压区和正压区,负压区的存在表明材料中可能产生微损伤.同时,统计结果表明微喷射体的速度随沟槽半角增加而增加的趋势,微喷射体的粒子数随沟槽半角增大而减少的趋势,当沟槽半角大于60°,微喷射效应消失.以上计算结果可以定性说明射流是沟槽型 关键词： 冲击波 沟槽型缺陷 微喷射 分子动力学     

基于"块-单元"数据结构的分子动力学并行计算

开发了一种基于"块-单元"数据结构的可扩展并行算法,以实现大规模、非均匀分子动力学模拟.它采用空间填充曲线将三维区域分解转换为-维负载平衡问题,然后用基于实测的多层均权法求解,以保持处理机间负载均衡.在一个MPP并行机的500个CPU上,模拟包含2.1×108个粒子的三维金属微喷射模型,该算法获得了420倍的加速比.     

基于跃迁系不可分辨方法的辐射不透明度

描述了基于跃迁系不可分辨(UTA)辐射不透明度的计算方法,给出了锗(Ge)和铝(Al)的计算结果,并同实验测量结果作了比较。     

Molecular dynamics simulation of thermodynamic properties of YAG

In this paper we study the thermodynamic properties of Y分子动力学;热力学;扩散;热学YAG, diffusion, elastic constant, molecular dynamicsProject supported by the National Natural Science Foundation of China (Grant No~10744002).2/2/2007 12:00:00 AMIn this paper we study the thermodynamic properties of Y3Al5O12 （YAG） by using molecular dynamic method combined with two- and three-body potentials. The dependences of melting process, elastic constant and diffusion coefficient on temperature of crystal YAG are simulated and compared with the experimental results. Our results show that anion O has the biggest self-diffusivity and cation Y has the smallest self-diffusivity in a crystal YAG. The calculated diffusion activation energies of ions O, Al and Y are 282.55, 439.46, 469.71kJ/mol, respectively. Comparing with experimental creep activation energy of YAG confirms that cation Y can restrict the diffusional creep rate of crystal YAG.