爆轰驱动铝飞层扰动增长的数值模拟（英文）

建立了研究炸药爆轰驱动条件下金属材料Rayleigh-Taylor不稳定性问题的实验技术和数值模拟方法。利用该实验技术和数值模拟方法研究了炸药爆轰驱动条件下,铝飞层界面Rayleigh-Taylor不稳定性增长规律,数值模拟显示界面扰动振幅以指数规律增长。数值模拟结果和实验定性相符,但是定量相比有较大差别,原因是高压高应变率加载条件下铝的强度增强,而数值模拟时所采用的SG本构模型在这样的加载条件下低估了铝的强度而导致对扰动增长致稳作用不足。然后在数值模拟中,通过改变材料的初始剪切模量和初始屈服强度,发现在一定范围内,初始剪切模量对材料动态屈服强度没有影响,而初始屈服强度增大可以明显提高材料的动态屈服强度,达到抑制扰动增长的目的,表明材料屈服强度主导界面扰动增长。     

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

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

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

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

炸药爆轰的连续介质本构模型和数值计算方法

假定炸药和爆轰产物处于局部热力学平衡状态, 即它们的压力和温度相同, 利用热力学基本关系建立炸药爆轰过程的连续介质本构模型的一般理论框架. 在此框架下, 炸药爆轰本构模型由一组常微分方程构成, 包括炸药和爆轰产物的状态方程、简单混合法则、化学反应速率方程和能量守恒方程, 易于由成熟的计算方法如梯形法等进行求解. 一组广义Maxwell型非线性固体本构形式的微分方程描述了压力和温度随时间的演化速率与应变率和化学反应速率的关系, 借助简单混合物理论, 其中的系数由炸药和爆轰产物的材料参数确定. 未反应的炸药和爆轰产物采用JWL状态方程, 化学反应率方程采用Lee-Tarver点火-燃烧二项式模型, 模拟PBX-9404炸药的一维冲击波起爆过程和爆轰波传播过程. 计算结果表明了本文给出的本构模型和相应计算方法的有效性. 关键词： 炸药爆轰 本构模型 化学反应率方程 数值模拟     

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

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

爆轰加载下弹塑性固体Richtmyer-Meshkov流动的扰动增长规律

研究了冲击波加载弹塑性材料扰动自由面的动力学演化过程,分析了高能炸药爆轰驱动时初始扰动与材料性质对扰动增长的影响.研究结果表明:初始扰动的振幅与波长之比越高,扰动越易增长,强度越高的材料扰动增长幅度越小;扰动增长被抑制时,尖钉的最大振幅与增长速度无量纲数之间存在线性近似关系,进一步理论分析表明尖钉的振幅增长因子与加载压力、初始扰动形态和材料强度有关,该理论关系作为扰动增长规律的线性近似在一定范围内适用于多种金属材料.     

烧蚀RT不稳定性X射线分幅诊断研究进展

在神光II装置上,利用高动态范围高性能X射线分幅相机开展了辐射驱动烧蚀RT不稳定性面背光实验研究.在神光II 8路2ns辐射源和第九路Mo背光条件下,利用二维时空照相获得了周期20 μupm、初始扰动1 μupm烧蚀RT样品清晰的增长过程,并通过掺Br比例1.1%样品观测到RT非线性增长的结果.实验为惯性约束聚变(inertial confinement fusion,ICF)RT不稳定性定量表征和数值模拟奠定了良好的基础.     

工字形太赫兹超材料吸波体的传感特性研究

利用超材料吸波体对材料参数的电磁响应, 可将其应用于传感. 本文设计了一种工字形单元结构的超材料吸波体, 基于频域算法对其在太赫兹频段的传感特性进行数值模拟, 研究了待测样品折射率、厚度及电介质隔层厚度对超材料吸波体传感器的频率灵敏度、振幅灵敏度及品质因数的影响. 研究结果表明:当待测样品厚度为40 μm时, 折射率频率灵敏度可达到153.17 GHz/RIU, 折射率振幅灵敏度可达到41.37%/RIU; 待测样品折射率一定时, 厚度频率灵敏度随其厚度的增大而线性减小; 随着待测样品厚度的增加, RFOM呈增大趋势, 但增大幅度在逐渐减小; TFOM随待测样品厚度的增加而减小.     

