直接驱动柱内爆流体力学不稳定性数值模拟

给出了LARED-S程序对直接驱动柱内爆流体力学不稳定性模拟的物理方程,介绍了一种直接驱动柱内爆靶的结构,给出了模拟中使用活动网格追踪和其他一些需要注意的问题。对OMEGA激光器直接驱动柱内爆实验物理模型进行了简化处理：使用理想气体状态方程,假设烧蚀层、标识层和内部填充材料为同一介质,只在密度上有差别;采用全电离和单温近似,在内爆过程中仅考虑热传导过程而忽略辐射和非局域电子热传导等过程的影响。模拟结果表明：LARED-S程序与LASNEX程序的计算结果基本上符合。     

柱面内爆驱动金属界面不稳定性的数值模拟研究

对柱面爆轰驱动内壁刻有正弦扰动的金属钢壳与内部硅橡胶界面产生不稳定性问题进行数值模拟,计算结果与实验结果定性符合.与不考虑金属强度情况对比分析认为,未熔化状态下金属强度对不稳定性具有较强抑制作用,在某些加载条件下会使扰动增长率随扰动模数增加而减小.之后,对强度因素影响下内爆压缩驱动金属不稳定性问题的扰动发展规律进行了总结.在聚心反射波到达壳体之前,造成初始界面反转的RM不稳定性起主导作用,随着扰动模数增加扰动由呈近似线性发展到基本不发展变化,基本不变化后的扰动振幅也随模数增加而减小.聚心反射波作用到壳体内界面后,减速RT不稳定性作用明显增强,与强度等因素共同作用造成扰动呈明显非线性发展.无论是前期RM不稳定性主导阶段还是之后以减速RT不稳定性为主的扰动发展阶段,强度因素均能造成未熔化状态下金属不稳定性截止波长存在.     

直接驱动柱内爆流体力学不稳定性数值模拟

 给出了LARED-S程序对直接驱动柱内爆流体力学不稳定性模拟的物理方程，介绍了一种直接驱动柱内爆靶的结构，给出了模拟中使用活动网格追踪和其他一些需要注意的问题。对OMEGA激光器直接驱动柱内爆实验物理模型进行了简化处理：使用理想气体状态方程，假设烧蚀层、标识层和内部填充材料为同一介质，只在密度上有差别；采用全电离和单温近似，在内爆过程中仅考虑热传导过程而忽略辐射和非局域电子热传导等过程的影响。模拟结果表明：LARED-S程序与LASNEX程序的计算结果基本上符合。     

柱几何Rayleigh-Taylor不稳定性的数值模拟

给出了柱几何中流体力学方程组及其在数值模拟中采用的计算方法。对二维柱几何Rayleigh-Taylor不稳定性进行数值模拟,在线性阶段与线性理论符合得很好;不稳定性增长进入非线性区域的阈值依赖于界面的位置,并且明显不同于平面情况。     

ICF流体力学不稳定性实验用柱状激波管的研制

在惯性约束聚变(ICF)实验中,点火靶丸表面(界面)的粗糙度和缺陷所产生的流体力学不稳定性是决定点火成功与否的关键因素之一,设计和研制流体力学不稳定性分解实验用靶是解决该问题的主要技术手段。结合国内外的研究现状和神光-Ⅱ激光装置的特点,设计并研制了一种新型柱状激波管。该靶型由三种介质组成,分别为调制聚苯乙烯(CH)圆片、柱状碳气凝胶(CRF)和CH微套管。调制CH圆片和柱状CRF通过微加工技术装配到CH微套管内,封装后形成柱状激波管。介绍了该靶型的设计原理和详细的制备工艺,并对相应的靶参数进行了测量。结果表明:柱状CRF气凝胶具有较好的成型性,长度、直径和密度分别为1000μm、730μm和250mg·cm-3;CH圆片的厚度和直径分别为15μm和730μm,表面调制图形的周期和峰谷差分别为100μm和4.3μm;实验得到的柱状激波管的轴向和径向最大装配误差分别为2μm和3μm。     

