极端冲击下激波诱导附加电场加速金属/气体界面的经验模型

张嘉炜; 黄生洪

doi:10.11858/gywlxb.20180607

极端冲击下激波诱导附加电场加速金属/气体界面的经验模型

doi: 10.11858/gywlxb.20180607

张嘉炜,
黄生洪^*, ,

中国科学技术大学近代力学系材料力学行为和设计重点实验室，安徽合肥　230026

基金项目: 国家自然科学基金(U1530125)；科学挑战专题(TZ2016001)

详细信息

作者简介:
张嘉炜(1994－)，男，硕士研究生，主要从事界面不稳定性分子动力学研究. E-mail: jwzhang9@mail.ustc.edu.cn

通讯作者:
黄生洪(1974－)，男，博士，副教授，主要从事极端冲击动力学研究. E-mail: hshnpu@ustc.edu.cn

中图分类号: O521.2; O357.4
计量
- 文章访问数: 6911
- HTML全文浏览量: 3041
- PDF下载量: 29
出版历程
- 收稿日期: 2018-07-31
- 修回日期: 2018-09-04

Acceleration Evaluation Model of Metal/Gas Interface by Extra Electric Field Induced by Shock under Extreme Impacting Conditions

ZHANG Jiawei,
HUANG Shenghong^{*
, ,}

CAS Key Laboratory of Mechanical Behavior and Design of Materials, Department of Modern Mechanics, University of Science and Technology of China, Hefei 230026, China

摘要

摘要: 采用基于电子力场势函数的分子动力学方法，模拟了极端冲击压缩（22.50～78.75 km/s）下Li/H₂界面的演化过程。计算结果表明，在极端冲击压缩条件下，激波压缩导致物质电离，由于电子和离子的扩散差异，激波附近出现电荷分离区，进而诱导出现附加电场。通过沿激波传播方向进行一维电荷统计和理论分析，发现在激波强度一定的情况下，随着激波的运动，电荷分离区的强度和宽度在整个过程中基本保持恒定，与激波强度正相关；进一步对电荷分布进行积分，得到由电荷分离引起的附加电场沿冲击波传播方向的变化。统计分析了附加电场影响下金属/气体界面附近金属侧物质产生的加速度随时间的变化，发现附加加速过程呈脉冲形态，满足Rayleigh模型。通过计算拟合，得到该模型的关键参数与初始冲击加载速度的关系，最终获得电离诱导附加电场引起界面金属侧物质附加加速度的经验计算公式。
- 极端条件 /
- 电荷分离 /
- Richtmyer-Meshkov不稳定性 /
- 加速度 /
- 电场 /
- 分子动力学
Abstract: The evolution processes of metal/gas (Li/H₂) interface at extreme impacting conditions (22.50–78.75 km/s) were numerically studied by molecular dynamics (MD) method incorporated with the electron force field (eFF) model. It was found that the strong shock compression leads to ionization and the electron/ion separation is produced due to different diffusivities of ions and electrons. Then a strong extra electric field was established adjacent to shock font. Through 1D statistic along shock propagating direction from MD results and theoretical analysis, it was found that the electron/ion separation is moving with shock and the intensity and width of electron/ion separation zone are kept to be constant during shock propagating process and determined by shock strength. Further integrating the extra electric field and extra acceleration of metal material adjacent to the interface, the time histories of material acceleration were obtained. It was found that the extra material acceleration curves were in accordance with Rayleigh model. The key parameters were fitted based on computation results. Finally, an empirical extra acceleration evaluation model of metal material on Li/H₂ interface under impact velocity range of 20–80 km/s was established.
- extreme condition /
- electron/ion separation /
- Richtmyer-Meshkov instability /
- acceleration /
- electric field /
- molecular dynamics

HTML全文

当冲击波通过两种不同物质的界面时，界面扰动在冲击加载作用下会增长演化，即Richtmyer-Meshkov（RM）不稳定性^[1]。RM不稳定性是目前惯性约束核聚变研究的关键基础问题，其挑战主要来源于极端热动力学条件（压力远大于100 GPa，温度远高于10⁴ K）^[2-3]。目前，尚未发现针对极端条件的新RM不稳定性理论模型，进行相关实验需要高昂的费用，并且数据测量极为困难，所获数据十分有限。