靶参数对收缩几何瑞利-泰勒不稳定性诊断的影响

面向背光照相技术是诊断惯性约束聚变(ICF)中瑞利-泰勒(RT)不稳定性的重要方法,讨论了利用该技术对球形靶丸扰动幅度测量时由于收缩几何效应带来的影响。通过计算,分析了一个简化模型下扰动振幅的实际值和测量值,讨论了偏移距离、靶丸外半径、扰动波长和扰动振幅等因素对实验结果相对误差的影响。计算表明,合理选择这些参数能使诊断的系统误差小于3%,而且可通过计算模型对测量结果进行修正。研究结果可为即将开展的神光III激光装置上的收缩几何烧蚀RT不稳定性实验参数设计和结果分析提供依据。     

爆轰作用下飞层的界面不稳定性

 计算了在平面一维爆轰驱动和滑移爆轰驱动下飞层的界面不稳定性问题。一维爆轰驱动的计算结果与实验符合较好；滑移爆轰驱动尚无实验可比，只是给出了一个算例，由于Helmholtz不稳定性，不仅导致振幅增长，还导致了位相移动。     

Ablation effects on weakly nonlinear Rayleigh-Taylor instability with a finite bandwidth

A weakly nonlinear but numerically tractable model (to third order in amplitude and including bandwidth effects) has been developed for the ablative Rayleigh-Taylor (RT) instability. Model results clearly show growth reduction from linear ablative RT values and even amplitude saturation in some realistic cases. For excitation of a band of wave numbers near the cutoff for growth, the behavior is dominated by the mode with the largest linear growth rate, and not by the mode with the largest initial amplitude. This type of model is likely to be important for the future assessment of the RT effects on specific target designs of the inertial confinement fusion.     

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.     

爆轰加载下高纯铜界面Rayleigh-Taylor不稳定性实验研究

金属界面不稳定性是内爆物理压缩过程中关注的重要问题,与传统流体界面不稳定性具有显著区别.由于相关理论和实验诊断技术的限制,目前该问题的研究还明显不足.为加深对金属界面不稳定性扰动增长行为的认识,本文建立了爆轰加载下高纯铜界面Rayleigh-Taylor不稳定性研究的实验诊断技术和数据处理方法,得到了扰动发展早期不同时刻界面扰动增长的X光图像.实验结果分析表明:在爆轰产物的无冲击加载条件下扰动波长基本保持不变,而初始扰动幅值越大,界面扰动增长的趋势就越明显;同时随着样品前界面扰动的不断发展,在样品的后自由面也出现了与前界面初始相位相反的扰动特征,即样品前界面扰动为波谷的位置所对应的后界面先运动而逐渐演变为波峰,而前界面扰动为波峰的位置所对应的后界面则演变为波谷;在5.26μs时刻,界面扰动幅值增长为初始值的700%左右,应变率达到了约105/s.结合数值模拟研究表明:在此情况下常用的Steinberg-Cochran-Guinan模型在一定程度上低估了高纯铜材料强度的强化特性,无法准确地描述强度对界面扰动增长的制稳作用,从而导致数值模拟结果要大于实验测量结果.     

Instability of a spherical drop in the field of an electrical dipole

Analytical calculations show that, as a field in which an initially spherical charged conducting incompressible drop is placed becomes more and more nonuniform, coupling between the drop’s oscillation modes grows and the threshold of stability against the electrical field pressure declines. When an electrostatic parameter characterizing the electric field pressure exceeds a value that is critical for a certain mode to be unstable, the amplitude of this mode exponentially grows in an aperiodic manner and the amplitudes of modes coupled with this mode build up in an oscillatory manner, each mode having its own instability growth rate. In all cases, there exists a threshold value of the dimensionless electric parameter above which all oscillation modes are unstable.     