储存环束流Robinson不稳定性的分析和测量

探讨了储存环束流的Robinson不稳定性问题,提出用“等效失谐角大于零”取代“失谐角大于零”作为束流稳定的基本条件。在合肥光源电子储存环200 MeV注入状态下,对束流不稳定性与高频腔失谐之间的关系进行了实验测量。结果表明：当束流稳定条件不满足时,如果高频腔大失谐,束流将全部丢失;小失谐时束流容易部分丢失;当高频腔处于负失谐状态,束流流强将限制在较低水平。     

激光间接驱动惯性约束聚变内爆物理实验研究

激光间接驱动惯性约束聚变利用辐射烧蚀驱动靶丸球形内爆,在减速阶段将内爆动能转化成热斑内能,同时压缩燃料,达到点火条件,实现聚变点火。根据目前认识,影响内爆压缩过程的主要因素包括内爆对称性、燃料熵增因子、内爆速度和混合。内爆物理实验研究的目的是发展对上述影响因素的实验表征方法,获取这些影响因素随靶设计参数的变化规律,建立相应的实验调控能力,最终达到不断提升内爆性能的目的。为此,在内爆对称性方面,开展了Bi球自发光实验,用于研究点火脉冲前2 ns驱动不对称性;在内爆速度方面,开展了球面弯晶单能流线实验,测量得到内爆速度和剩余质量随时间的变化;在混合方面,开展了内壳层示踪涂层内爆混合实验,测量得到环形发光图像。为考察综合内爆性能,在神光Ⅱ和神光Ⅲ原型装置上开展了DT内爆实验,获得了中子产额随初始靶参数的变化规律。     

Effects of the initial perturbations on the Rayleigh–Taylor–Kelvin–Helmholtz instability system

The effects of initial perturbations on the Rayleigh–Taylor instability (RTI), Kelvin–Helmholtz instability (KHI), and the coupled Rayleigh–Taylor–Kelvin–Helmholtz instability (RTKHI) systems are investigated using a multiple-relaxation-time discrete Boltzmann model. Six different perturbation interfaces are designed to study the effects of the initial perturbations on the instability systems. It is found that the initial perturbation has a significant influence on the evolution of RTI. The sharper the interface, the faster the growth of bubble or spike. While the influence of initial interface shape on KHI evolution can be ignored. Based on the mean heat flux strength D_3,1, the effects of initial interfaces on the coupled RTKHI are examined in detail. The research is focused on two aspects: (i) the main mechanism in the early stage of the RTKHI, (ii) the transition point from KHI-like to RTI-like for the case where the KHI dominates at earlier time and the RTI dominates at later time. It is found that the early main mechanism is related to the shape of the initial interface, which is represented by both the bilateral contact angle θ₁ and the middle contact angle θ₂. The increase of θ₁ and the decrease of θ₂ have opposite effects on the critical velocity. When θ₂ remains roughly unchanged at 90 degrees, if θ₁ is greater than 90 degrees (such as the parabolic interface), the critical shear velocity increases with the increase of θ₁, and the ellipse perturbation is its limiting case; If θ₁ is less than 90 degrees (such as the inverted parabolic and the inverted ellipse disturbances), the critical shear velocities are basically the same, which is less than that of the sinusoidal and sawtooth disturbances. The influence of inverted parabolic and inverted ellipse perturbations on the transition point of the RTKHI system is greater than that of other interfaces: (i) For the same amplitude, the smaller the contact angle θ₁, the later the transition point appears; (ii) For the same interface morphology, the disturbance amplitude increases, resulting in a shorter duration of the linear growth stage, so the transition point is greatly advanced.     

ICF流体力学不稳定性实验用柱状激波管的研制

