微型磁探针在Z箍缩负载电流诊断中的应用

介绍了磁探针测量等离子体电流的设计原理,针对Z箍缩实验负载的实际结构特点和现场布局,制作出应用于诊断脉冲功率装置Z箍缩实验负载通过电流的微型磁探针,并通过建立相同结构尺寸的模拟负载装置的方法,实现了对其灵敏度的标定。实验结果显示:在脉冲功率装置峰值电流1.2MA、电流上升时间60ns时,由微型磁探针测得的负载电流与加速器监测电流存在12%的幅度差异,电流峰值时刻存在5ns的差异,说明微型磁探针技术测量Z箍缩负载电流的结果是可靠的。     

磁驱动飞片发射实验结构系数初步研究

为了确定磁驱动飞片发射实验结构系数的范围、影响因素、结构系数与影响因素的关系，对聚龙一号装置上的磁驱动飞片发射实验进行了数值模拟和分析。数值模拟表明，磁流体力学程序能正确模拟聚龙一号装置上各个磁驱动飞片发射实验；磁驱动双侧飞片发射实验的结构系数为0.7～0.8；磁驱动单侧飞片发射实验的结构系数为0.80～0.85。磁驱动飞片发射实验的结构系数与实验加载电流无关，仅由磁驱动飞片发射实验的负载结构决定。磁驱动飞片发射实验的结构系数取决于阴阳电极极板的初始宽度、阴阳电极之间的初始间隙以及阴阳电极上飞片厚度之和等三个因素。在磁驱动飞片发射实验中，电极初始宽度、阴阳电极之间的初始间隙不变的情况下，结构系数由阴阳电极上飞片厚度之和确定，阴阳电极上飞片厚度之和越大，结构系数越大。     

Z箍缩装置外磁绝缘传输线全尺寸粒子模拟研究

采用电路模拟得到的柱孔结构处向外磁绝缘传输线传输的电压反射波来等效Z箍缩装置中的柱孔结构至丝阵负载部分,实现了Z箍缩装置四层外磁绝缘传输线的全尺寸粒子模拟.为了进一步提高柱孔结构和丝阵负载等效电路模型的精度,通过粒子模拟,对电路模拟得到的电压反射波进行了修正.PBFA Z装置四层外磁绝缘传输线部分的全尺寸粒子模拟结果表明,修正电压反射波后得到的绝缘堆处电压波形和电流波形比原有电路模拟结果更接近实验结果.另外,利用粒子模拟结果分析和解释了丝阵负载内爆对外磁绝缘传输线脉冲功率传输物理过程的影响. 关键词： Z箍缩 外磁绝缘传输线 粒子模拟 等效电路模型     

PTS装置Z箍缩负载设计分析

介绍了PTS装置丝阵负载设计思路,用简化的薄壳模型(零维)和辐射磁流体力学模型(一维)数值模拟了丝阵Z箍缩内爆过程,给出了PTS丝阵负载设计参数,分析了初步优化结果,确定了优化选择的箍缩常数.经脉冲功率驱动器与负载优化匹配设计,PTS实验平台能够获得100TW量级的X射线功率,最佳插头能量转换效率10%.预计PTS实验平台的X射线输出能力基本达到Saturn装置长脉冲模式放电的水平.结果将为PTS开展Z箍缩研究提供有价值的参考. 关键词： PTS装置 Z箍缩丝阵负载 零维薄壳模型 一维辐射磁流体力学数值模拟     

Z箍缩内爆的MARED程序数值模拟分析

MARED程序是模拟Z箍缩内爆过程的二维三温辐射磁流体力学程序,它适用于不同装置条件和不同负载参数.利用MARED程序对Z箍缩内爆进行模拟,结合丝阵Z箍缩实验分析表明:相同负载质量条件下,钨丝阵内爆产生的X射线辐射功率远大于铝丝阵产生的X射线功率;相同负载电流条件下,负载质量越大,计算得到的X射线功率越低;X射线功率随着负载电流增加而增加.MARED程序能够较好地反映Z箍缩内爆动力学过程,特别是不稳定性发展的二维图像,能够给出与不稳定性简化模型的理论分析及实验结果定性一致的演化规律.MARED程序模拟丝阵填充泡沫形成辐射场的初步计算得到了与Sandia实验室模拟Z装置上丝阵填充泡沫定性一致的结果.     

