基于MULTI2D-Z程序的Z箍缩动态黑腔形成过程模拟

对辐射流体力学程序MULTI-2D进行改造,增加磁场演化方程程序模块,自洽地在运动方程模块中增加洛伦兹力,在能量方程模块中增加欧姆加热,将它改造成辐射磁流体力学程序MULTI2D-Z.验证了新增磁场程序模块的可靠性,并发现温度和密度的增大会抑制磁场的扩散,负径向速度梯度的流体对流也会抑制磁场的扩散.利用改造好的MULTI2D-Z程序模拟了峰值为8 MA的脉冲电流驱动的钨丝阵Z箍缩动态黑腔形成过程.得到了X光功率(约30 TW)和能量(约300 k J)、泡沫辐射温度(约120 eV)、箍缩轨迹等模拟结果.在动态黑腔形成过程中,磁场主要分布在钨主体等离子体中;辐射向内传播,烧蚀泡沫柱而使它膨胀;辐射热波在被撞击的泡沫柱中传播,其传播速度比物质温度传播得快,当辐射热波传播到中心轴时泡沫柱中的辐射场变得比较均匀,并且除了冲击波处外辐射温度与物质温度基本上没有分离.这些模拟结果可增强人们对磁场扩散和对流规律以及动态黑腔形成机制的理解,同时表明了MULTI2D-Z程序可成为Z箍缩及其应用的新的程序模拟工具.

Simulation of forming process of Z-pinch dynamic hohlraum based on the program MULTI2D-Z

The radiation hydrodynamics code MULTI-2D, which was developed by Ramis et al. in 2009 (2009 Comput. Phys. Commun. 180 977) and adopted the single temperature fluid and unstructured lagrangian mesh, is modified into a radiation magnetohydrodynamics code MULTI2D-Z by adding the program module of evolution equation of magnetic field, and self-consistently considering the Lorentz force in the module of motion equation and the Ohmic heating in the module of energy equation. The newly developed module for magnetic field was validated to be reliable. The module is used to study the magnetic field diffusion process, and it is found that the diffusion is weakened due to the increasing of plasma temperature and density and the fluid convection, in which there is minus grads of velocity in radial direction. The new code MULTI2D-Z is used to simulate the formation process of dynamic hohlraum driven by tungsten wire-array Z-pinch at an 8 MA current level. The obtained results are that X-ray power and energy are, respectively, ～30 TW and ～300 kJ, radiation temperature in foam is ～120 eV, and the implosion trajectory of wire-array is also obtained. The calculated results reveal that the magnetic field is mainly distributed in the outside of tungsten plasma during the hohlraum formation. The foam expands due to the radiation heating from the shock wave created by the collision between wire-array plasma and the foam. The thermal radiation wave, which is characterized by radiation temperature, spreads towards the central axis faster than the plasma temperature. When the thermal radiation wave spreads to the central axis, the radiation temperature becomes comparatively uniform in space, and is almost equal to the plasma temperature except at the place of the shock wave. These results help the people to better understand the magnetic field diffusion and convection in Z-pinch, as well as the formation mechanism of dynamic hohlraum driven by wire-array Z-pinch. It is also indicated that the newly developed code MULTI2D-Z can be considered as a new tool for simulating Z-pinch and its applications, such as inertial confinement fusion and magnetically accelerated flyer plates.