Z-pinch内爆等离子体时空分辨诊断技术

 为了深入认识和理解Z-pinch物理机制并对物理实验进行优化设计，需要研究内爆等离子体及其辐射的时空分布特性诊断技术，建立诊断装置与设备，通过实验研究等离子体内爆时间过程及其空间分布特性，利用所获得的实验诊断数据校核数值模拟程序。简要介绍近几年开展的X光空间积分功率谱测量、X光时间积分针孔照相、X光时分幅针孔照相、X光时间分辨1维空间分布、紫外探针阴影照相等时空分辨诊断技术研究工作，并给出在1~3 MA电流驱动装置上进行的喷气和丝阵内爆实验中获得的部分时空分辨诊断结果，包括辐射功率、拉链速度、内爆速度、收缩比等。     

斜界面Rayleigh-Taylor不稳定性混合实验研究

 实验研究了不相溶流体斜界面Rayleigh-Taylor（R-T）不稳定性湍流混合区的混合不对称性特征。利用高压气体加速装有不同液体的箱体,加速度方向由轻液体指向重液体,此时界面是R-T不稳定性的。利用阴影测试技术,研究了初始倾角9°的ZnCl₂溶液/正己烷斜界面的演化规律,得出混合区内混合不对称、斜界面倾角渐增的规律。     

液-液斜界面R-T失稳特征的实验研究

实验研究了不同加速度下不相溶流体斜界面Rayleigh-Taylor(R-T)不稳定性诱导湍流混合区的发展特征。利用高速阴影测试技术研究了硅油/碘化钾溶液界面R-T不稳定性诱导界面的演化规律,定量分析了混合区发展宽度、界面倾角与时间的关系。通过实验结果分析得出:水平界面与斜界面混合区宽度演化的差别主要体现在后期,前、中期两者演化规律基本一致;界面倾角的演化呈现出两种不同的趋势,随时间先呈抛物线增长,随后线性增长,演化后期流体界面发生翻转;由R-T不稳定性诱导的界面混合与K-H不稳定性诱导的界面混合呈相互竞争的关系。     

神光Ⅱ升级装置激光加速质子照相初步实验研究

在神光Ⅱ升级装置上开展了首轮激光加速质子对间接驱动快点火靶内爆过程的照相实验研究。通过激光与靶参数的优化，获得了能量高于18 MeV的质子束。通过静态客体的照相，获得了优于20 m的高空间分辨网格图像，为开展时间分辨的啁啾质子照相奠定了基础。开展了质子动态照相实验，获得了内爆压缩晚期的质子照相图像。实验发现内爆区域质子照相图像存在大量排空现象。内爆压缩区域不足以阻挡如此大范围质子束，证明了其中存在电磁场使得质子向外排开。动态照相的质子能量较低，分析是ns激光打靶过程产生的X射线及等离子体对质子加速存在影响。后续实验中需要进一步优化靶的屏蔽设计。     

含尘气体Rayleigh-Taylor不稳定性诱导混合的多相浮阻力模型研究

文章首先采用单相浮阻力模型对不同加速度下Rayleigh-Taylor不稳定性诱发的物质渗透边界的演化过程进行了计算, 揭示了该混合在常加速度和变加速度情况下不同的发展规律, 并通过与实验结果的比较分析, 验证了该模型的适用性.在此基础上, 发展了多相浮阻力模型, 采用该模型对常加速度情况下含尘气体中的Rayleigh-Taylor不稳定性诱导混合进行了研究, 发现混合区宽度随着颗粒数密度和颗粒尺寸的增大而减小, 揭示了气体中所含杂质抑制混合发展的规律.      

示踪材料在内爆X光诊断过程中的作用

 激光聚变内爆实验中,在燃料中掺杂少量比例的中高Z材料,用X光光谱和X光成像测量掺杂元素的发射信息,诊断燃料的温度、密度和压缩形状。用辐射流体力学数值计算和X光成像后处理程序综合分析方法,给出了内爆靶丸优化设计,并讨论示踪材料在X光诊断中的作用。结果表明:在靶丸燃料D2中掺原子分数约1.0％的氩,内爆压缩中子产额下降约15％。由于氩线发射使整个燃料区X光发射强度提高约50倍,X光成像区域增大约30％,有利于实验诊断测量燃料芯部。为了测量燃料区的边界,在CH内壳层涂厚度0.05 μm的硫,分析表明硫Ly-α单能成像大小与流体力学计算的燃料区大小一致,可用于诊断燃料最终压缩界面。数值分析结果得到了神光Ⅱ间接驱动内爆物理相关实验的验证。     

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

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

内表面掺硅聚苯乙烯空心微球初步研究