Z箍缩内爆的MARED程序1维模拟分析

 研制可靠的数值模拟工具对Z箍缩内爆产生X光辐射过程进行理论研究、实验分析以及负载设计至关重要。介绍了2维三温辐射磁流体力学程序（MARED）的物理方案,给出了MARED程序的1维检验结果,验证表明MARED程序适用不同装置条件、不同负载参数。结合丝阵Z箍缩实验的数值模拟和分析表明：相同负载质量条件下,钨丝阵内爆产生的X光辐射功率远大于铝丝阵产生的X光功率;相同负载电流条件下,负载质量越大,计算得到X光功率越低;X光功率随着负载电流增加而增加。     

X箍缩三维数值模拟

利用FOI-PERFECT程序对X箍缩进行了3D数值模拟研究,给出了X箍缩的物理图像和动力学过程,探讨了Z箍缩中出现磁重联的可能性,并指出如果双丝Z箍缩中能够出现磁重联,那么X箍缩是更有利于磁重联出现的位形,并且,X箍缩中出现多重X射线暴的一个可能原因是z轴上多个位置出现磁重联。     

二维弹塑性磁流体力学数值模拟

为了研究物质弹塑性对磁驱动实验运动过程、不稳定性发展等的影响,在MDSC2程序的基础上,增加了弹塑性模块,研制了包括弹塑性的磁流体力学程序,并进行了弹塑性项影响的数值模拟和分析。数值模拟表明:没有初始扰动时,弹塑性项几乎不影响套筒内外半径的运动轨迹;有初始扰动时,弹塑性项对磁驱动固体套筒的Rayleigh-Tayor不稳定性有明显的抑制作用。     

二维磁驱动数值模拟程序MDSC2的验证与确认

为了对磁驱动实验提供高置信度的数值模拟,需要开展磁流体力学程序的验证与确认。采用人为解比较法、网格收敛性研究和与成熟程序比较等方法,对二维磁驱动数值模拟程序MDSC2进行了程序验证。数值模拟表明:MDSC2程序正确地表示了磁流体力学模型,其中热扩散、磁扩散的离散格式具有二阶收敛精度。采用与磁驱动实验相比较的方法,进行了MDSC2程序的确认。对聚龙一号装置上的PTS-061发次磁驱动单侧飞片发射和PTS-122发次磁驱动双侧飞片发射实验进行了模拟,模拟的飞片自由面速度与实验测量的飞片自由面速度相一致;对FP-1装置上的固体套筒实验进行了模拟,模拟的套筒内外半径与实验测量结果相一致。MDSC2程序能正确模拟磁驱动单侧飞片发射、磁驱动双侧飞片发射和磁驱动固体套筒等磁驱动实验。     

紧凑型X箍缩脉冲功率发生器

为研究丝阵负载对X箍缩的影响,设计了一台紧凑型脉冲功率发生器。介绍了该脉冲功率发生器的设计及初步实验。根据装置的结构特点,研制了一种利用金属膜连接传输线外筒与负载内筒,构成回路测量负载电流的探头。结果表明:在对Marx电容器充电60 kV时,输出电流峰值117 kA,脉宽 60 ns,上升沿22 ns,该装置已经成功箍缩出X射线。     

Simulation of Z-Pinch Processes of Nested Tungsten Wire-Array on Angara-5-1 Facility

Based on the hydrodynamic shell-on-shell collision model, the Z-pinch processes of nested tungsten wire-array in Sino-Russian joint experiments on Angara-5-1 facility are simulated by means of our one-dimensional three-temperature radiation magneto-hydrodynamic code. The results show the evolutions of x-ray radiation burst, implosion trajectories of interfaces, current transfer in inner and outer wire-array plasmas, and the temporal and spatial changes of magnetic field and current density in the process. About 20% of the total driven current is transferred into the inner wire-array plasma by convection and diffusion of magnetic field when the two shells are pinched closest. Compared to the measured x-ray power, the simulated full width at half maximum and time at the strongest radiation agree approximately with the measured values. It is also demonstrated in our simulation that the radiation of nested wire-array Z-pinch is enhanced. The effects of fluctuations of driven current on yields of x-ray are also investigated.     

