Numerical Simulation of Plasma Flows in the Interior Cavities of Laser Fusion Targets

This paper deals with a numerical simulation and theoretical investigation of the behavior of plasma flows produced when injecting a laser beam into a target for inertial thermonuclear fusion and absorbing laser radiation in the target internal cavity. Results of two-dimensional numerical calculation are presented for laser beam interaction with the sub critical plasma contained in the gap between two coaxial shells of a laser greenhouse target. The calculations confirm the possibility of initiating a supersonic wave of electron heat conduction in the absorber of the laser greenhouse target.     

Numerical Modelling of Impurity Retention in the ITER Tokamak Scrape-off Layer Plasma

The hydrogen plasma and carbon ions flows generated at the target plates in the ITER tokamak scrape-off layer are numerically investigated. A 2D model of hydrogen plasma and impurity ions is presented. The model is based on the electron-ion Braginski fluid equations [1] for hydrogen plasma and rate equations for impurities. Arbitrary level of impurity ions concentration is assumed. Recycling of hydrogen, sputtering and self sputtering of carbon atoms at the target plates are taken into account. Equations of the model are solved in the slab geometry using 2D-multifluid numerical code EPIT [2]. Problems of impurity ions retention and radiation in the divertor volume are analyzed. Results of calculations for ITER tokamak boundary plasma are presented, showing that poor retention is likely at high impurity concentration in the divertor volume. The radiation power can be a significant part of ingoing energy.     

直流电弧等离子体炬的特性研究

介绍了局部热力学平衡下描述直流电弧等离子体炬的磁流体力学模型，采用数值模拟方法研 究了等离子体炬内等离子体的传热和流动特性、湍流对等离子体炬内等离子体特性的影响， 以及等离子体炬运行参数对等离子体特性的影响.将数值模拟结果与已有的实验结果进行了 比较. 关键词： 等离子体炬 磁流体力学 数值模拟     

Modelling of the Effect of the Sheared Poloidal Flow on the Electrostatic Turbulence on the CASTOR Tokamak

Electrostatic drift turbulence of the edge plasma in the CASTOR tokamak is studied numerically by using the Hasegawa-Wakatani equations. The fluctuations of plasma density and potential as well as the corresponding fluctuation-induced particle flux are calculated for regimes with various plasma poloidal flows. Results of numerical simulations are in a reasonable agreement with experimental results.     

Change in the relief of a target surface treated by compression plasma flows

The change in the surface relief of a steel-3 target (GOST 380) treated by compression plasma flows is experimentally studied. The energy density absorbed by the target varies in the range of 10–35 J/cm² and the pulse duration is 100 μs. It is shown experimentally and numerically that the development of KelvinHelmholtz instability strongly affects the formation of the target surface treated with compression plasma flows: a large-scale wave-like relief with characteristic sizes of 200 × 1000 μm is formed on the target surface and, as a result, the roughness of the surface increases. However, the microrelief at the scale of individual elements is smoothed to a maximum roughness of about 0.5 μm.     

Theoretical analysis and numerical simulation of the impulse delivering from laser-produced plasma to solid target

A plasma is produced in air by using a high-intensity Q-switch Nd：YAG pulsed laser to irradiate a solid target, and the impulses delivering from the plasma to the target are measured at different laser power densities. Analysing the formation process of laser plasma and the laser supported detonation wave （LSDW） and using fluid mechanics theory and Pirri＇s methods, an approximately theoretical solution of the impulse delivering from the plasma to the target under our experimental condition is found. Furthermore, according to the formation time of plasma and the variation of pressure in plasma in a non-equilibrium state, a physical model of the interaction between the pulse laser and the solid target is developed. The plasma evolutions with time during and after the laser pulse irradiating the target are simulated numerically by using a three-dimensional difference scheme. And the numerical solutions of the impulse delivering from the plasma to the target are obtained. A comparison among the theoretical, numerical and experimental results and their analyses are performed. The experimental results are explained reasonably. The consistency between numerical results and experimental results implies that the numerical calculation model used in this paper can well describe the mechanical action of the laser on the target.     

Experimental study and computer simulations of the implantation of laser plasma ions in pulsed electric fields

Results are presented from experimental studies of pulsed plasma flows generated by nanosecond laser pulses with an intensity of 7 × 10⁸ W/cm² from a solid-state target in a strong electric field. The current pulses through the laser target and the depth distributions of the iron ions implanted in a silicon substrate to which a negative high-voltage pulse was applied are measured. The physical processes occurring in laser plasma with an initial iron ion density of 6 × 10¹⁰ cm⁻³ are simulated numerically by the particle-in-cell method for different delay times and different shapes of the accelerating high-voltage pulse. The model developed allows one to calculate the ion flows onto the processed substrate, the electron flows onto the target, and the energy spectra of the implanted ions. The results from computer simulations are found to be in good agreement the experimental data.     