靶丸内表面掺杂诊断元素可以为内爆压缩界面的研究提供必要的手段。以二甲基二乙氧基硅烷为原料，利用乳液微封装技术和界面聚合技术可以实现单一内表面掺硅的空心聚苯乙烯微球的制备。X光照相和X能谱分析表明:内表面掺硅量比外表面高出2~3个数量级。均匀内表面掺硅层的厚度小于0.3μm，掺杂层厚度越小，均匀性越好。     

内表面掺硅聚苯乙烯空心微球初步研究

 靶丸内表面掺杂诊断元素可以为内爆压缩界面的研究提供必要的手段。以二甲基二乙氧基硅烷为原料，利用乳液微封装技术和界面聚合技术可以实现单一内表面掺硅的空心聚苯乙烯微球的制备。X光照相和X能谱分析表明：内表面掺硅量比外表面高出2~3个数量级。均匀内表面掺硅层的厚度小于0.3μm，掺杂层厚度越小，均匀性越好。     

不同界面组分分布对Richtmyer-Meshkov不稳定性的影响

基于二维非定常Euler方程,对平面激波与不同界面组分分布下氦气气柱作用过程所引起的Richtmyer-Meshkov不稳定性现象进行了数值模拟,探讨了激波冲击轻质气柱后气柱界面形态的演变及流场波系结构,定量分析了气柱特征尺度(气柱长度、高度和中轴宽度)和气柱体积压缩率随时间变化.此外,结合流场压强、速度、环量和气体混合率,多角度分析了激波驱动界面气体混合的流动机制,获得了不同界面组分分布对界面不稳定性的影响.结果表明,随着气柱界面从完全扩散界面向间断界面的过渡,界面两侧的声反射系数随之增大,使入射激波与气柱界面的作用由常规透射转变为非常规透射,反射激波逐渐加强,透射激波逐渐减弱,使得Richtmyer-Meshkov不稳定性随之增强;同时,界面两侧阿特伍德数的增大,加强了Rayleigh-Taylor不稳定性和Kelvin-Helmholtz不稳定性的发展.此外,界面不稳定性的加强使得流场环量增大,导致气体混合率的增长速率随之升高.     

激光加载下金属锡材料微喷颗粒与低密度泡沫混合实验研究

金属材料的微喷是冲击加载下金属表面发生的一种动态破碎现象,微喷研究在很多领域都具有重要意义,包括惯性约束聚变(ICF)和烟火制造等.由于激光实验特有的优势,近几年国内外开展了很多利用强激光驱动冲击加载研究材料微喷过程的实验.利用泡沫材料对微喷颗粒进行静态软回收虽然可以获得颗粒的形态分布、颗粒尺寸及颗粒质量等定量结果,但并不能反演微喷颗粒从进入泡沫到停滞过程中的动态混合过程.为此,在神光Ⅱ升级装置上利用皮秒脉冲激光照射金丝产生高能X射线,实现了对锡微喷颗粒与低密度泡沫混合过程的高时间分辨和高空间分辨背光照相.背光图像面密度结果证实微喷颗粒在泡沫中并没有发生二次破碎.静态回收结果表明,在锡材料与泡沫紧贴放置的情况下,微喷颗粒在泡沫中的穿透深度随着加载压强升高呈现先增大后减小的规律,与非紧贴放置的实验结果有明显的差别.     

Analysis of the flame thickness of turbulent flamelets in the thin reaction zones regime

The thickness of the instantaneous flamelets in a turbulent flame brush on a weak-swirl burner burning in the thin reaction zones regime has been analysed experimentally, theoretically, and numerically. The experimental flame thickness has been measured correlating two simultaneous Rayleigh images and one OH-image from two closely spaced cross sections in the flame. It appears that the low temperature edge of the flame is thickened by turbulent eddies but that these structures cannot penetrate far enough into the flame front to distort the inner layer for the moderate Karlovitz numbers used. The flame front based on the temperature gradient at the inner layer becomes thinner for lean flames and thicker for rich methane–air flames. This has been explained theoretically and numerically by studying the influence of flame stretch and preferential diffusion on the flame thickness. It appears that the flame front thickness at the inner layer (and mass burning rate) is not influenced by turbulent mixing processes, and it seems that eddies of the size of the inner layer have to be used to change this picture. Experiments closer to the boundary of the broken reaction zones regime have to confirm this in the future.     

闪光照相中FXRMC和MCNP4B的散射比较研究

 散射问题是高能辐射成像研究中的一个重要问题，采用蒙特卡罗模拟来确定散射对提取客体信息的影响是一种重要的研究手段。简单介绍了FXRMC和MCNP4B程序的特点及其记录方式；在确保相同输入参数的条件下，针对不同的照相模型进行了对比计算。结果表明两个程序计算的散射照射量相对差别小于5%，说明这两个程序具有较高的符合程度。通过与实验结果的比较发现，这两个程序模拟的散射分布与实验结果基本一致，均可用于高能闪光照相的模拟研究。还给出了在散射检验方面的一些建议。     

滚边法测量高能X射线源焦斑大小

