Formation of compression shocks in a nonequilibrium plasma flow in a magnetic field

Plasma flow in a linearly widening, ideally sectioned, short-circuited magnetohydrodynamic (MHD) channel is studied. MHD flows are classified into two types: continuous flows and flows with a compressional MHD shock in plasmas that are stable and unstable against the onset of ionization instability. Specific features in the evolution of a stationary compression MHD shock are investigated, and its position as a function of the Stewart number is determined. It is found that, in a plasma flow in which ionization instability develops, a compression MHD shock arises at lower values of the MHD interaction parameter than in a stable plasma flow. An unidentified type of instability of MHD discontinuities is revealed.     

高超声速飞行器磁控热防护系统建模分析

针对高超声速飞行器防热, 搭建了螺线管磁控热防护系统的物理模型. 采用低磁雷诺数磁流体数学模型, 分析了外加磁场强度及磁场形态对磁控热防护效果的影响. 对比了三种磁场类型(磁偶极子、螺线管、均布磁场)下磁控热防护效果的差异, 分析了螺线管几何参数对磁控热防护效果的影响. 研究表明, 磁场降低表面热流作用存在“饱和现象”; 三种磁场形态的磁控热防护能力从小到大依次为磁偶极子、螺线管、均布磁场; 相同驻点磁感应强度条件下, 增大螺线管半径有利于提高磁控热防护效果, 缩短螺线管与驻点距离不利于驻点和肩部防热, 螺线管长度对磁控热防护效果影响相对较小.     

Experimental investigation of the magnetohydrodynamic parachute effect in a hypersonic air flow

New data on experimental implementation of the magnetohydrodynamic (MHD) parachute configuration in an air flow with Mach number M = 6 about a flat plate are considered. It is shown that MHD interaction near a flat plate may transform an attached oblique shock wave into a normal detached one, which considerably extends the area of body-incoming flow interaction. This effect can be employed in optimizing return space vehicle deceleration conditions in the upper atmosphere.     

Numerical Simulation of the 12 May 1997 CME Event

Our newly developed CESE MHD model is used to simulate sun-earth connection event with the well-studied 12 May 1997 CME event as an example. The main features and approximations of our numerical model are as follows： （1） The modifed conservation element and solution element （CESE） numerical scheme in spherical geometry is implemented in our code. （2） The background solar wind is derived from a 3D time-dependent numerical MHD model by input measured photospheric magnetic fields. （3） Transient disturbances are derived from solar surface by introducing a mass flow of hot plasma. The numerical simulation has enabled us to predict the arrival of the interplanetary shock and provided us with a relatively satisfactory comparison with the WIND spacecraft observations.     

An adaptive moving mesh method for two-dimensional ideal magnetohydrodynamics

This paper presents an adaptive moving mesh algorithm for two-dimensional (2D) ideal magnetohydrodynamics (MHD) that utilizes a staggered constrained transport technique to keep the magnetic field divergence-free. The algorithm consists of two independent parts: MHD evolution and mesh-redistribution. The first part is a high-resolution, divergence-free, shock-capturing scheme on a fixed quadrangular mesh, while the second part is an iterative procedure. In each iteration, mesh points are first redistributed, and then a conservative-interpolation formula is used to calculate the remapped cell-averages of the mass, momentum, and total energy on the resulting new mesh; the magnetic potential is remapped to the new mesh in a non-conservative way and is reconstructed to give a divergence-free magnetic field on the new mesh. Several numerical examples are given to demonstrate that the proposed method can achieve high numerical accuracy, track and resolve strong shock waves in ideal MHD problems, and preserve divergence-free property of the magnetic field. Numerical examples include the smooth Alfvén wave problem, 2D and 2.5D shock tube problems, two rotor problems, the stringent blast problem, and the cloud–shock interaction problem.     

Supersonic flow about a body exposed to electric and magnetic fields