数值模拟方面，对于温稠密物质及高密度等离子体，常规的流体动力学模拟缺少有效的模型和输运参数。不过，从原子尺度，采用分子动力学的建模和模拟方法正成为一种有效的技术手段。这类方法从描述原子或核子的基本相互作用势出发，不要求预设任何宏观模型、本构关系、宏观输运系数等，通过计算大量原子或核子的运动，直接模拟各种能量密度范围内的物质运动，可以涉及非平衡态、非连续效应、多组分/多相效应、电离以及与材料混合相关的随机反应等各种复杂物理机制。目前，在极端热力学条件下材料混合问题的理论研究中已经普遍采用分子动力学方法进行计算^[4-8]，特别是基于量子力学方程的第一性原理分子动力学方法（First Principles Molecular Dynamics，FPMD）^[9-19]，且正有计划地发展大规模分子动力学程序^[4]。然而，利用经典势函数的分子动力学模拟受到势函数的限制，一般只能模拟常规条件下或非电离状态下的RM不稳定性；基于密度泛函理论的QMD（Quantum Molecular Dynamics）方法可以计算极端条件下的物性参数和状态方程，但是由于计算量巨大，无法直接模拟RM不稳定过程。Su等^[20-21]提出将单电子显式处理成波包（Wave Packet，有位置，大小可变）的模型，从而发展出一种新的具有第一性原理特质的分子动力学方法—电子力场（Electron Force Field，eFF）分子动力学。近年来，该方法得到不断发展，已经具备处理10⁶量级离子和电子及长时间尺度的能力^[22]。

近期本研究团队利用该方法^[23-24]，从微观角度开展了不同热动力学条件下（3～30 km/s高速冲击加载）金属锂和气态氢材料混合的分子动力学数值模拟，结果发现极端加载条件下RM不稳定性的演化呈现新的特征，即振幅增长的线性段加长以及非线性段振幅增长速度高于无电离效应工况，这一结果得到了Nova实验结果^[25-26]的印证。其机理在于：极端条件下激波压缩导致物质电离，由于电子和离子扩散存在差异，电子将扩散到激波之前，形成电荷分离区，进而诱导产生附加电场，由附加电场引起的界面物质附加加速度会显著影响界面演化。为深入理解极端条件下RM不稳定性的演化特征，特别是附加电场对界面物质的加速效应，本研究继续采用eFF MD方法，模拟强激波作用下电荷分离情况，对电荷沿激波和界面运动方向的分布进行统计分析，得到电荷分离区强度和宽度随激波强度变化的规律，并进一步获得附加电场影响下界面附近金属侧物质产生的加速度，通过拟合，获得附加电场引起的界面附加加速度的计算模型。

1. 方法与模型

1.1 MD方法

MD方法的本质是通过牛顿第二定律求解每个原子（或分子）位置、速度、受力等基本物理量随时间的变化，其控制方程为

${{{F}}_i} = - {\nabla _i}U = - \left( {\frac{\partial }{{\partial {x_i}}}{{i}} + \frac{\partial }{{\partial {y_i}}}{{j}} + \frac{\partial }{{\partial {z_i}}}{{k}}} \right)U$

(1)

式中：i代表粒子编号；F_i为系统中其他粒子施加在粒子i上的力的总和，可以由粒子间的势函数 $U$ 求梯度得到。势函数的表达形式为

$U = U\left( {{{{r}}_1},{{{r}}_2}, \cdot \cdot \cdot ,{{{r}}_i}, \cdot \cdot \cdot ,{{{r}}_n}} \right)$

(2)

式中：r_i为第i个原子的坐标。获得施加在粒子i上的力之后，它相应的加速度a_i、速度v_i以及位置都可以由下列力学式计算得到

