神光Ⅱ装置上辐射驱动瑞利-泰勒不稳定性实验

在神光Ⅱ装置上开展了辐射驱动RT不稳定性的一系列实验,获得了不同初始扰动幅度、不同扰动波长、不同材料样品等条件下辐射烧蚀RT不稳定性增长的高时空分辨背光图像,特别是在大初始扰动幅度样品实验中获得了扰动增长的清晰图像,观察到了扰动增长从线性区到非线性区的过渡过程,二次和三次谐波的产生和发展清楚可见。充实了数值模拟程序考核的实验数据库,对间接驱动ICF点火靶设计和研究具有重要作用。     

神光Ⅱ装置上辐射驱动瑞利-泰勒不稳定性实验北大核心CSCD

在神光Ⅱ装置上开展了辐射驱动RT不稳定性的一系列实验,获得了不同初始扰动幅度、不同扰动波长、不同材料样品等条件下辐射烧蚀RT不稳定性增长的高时空分辨背光图像,特别是在大初始扰动幅度样品实验中获得了扰动增长的清晰图像,观察到了扰动增长从线性区到非线性区的过渡过程,二次和三次谐波的产生和发展清楚可见。充实了数值模拟程序考核的实验数据库,对间接驱动ICF点火靶设计和研究具有重要作用。     

远场来流条件下过冷熔体球晶生长的稳定性

建立了远场来流条件下过冷熔体球晶生长的温度场和浓度场稳态模型,分析了对流对球晶周围温度场和浓度场的影响,并以Trivedi的纯扩散球晶稳定性判据为基础,推导出远场来流条件下过冷熔体球晶生长的临界稳定性判据. 研究表明:远场来流条件下,迎流面的扰动振幅增加速率明显大于背流面的扰动振幅增加速率. 振幅增加速率最大值对应的扰动阶次从迎流面到背流面逐渐减小,随着球晶半径增加而增大. 对流使迎流面的稳定性降低,背流面的稳定性增大. 随着流速的增加,球晶的临界稳定半径减小. 关键词： 球晶 远场来流 界面形态稳定性 Trivedi判据     

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

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

非理想流体中Rayleigh-Taylor和Richtmyer-Meshkov不稳定性气泡速度研究

在随气泡顶端运动的坐标系中, 通过将理想流体模型推广到非理想流体的情况, 研究了流体黏性和表面张力对Rayleigh-Taylor (RT)和Richtmyer-Meshkov (RM)不稳定性气泡速度的影响. 首先得到了RT和RM不稳定性气泡运动的控制方程 (自洽的微分方程组); 其次给出了二维平面坐标和三维柱坐标中气泡速度的数值解和渐近解, 并定量分析了流体黏性和表面张力对RT和RM气泡速度和振幅的影响. 结果表明: 从线性阶段到非线性阶段的全过程中, 非理想流体中的气泡速度和振幅小于理想流体中的气泡速度和振幅. 也就是说, 流体黏性和表面张力对RT和RM不稳定性的发展都具有致稳作用. 关键词： Rayleigh-Taylor不稳定性 Richtmyer-Meshkov不稳定性 气泡速度 非理想流体     

偏心对汇聚激波诱导的RM不稳定性影响的数值研究

数值研究汇聚激波与四种形状（圆形、小振幅单模、大振幅单模和正方形）的二维气柱界面相互作用,激波汇聚中心与界面同心和不同心（即偏心）时Richtmyer-Meshkov（RM）不稳定性的发展规律,重点考察界面中心的压力及混合区面积在两种情况下随时间的变化.数值方法使用VAS2D程序,该方法采用有限体积法结合网格自适应技术,能够达到时间和空间的二阶精度.结果表明,偏心情况下RM不稳定性是其在同心情况下的扰动和偏心小扰动叠加的结果.在本文采用的偏心程度下（20%）,偏心对于圆形无扰动界面发展的影响主要表现在后期界面出现微小扰动结构;而对于单模和正方形这种原本有扰动的界面,偏心使扰动结构呈现不对称及扭曲,同时也影响了界面中心压力和混合区面积,因而加剧了不稳定性的发展.     

涡旋中双臂形成的物理机制及其演变

利用柱坐标下的流体涡度方程,研究了涡旋中双臂的形成机制问题.结果表明：涡旋中切向波数为2（m=2）的扰动在基本流场中可出现不稳定快速增长,但其增长率受时间限制;扰动可演变为螺旋双臂并向外扩展;扰动发展的物理机制是扰动的速度场对基本涡度场强平流效应所致. 关键词： 旋涡星系 不稳定波动 螺旋云带 台风     