The feasibility of using nonmechanical (electrogasdynamic, EGD, and magnetohydrodynamic, MHD) methods to control shock-wave configurations emerging in supersonic flows is investigated. In the EGD method, the flow is heated by a gas discharge; in the MHD one, the flow is influenced by a Lorentz force arising in a gas discharge upon applying a magnetic field. The influence of the gas discharge and MHD interaction on the position of a detached shock wave appearing in a supersonic xenon flow about a semicylindrical body is studied. A discharge is initiated in the immediate vicinity of the leading edge of the body, and the variation of the shock wave position with the intensity of the discharge (discharge current density) is traced when the influence of the EGD action increases and/or an external magnetic field is applied (the influence of the MHD action increases). Preliminary data for a supersonic air flow about a body are presented.     

On the shock induced instabilities in collisionless plasma

Firehose, mirror and ion-acoustic instabilities behind MHD shock waves are discussed for collisionless anisotropic plasma with heat fluxes. For the parallel shock wave it has been demonstrated that initially stable plasma can be destabilized by the shock which leads to turbulence generation in the downstream flow. Upstream parameter domains are determined where such destabilization occurs.     

Anomalous flow deflection at earth's low-Alfvén-Mach-Number bow shock

Earth's magnetosphere is an obstacle to the supersonic solar wind and the bow shock is formed in the front side of it. In ordinary hydrodynamics, the flow decelerated at the shock is diverted around the obstacle symmetrically about the Earth-Sun line, which is indeed observed in the magnetosheath most of the time. Here we show a case under a very low-density solar wind in which duskward flow was observed in the dawnside magnetosheath. A Rankine-Hugoniot test shows that the magnetic effect is crucial for this "wrong flow" to appear. A full three-dimensional magnetohydrodynamics (MHD) simulation of the situation confirming this interpretation and earlier simulations is also performed. It is illustrated that in addition to the "wrong flow" feature, various peculiar characteristics appear in the global picture of the MHD flow interaction with the obstacle.     

Shock waves in three-dimensional transonic flow with weak spatial heat release

The presence or absence of shock waves in the case of weak 3D heat release to homogeneous uniform transonic gas flow is analyzed. It is well known that shock waves can exist in a steady-state flow in the 2D case for a small transonic similarity number. It is shown that in the main approximation, shock waves are absent in the axisymmetric and 3D versions. In the next approximation, in which the logarithm factor of a small parameter is absent, shock waves can exist.     

Local effect of electric and magnetic fields on the position of an attached shock in a supersonic diffuser

Effective ways for controlling shock wave configurations by means of external actions are sought. One such way is a local effect of electric and magnetic fields. In this paper, the local effect of external fields is implemented by current localization in a limited region of a diffuser. The experiment is carried out in a diffuser providing the complete internal compression of the gas with a Mach number at the inlet M=4.3. As a working medium, a xenon plasma is used. The plasma flow is formed in a shock tube equipped with an accelerating nozzle. Two ways of current localization are tested. In the first one, the diffuser inlet is a short channel of Faraday generator type. In this case, the ponderomotive force basically decelerates or accelerates the flow depending on the direction of the electric field. In the second way, the current flows through a narrow near-wall region between adjacent electrodes. In this case, the ponderomotive force compresses or expands the gas. In both cases, it is shown that the angle of an attached shock due to MHD interaction can be both decreased and increased. The central problem with the MHD control of shock waves is near-electrode and near-wall phenomena.     

FCI结构及其电导率对MHD效应影响实验研究

实验研究了通道插件(FCI)压力平衡孔(PEH)和压力平衡狭缝(PES)结构、绝缘和非绝缘材质对磁流体动力学(MHD)效应的影响,结果表明：带PEH或PES、绝缘或非绝缘材质的FCI产生的流速分布有很大的差异;对MHD压降而言,FCI结构PEHs或PES差异的影响大于FCI绝缘或非绝缘材质差异引起的影响;对流速分布而则刚好相反。二次流引起的局部流速骤增的MHD效应(简称为S-MHD效应)对控制管道内流速分布和降低MHD压降有作用,这些实验结果为弄清FCI的MHD效应机制和液态金属包层设计提供了实验认知与有价值的参考。     