毛细管放电类氖氩69.8nm激光增益特性研究

建立了计算69.8nm激光增益系数的理论模型,根据实验参数,计算了在主脉冲电流为12 kA时,69.8nm激光增益系数最大值为0.32 cm~(-1).理论模拟了不同初始气压下增益系数在毛细管径向上的分布情况.对理论结果的分析表明,最佳的初始气压在12—14 Pa范围内,此时69.8nm激光增益系数的极值最大.实验上,利用毛细管放电装置和罗兰光谱仪,测量了不同气压下的69.8nm激光强度,实验确定的最佳气压为16 Pa,与理论结果相近.此外,实验测量的增益系数(0.4 cm~(-1))略高于理论计算的增益系数(0.32 cm~(-1)).     

带窗口磁驱动准等熵压缩实验模拟

在二维磁驱动数值模拟程序MDSC2中增加了LiF材料的材料参数和功能模块,使MDSC2程序具有了求解带窗口磁驱动准等熵压缩实验的能力。采用MDSC2程序,对大电流脉冲功率装置上的exp-3-window、exp-6-window带窗口磁驱动准等熵压缩实验进行了模拟。数值模拟结果表明,二维磁驱动数值模拟程序MDSC2能正确模拟带窗口磁驱动准等熵压缩实验exp-3-window和exp-6-window的全过程,模拟的飞片/窗口界面速度在飞片/窗口界面速度的上升阶段、峰值附近和卸载阶段与实验测量基本一致,验证了新程序的计算有效性。MDSC2程序对带窗口磁驱动准等熵压缩实验的正确模拟有助于磁驱动样品物性实验的研究。     

基于等离子体粒子模拟的喷气Z箍缩过程物理研究

给出了喷气Z箍缩动力学过程在二维柱坐标系下的等离子体粒子模拟物理模型,编写了相应的程序.对低电流驱动下的稀薄喷气Z箍缩动力学过程进行了验证性的等离子体粒子模拟,得到了许多微观的Z箍缩物理信息,如负载中的电流(密度)、电磁场、粒子位置和密度的时空演化,以及总的z箍缩拖尾质量和拖尾电流等信息.发现在Z箍缩过程中,模拟得到的等离子体电流随时间的变化反映出了等离子体箍缩到心和反弹的过程特征,磁场随径向的变化与长直导线电流给出的磁场很接近;电子所受到的电场力和磁场力(洛伦兹力)是相当的,而离子所受到的力主要是电场力;电子首先在z方向加速,然后在自身运动产生电流的磁场的作用下向轴心箍缩,而离子是在电子和离子电荷分离所产生的电场力的作用下向轴心运动;在压缩到轴心附近时,电子首先因静电排斥而飞散,而离子则在惯性的作用下继续向轴心箍缩,而后滞止飞散.Z箍缩等离子体的拖尾质量在20%左右,拖尾电流最大时在7%左右.     

Reversed current structure in a Z-pinch plasma

The current profile of a Z-pinch plasma is investigated using a one-dimensional magnetohydrodynamic code. Simulation results reveal the formation of a reversed current profile, its propagation, and an ejection of plasma at boundary region, which have been observed in previous experiments. A new physical mechanism is proposed to account for such phenomena. The physical mechanism involves the propagation of a shock wave. It is found that a reversed current profile appears when a shock wave reflected at axis expands in a compressing plasma column.     

Implosion dynamics of a radiative composite Z-pinch