在惯性约束聚变（ICF）实验中,点火靶丸表面（界面）的粗糙度和缺陷所产生的流体力学不稳定性是决定点火成功与否的关键因素之一,设计和研制流体力学不稳定性分解实验用靶是解决该问题的主要技术手段。结合国内外的研究现状和神光-Ⅱ激光装置的特点,设计并研制了一种新型柱状激波管。该靶型由三种介质组成,分别为调制聚苯乙烯（CH）圆片、柱状碳气凝胶（CRF）和CH微套管。调制CH圆片和柱状CRF通过微加工技术装配到CH微套管内,封装后形成柱状激波管。介绍了该靶型的设计原理和详细的制备工艺,并对相应的靶参数进行了测量。结果表明：柱状CRF气凝胶具有较好的成型性,长度、直径和密度分别为1000 m、730 m和250 mgcm-3;CH圆片的厚度和直径分别为15 m和730 m,表面调制图形的周期和峰谷差分别为100 m和4.3 m;实验得到的柱状激波管的轴向和径向最大装配误差分别为2 m和3 m。     

Electron proton instability in the CSNS ring

The electron proton (e-p) instability has been observed in many proton accelerators. It will induce transverse beam size blow-up, cause beam loss and restrict the machine performance. Much research work has been done on the causes, dynamics and cures of this instability. A simulation code is developed to study the e-p instability in the ring of the China Spallation Neutron Source (CSNS).     

单层钨丝阵Z箍缩内爆不稳定性特征实验研究

基于激光阴影像,研究了直径为8 mm的钨丝阵在滞止阶段的不稳定性发展。实验结果表明,在消融阶段早期,不稳定性扰动的主导波长为0.2 mm。在内爆后期和滞止阶段,不稳定主导波长迅速向长波长方向发展。在功率峰值前后,扭曲不稳定性特征开始出现,在滞止等离子体柱达到最小收缩直径后,辐射功率出现了明显下降。     

Electron proton instability in the CSNS ring

The electron proton (e-p) instability has been observed in many proton accelerators. It will induce transverse beam size blow-up, cause beam loss and restrict the machine performance. Much research work has been done on the causes, dynamics and cures of this instability. A simulation code is developed to study the e-p instability in the ring of the China Spallation Neutron Source (CSNS).     

Observation of beam instability in the SSRF storage ring

During the SSRF(Shanghai Synchrotron Radiation Facility) storage ring phase-I commissioning, some instabilities have been observed, and the broadband impedance has also been measured. The primary instabilities at present stage are vertical beam blow up and resistive wall instability.     

Rayleigh-Taylor instability of multi-fluid layers in cylindrical geometry

Rayleigh-Taylor instability of three fluid layers with two interfaces in cylindrical geometry is investigated analytically.The growth rates and the amplitudes of perturbation on the two interfaces are obtained. The feedback factor from outer to inner interface is larger than that from inner to outer interface under the same conditions. The growth rate on the initially unstable interface is larger than the corresponding result in planar geometry for low mode perturbation. The two interfaces are decoupled for a larger mode number perturbation. The dependencies of the amplitudes of perturbation on different initial conditions are analyzed. The negative feedback effect from initially stable interface to another unstable interface is observed. In the limit of infinity inner radius and finite shell thickness, the results in planar geometry are recovered.     

Observation of beam instability in the SSRF storage ring

During the SSRF(Shanghai Synchrotron Radiation Facility) storage ring phase-I commissioning, some instabilities have been observed, and the broadband impedance has also been measured. The primary instabilities at present stage are vertical beam blow up and resistive wall instability.     

柱几何Rayleigh-Taylor不稳定性的数值模拟

 给出了柱几何中流体力学方程组及其在数值模拟中采用的计算方法。对二维柱几何Rayleigh-Taylor不稳定性进行数值模拟，在线性阶段与线性理论符合得很好；不稳定性增长进入非线性区域的阈值依赖于界面的位置，并且明显不同于平面情况。     