FCI�ṹ����絼�ʶ�MHDЧӦӰ��ʵ���о�

实验研究了通道插件(FCI)压力平衡孔(PEH)和压力平衡狭缝(PES)结构、绝缘和非绝缘材质对磁流体动力学(MHD)效应的影响,结果表明：带PEH或PES、绝缘或非绝缘材质的FCI产生的流速分布有很大的差异;对MHD压降而言,FCI结构PEHs或PES差异的影响大于FCI绝缘或非绝缘材质差异引起的影响;对流速分布而则刚好相反。二次流引起的局部流速骤增的MHD效应(简称为S-MHD效应)对控制管道内流速分布和降低MHD压降有作用,这些实验结果为弄清FCI的MHD效应机制和液态金属包层设计提供了实验认知与有价值的参考。     

Numerical modelling of magnetogasdynamic processes in the pulsed setup duct

A computational model for processes in the duct of an experimental setup consisting of a shock tube and an MHD duct is presented. The one-dimensional model is used for determining the flow characteristics in the entire setup duct, and the three-dimensional model is used for studying the current layer dynamics in the MHD channel. Computations have enabled the elucidation of flow structure and of the peculiarities of current layer formation.     

An exponential compact difference scheme for solving 2D steady magnetohydrodynamic (MHD) duct flow problems

In this article, an exponential high-order compact (EHOC) difference scheme on the nine-point stencil is developed for the solution of the coupled equations representing the steady incompressible, viscous magnetohydrodynamic (MHD) flow through a straight channel of rectangular section. A key property of the EHOC scheme is that it has excellent stability and higher accuracy so that the high gradients near the boundary layer areas can be effectively resolved without refining the mesh. Numerical experiments are carried out to validate the performance of the currently proposed scheme. Computation results of the MHD flow in the 2D square-channel problems with different wall conductivities are presented for Hartmann numbers ranging from 10 to 10⁶. The numerical solutions obtained with the newly developed EHOC scheme are also compared with analytic solutions and numerical results by other available methods in the literature.     

Use of MHD systems in hypersonic aircraft

The possibilities of using magnetohydrodynamic (MHD) systems on hypersonic aircraft are discussed. The distinctive features of using MHD systems in the flow path of ramjet engines are examined. A quasi-one-dimensional mathematical model for the engine is presented which includes the MHD interaction with the flow. It is shown that the specific impulse of an engine system can be raised by using MHD systems. Zh. Tekh. Fiz. 68, 43–47 (November 1998)     

一种适用于磁流体力学切向间断的HLLC黎曼算子

磁场的存在使得磁流体力学特征波模不同于流体力学, 因此直接由流体力学HLLC黎曼算子导出的HLLC双中间态在交界的间断面会出现不守恒的问题。通常降级采用HLL单磁场中间态代替HLLC双磁场中间态以实现守恒和计算稳定, 代价是切向间断的模拟精度不足。本文对此进行改进, 在模拟切向间断时仍然保留原有的HLLC双磁场中间态, 同时各守恒量仍然能够满足Toro相容条件; 改进型HLLC算子在间断两侧的磁场分量存在差异, 因此能够更精确还原切向间断面。基于数值测试, 包括一维激波管和切向间断的时变模拟, 以及地球磁层三维数值模拟, 将模拟结果进行对比, 结果表明: 相比于已发展的HLLC算子, 改进型HLLC算子对切向间断具有更好的捕捉精度, 能够达到或接近耗时更多的HLLD算子的模拟精度。     

