飞片温升对Hugoniot参数测量的影响

 在格临爱森物态方程框架内建立了一种估算高温材料Hugoniot参数的理论方法,这种方法计算简单,原则上适用于各种固体材料,弥补了BACT方法的缺点。结合一些Hugoniot实验数据研究了实验中飞片的温升对粒子速度的影响。数值计算表明：飞片温升对粒子速度造成的不确定度非常小,因此在数据处理中,假定飞片的Hugoniot参数不变是合理的。     

Shock compressibility of polycrystalline nickel aluminide

ABSTRACT

We present experimental data on Hugoniot and release isentropes of nickel aluminide obtained for shock pressures up to 83?GPa and make comparison with previously published calculation data. We estimate a Hugoniot elastic limit of samples, construct the EOS of nickel aluminide and investigate the Hugoniot position relative to the melting curve.     

200 GPa压力范围内铝和铜的等熵压缩线计算

 将冲击Hugoniot线作为Grüneisen物态方程的参考线,以冲击的初始状态为参考状态,推导得到线性和二次曲线表示的冲击绝热线所对应的等熵压缩线方程,计算了200 GPa压力范围内铝和铜两种材料的等熵压缩线,并且计算了以Hugoniot关系为基础的Appy经验物态方程导出的等熵压缩线。计算结果表明,以Appy经验物态方程导出的等熵压缩线与以线性冲击绝热线导出的等熵压缩线接近,在200 GPa压力范围内两者相差不到1.5%。将计算得到的铝的等熵压缩线与美国Sandia实验室ICE实验Z864数据进行了比较,由线性Hugoniot得到的等熵压缩线与实验数据相差不到1%,由Appy经验物态方程得到的等熵线与实验数据几乎重合,说明在200 GPa压力范围内,以Appy物态方程和以线性Hugoniot为参考来计算的等熵压缩线有较高的精度。     

液氮与液态一氧化碳混合物的冲击压缩特性研究

 利用二级轻气炮研究了等体积的液态一氧化碳（LCO）与液氮（LN₂）混合物的冲击压缩特性，在10～25 GPa压力范围内获得了5个Hugoniot数据点。混合物冲击绝热线位于单质一氧化碳和液氮Hugoniot拟合线之间，并靠近前者。这意味着混合物的冲击压缩特性不是其中各单质行为的简单平均。     

凝聚态物质状态方程的一个数值模型

建立了凝聚态物质的一个三项式状态方程：以Faussurier平均原子模型为基础计算电子热压和电子热能;以Cowan模型为基础计算离子热压和离子热能;用基于实验数据的半经验拟合公式计算物质的冷压和冷能。用实验数据检验了用平均原子模型计算的平均电离度。将状态方程与Hugoniot关系式相结合,计算了Be和Al的冲击绝热曲线,结果充分地展现出电子在高温、高密度条件下的壳层结构效应。     

泡沫态硅橡胶冲击绝热线的近似计算

在孔穴塌缩模型的基础上,建立了计算泡沫态硅橡胶在不同初始密度下的冲击绝热线参数的理论方法。获得了泡沫态硅橡胶在两种不同初始密度下(0.52 g/cm3和0.92 g/cm3)的冲击绝热线,并与已有的实验数据进行对比分析。结果表明,建立的理论方法对于泡沫态硅橡胶的冲击绝热线计算是适用的。     

Fe/FeO/FeS混合物的Hugoniot线研究

 在50～210 GPa的压力范围内,用二级轻气炮和电探针技术对平均密度ρ₀＝(6.69±0.06) g/cm³的三组元Fe/FeO/FeS（质量分数分别为58.96%、35.83%、5.21%）混合物的Hugoniot线进行了实验测量,所得的Hugoniot参数为：C₀＝(3.97±0.07) km/s,λ＝1.58±0.03。该混合物Hugoniot线的实测结果与用体积可加性原理计算得到的相同组分混合物的Hugoniot线的符合性很好;根据实验数据还计算了混合物的0 K等温压缩线,发现它与体积可加性原理对单质Fe、FeO和FeS的0 K等温压缩线的计算结果相一致,证明了实验结果的合理性与体积可加性原理的适用性,也表明了此混合物在冲击压缩过程中没有发生过可察觉的化学反应。研究结果亦为今后对外地核各种浓度比Fe-O-S体系候选组分高温高压物态方程及物性的进一步研究奠定了基础。     

C30混凝土冲击绝热关系和Grüneisen型状态方程的实验研究

采用一级气炮加载技术和锰铜压力计测试技术,对含初始空隙的C30混凝土在一维应变条件下的冲击特性进行了实验测量和分析.基于锰铜压力计测量的压力波形,确定了C30混凝土材料的冲击绝热关系,即冲击波速度D与波后粒子速度u之间满足线性关系.再从C30混凝土的冲击绝热数据出发,获得了计及初始空隙度₀影响的多项式形式Grüneisen型状态方程中的各项系数.实测压力波形还显示：不同位置处的压力波形在迅速上升至峰值后均随时间逐渐衰减,而冲击波峰值又随传播 关键词： Grüneisen型状态方程 冲击绝热关系 混凝土     