Simulation of Processes on a Target Induced by a Jet System and Laser Radiation

The mathematical model and the results of computer simulation of the processes of compression, heating, and energy release in a target under the hybrid influence of a system of pulsed plasma jets and intense laser radiation in magnetic inertial thermonuclear fusion are presented. The hybrid effects of intense energy flows on a single-layer cylindrical target are compared in the experimental and simulation conditions. The main radiation magnetic and gas dynamic parameters of the compressed and heated target plasma are evaluated as well.     

激光聚变内爆流体不稳定性基础问题研究进展

激光聚变有望一劳永逸地解决人类的能源问题,因而受到国际社会的普遍重视,一直是国际研究的前沿热点。目前实现激光惯性约束聚变所面临的最大科学障碍（属于内禀困难）是对内爆过程中高能量密度流体力学不稳定性引起的非线性流动的有效控制,对其研究涵盖高能量密度物理、等离子体物理、流体力学、计算科学、强冲击物理和高压原子物理等多个学科,同时还要具备大规模多物理多尺度多介质流动的数值模拟能力和高功率大型激光装置等研究条件。作为新兴研究课题,高能量密度非线性流动问题充满了各种新奇的现象亟待探索。此外,流体力学不稳定性及其引起的湍流混合,还是天体物理现象（如星系碰撞与合并、恒星演化、原始恒星的形成以及超新星爆炸）中的重要过程,涉及天体物理的一些核心研究内容。本文首先综述了高能量密度非线性流动研究的现状和进展,梳理了其中的挑战和机遇。然后介绍了传统中心点火激光聚变内爆过程发生的主要流体力学不稳定性,在大量分解和综合物理研究基础上,凝练出了目前制约美国国家点火装置（NIF）内爆性能的主要流体不稳定性问题。接下来,总结了国外激光聚变流体不稳定性实验物理的研究概况。最后,展示了内爆物理团队近些年在激光聚变内爆流体不稳定性基础性问题方面的主要研究进展。该团队一直从事激光聚变内爆非线性流动研究与控制,以及聚变靶物理研究与设计,注重理论探索和实验研究相结合,近年来在内爆重要流体力学不稳定性问题的解析理论、数值模拟和激光装置实验设计与数据分析等方面取得了一系列重要成果,有力地推动了该研究方向在国内的发展。     

等离子体湍流理论研究的最新进展

评述和分析了磁约束核聚变理论研究、数值模拟和实验研究等方面最近几年共同关心的一个重要问题——寻找托卡马克等离子体湍流中的带状剪切流（zonal flows）。简要介绍了作者最近对电阻性-重力模湍流中的带状剪切流的研究结果。 Progress of the research on the zonal flows in tokamak plasma turbulence is surveyed, especially it is reviewed that the zonal flows observed in the experiments and numerical simulations on atmosphere turbulence and ocean turbulence and the discovery of H-mode in tokamak experiments how lead the researchers in magnetic confinement fusion community to find out the existence of the zonal flow in tokamak plasma turbulence and subsequently give the experimental verification of its existence. Finally, the results of our research on zonal flow generation and evolution in resistive-g mode turbulence are presented in brief.     

Thermal and reactive nonequilibrium of magnetoactive plasma flows

Based on multifluid equations derived from the Boltzman equation under inclusion of ionization and recombination reactions, the following convective processes in magnetoactive plasma flows transverse to homogeneous magnetic fields are investigated:The build-up of the nonuniform plasma state, the intercomponent thermal nonequilibrium and the reactive nonequilibrium in dependence of the flow coordinate. Due to the complexity of the problem, a numerical approach is used, parallel to which an approximate analytical theory is developed.     

Generation of Broadband High Harmonics through Linear Mode Conversion in Inhomogeneous Plasmas

The generation of high order harmonics from an inhomogeneous ovderdense plasma target irradiated by an ultrashort intense laser pulse is studied by numerical simulation. During such interaction, ultrafast electron bunches are generated and excite electron plasma oscillations as they pass through the overdense target. These plasma oscillations will emit high-frequency electromagnetic emission by linear mode conversion. Instead of the integer harmonies generation, the emission appears with a broadband and even continuous spectrum corresponding to the electron plasma frequency range of the inhomogeneous plasma density.     