烧蚀RT不稳定性非线性过程模拟

重点研究了不同预热对烧蚀RT弱非线性模耦合的影响。首先不同预热情况烧蚀RT不稳定性的线性增长率曲线明显不同。不同预热情况的模拟结果都表明：烧蚀RT的模耦合系数是k／kc的函数，二次和三次谐波的产生系数C（k）和D（k）遵从定标关系C（k）：c1（1-c2k/kc），D（k）=d1-d2k/kc＋d3（k/kc）^2。弱预热和中等预热情况，二次谐波的产生系数约只有经典RT的50％，三次谐波的产生系数约只有经典RT的20％。两体模耦合起主要作用，n次高次谐波可近似看作是由（”一1）次两体模耦合产生的。由于烧蚀RT不稳定性两体模耦合系数小于1，所以高次谐波的产生系数很快衰减。由于非线性作用变弱，单模扰动的非线性饱和阈值明显增大，意味着基模有较长时间的线性增长，另外，向长波长方向的模耦合系数明显大于经典尺丁不稳定性数值，因此，烧蚀RT非线性模耦合容易形成长的尖顶，从而对点火构成严重威胁。     

二维不可压流体瑞利-泰勒不稳定性的非线性阈值公式

分析了二维情况下平面几何、柱几何和球几何中瑞利-泰勒不稳定性发生非线性偏离的阈值问题，给出了三种几何中密度扰动振幅的定义，以及与界面扰动振幅的关系.由此得到了三种几何中密度扰动的非线性阈值公式，用高精度流体程序对三种几何中的不稳定性进行了数值模拟，验证了得到的非线性阈值公式. 关键词： 瑞利-泰勒不稳定性 非线性阈值 密度振幅     

离散效应对高强度激光三次谐波转换的影响

使用分步傅里叶变换和四阶龙格库塔法(R-K) ，对高强度激光以Ⅰ/Ⅱ类角度失谐方式，在KDP晶体中的谐波转换进行了研究，详细讨论了离散效应对三次谐波转换的影响。结果表明，离散效应不但降低了三次谐波转换效率，而且使谐波光束质量显著降低；初始入射基频光束腰半径较小时，离散效应是二次谐波转换和三次谐波转换效率降低的主要因素，失谐角对三次谐波转换的影响较小；随着束腰半径的增加，离散效应的影响逐渐减小，失谐角对三次谐波转换的影响逐渐增加。     

表面张力对Rayleigh-Taylor不稳定性谐波的影响

Using the method of the parameter expansion up to the third order, explicitly investigates surface tension effect on harmonics at weakly nonlinear stage in Rayleigh-Taylor instability (RTI) for arbitrary Atwood numbers and compares the results with those of classical RTI within the framework of the third-order weakly nonlinear theory. It is found that surface tension strongly reduces the linear growth rate of time, resulting in mild growth of the amplitude of the fundamental mode, and changes amplitudes of the second and third harmonics, as is expressed as a tension factor coupling in amplitudes of the harmonics. On the one hand, surface tension can either decrease or increase the space amplitude; on the other hand, surface tension can also change their phases for some conditions which are explicitly determined.     

Nonlinear Saturation Amplitude in Classical Planar Richtmyer–Meshkov Instability

The classical planar Richtmyer–Meshkov instability(RMI) at a fluid interface supported by a constant pressure is investigated by a formal perturbation expansion up to the third order,and then according to definition of nonlinear saturation amplitude(NSA) in Rayleigh–Taylor instability(RTI),the NSA in planar RMI is obtained explicitly.It is found that the NSA in planar RMI is affected by the initial perturbation wavelength and the initial amplitude of the interface,while the effect of the initial amplitude of the interface on the NSA is less than that of the initial perturbation wavelength.Without marginal influence of the initial amplitude,the NSA increases linearly with wavelength.The NSA normalized by the wavelength in planar RMI is about 0.11,larger than that corresponding to RTI.     

Simulation of the Weakly Nonlinear Rayleigh–Taylor Instability in Spherical Geometry

The Rayleigh-Taylor instability at the weakly nonlinear(WN) stage in spherical geometry is studied by numerical simulation.The mode coupling processes are revealed.The results are consistent with the WN model based on parameter expansion,while higher order effects are found to be non-negligible.For Legendre mode perturbation Pn(cos B),the nonlinear saturation amplitude(NS A) of the fundamental mode decreases with the mode number n.When n is large,the spherical NSA is lower than the corresponding planar one.However,for large n,the planar NSA can be recovered by applying Fourier transformation to the bubble/spike near the equator and calculating the NSA of the converted trigonometric harmonic.     

