液H2+D2混合冲击压缩特性的理论计算

在液体变分微扰理论加量子力学修正基础上,用指数6有效作用势对液H2及其同位素液D2的混合物态方程进行了计算,首次给出了液H2+D2混合冲击压缩Hugoniot曲线和冲击温度.     

固氢固氦等温压缩线的理论计算

用原子团簇理论和完全量子力学计算,研究了固氦和固氢的等温压缩线.在固氢的分子间相互作用势计算中,采用了球形势近似和等效氦原子模型.计算结果与100GPa以下的实验等温线相比较,发现两者有很好的符合. 关键词：     

高密度氦状态方程的分子动力学模拟研究

 用分子动力学模拟和与固氦实验等温线比较的方法,确定了高密度氦指数6势函数中劲度系数α的优选值为12.7。然后我们用这个新的α值计算了fcc固氦的状态方程及其径向分布函数,发现当把ρ限定为1.60 g/cm³时,其径向分布函数的第二个峰将在2 000 K到3 040 K区间消失,表明此时固氦的长程有序性降低,或者说发生了固-液相变过程。     

中子辐照金属材料的脆化模型研究

金属材料的辐照脆化问题一直以来都是核能安全领域亟待解决的关键问题之一.为了更准确地预测金属材料的辐照脆化行为,基于Johnson-Cook本构模型,将未辐照金属材料的断裂真应力取作辐照材料的断裂真应力,建立了通过辐照退火态金属材料屈服强度就能够预测其整个真应力-应变曲线,以及断裂真应变的辐照脆化模型.实验研究了不同中子剂量辐照退火态高纯铝的准静态拉伸真应力-应变曲线、断裂真应力和断裂真应变随辐照剂量的变化规律.结果表明,辐照剂量越高,高纯铝的屈服强度越高,断裂真应变越低,但断裂真应力几乎不变.通过TEM显微分析获得了高纯铝内部辐照缺陷的尺寸和数密度随辐照剂量的变化规律,结果表明,辐照剂量越高,孔洞的尺寸和数密度越高,但位错环尺寸和数密度始终很小,难以准确统计.由辐照高纯铝实验数据拟合得到了辐照脆化模型所需参数,并检验了该模型的预测效果.结果表明,无论是通过实验还是显微分析得到辐照高纯铝的屈服强度,模型的预测结果均能够与实验结果较好地吻合,且模型对退火态高纯铝临界中子剂量的预测值也与文献结果一致.     

面心立方金属铈的高压电子结构研究

本文简要介绍了我们开发的基于密度泛函理论和动力学平均场理论（DFT+DMFT）的全自洽的的计算程序,并运用该程序研究了面心立方金属铈的等结构相变. 动力学平均场采用的是目前公认的最为准确和高效的连续时间蒙特卡罗方法作为杂质求解器. 我们首先分别计算了 和 相的电子态密度,该结果可以正确的描述两个相的电子结构特征. 然后计算了不同温度下的相变前后的电子结构,从中可以得到电子结构在不同温度下的相变过程,与之前的实验结果一致. 这一工作为以后进一步的金属铈的等结构相变研究打下了基础.     

预应变对中子辐照高纯铝拉伸性能的影响

金属材料的中子辐照硬化和脆化一直都是核能安全领域十分关注的重要问题之一. 为了进一步认识预应变对中子辐照金属材料塑性形变和最终断裂特性的影响规律, 及其微观机理, 本文研究了10%拉伸预应变高纯铝的拉伸应力-应变曲线、失稳应力和失稳应变等随辐照剂量的变化规律. 结果表明, 辐照剂量越高, 预应变高纯铝内部孔洞的尺寸和数密度越高, 导致屈服强度和极限拉伸强度越高, 均匀延伸率和失稳应变越小, 表现出典型的辐照硬化和脆化效应, 但失稳应力与辐照剂量几乎无关. 相同辐照剂量条件下, 预应变引入的高密度位错能够显著降低辐照孔洞的尺寸和数密度, 加之辐照退火效应的综合影响, 导致预应变能够降低高纯铝屈服强度的增长率和失稳应变的下降率, 从而表现出一定的抑制辐照硬化和脆化的能力, 预应变还能够提高高纯铝的失稳应力, 但整体而言预应变并不能提高高纯铝的延性. 最后, 基于J-C本构模型的中子辐照退火态金属材料的脆化模型能够直接应用于预应变金属材料, 且模型预测结果与实验结果吻合较好.      

Extension of Compressibility-Route Cubic Equations of State and the Radial Distribution Functions at Contact to Multi-Component Hard-Sphere Mixtures

The equation of state(EOS) for hard-sphere fluid derived from compressibility routes of Percus-Yevick theory(PYC) is extended. The two parameters are determined by fitting well-known virial coefficients of pure fluid.The extended cubic EOS can be directly extended to multi-component mixtures, merely demanding the EOS of mixtures also is cubic and combining two physical conditions for the radial distribution functions at contact(RDFC) of mixtures.The calculated virial coefficients of pure fluid and predicted compressibility factors and RDFC for both pure fluid and mixtures are excellent as compared with the simulation data. The values of RDFC for mixtures with extremely large size ratio 10 are far better than the BGHLL expressions in literature.     

稠密氦气辐射不透明度的实验研究

用二级氢气炮作为冲击压缩加载工具和多通道瞬态辐射高温计作为主要测量系统,测量了冲击压缩氦气等离子体的光辐亮度历史(初始温度293 K,初始压力为0.6 MPa和1.2 MPa两种).根据实测记录信号波形的有关特征量,计算得到了氦等离子体的光谱吸收系数κ(λ)和冲击压缩铝基板表面的光反射率R.结果发现:受冲击铝基板表面的光反射率～0.4(比其初始反射率0.8约降低一半),与Erskine的数据相同;对于光谱吸收系数的实测数据,本文目前未能给出合理的物理解释.     

Simulated Melting Curve of NaC1 up to 200 kbar

The melting curve of NaCI is studied up to 200 kbar by means of the shell-model molecular dynamics method, usiug massive shell core interaction potentials. The model for the interatomic interaction is shown to produce reasonable results at a wide range of pressures in bulk transitions. The pressure dependence of the melting curve of NaCl was calculated and the result was modified on the assumption of overheating due to the small system size and small time scale simulation. The final result is in good agreement with the corrected experimental values, accounting for melting mechanisms such as surface heating or superheating. Therefore, it is believed that bulk transition simulation at constant pressure indeed provides a useful tool for studying the melting transition.     

First-Principle-Based Calculations of the Hugoniot of Cu

The equation of state of face-centred-cubic (fcc) copper crystals at pressures up to 500 GPa and relative volume to 0.55 have been evaluated by using the full-potential linear muffin-tin orbital (FPLMTO) total-energy method combining with a mean-field model of the vibrational partition function. The mean-field is constructed from the sum of all the pair potentials between the reference atom and the others of the system. The calculated properties are in good agreement with the available shock-wave experimental measurements.