在200 GPa冲击压强下铁是否会发生固-固相变？

 根据Grüneisen物态方程、Hugoniot内能守恒方程和最新发表的铁的热物理参数，计算了ε-(hcp)铁和液态铁的理论Hugoniot曲线。计算的Hugoniot曲线与最新修正的铁的实验数据[J Appl Phys, 2000, 88: 5496]在总体上符合很好，并且可以细分为两个部分：在约低于200 GPa的压强区，用ε-铁模型的计算结果与实验结果符合很好；在约高于260 GPa的压强区，用液态铁模型的计算结果与实验数据也符合得比较好；对介于200~260 GPa之间的压强区，则归属于由ε-铁向液态铁转变的混合相区。这意味着铁的Hugoniot曲线在约200 GPa处出现的微小偏折是由固-液相变引起的，从而否定了Brown [Geophy Res Lett, 2001, 28: 4339]提出的它是一次由ε-铁向另一个未知结构的某个固相铁的相转变的见解。     

Convergence acceleration of molecular dynamics methods for shocked materials using velocity scaling

ABSTRACT

In this work, a convergence acceleration method applicable to extended system molecular dynamics techniques for shock simulations of materials is presented. The method uses velocity scaling to reduce the instantaneous value of the Rankine–Hugoniot conservation of energy constraint used in extended system molecular dynamics methods to more rapidly drive the system towards a converged Hugoniot state. When used in conjunction with the constant stress Hugoniostat method, the velocity scaled trajectories show faster convergence to the final Hugoniot state with little difference observed in the converged Hugoniot energy, pressure, volume and temperature. A derivation of the scale factor is presented and the performance of the technique is demonstrated using the boron carbide armour ceramic as a test material. It is shown that simulation of boron carbide Hugoniot states, from 5 to 20 GPa, using both a classical Tersoff potential and an ab initio density functional, are more rapidly convergent when the velocity scaling algorithm is applied. The accelerated convergence afforded by the current algorithm enables more rapid determination of Hugoniot states thus reducing the computational demand of such studies when using expensive ab initio or classical potentials.     

Laser-shock-driven laboratory measurements of the equation of state of hydrogen isotopes in the megabar regime

Abstract

A high intensity laser has been used to shock-compress liquid deuterium to pressures from 0.22 Mbar, near the upper end of gas-gun shock data, to 3.4Mbar. Using a variety of diagnostics, the density, pressure, and reflectivity of the shocked deuterium were determined, culminating in a knowledge of the behavior of hydrogen isotopes along the principal Hugoniot into the Mbar regime. It was found that hydrogen isotopes become highly compressible near 1 Mbar as a result of molecular dissociation and simultaneous ionization. Reflectivity measurements confirmed that deuterium transitions from an insulating state to a conducting phase at pressures near 0.5 Mbar on the Hugoniot and appears to do so in a continuous manner. At pressures above 0.5 Mbar, the measured reflectivity of the shocked deuterium is characteristic of a liquid metal.     

Molecular Dynamics Simulations of the Elastic Anisotropy of Pd at Extreme Conditions

It is very interesting to discover the elastic properties of engineering material palladium, especially its elastic anisotropy along Hugoniot states. We here investigate the evolution of its high pressure and temperature(PT) elastic ansotropy along Hugoniot using molecular dynamics simulations based on accurate classical interatomic potential. In order to testify the validity of the interatomic potential of Pd in describing the high PT elastic properties, we calculate its isothermal and adiabatic elastic moduli using molecular dynamics method. The obtained data are in good agreement with experimental data. From the isothermal elastic constants, we deduce the Hugoniot acoustic velocities and find that the resulting data are in good agreement with experimental acoustic velocity data. Based on the reliable elastic constants, we further investigate the spacial elastic ansotropy along Hugoniot PT states. It is found that the spacial elastic anisotropy of Pd increases along Hugoniot states.     

Aluminum Hugoniot from model calculations including semicore electrons based on density functional theory

We present a method to obtain Hugoniot from model calculations based on density functional theory, and apply the method to aluminum Hugoniot. Technological advances have extended the experimental research of high energy density physics, and call for quantitative theoretical analysis. However, direct computation of Hugoniot from density functional theory is very difficult. We propose two step calculations of Hugoniot from density functional theory. The first step is molecular dynamics simulations with an ambient temperature for electrons. The second step is total energy calculations of a crystal with desired high temperatures for electrons and with the ambient temperature for electrons. We treated the semicore 2s and 2p electrons of aluminum as valence electrons only for the total energy calculations of the aluminum crystal. The aluminum Hugoniot from our model calculations is in excellent agreement with available experimental data and the previous density functional theory calculations in the literature.     

Shock-wave compression of x-cut quartz as determined by electrical response measurements