Theoretical and simulation research of hydrodynamic instabilities in inertial-confinement fusion implosions

Inertial fusion energy(IFE)has been considered a promising,nearly inexhaustible source of sustainable carbon-free power for the world's energy future.It has long been recognized that the control of hydrodynamic instabilities is of critical importance for ignition and high-gain in the inertial-confinement fusion(ICF)hot-spot ignition scheme.In this mini-review,we summarize the progress of theoretical and simulation research of hydrodynamic instabilities in the ICF central hot-spot implosion in our group over the past decade.In order to obtain sufficient understanding of the growth of hydrodynamic instabilities in ICF,we first decompose the problem into different stages according to the implosion physics processes.The decomposed essential physics processes that are associated with ICF implosions,such as Rayleigh-Taylor instability(RTI),Richtmyer-Meshkov instability(RMI),Kelvin-Helmholtz instability(KHI),convergent geometry effects,as well as perturbation feed-through are reviewed.Analytical models in planar,cylindrical,and spherical geometries have been established to study different physical aspects,including density-gradient,interface-coupling,geometry,and convergent effects.The influence of ablation in the presence of preheating on the RTI has been extensively studied by numerical simulations.The KHI considering the ablation effect has been discussed in detail for the first time.A series of single-mode ablative RTI experiments has been performed on the Shenguang-Ⅱlaser facility.The theoretical and simulation research provides us the physical insights of linear and weakly nonlinear growths,and nonlinear evolutions of the hydrodynamic instabilities in ICF implosions,which has directly supported the research of ICF ignition target design.The ICF hot-spot ignition implosion design that uses several controlling features,based on our current understanding of hydrodynamic instabilities,to address shell implosion stability,has been briefly described,several of which are novel.     

掺杂对CH样品Rayleigh-Taylor不稳定性增长的影响

在神光II激光装置上进行了辐射驱动不同掺杂样品的单模Rayleigh-Taylor(RT)不稳定性实验.结果显示:与纯碳氢(CH)样品相比,掺Br的CH样品的扰动更早、更快地进入非线性区,产生二次谐波,并且掺Br比例越高,CH样品扰动进入非线性区的时间越早,相同时刻扰动的二次谐波的幅度越高.这是因为密度梯度效应抑制了二次谐波的产生,掺Br比例越高,密度梯度标长越小;同时密度梯度效应还抑制三次谐波对基模增长的负反馈,造成基模具有更大的线性增长,导致线性饱和幅值大于经典值0.1λ.     

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.     

Weak Nonlinearity of Ablative Rayleigh--Taylor Instability

The weakly nonlinear regime of single mode ablative Rayleigh-Taylor instability is studied, with consideration of preheat effect and the width of the ablation front. The Rayleigh-Taylor linear growth rate agrees well with the direct numerical simulation. For the density perturbation, the amplitude distribution of the fundamental mode has one peak value whereas those of the second and third harmonics have two and three peak values, respectively. Harmonics generation versus wave number is also given and it is close to the result of direct numerical simulation.     

Linear Growth of Rayleigh-Taylor Instability of Two Finite-Thickness Fluid Layers

The linear growth of Rayleigh-Taylor instability(RTI) of two superimposed finite-thickness fluids in a gravitational field is investigated analytically. Coupling evolution equations for perturbation on the upper, middle and lower interfaces of the two stratified fluids are derived. The growth rate of the RTI and the evolution of the amplitudes of perturbation on the three interfaces are obtained by solving the coupling equations. It is found that the finite-thickness fluids reduce the growth rate of perturbation on the middle interface. However, the finite-thickness effect plays an important role in perturbation growth even for the thin layers which will cause more severe RTI growth. Finally, the dependence of the interface position under different initial conditions are discussed in some detail.     

Effect of initial phase on the Rayleigh—Taylor instability of a finite-thickness fluid shell

Rayleigh—Taylor instability (RTI) of finite-thickness shell plays an important role in deep understanding the characteristics of shell deformation and material mixing. The RTI of a finite-thickness fluid layer is studied analytically considering an arbitrary perturbation phase difference on the two interfaces of the shell. The third-order weakly nonlinear (WN) solutions for RTI are derived. It is found the main feature (bubble-spike structure) of the interface is not affected by phase difference. However, the positions of bubble and spike are sensitive to the initial phase difference, especially for a thin shell (kd<1), which will be detrimental to the integrity of the shell. Furthermore, the larger phase difference results in much more serious RTI growth, significant shell deformation can be obtained in the WN stage for perturbations with large phase difference. Therefore, it should be considered in applications where the interface coupling and perturbation phase effects are important, such as inertial confinement fusion.