应用于高能闪光X光照相技术的X射线源焦斑大小是闪光照相装置的关键参数, 直接影响成像的分辨能力。由于高能X射线的强穿透性和强辐射环境, 给焦斑测量带来一定困难。介绍了一种间接测量方法, 采用滚边装置(rollbar)成像得到X射线源的边扩展函数, 微分后得到光源的线扩展函数并计算调制传递函数(MTF), 而后从MTF为0.5所对应的空间频率之值确定出光源的光斑大小。给出了神龙二号加速器电子束聚焦调试实验中得到的X射线焦斑测量结果, 分析影响测量结果的因素并提出了解决方法。     

针对高压物质密度诊断的激光X射线照相优化设计

激光驱动冲击波能提供高压的物质状态, 是状态方程研究的重要工具. 超短超强激光与固体靶相互作用产生的X射线源, 具有短脉冲、微焦点、高产额、能量可调的特点, 是高压物质密度测量的首选背光源. 本文基于蒙特卡罗程序Geant4建立了X射线照相模型, 客体密度分布由流体力学程序Multi-1D模拟激光冲击加载高压物质获得. 在一维长方体形密度客体情况下, 定义了均方根、峰值偏差与上升沿斜率比三种指标, 对照相图像求解的密度结果进行评价, 开展了照相结果信噪比、分辨率与对比度等参数优化. 并开展了一维圆柱形密度客体的照相模拟, 建立了基于Radon逆变换法的Abel反演算法. 反演结果与模拟设计密度分布符合良好, 要求X射线源半径在5 μm以内; 反演结果与模拟设计密度分布较为一致, 要求X射线源半径在15 μm以内.     

Observation of laser light scattering by a turbulent gas without seeding

An experiment of laser light scattering by a turbulent highly subsonic free jet of air is presented. The originality of the method lies in the fact that no particle seeding is needed. The experimental scattering cross section, which can be interpreted from the three-dimensional spectral density of temperature fluctuations, is compared to the theoretical one obtained by turbulent mixing theories. The scattered field time correlations measured in the heterodyne mode are interpreted by considering the Doppler effect on moving turbulent eddies and compared with theory. The intensity time correlations measured in the homodyne mode are interpreted in terms of local properties of the turbulence.     

Aerooptic effects in a turbulent flow and their simulation

Large-scale eddies and related aerooptic effects in a turbulent boundary layer, free mixing layer, and a jet flowing out of a round nozzle into a submerged space are simulated. The results obtained are used to analyze the distortions produced in the phase function of a coherent light beam by the turbulent fluctuations of the parameters of the medium. The results of numerical calculation are compared with the data of natural experiment and data obtained by solving Reynolds-averaged Navier-Stokes equations.     

Use of Lagrangian statistics for the analysis of the scale separation hypothesis in turbulent channel flow

Turbulence models often involve Reynolds averaging, with a closure providing the Reynolds stress tensor as function of mean velocity gradients, through a turbulence constitutive equation. The main limitation of this linear closure is that it rests on an analogy with kinetic theory. For this analogy to be valid there has to be a scale separation between the mean velocity variations and the turbulent Lagrangian free path whose mean value is the turbulent mixing length. The aim of this work is to better understand this hypothesis from a microscopic point of view. Therefore, fluid elements are tracked in a turbulent channel flow. The flow is resolved by direct numerical simulation (DNS). Statistics on particle trajectories ending on a certain distance y₀ from the wall are computed, leading to estimations of the turbulent mixing length scale and the Knudsen number. Comparing the computed values to the Knudsen number in the case of scale separation, we may know in which region of the flow and to what extent the turbulence constitutive equation is not verified. Finally, a new non-local formulation for predicting the Reynolds stress is proposed.     

Characterisation of mixing rate due to high power ultrasound

For meaningful assessment of results from laboratory and pilot plant trials, it is often necessary to know the mixing characteristics within the ultrasonic reactors. Previous workers have used conductivity measurements in an attempt to characterize the residence time distribution in ultrasonic reactors, but these results do not provide direct data on the mixing within the high power region adjacent to the ultrasonic probe tip. We overcome this difficulty through direct visualization of the mixing process within the high energy region close to the tip of the ultrasonic probe. Our analysis proceeds by determining an approximate turbulent diffusivity in a batch reactor arrangement for different values of ultrasonic energy input. For input electrical power levels between 70 and 120 W and a processing volume of 30 ml, the effective turbulent diffusivity varied from about 0.2 x 10(-3) to 0.7 x 10(-3)m(2)/s. We demonstrate that such results can be coupled to a suitable dispersion model to estimate the actual residence time distribution in flow-though arrangements when the through-put adds little additional mixing energy. Therefore, coupling the effective turbulent diffusivity identified in a batch reactor with a suitable dispersion model for the reactor offers an alternative approach to the deduction of RTD when determining the actual RTD in the high intensity zone of steady flow sonochemical reactors is problematic.