Nonlinear saturation amplitude of cylindrical Rayleigh Taylor instability

The nonlinear saturation amplitude （NSA） of the fundamental mode in the classical Rayleigh-Taylor instability with a cylindrical geometry for an arbitrary Atwood number is analytically investigated by considering the nonlinear corrections up to the third order. The analytic results indicate that the effects of the initial radius of the interface （r0） and the Atwood number （A） play an important role in the NSA of the fundamental mode. The NSA of the fundamental mode first increases gently and then decreases quickly with increasing A. For a given A, the smaller the ro/λ（λ is the perturbation wavelength）, the larger the NSA of the fundamental mode. When ro/λ is large enough （r0 〉〉 λ）, the NSA of the fundamental mode is reduced to the prediction in the previous literatures within the framework of the third-order perturbation theory.     

Directly driven Rayleigh—Taylor instability of modulated CH targets

Directly driven ablative Rayleigh-Taylor (R-T) instability of modulated CH targets was studied using the face- on X-ray radiography on the Shen-Guang II device. We obtained temporal evolution images of the R-T instability perturbation. The R-T instability growth factor has been obtained by using the methods of fast Fourier transform and seeking the difference of light intensity between the peak and the valley of the targets. Through comparison with the the theoretical simulation, we found that the experimental data had a good agreement with the theoretical simulation results before 1.8 ns, and was lower than the theoretical simulation results after that.     

双锥对撞点火机制2020年冬季实验中的瑞利-泰勒不稳定性分析

直接驱动激光聚变通过整形后的纳秒脉冲激光辐照氘氚(DT)球壳靶,经球对称压缩加速后,在中心转滞获得高温等离子体热斑,实现聚变点火.在球壳靶受到压缩和加速过程中等离子体界面的流体力学不稳定性,特别是瑞利-泰勒不稳定性的增长有可能会对压缩壳层造成破坏,导致点火的失败.本文通过理论解析和数值模拟,对基于Zhang等提出的双锥对撞点火方案(2020 Philos.Trans.A Math.Phys.Eng.Sci.378 20200015)在2020年冬季实验条件下的流体力学不稳定性增长进行了分析.结果显示理论模型与一维数值模拟中对整体压缩和加速过程的描述基本一致,在当前的近等熵波形下金锥中的壳层靶实现了低熵压缩,同时瑞利-泰勒不稳定性增长导致的最危险时刻扰动振幅和壳层厚度比可以达到约0.25,壳层依然处于安全状态,但当初始壳层表面扰动均方根振幅大于22 nm时,则可能出现壳层的破裂.因此,未来实验中的靶设计与驱动激光脉冲波形设计中可以通过增加靶壳层厚度、提高预脉冲强度、减小靶表面的粗糙度和提高激光辐照的匀滑度等方式来抑制不稳定性增长.     

The numerical study of shock-induced hydrodynamic instability and mixing

Based on multi-fluid volume fraction and piecewise parabolic method (PPM), a multi-viscosity-fluid hydrodynamic code MVPPM (Multi-Viscosity-Fluid Piecewise Parabolic Method) is developed and applied to the problems of shock-induced hydrodynamic interfacial instability and mixing. Simulations of gas/liquid interface instability show that the influences of initial perturbations on the fluid mixing zone (FMZ) growth are significant, especially at the late stages, while grids have only a slight effect on the FMZ width, when the interface is impulsively accelerated by a shock wave passing through it. A numerical study of the hydrodynamic interfacial instability and mixing of gaseous flows impacted by re-shocks is presented. It reveals that the numerical results are in good agreement with the experimental results and the mixing growth rate strongly depends on initial conditions. Ultimately, the jelly layer experiment relevant to the instability impacted by exploding is simulated. The shape of jelly interface, position of front face of jelly layer, crest and trough of perturbation versus time are given; their simulated results are in good agreement with experimental results.