Supersonic air-scoop flows of a weakly ionized gas in external electromagnetic field

Steady supersonic flows of a weakly ionized gas in external electric and magnetic fields computed for the problem of flow control in the air scoop of a supersonic vehicle. The numerical simulation was based on a second-order accurate marching scheme for the Euler equations. Effects of plasma and electromagnetic-field parameters on flow structure and distributions of flow variables are analyzed.     

Features of the generation of a collisionless electrostatic shock wave in a laser-ablation plasma

The appearance of a density bump is experimentally revealed in an electrostatic shock wave during the ablation of an aluminum foil by a femtosecond laser pulse. The numerical simulation shows that this phenomenon can be explained by the generation of a packet of ion acoustic waves under the action of high-energy electron flows in a collisionless plasma. It is found that, for the formation and maintenance of the dense plasma layer in the shock wave, the contributions of accelerated ions overtaking it and wave-captured ions of the background plasma formed by a nanosecond laser prepulse in the process of ablation are significant.     

Numerical Simulation on Expansion Process of Ablation Plasma Induced by Intense Pulsed Ion Beam

We present a one-dimensional time-dependent numerical model for the expansion process of ablation plasma induced by intense pulsed ion beam （IPIB）. The evolutions of density, velocity, temperature, and pressure of the ablation plasma of the aluminium target are obtained. The numerical results are well in agreement with the relative experimental data. It is shown that the expansion process of ablation plasma induced by IPIB includes strongly nonlinear effects and that shock waves appear during the propagation of the ablation plasma.     

Melting and crystallization dynamics of single-crystal silicon exposed to compression plasma flows

The melting and crystallization of single-crystal silicon wafers exposed to compression plasma flows generated by quasi-stationary plasma accelerators are studied by numerical simulation. The results include the phase transition kinetics described on the basis of the Kolmogorov equation. The space-time characteristics of the melting and crystallization processes for variously shaped plasma pulses are discussed. Based on the experimental data and estimates, it is concluded that thermoelectric instability plays an essential role in the formation of three-dimensional periodic structures.     

等离子体源离子注入鞘层随时空演化的数值模拟

在等离子体源离子注入装置中(PSⅡ),靶被直接放入均匀等离子体源中,并在其上加负高压脉冲,从而使靶的周围形成了一个只存在离子的等离子体鞘层,离子以较高速度注入靶中。利用冷流体方程和泊松方程组成的非线性方程组对一维平面几何鞘层时空演化过程进行数值模拟,所采用的负高电势脉冲波形为方波和梯形波,得到了电子密度、离子密度、电势分布的时空演化曲线,并首次对梯形波下降沿时间内鞘层时空演化作了详细讨论。     

Orbital motion of dust particles in an rf magnetron discharge. Ion drag force or neutral atom wind force

Microparticles with sizes up to 130 μm have been confined and the velocity and diameter of particles in a plasma trap of an rf magnetron discharge with an arc magnetic field have been simultaneously measured. The motion of the gas induced by electron and ion cyclotron currents has been numerically simulated using the Navier-Stokes equation. The experimental and numerical results confirm the mechanism of the orbital motion of dust particles in the magnetron discharge plasma that is associated with the orbital motion of the neutral gas accelerated by electron and ion drift flows in crossed electric and magnetic fields.     

Experimental investigation of magnetodynamic processes in two-phase flows

Power plants for use in investigation of flows of products of combustion of metallized plasma-forming fuels in a combustion chamber, Laval nozzle, and MHD channel and for determining the physical characteristics of plasma are described. The results of experimental investigations and numerical calculations of the combustion of fuel in a combustion chamber and flows of combustion products through gas-dynamic channels were employed in the design of the power plants.Scientific-Research Institute of Applied Mathematics and Mechanics at Tomsk State University. Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 4, pp. 92–100, April, 1993.     

激光与固体靶面烧蚀等离子体的能量耦合计算

 强激光辐照下固体靶表面迅速汽化产生靶蒸气等离子体,激光穿过等离子体区到达固体靶表面的过程就是激光束与等离子体的能量耦合与交换过程。采用具有五阶精度的WENO差分格式和简易等离子体状态方程模型对激光与等离子体相互作用的复杂物理过程进行了数值计算,分析了激光束能量在等离子体区中的吸收、屏蔽效应等动态耦合规律以及激光支持等离子体前驱冲击波传播。数值模拟结果表明：激光能量是支持靶面等离子体运动的唯一原因,能量屏蔽效应对激光与等离子体能量耦合有很大影响,通过控制激光脉冲宽度,可以合理调节屏蔽效应的影响。