Two-dimensional simulation of a composite Z-pinch was performed by the complete radiative magnetohydrodynamic (MHD) code ZETA including detailed calculation of equations of state, spectral properties of materials, and radiation transport in non-local thermodynamic equilibrium multicharged ions plasma. The initial geometry, the substance components, and the electric current through the Z-pinch were similar to the joint experiment set up JEX-94 at Angara-5 facility. The geometry was: annular argon gas puff with inner and outer diameter of a nozzle cross section, 3 cm and 3.4 cm, respectively, and specific mass of 80 μg/cm. Inside it along the axis a foam cylinder 30% KCl in agar-agar with total mass 60 μg and diameter 1 mm was put. The initial gas distribution was modelled with a divergence of a jet along the Z-axis. A coupling of the Z-pinch with the electric current generator was modelled by an electrical circuit with a given electromagnetic wave, a resistance, an inductance, and a variable load (Z-pinch) similar to the Angara-5-1 output parameters. During the plasma implosion the total current reached the value of 3 MA at a time of 85 ns from the voltage start. Such current amplitude is much less than through the matched load (up to 4 MA as a rule): The plasma implosion is accompanied by the development of different types of short and long wave instabilities (thermal, radiative, nonisothermal, and MHD Rayleigh-Taylor modes). In this report, the detailed plasma implosion dynamics, the influence of instabilities, and the spectral radiation yield are discussed, and a comparison with the experimental results is done     

“聚龙一号”装置上磁驱动铝飞片实验的数值模拟

采用二维磁驱动数值模拟程序MDSC2,对大电流脉冲装置"聚龙一号"上的151发次磁驱动370μm厚铝飞片实验、164发次磁驱动330μm厚铝飞片实验进行了数值模拟和分析。数值模拟表明:370μm厚铝飞片和330μm厚铝飞片的磁驱动过程中,VISAR测量的速度不是飞片自由面的速度,而是飞片中邻近自由面最近的固体反射面的速度。这是由于磁驱动飞片发射过程中,飞片自由面部分被烧蚀,密度低于固体密度状态,而飞片自由面和加载面中间的飞片还保持固体密度状态。VISAR测量的激光将穿过自由面的低于固体密度状态的飞片部分,到达飞片自由面最近的固体密度位置再反射回去,获得这一位置的飞片速度。数值模拟的飞片固体反射面速度历史和VISAR测量的速度历史相吻合。     

Characteristics of high-power radiating imploding discharge with cold start

A qualitative model of the dynamics of a multiterawatt radiating Z-pinch with cold start and high rate of current rise is proposed. The model is used to analyze discharges with currents I ～ 2–5 MA (with dI/dt > 10¹³ A/s) through uniform or structured plasma-producing loads, including wire arrays. The most important consequence of cold start is that spatially nonuniform plasma production is prolonged to almost the entire current rise time. Under these conditions, the Ampére force begins to play a dominant role in the plasma dynamics before the plasma-producing load is completely transformed into an accelerated plasma. The results of computations of wire-array vaporization are presented. A formula is proposed for estimating the highest attainable velocity of plasma flow into a heterogeneous liner driven by the Ampére force. It is shown that local imbalance between radial motion of the produced plasma and supply of the plasma-producing substance to be ionized leads to axially nonuniform breakthrough of magnetic flux into the liner, which precedes plasma collapse. The magnetic-flux breakthrough gives rise to a chaotic azimuthal-axial plasma structure consisting of radial plasma jets of relatively small diameter, which is called a radial plasma rainstorm. The breaking-through azimuthal magnetic flux obstructs further current flow in the breakthrough region. Analyses of Z-pinch implosion based on the theory of Rayleigh-Taylor instability or the snowplow model are incorrect under the plasma-rainstorm conditions. The processes taking place in a stagnant Z-pinch include conversion of the energy carried by the current-generated magnetic field into turbulent MHD flow of the ion component of the plasma, its convective mixing with magnetic field, heating, energy transfer from ions to electrons, and emission from the plasma. Under typical experimental conditions, emission plays a key role in the energy balance in an imploding pinch. Z-pinch is modeled by an electric-circuit component that has a time-dependent nonlinear impedance and consumes the magnetic energy supplied by a generator through a magnetically insulated transmission line (MITL). The peak power reached in the circuit is comparable to the peak soft X-ray power output emitted by the pinch in terms of magnitude and timing. Optimum matching conditions are formulated for the generator-MITL-pinch circuit.