The mechanical response of x-cut quartz in the vicinity of the Hugoniot elastic limit is determined from measurements of the piezoelectric current from samples impact loaded from 26 to 130 kbar. The Hugoniot elastic limit is determined to be 60_?1·5⁺³ kbar at a compression of 0·066_?0·002^+0·004 This Hugoniot elastic limit corresponds to a shear strength of 5·5 per cent of the C₄₄ shear modulus. For stresses well above the Hugoniot elastic limit the electrical current measurements show that the material exhibits a substantial reduction of shear strength. The pressure derivative of the bulk modulus is determined to be 4·5, substantially less than the ultrasonic value. The experimental records show evidence for a time delay for reduction of shear strength which varies from about 10^?7 sec immediately above the 60 kbar Hugoniot elastic limit to about 10^?8 sec for stresses well above the Hugoniot elastic limit. The measurements also show stress relaxation below the Hugoniot elastic limit between 40 and 60 kbar.     

Hugoniot线上的一个经验常数——来自多孔材料Hugoniot数据的对比研究

 对不同初始密度样品冲击压缩Hugoniot数据的系统分析显示，对每种固体物质而言存在一个表征材料特性的沿其Hugoniot线不随冲击压力而变的经验常数β。β的数值随材料而异。对金属元素而言，该经验特性常数β与初始密实密度ρ₀之间遵循同一个幂指数关系。利用该常数β，可很方便地对不同初始密度的Hugoniot数据进行换算，并可以很方便地获得Grüneisen参数。还对β值为常数的适用范围作了讨论。     

用于预测疏松材料冲击压缩特性的热力学模型

 在一定的假设条件下，从基本的热力学关系出发，导出了一个新的物态方程表达式。利用这一新形式特态方程所建立的热力学模型，可预测疏松材料的冲击压缩行为。以低、中、高冲击阻抗的铝、铜和钨为示例，将本模型的理论计算曲线与实验冲击压缩数据比较，表明本模型的在宽广的压力范围内具有良好的普活性；对不同初妈密度，疏松材料的冲击压缩行为均显示出较为满意的理论经历测能力；与现有的其它疏松材料物态方程模型相比，本模型在理论上和实际应用中均具有较大的优越性。     

A caloric equation of state of a condensed medium derived based on generalized Hugoniot curves

The caloric equation of state, which expresses the internal energy through the volume and pressure is required in analyzing and calculating shock-wave processes in reacting media. The present work develops a method for constructing such an equation for condensed media on the basis of generalized Hugoniot curves. Expressions approximating the dependence of the Grüneisen coefficient and pressure on the volume along the Hugoniot curve. Formulas for calculating detonation Hugoniot curves from the shock Hugoniot curve of the condensed reaction products. The expressions obtained can be applied to approximate calculations of shockwave processes (for example, gasless detonation) in reacting condensed media. If necessary, such calculations can be performed with the use of physical quantities from standard handbooks.     

纳米结构泡沫金冲击响应的分子动力学模拟

采用分子动力学计算程序对纳米结构泡沫金(Au)的冲击响应进行了模拟,得到了不同疏松度条件下泡沫Au的冲击压缩特性。通过获取不同势函数条件下实密Au的冲击Hugoniot关系以及泡沫结构稳定性测试选取适合描述Au泡沫冲击过程中原子的相互作用势。采用密堆积球壳的方式建立泡沫Au的初始构型。通过改变空心球壳的尺寸得到不同疏松度的稳定的泡沫Au结构。对泡沫Au的冲击过程进行分子动力学模拟,获得了不同疏松度泡沫Au在不同冲击压缩强度下的热力学状态参数。将模拟结果与已有的状态方程数据库以及疏松物质冲击压缩模型进行比较,结果表明,计算和理论模型给出的结果仍然存在明显的差异性,亟需通过进一步实验研究来验证模拟计算和理论模型结果的可靠性。     

氢及其同位素宽区状态方程的核量子效应

应用大规模第一性原理分子动力学数值模拟方法,系统研究了氢及其同位素在极端条件下的分子解离规律及状态方程,并给出了参数化的解函数拟合公式。结合分子、原子流体量子振动的一阶修正,揭示了氢及其同位素的核量子效应,解析获得了氢、氘分子在温稠密区域分子解离规律的差别。由第一性原理状态方程给出了雨贡纽曲线,与气炮、化学炸药、磁驱动、高能激光等实验数据相符合,并详细讨论了由核量子效应导致的氢、氘雨贡纽曲线的同位素效应。     

派热克斯玻璃的固态相变区物态方程

派热克斯玻璃的Hugoniot实验数据显示了在20GPa以下存在相变,我们认为这是一种固态-固态的一级相变,根据Gibbs相平衡条件并应用了相变前的低压物态方程,计算了这种固态-固态相变的熵增。计算表明熵增为负值,即相变为放热过程。在此基础上给出了相变区的物态方程,其中冷压为平台,内能计入了相变熵的影响。由相变区的物态方程得到的理论的与实验的Hugoniot曲线基本一致。