Measurements of Rayleigh–Taylor instability growth of laser-shocked iron–silicon alloy

The Rayleigh–Taylor (RT) instability of liquid iron alloys is important for understanding the core formation mechanism in the Earth. Here we first report the measurement of RT instability growth for a liquid iron–silicon (Fe–Si) alloy, which is one of the major candidate for the material of the Earth’s core, using a high power laser. We optimized the measurement setup and analytical technique to observe the growth of perturbation on an Fe–Si sample surface. The growth of perturbation amplitude on the Fe–Si alloy under high pressure and temperature was successfully observed using in situ X-ray radiography. The growth rate of the RT instability for the Fe–Si alloy on about 1000?GPa was estimated to be 0.3 ns^?1.     

金属镀层两层膜系长周期光纤光栅谐振特性研究

采用严格的耦合模理论,建立了镀会属膜和敏感膜的两层膜系长周期光纤光栅复特征方程,用微扰法对复数超越特征方程进行求解.结果和文献[5]给出的数据相符.理论模拟数据表明.金属的消光系数对功率密度和透射谱影响很大.低阶偶次模式在纤芯内的功率密度所占比重远大于奇次模式,其衰减峰幅度也比奇次模高,因此,金属镀层长周期光纤光栅实际应用时,需要选择偶次模作为研究对象.进一步理论分析了金属和敏感膜层参数对长周期光纤光栅谐振特性的影响,研究发现,金属膜厚、敏感膜厚和折射率对谐振波长有较大的影响,谐振波长随金属膜厚、敏感膜厚和折射率的增大向短波方向漂移.在某些区域谐振波长存在跳变现象,长周期光纤光栅实际应用时应避开此区域.     

Oscillating two-stream instability of laser wakefield-driven plasma wave

The laser wakefield-driven plasma wave in a low-density plasma is seen to be susceptible to the oscillating two-stream instability (OTSI). The plasma wave couples to two short wavelength plasma wave sidebands. The pump plasma wave and sidebands exert a ponderomotive force on the electrons driving a low-frequency quasimode. The electron density perturbation associated with this mode couples with the pump-driven electron oscillatory velocity to produce nonlinear currents driving the sidebands. At large pump amplitude, the instability grows faster than the ion plasma frequency and ions do not play a significant role. The growth rate of the quasimode, at large pump amplitude scales faster than linear. The growth rate is maximum for an optimum wave number of the quasimode and also increases with pump amplitude. Nonlocal effects, however reduce the growth rate by about half.     

Behavior of rigid-porous layers at high levels of continuous acoustic excitation: theory and experiment

A model for the propagation of high amplitude continuous sound through hard-backed rigid-porous layers has been developed which allows for Forchheimer's correction to Darcy's law. The nonlinearity associated with this is shown to be particularly important in the range of frequencies around layer resonance. The model is based on the introduction of particle velocity dependent flow resistivity into the equivalent fluid model expression for complex tortuosity. Thermal effects are accounted for by means of a linear complex compressibility function. The model has been used to derive analytical expressions for surface impedance and reflection coefficient as a function of incident pressure amplitude. Depending on the material parameters, sample thickness, and frequency range the model predicts either growth or decrease of reflection coefficient with sound amplitude. Good agreement between model predictions and data for three rigid-porous materials is demonstrated.     

金属纳米结构对光谱响应及折射率灵敏度的影响

推导了任意形状金属纳米结构对复杂环境折射率的局域表面等离子体共振(LSPR)光谱响应的解析解,给出了体折射率和局部折射率灵敏度公式,并进行了详细的理论分析和数值验证。结果表明,金属纳米结构的光谱响应与被分析物厚度呈指数增长关系,并与被分析物的折射率、敏感层厚度、金属纳米结构的材料和形状等因素有关。被分析物折射率越高,光谱响应越大。敏感层越厚,局部折射率灵敏度越差。同时,尖锐的金属纳米结构形状可以获得更大的体折射率灵敏度。对金属纳米结构光谱响应及折射率灵敏度的研究对制作高灵敏度LSPR传感器具有重要的指导意义。