$\left\{ \begin{array}{l} {{{F}}_i} = {m_i}{{{a}}_i}\\ {{{v}}_i} = {{{v}}_i^0} + {{{a}}_i}t\\ {{{r}}_i} = {{{r}}_i^0} + {{{v}}_i^0}t + \dfrac{1}{2}{{{a}}_i}{t^2} \end{array} \right.$

(3)

式中：m_i为粒子质量，t为时间。物质其他相关的物理量，如体积、温度、压力等，均可以通过数值统计方法得到。

1.2 电子力场势函数

在目前的研究中，包含电子力场势函数的模型将原子核处理成点电荷^[20-21]，而将电子处理成满足高斯分布的波包。整个系统的能量是由原子核位置R、电子位置r和相对距离s构成的函数，表达形式如下

$U\left( {R,r,s} \right) = E{}_{{\rm{NN}}}\left( R \right) + {E_{{\rm{Ne}}}}\left( {R,r,s} \right) + {E_{{\rm{ee}}}}\left( {r,s} \right) + {E_{{\rm{KE}}}}\left( {r,s} \right) + {E_{{\rm{pauli}}}}\left( { \uparrow \downarrow ,s} \right)$

(4)

${E_{{\rm{NN}}}} = \frac{1}{{4{\text{π}} {\varepsilon _0}}}\sum\limits_{i < j} {\frac{{{Z_i}{Z_j}}}{{{R_{ij}}}}}$

(5)

${E_{{\rm{Ne}}}} = - \frac{1}{{4{\text{π}}{\varepsilon _0}}}\sum\limits_{i,j} {\frac{{{Z_i}}}{{{R_{ij}}}}} {\rm{erf}}\left( {\frac{{\sqrt 2 {R_{ij}}}}{{{s_j}}}} \right)$

(6)

${E_{{\rm{ee}}}} = \frac{1}{{4{\text{π}}{\varepsilon _0}}}\sum\limits_{i < j} {\frac{1}{{{r_{ij}}}}} {\rm{erf}}\left( {\frac{{\sqrt 2 {r_{ij}}}}{{\sqrt {{s_i}^2 + {s_j}^2} }}} \right)$

(7)

${E_{{\rm{KE}}}} = \frac{{{h^2}}}{{{m_\varepsilon }}}\sum\limits_i {\frac{3}{{2{s_i}^2}}}$

(8)

${E_{{\rm{pauli}}}} = \sum\limits_{{\sigma _i} = {\sigma _j}} {E{{\left( { \uparrow \downarrow } \right)}_{ij}}} + \sum\limits_{{\sigma _i} \ne {\sigma _j}} {E{{\left( { \uparrow \downarrow } \right)}_{ij}}}$

(9)

式中：E_NN表示原子核之间的库仑作用力，E_Ne表示原子核与电子之间的库仑力，E_ee表示电子之间的库仑力，这3项能量与经典分子动力学下的相一致；E_KE表示电子的动能，E_pauli表示Pauli排斥能，这两项能量反映了量子影响； ${\varepsilon _0}$ 和 $\sigma$ 是经典Lennard-Jones势函数下的参量； $Z$ 表示原子核的带电量； ${R_{ij}}$ 和 ${r_{ij}}$ 表示相互作用的距离；m_ε为电子质量，h为普朗克常数。单个电子被处理成满足高斯分布的波包，因而波包的位置x和大小s会随时间变化。按照运动和静止体系，其波函数形式为

$\psi \left( r \right) \propto \exp \left[ { - \left( {\frac{1}{{{s^2}}} - \frac{{2{P_s}}}{s}\frac{i}{h}} \right){{\left( { { r} - { x} } \right)}^2}} \right]\exp \left( {\frac{i}{h} {{{ P}_x}} \cdot { r}} \right)$

(10)

式中：P_x和P_s分别表示高斯波包的平动动量和径向伸缩动量。将（10）式代入时间相关的薛定谔方程，得到简化Hamiltonian的运动方程

$\dot x = \frac{{{\rm{d}}x}}{{{\rm{d}}t}} = \frac{{{P_x}}}{{{m_{\rm{e}}}}}$

(11)

${\dot P_x} = \frac{{{\rm{d}}{P_x}}}{{{\rm{d}}t}} = - \frac{{\partial V}}{{\partial x}}$

(12)

$\dot s = \frac{{{\rm{d}}s}}{{{\rm{d}}t}} = \frac{4}{3}\frac{{{P_{\rm s}}}}{{{m_{\rm{e}}}}}$

(13)

${\dot P_{\rm s}} = \frac{{{\rm{d}}{P_{\rm s}}}}{{{\rm{d}}t}} = - \frac{{\partial V}}{{\partial s}}$

(14)

式中：m_e为电子质量，V为系统势能。这些方程说明电子的位置和大小遵循经典分子动力学，并且独立地随时间变化。因此，eFF方法可以用来模拟长时间、大空间尺度各种极端条件下材料的非绝热动力学过程。

1.3 计算模型

如图1所示，计算模型的大小为123.760 nm×3.536 nm×1.768 nm (x×y×z)，由金属Li和气体H₂组成。Li和H₂之间的界面为光滑平面，无初始扰动。于是，冲击压缩过程可近似成一维流动。e⁻表示电子，±1代表电子的自旋方向。

图 1 计算模型

Figure 1. Computational model

下载: 全尺寸图片幻灯片

界面冲击激波首先由Li区域右侧刚性壁模拟活塞向−x方向以恒速v_p冲击压缩产生。当它经过界面位置时，会形成向左的透射激波和向右的反射激波，电荷分离将在此冲击压缩过程中产生。

2. 结果与讨论

2.1 电荷分离

为了清晰地看到电荷分离情况，选取初始压缩速度v_p=22.50 km/s的工况，如图2（a）所示。可见，强激波冲击运动过程中，激波压缩产生的高温导致电离，电离出自由电子和离子。由于激波压缩区的密度高，电子和离子浓度均高于激波未到达区。根据等离子体动力学，电子和离子将向低浓度区域扩散。由于电子质量极小，碰撞截面小，其扩散速度极高，且未电离原子对其扩散阻力小，因此电子很快跨越激波，扩散到激波未到达区；而离子与原子的质量基本相同，且碰撞截面大，扩散速度低，受其他原子碰撞后扩散阻力大，无法跨越激波。因此，在激波附近必然出现电子-离子分离区。一旦出现电子-离子分离区，等离子体的空间局部中性电场将不复存在，而随着空间电荷的分布出现局部附加电场。附加电场的建立将阻止电子进一步向低电荷密度区扩散，直至达到动态平衡为止，最终形成如图2（a）所示的台阶状电荷分布状态。

图2（b）进一步显示了电荷密度随时间的变化。可见，电荷分离现象伴随整个激波压缩过程，且随着激波向前运动而移动。这是因为运动激波不断压缩导致其所经过区域物质电离，电子和离子不断扩散，建立附加电场，并形成新的扩散平衡。可以确定，电荷分离区将伴随运动激波，以相同的速度运动。

2.2 电荷分离区的强度与宽度

根据图2所描述的电荷一维分布曲线，构建了如图3（a）所示的一维电荷分布模型。它由3部分构成：区域I₁电荷量为负，表示电子分布区域；区域I₃电荷量为正，表示离子区域；中间的区域I₂是从负（电子区域）到正（离子区域）的过渡区，可近似地看成线性分布，它随着激波不断向激波运动方向（x负方向）移动。

图 3 电荷密度分布一维模型以及电荷分离区宽度和强度的时间历程 (v_p=22.50 km/s)

Figure 3. 1D model of charge density distribution and time histories of the width and intensity of the electron-ion separation zone for v_p =22.50 km/s

下载: 全尺寸图片幻灯片

图 2 v_p=22.50 km/s时电荷密度分布一维统计及其随时间的变化

Figure 2. 1D statistic of charge density distribution for v_p=22.50 km/s case and its patterns at different moments

下载: 全尺寸图片幻灯片

过渡区宽度定义为电荷分离区的宽度，过渡区两端的电荷密度差（q₃−q₁）定义为电荷分离区的强度。图3（b）和图3（c）分别显示v_p=22.50 km/s的工况下电荷分离区宽度和强度随时间变化的曲线。可以看出，这两个新定义的量在激波运动过程中基本保持恒定，宽度维持在4.65 nm左右，强度在4.55×10⁸ C/m³附近波动。这是由于：分离区诱导的附加电场抑制了电子和离子的进一步扩散，达到一种局部平衡状态。而扩散动力来自于电离强度，当压缩激波强度一定时，其波后温度和压力一定，电离强度趋于恒定，则扩散动力恒定，因而分离区的强度和宽度必然趋于恒定。

表1进一步统计了其他加载速度下电荷分离区的强度和宽度数据。可见，随着冲击速度的增加，电荷分离区的宽度和强度持续增加。这是因为随着冲击速度的增加，激波压缩强度增加，电离强度增加，分离区宽度和强度自然相应增加。不过从增加的程度来看，宽度增加程度较小，而强度增加程度较高，且均呈现非线性增加特征。

表 1 不同冲击加载速度下分离区的宽度和强度统计

Table 1. Statistics of the width and intensity of the electron-ion separation zone at different impact velocities

v_p/(km·s^–1)	Width/nm	Intensity/(10⁸ C·m^–3)
22.50	4.65	4.55
45.00	6.40	6.11
67.50	8.65	10.25

下载: 导出CSV

| 显示表格

2.3 附加电场对界面附近金属物质的附加加速度模型

2.3.1 附加电场和附加加速度

如前所述，在极端加载条件下，电荷分离会诱导附加电场的产生。本研究根据模拟结果构建了如图3（a）所示的一维电荷分布模型。为验证该模型的合理性，对模型进行电势积分获得一维分布的电场，如图4中红实线所示。

由于电荷分离，激波前后的物质不再呈中性状态，而是分别带部分负电和正电。在附加电场中必然产生加速度，其计算表达式如下

$a = \frac{{{\rho _{\rm{q}}}E}}{{{\rho _{\rm{m}}}}}\quad\quad\quad\quad\quad\quad\quad\quad\quad(15)\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad$

(15)

式中：E为附加电场强度，ρ_q和ρ_m分别为单位体积下的物质带电量和物质密度，可通过分子动力学模拟结果的统计得到。

同时根据积分获得的一维分布的电场，计算各处带电物质获得的一维加速度分布，如图4中蓝色实线所示。图4还给出了根据eFF MD模拟结果一维统计获得的电场分布（红色虚线）和物质加速度分布（蓝色虚线），可见二者基本吻合。说明本研究建立的一维电荷分布模型基本合理。

图 4 一维附加电场和附加加速度(t=2.75 ps，v_p=22.50 km/s)的eFF MD统计分析结果和基于1D电荷密度模型的计算结果对比

Figure 4. Comparison of 1D extra electric field and extra acceleration of materials between results of eFF MD statistics and 1D charge density model (t=2.75 ps of v_p =22.50 km/s case)

下载: 全尺寸图片幻灯片

2.3.2 金属/气体界面附近金属物质的加速度模型

由于RM不稳定性的主要演化特征表现为界面物质扰动的振幅增长及物质混合，因此在极端条件下，激波诱导电荷分离导致附加电场的产生，进而引起界面物质的附加加速效应，可能显著影响界面演化。为深入理解极端条件下RM不稳定性的演化特征，特别是附加电场对界面物质的加速效应，对金属/气体界面附近金属侧物质（Li和H₂物质的量之比0.6～0.9，厚度在0～2.0 nm区间）的加速度进行跟踪统计（以激波到达界面位置为统计起点），获得了不同冲击加载速度v_p下界面附近金属物质的附加加速度随时间变化曲线，如图5所示。图5中虚线为根据eFF MD模拟结果统计分析获得，由于电子和离子的动态运动效应，统计值有一定涨幅落差，因此虚线呈现一定的振荡特征。对结果采用Rayleigh模型（（16）式）进行拟合，如图5中红实线所示。综合来看，可以获得以下几点认识。

图 5 不同冲击加载速度v_p下界面金属侧物质附加加速度的变化历程

Figure 5. Acceleration time histories of material located at the metal side of the interface at different v_p

下载: 全尺寸图片幻灯片

（1）界面附近金属物质获得的加速度呈脉冲形式。当激波刚到达界面位置时，电荷分离区跨越气/固界面，界面附近金属侧物质处于分离区中心区域，净电荷接近零；尽管此时分离区电场强度较高，但目标区受到的电场力较弱，加速度较低。随着激波往前推进，电荷分离区亦往前推进，分离区的电荷扩散处于动态平衡，目标区的净电荷逐渐偏正，并迅速向峰值增加；此时目标区受到的电场力开始增强，加速度迅速向峰值挺近，形成脉冲尖峰。随后，由于激波速度大于界面运动速度，随着激波继续向前推进，电荷分离区逐渐远离界面；此时目标区还带有较强的正电荷，但该区电场强度迅速衰减，对应的加速度也迅速衰减，形成完整的加速脉冲效应。

（2）随着冲击加载速度的增加（从22.50 km/s增加至78.75 km/s），目标区域的峰值加速度迅速提高（从1.7×10¹⁷ m/s²增加至1.2×10¹⁸ m/s²），同时脉冲持续时间逐渐减小（从1.5 ps降至1.0 ps）。这是因为随着冲击加载速度的提高，压缩程度增强，电离强度增加，电荷分离区强度相应提高，导致峰值加速度提升；但同时激波运动速度加快，激波跨越界面的时间相应缩短，脉冲持续时间自然有所降低。

（3）需要指出的是，曲线后期呈现较大幅度的振荡，这是因为激波到达左侧边界发生了Re-shock现象，Re-shock导致电荷分离呈现较大幅度的振荡。不过这一过程不是本研究关注内容，在此不作讨论。

本研究拟合所采用的模型公式为

$a\left( t \right) = {a_{\max }}{{\rm{e}}^{{1 / 2}}}\frac{t}{{{\tau _{{\rm{rise}}}}}}\exp \left( { - \frac{{{t^2}}}{{2{\tau _{{\rm{rise}}}}}}} \right)$

(16)

式中：a_max和τ_rise分别代表加速度峰值和上升时间，拟合参数如表2所示，拟合曲线如图6所示。

表 2 不同加载速度下Rayleigh模型的拟合参数

Table 2. The fitted parameters of Rayleigh model at different v_p

v_p/(km·s^–1)	a_max/(10¹⁷ m·s^–2)	τ_rise/ps
22.50	1.7	0.50
33.75	1.9	0.38
45.00	3.5	0.30
56.25	5.3	0.24
67.50	7.5	0.20
78.75	12.0	0.17

下载: 导出CSV

| 显示表格

图 6 Rayleigh模型的a_max和τ_rise与v_p的拟合关系

Figure 6. Fitted correlations of a_max and τ_rise with v_p in Rayleigh model

下载: 全尺寸图片幻灯片

从图6中可以明显看出a_max和τ_rise分别随初始加载速度呈指数上升和指数下降趋势。这反映了拟合模型的合理性，即峰值加速度及脉冲上升时间与冲击强度的非线性关系。进一步对两个参数进行拟合，得到如下表达式

${a_{\max }} = 6.233\;20 \times {10^{16}} \times {{\rm{e}}^{0.037\; 44{v_{\rm{p}}}}}\quad\quad\quad(17)\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\;\quad\quad\quad\quad\quad$

(17)

${\tau _{{\rm{rise}}}} = 0.773\;69 \times {{\rm{e}}^{ - 0.020\;35{v_{\rm{p}}}}}\quad\quad\quad\quad(18)\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\;\quad\quad\quad\;\;\;$

(18)

将（17）式和（18）式代入(16)式，可以获得界面附近金属侧物质受到激波压缩电离诱导形成附加电场的加速度经验公式

$a\left( {t,{v_{\rm{p}}}} \right) = 6.233\;20 \times {10^{16}} \times {{\rm{e}}^{0.037\;44{v_{\rm{p}}} + 0.5}} \times \frac{t}{{0.773\;69 \times {{\rm{e}}^{ - 0.020\;35{v_{\rm{p}}}}}}} \times \exp \left( { - \frac{{{t^2}}}{{1.547\;38 \times {{\rm{e}}^{ - 0.020\;35{v_{\rm{p}}}}}}}} \right)$

(19)

这是在本计算模型条件下冲击加载速度在20～80 km/s范围内激波压缩电离诱导附加电场引起的界面金属侧物质附加加速度的经验模型，揭示了极端加载条件下物质电离引起的附加电场对界面物质运动影响的基本特征。从积分效应来看，虽然其脉冲时间在皮秒量级，但所导致的界面物质的附加速度增量可达百千米每秒量级，显然不可忽视。当然，这只是界面物质加速的一方面，其他作用力耦合后，最终效应取决于综合效果。

3. 结　论

采用基于电子力场势函数的分子动力学方法，模拟了极端加载条件下金属/气体(Li/H₂)平滑界面运动过程，对激波压缩电离过程及电荷扩散在激波附近诱导的电荷分离效应进行了理论分析和结果统计，发现了电荷分离诱导附加电场及其对界面物质加速作用的现象，进一步开展了理论建模和经验拟合分析，获得了电离引起界面附近金属侧物质附加加速度的经验模型，主要获得如下结论：

（1）极端冲击加载下，激波诱导附加电场的基本机理是，极端冲击使材料电离出较轻的电子和较重的离子，电子和离子的扩散差异导致在激波附近形成电荷分离区并达到动态扩散平衡，诱导出附加电场；

（2）统计分析发现，电荷分离区的宽度和强度在激波运动过程中基本保持稳定，它们仅与初始冲击加载速度相关，即与激波强度相关，呈非线性关系；

（3）由附加电场引起的界面物质的附加加速效应可能显著影响界面演化，统计分析发现，界面处金属侧的加速度满足Rayleigh模型，呈脉冲形态，揭示了极端加载条件下物质电离引起的附加电场对界面物质运动影响的基本特征，从积分效应来看，虽然其脉冲时间在皮秒量级，但其影响却不可忽视。

需要指出的是，本研究结果基于微观分子动力学模拟，其尺度在百纳米量级，相关结论尽管揭示了内在机理，但在一定程度上还受尺度效应的限制。后续工作将进一步加大模拟尺度并开展理论建模，进而推广到宏观尺度范围，使之更贴合工程实际工况。

图 1 计算模型

Figure 1. Computational model

下载: 全尺寸图片幻灯片

图 3 电荷密度分布一维模型以及电荷分离区宽度和强度的时间历程 (v_p=22.50 km/s)

Figure 3. 1D model of charge density distribution and time histories of the width and intensity of the electron-ion separation zone for v_p =22.50 km/s

下载: 全尺寸图片幻灯片

图 2 v_p=22.50 km/s时电荷密度分布一维统计及其随时间的变化

Figure 2. 1D statistic of charge density distribution for v_p=22.50 km/s case and its patterns at different moments

下载: 全尺寸图片幻灯片

图 4 一维附加电场和附加加速度(t=2.75 ps，v_p=22.50 km/s)的eFF MD统计分析结果和基于1D电荷密度模型的计算结果对比

Figure 4. Comparison of 1D extra electric field and extra acceleration of materials between results of eFF MD statistics and 1D charge density model (t=2.75 ps of v_p =22.50 km/s case)

下载: 全尺寸图片幻灯片

图 5 不同冲击加载速度v_p下界面金属侧物质附加加速度的变化历程

Figure 5. Acceleration time histories of material located at the metal side of the interface at different v_p

下载: 全尺寸图片幻灯片

图 6 Rayleigh模型的a_max和τ_rise与v_p的拟合关系

Figure 6. Fitted correlations of a_max and τ_rise with v_p in Rayleigh model

下载: 全尺寸图片幻灯片

表 1 不同冲击加载速度下分离区的宽度和强度统计

Table 1. Statistics of the width and intensity of the electron-ion separation zone at different impact velocities

v_p/(km·s^–1)	Width/nm	Intensity/(10⁸ C·m^–3)
22.50	4.65	4.55
45.00	6.40	6.11
67.50	8.65	10.25

下载: 导出CSV

表 2 不同加载速度下Rayleigh模型的拟合参数

Table 2. The fitted parameters of Rayleigh model at different v_p

v_p/(km·s^–1)	a_max/(10¹⁷ m·s^–2)	τ_rise/ps
22.50	1.7	0.50
33.75	1.9	0.38
45.00	3.5	0.30
56.25	5.3	0.24
67.50	7.5	0.20
78.75	12.0	0.17

下载: 导出CSV

参考文献(26)

[1]	VELIKOVICH A L, DIMONTE G. Nonlinear perturbation theory of the incompressible Richtmyer-Meshkovinstability [J]. Physical Review Letters, 1996, 76(17): 3112. doi: 10.1103/PhysRevLett.76.3112
[2]	NAKAI S, TAKABE H. Principles of inertial confinement fusion-physics of implosion and the concept of inertial fusion energy [J]. Reports on Progress in Physics, 1996, 59(9): 1071. doi: 10.1088/0034-4885/59/9/002
[3]	ANDREWS M J. Workshop: research needs for material mixing at extremes: LA-UR-11-02565 [R]. Los Alamos: Los Alamos National Laboratory, 2011.
[4]	GRAZIANI F R, BATISTA V S, BENEDICT L X, et al. Large-scale molecular dynamics simulations of dense plasmas: the Cimarron Project [J]. High Energy Density Physics, 2012, 8(1): 105–131. doi: 10.1016/j.hedp.2011.06.010
[5]	ZHAKHOVSKII V, NISHIHARA K, ABE M. Molecular dynamics simulation on stability of converging shocks [C]//TANAKA K A, MEYERHOFER D D, MEYER-TER-VEHN J. Proceedings of the 2nd International Conference on Inertial Fusion Science and Applications. Elsevier, 2002: 106-109.
[6]	KADAU K, GERMANN T C, HADJICONSTANTINOU N G, et al. Nanohydrodynamics simulations: an atomistic view of the Rayleigh–Taylor instability [J]. Proceedings of the National Academy of Sciences, 2004, 101(16): 5851–5855. doi: 10.1073/pnas.0401228101
[7]	ZYBIN S V, ZHAKHOVSKII V V, BRINGA E M, et al. Molecular dynamics simulations of the Richtmyer-Meshkov instability in shock loaded solids [J]. AIP Conference Proceedings, 2006, 845(1): 437–441.
[8]	CHERNE F J, DIMONTE G, GERMANN T C. Richtmyer-Meshkov instability examinedwith large-scalemolecular dynamics simulations [J]. AIP Conference Proceedings, 2012, 1426(1): 1307–1310.
[9]	KOHANOFF J, HANSEN J P. Ab initio molecular dynamics of metallic hydrogen at high densities [J]. Physical Review Letters, 1995, 74(5): 626. doi: 10.1103/PhysRevLett.74.626
[10]	KOHANOFF J, HANSEN J P. Statistical properties of the dense hydrogen plasma: an ab initio molecular dynamics investigation [J]. Physical Review E, 1996, 54(1): 768. doi: 10.1103/PhysRevE.54.768
[11]	KRESSE G, HAFNER J. Ab initio molecular dynamics for liquid metals [J]. Physical Review B, 1993, 47(1): 558. doi: 10.1103/PhysRevB.47.558
[12]	COLLINS L, KWON I, KRESS J, et al. Quantum molecular dynamics simulations of hot, dense hydrogen [J]. Physical Review E, 1995, 52(6): 6202. doi: 10.1103/PhysRevE.52.6202
[13]	龚新高. 高温及高压下液体镓的结构——第一性原理分子动力学方法研究 [J]. 物理学报, 1995, 44(6): 885–896 GONG X G. Structural properties of liquid gallium at high temperature and high pressure: an ab initio molecular dynamics study [J]. Acta Physica Sinica, 1995, 44(6): 885–896
[14]	何以广. 氢和氦高压物性的第一原理分子动力学研究及实验探索 [D]. 北京: 清华大学, 2010.
[15]	张玉娟. 温稠密乙烷等流体物性的第一性原理分子动力学研究 [D]. 北京: 中国工程物理研究院, 2013.
[16]	刘海, 李启楷, 何远航. 高速冲击压缩梯恩梯的分子动力学模拟 [J]. 力学学报, 2015, 47(1): 174–179 LIU H, LI Q K, HE Y H. Molecular dynamics simulations of high velocity shock compressed TNT [J]. Chinese Journal of Theoretical and Applied Mechanics, 2015, 47(1): 174–179
[17]	宋海峰, 刘海风. 金属铍热力学性质的理论研究 [J]. 物理学报, 2007, 56(5): 2833–2837 doi: 10.3321/j.issn:1000-3290.2007.05.060 SONG H F, LIU H F. Theoretical study of thermodynamic properties of metal Be [J]. Acta Physica Sinica, 2007, 56(5): 2833–2837 doi: 10.3321/j.issn:1000-3290.2007.05.060
[18]	王聪, 贺贤土, 张平. 温稠密物质的量子及半经典分子动力学研究 [C]//中国力学大会2013论文摘要集. 西安, 2013.
[19]	戴佳钰, 康冬冬, 侯永, 等. 高温稠密物质的多尺度动力学研究 [C]//中国力学大会2013论文摘要集. 西安, 2013.
[20]	SU J T, GODDARD III W A. Excited electron dynamics modeling of warm dense matter [J]. Physical Review Letters, 2007, 99(18): 185003. doi: 10.1103/PhysRevLett.99.185003
[21]	SU J T, GODDARD III W A. The dynamics of highly excited electronic systems: applications of the electron force field [J]. The Journal of Chemical Physics, 2009, 131(24): 244501. doi: 10.1063/1.3272671
[22]	JARAMILLO-BOTERO A, SU J, QI A, et al. Large-scale, long-term nonadiabatic electron molecular dynamics for describing material properties and phenomena in extreme environments [J]. Journal of Computational Chemistry, 2011, 32(3): 497–512. doi: 10.1002/jcc.21637
[23]	王维荣. 极端条件下单模界面不稳定性的分子动力学研究 [D]. 合肥: 中国科学技术大学, 2017.
[24]	HUANG S, WANG W, LUO X. Molecular-dynamics simulation of Richtmyer-Meshkov instability on a Li-H₂ interface at extreme compressing conditions [J]. Physics of Plasmas, 2018, 25(6): 062705. doi: 10.1063/1.5018845
[25]	DIMONTE G, REMINGTON B. Richtmyer-Meshkov experiments on the Nova laser at high compression [J]. Physical Review Letters, 1993, 70(12): 1806. doi: 10.1103/PhysRevLett.70.1806
[26]	DIMONTE G, FRERKING C E, SCHNEIDER M, et al. Richtmyer–Meshkov instability with strong radiatively driven shocks [J]. Physics of Plasmas, 1996, 3(2): 614–630. doi: 10.1063/1.871889

施引文献

资源附件(0)

访问统计

点击查看大图

图(6) / 表(2)

计量

文章访问数: 6911
HTML全文浏览量: 3041
PDF下载量: 29

1. 方法与模型
1.1 MD方法
1.2 电子力场势函数
1.3 计算模型
2. 结果与讨论
2.1 电荷分离
2.2 电荷分离区的强度与宽度
2.3 附加电场对界面附近金属物质的附加加速度模型
3. 结　论

极端冲击下激波诱导附加电场加速金属/气体界面的经验模型

doi: 10.11858/gywlxb.20180607

作者简介: 张嘉炜(1994－)，男，硕士研究生，主要从事界面不稳定性分子动力学研究. E-mail: jwzhang9@mail.ustc.edu.cn

通讯作者: 黄生洪(1974－)，男，博士，副教授，主要从事极端冲击动力学研究. E-mail: hshnpu@ustc.edu.cn

计量

出版历程