A Complete 2D Stability Analysis of Fast MHD Shocks in an Ideal Gas

 An algorithm of numerical testing of the uniform Lopatinski condition for linearized stability problems for 1-shocks is suggested. The algorithm is used for finding the domains of uniform stability, neutral stability, and instability of planar fast MHD shocks. A complete stability analysis of fast MHD shock waves is first carried out in two space dimensions for the case of an ideal gas. Main results are given for the adiabatic constant γ=5/3 (mono-atomic gas), that is most natural for the MHD model. The cases γ=7/5 (two-atomic gas) and γ>5/3 are briefly discussed. Not only the domains of instability and linear (in the usual sense) stability, but also the domains of uniform stability, for which a corresponding linearized stability problem satisfies the uniform Lopatinski condition, are numerically found for different given angles of inclination of the magnetic field behind the shock to the planar shock front. As is known, uniform linearized stability implies the nonlinear stability, that is local existence of discontinuous shock front solutions of a quasilinear system of hyperbolic conservation laws. Received: 12 March 2002 / Accepted: 8 November 2002 Published online: 10 February 2003 Communicated by P. Constantin     

Calculation of the electric potential and Lorentz force acting on a locally ionized region in a magnetohydrodynamic flow placed in a nonuniform magnetic field

An analytical solution to electrodynamic equations for the electric potential in a locally ionized magnetohydrodynamic (MHD) flow in the nonuniform magnetic field produced by a straight-line conductor is found. Analytical formulas are obtained to evaluate the volume density of the Lorentz force and the integral Lorentz force acting on the locally ionized region of the MHD flow. It is shown that the MHD action on the locally ionized flow in the nonuniform magnetic field can be used to control the elevating force as well as the ratio of the elevating force to the drag force.     

高速飞行器磁控阻力特性

采用低磁雷诺数磁流体数学模型,对外加磁场下的高超声速半球体流场进行数值模拟.选取三种简单理想磁场(轴向、径向、周向均布磁场),分析了不同磁场类型对流场结构、气动阻力与洛伦兹阻力的影响及作用机理.研究发现,轴向磁场径向"挤压"效应使得激波外形凸出,且壁面静压存在"饱和现象";径向磁场存在轴向"外推"效应,较大的磁场强度会导致肩部形成高温区;周向磁场下感应电场的存在导致增阻效果很差.进而对比了两种相同驻点磁感应强度特殊分布磁场(偶极子磁场、螺线管磁场)下的流场,发现了不同于理想磁场的径向"扩张"效应.按增阻效果从大到小依次为径向磁场、螺线管磁场、轴向磁场、偶极子磁场、周向磁场.     

Analytical studies on the evolution processes of rarefied deuterium plasma shell Z-pinch by PIC and MHD simulations

In this paper, we analytically explore the magnetic field and mass density evolutions obtained in particle-in-cell(PIC)and magnetohydrodynamics(MHD) simulations of a rarefied deuterium shell Z-pinch and compare those results, and also we study the effects of artificially increased Spitzer resistivity on the magnetic field evolution and Z-pinch dynamic process in the MHD simulation. There are significant differences between the profiles of mass density in the PIC and MHD simulations before 45 ns of the Z-pinch in this study. However, after the shock formation in the PIC simulation,the mass density profile is similar to that in the MHD simulation in the case of using multiplier 2 to modify the Spitzer resistivity. Compared with the magnetic field profiles of the PIC simulation of the shell, the magnetic field diffusion has still not been sufficiently revealed in the MHD simulation even though their convergence ratios become the same by using larger multipliers in the resistivity. The MHD simulation results suggest that the magnetic field diffusion is greatly enhanced by increasing the Spitzer resistivity used, which, however, causes the implosion characteristic to change from shock compression to weak shock, even shockless evolution, and expedites the expansion of the shell. Too large a multiplier is not suggested to be used to modify the resistivity in some Z-pinch applications, such as the Z-pinch driven inertial confinement fusion(ICF) in a dynamic hohlraum. Two-fluid or Hall MHD model, even the PIC/fluid hybrid simulation would be considered as a suitable physical model when there exist the plasma regions with very low density in the simulated domain.