Computation of the electric potential and the Lorentz force in a locally ionized magnetohydrodynamic flow in a nonuniform magnetic field for a transverse flow past a circular cylinder

An analytical solution to electrodynamic equations is obtained for the electric potential in a locally ionized magnetohydrodynamic (MHD) flow for a transverse flow past a circular cylinder in the non-uniform magnetic field of a rectilinear conductor. Analytical formulas for computing the volume density of the Lorentz force acting on the flow in a locally ionized MHD flow are obtained for the case of the conducting and nonconducting surfaces of the cylinder. The influence of the Hall parameter and width of the MHD interaction region on the value of the Lorentz force is analyzed. It is demonstrated that the Lorentz force, which accelerates and not decelerates the flow, appears under certain conditions near the surface of the cylinder in the neighborhood of the critical point.     

Calculation of the electric potential and the Lorentz force in a transverse flow past a circular cylinder in a nonuniform magnetic field for various configurations of a locally ionized region at the cylinder surface

We obtain a solution to the equation for the electric potential in a locally ionized transverse magnetohydrodynamic flow past a circular cylinder in a nonuniform magnetic field produced by a linear conductor for various configurations of the ionization region. Analytical formulas are derived for the volume density of the Lorentz force acting on the flow in a locally ionized region. The effect of the Hall parameter and of the configuration of the region of the magnetohydrodynamic interaction on the Lorentz force is analyzed. It is shown that an increase in the Hall parameter leads to a decrease in the Lorentz force acting on the flow, and a change in the configuration of the locally ionized region makes it possible to suppress the effect of the Hall parameter on the Lorentz force.     

高速飞行器磁控阻力特性

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

A current density conservative scheme for incompressible MHD flows at a low magnetic Reynolds number. Part II: On an arbitrary collocated mesh

A conservative formulation of the Lorentz force is given here for magnetohydrodynamic (MHD) flows at a low magnetic Reynolds number with the current density calculated based on Ohm’s law and the electrical potential formula. This conservative formula shows that the total momentum contributed from the Lorentz force is conservative when the applied magnetic field is constant. For the case with a non-constant applied magnetic field, the Lorentz force has been divided into two parts: a strong globally conservative part and a weak locally conservative part.The conservative formula has been employed to develop a conservative scheme for the calculation of the Lorentz force on an unstructured collocated mesh. Only the current density fluxes on the cell faces, which are calculated using a consistent scheme with good conservation, are needed for the calculation of the Lorentz force. Meanwhile, a conservative interpolation technique is designed to get the current density at the cell center from the current density fluxes on the cell faces. This conservative interpolation can keep the current density at the cell center conservative, which can be used to calculate the Lorentz force at the cell center with good accuracy. The Lorentz force calculated from the conservative current at the cell center is equivalent to the Lorentz force from the conservative formula when the applied magnetic field is constant, which can conserve the total momentum. We will further prove that the simple interpolation scheme used in the Part I [M.-J. Ni, R. Munipalli, N.B. Morley, P.Y. Huang, M. Abdou, A current density conservative scheme for MHD flows at a low magnetic Reynolds number. Part I. On a rectangular collocated grid system, Journal of Computational Physics, in press, doi:10.1016/j.jcp.2007.07.025] of this series of papers is conservative on a rectangular grid and can keep the total momentum conservative in a rectangular grid.     

A consistent and conservative scheme for incompressible MHD flows at a low magnetic Reynolds number. Part III: On a staggered mesh

The consistent and conservative scheme developed on a rectangular collocated mesh [M.-J. Ni, R. Munipalli, N.B. Morley, P. Huang, M.A. Abdou, A current density conservative scheme for incompressible MHD flows at a low magnetic Reynolds number. Part I: on a rectangular collocated grid system, Journal of Computational Physics 227 (2007) 174–204] and on an arbitrary collocated mesh [M.-J. Ni, R. Munipalli, P. Huang, N.B. Morley, M.A. Abdou, A current density conservative scheme for incompressible MHD flows at a low magnetic Reynolds number. Part II: on an arbitrary collocated mesh, Journal of Computational Physics 227 (2007) 205–228] has been extended and specially designed for calculation of the Lorentz force on a staggered grid system (Part III) by solving the electrical potential equation for magnetohydrodynamics (MHD) at a low magnetic Reynolds number. In a staggered mesh, pressure (p) and electrical potential (φ) are located in the cell center, while velocities and current fluxes are located on the cell faces of a main control volume. The scheme numerically meets the physical conservation laws, charge conservation law and momentum conservation law. Physically, the Lorentz force conserves the momentum when the magnetic field is constant or spatial coordinate independent. The calculation of current density fluxes on cell faces is conducted using a scheme consistent with the discretization for solution of the electrical potential Poisson equation, which can ensure the calculated current density conserves the charge. A divergence formula of the Lorentz force is used to calculate the Lorentz force at the cell center of a main control volume, which can numerically conserve the momentum at constant or spatial coordinate independent magnetic field. The calculated cell-center Lorentz forces are then interpolated to the cell faces, which are used to obtain the corresponding velocity fluxes by solving the momentum equations. The “conservative” is an important property of the scheme, which can guarantee computational accuracy of MHD flows at high Hartmann number with a strongly non-uniform mesh employed to resolve the Hartmann layers and side layers. 2D fully developed MHD flows with analytical solutions available have been conducted to validate the scheme at a staggered mesh. 3D MHD flows, with the experimental data available, at a constant magnetic field in a rectangular duct with sudden expansion and at a varying magnetic field in a rectangular duct are conducted on a staggered mesh to verify the computational accuracy of the scheme. It is expected that the scheme for the Lorentz force can be employed together with a fully conservative scheme for the convective term and the pressure term [Y. Morinishi, T.S. Lund, O.V. Vasilyev, P. Moin, Fully conservative higher order finite difference schemes for incompressible flow, Journal of Computational Physics 143 (1998) 90–124] for direct simulation of MHD turbulence and MHD instability with good accuracy at a staggered mesh.     

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

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

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.     

Light-induced drift of ions in a magnetic field

The effect of a magnetic field on the ion drift in a weakly ionized gas under the combined action of the light-induced drift and light pressure is studied theoretically. It is shown that, under the action of light, a component of ion drift velocity transverse to the direction of propagation of radiation may appear in a weakly ionized gas upon the application of an external magnetic field. It is shown that the Lorentz force acting on ions in the magnetic field radically changes the dependence of the ion drift velocity on the radiation frequency detuning. It is predicted that the ion drift velocity component along the direction of radiation must reverse its sign upon an increase in the magnetic field and an anomalous light-induced drift may be observed.     

绕圆柱体自由表面磁流体流动和传热的研究

本文对在不同雷诺数下,绕圆柱体的磁流体自由表面流动及传热进行了模拟,分析了磁场对绕流圆柱尾迹和涡分离的影响,获得了两种雷诺数下的电磁力密度、流场和温度场分布。结果表明,磁场不仅影响了流动的形态,而且对湍流有抑制作用,降低了自由表面的更新机制,从而减少了传热能力;在相同的Hartmann数下,相比低雷诺数下的流动换热情况,高雷诺数下的湍流不能被完全抑制,自由表面与尾迹的相互作用也较强,因而自由表面换热也较强。     

基于OpenFOAM的投影法磁流体求解器开发与验证

在开源计算流体力学C++工具包OpenFOAM环境下开发了低磁雷诺数条件下的磁流体求解器,并进行了验证。采用投影算法求解动量方程和压力泊松方程;采用非结构网格同位相容守恒算法求解电势泊松方程、感应电流和洛伦兹力;采用边界耦合方法求解流固耦合电势场。通过对均匀磁场下导电方管和导电圆管内的完全发展磁流体层流的数值模拟和解析解的对比,对求解器进行了验证。进一步对非均匀强磁场作用下导电方管和导电圆管内完全发展磁流体层流进行了数值模拟,并与ALEX实验结果进行了比较。数值解和实验结果吻合良好。所开发的求解器可用于复杂结构强磁场作用下磁流体的数值模拟研究。     

带引流条导电矩形管磁流体湍流大涡数值模拟

为研究引流条对磁流体湍流的影响，采用自主开发的低磁雷诺数流固耦合磁流体相干结构模型大涡模拟求解器，对均匀磁场作用下平行层内带引流条导电矩形管和标准导电矩形管中液态金属湍流进行了数值模拟研究。结果表明，外加垂直流动方向的均匀磁场与流动的导电流体相互作用产生与流动方向相反的洛伦兹力，能够抑制磁流体的湍流脉动，这种抑制作用随着哈特曼数增大而增强。在弱导电率条件下，当Re=16350、Ha=212 时，两种管道中的流动均转换为层流流动状态。管道内壁面摩擦系数随着哈特曼数的增大而增大。引流条能在其近壁局部区域增强横向速度，有效激发湍流，但在弱壁面导电率条件下，带引流条导电矩形管壁面摩擦系数较标准矩形管大。     

Direct numerical simulation of the turbulent MHD channel flow at low magnetic Reynolds number for electric correlation characteristics

Direct numerical simulation (DNS) of incompressible magnetohydrodynamic (MHD) turbulent channel flow has been performed under the low magnetic Reynolds number assumption.The velocity-electric field and electric-electric field correlations were studied in the present work for different magnetic field orientations.The Kenjeres-Hanjalic (K-H) model was validated with the DNS data in a term by term manner.The numerical results showed that the K-H model makes good predictions for most components of the velocity-electric field correlations.The mechanisms of turbulence suppression were also analyzed for different magnetic field orientations utilizing the DNS data and the K-H model.The results revealed that the dissipative MHD source term is responsible for the turbulence suppression for the case of streamwise and spanwise magnetic orientation,while the Lorentz force which speeds up the near-wall fluid and decreases the production term is responsible for the turbulence suppression for the case of the wall normal magnetic orientation.     

Magnetohydrodynamic Electroosmotic Flow with Patterned Charged Surface and Modulated Wettability in a Parallel Plate Microchannel*

This paper investigates the magnetohydrodynamic (MHD) electroosmotic flow (EOF) of Newtonian fluid through a zeta potential modulated parallel plate microchannel with patterned hydrodynamic slippage. The driven mechanism of the flow originates from the Lorentz force generated by the interaction of externally imposed lateral electric field $E_y$ and vertical magnetic field $B_z$ and electric field force produced by an externally applied electric field $E_x$. It is assumed that the wall zeta potential and the slip length are periodic functions of axial coordinate $x$, an analytical solution of the stream function is achieved by utilizing the method of separation of variables and perturbation expansion. The pictures of streamlines are plotted and the vortex configurations produced in flow field due to patterned wall potential and hydrodynamic slippage are discussed. Based on the stream function, the velocity field and volume flow rate are obtained, which are greatly depend on some dimensionless parameters, such as slip length $l_s$, electrokinetic width $\lambda$, the amplitude $\delta$ of the patterned slip length, the amplitude $m$ of the modulated zeta potential and Hartmann number $Ha$. The variations of velocity and volume flow rate with these dimensionless parameters are discussed in details. These theoretical results may provide some guidance effectively operating micropump in practical nanofluidic applications.     

A current density conservative scheme for incompressible MHD flows at a low magnetic Reynolds number. Part I: On a rectangular collocated grid system

A consistent, conservative and accurate scheme has been designed to calculate the current density and the Lorentz force by solving the electrical potential equation for magnetohydrodynamics (MHD) at low magnetic Reynolds numbers and high Hartmann numbers on a finite-volume structured collocated grid. In this collocated grid, velocity (u), pressure (p), and electrical potential (φ) are located in the grid center, while current fluxes are located on the cell faces. The calculation of current fluxes on the cell faces is conducted using a conservative scheme, which is consistent with the discretization scheme for the solution of electrical potential Poisson equation. A conservative interpolation is used to get the current density at the cell center, which is used to conduct the calculation of Lorentz force at the cell center for momentum equations. We will show that both “conservative” and “consistent” are important properties of the scheme to get an accurate result for high Hartmann number MHD flows with a strongly non-uniform mesh employed to resolve the Hartmann layers and side layers of Hunt’s conductive walls and Shercliff’s insulated walls. A general second-order projection method has been developed for the incompressible Navier–Stokes equations with the Lorentz force included. This projection method can accurately balance the pressure term and the Lorentz force for a fully developed core flow. This method can also simplify the pressure boundary conditions for MHD flows.     

交错网格下MHD相容守恒格式的发展

在低磁场雷诺数条件下,基于电势泊松方程,发展了交错网格下可以精确计算电流和洛伦兹力(电磁力)的相容守恒格式。采用压力为变量的原始变量法求解不可压缩Navier-Stokes方程,所计算的电流满足电荷守恒定律,所计算的电磁力满足动量守恒定律。对金属流体在Hartmann数50~5000范围内验证了格式的精确性。交错网格下相容守恒格式的发展为后续MHD稳定性分析、湍流的大涡模拟及直接数值模拟提供很好的选择。     

Performance Analysis of a Diagonal MHD Accelerator for Space Propulsion

The performance of a diagonal magnetohydrodynamic (MHD) accelerator has been numerically investigated. Studies were carried out using air plasma as a working gas in an equilibrium condition based on the MHD Augmented Propulsion Experiment channel designed by NASA. The MacCormack scheme is employed in order to solve the set of differential equations with MHD approximations. The fundamental performance of a diagonal MHD accelerator considering both flow performance along the channel and propulsion performance has been evaluated under various applied input currents and magnetic fields. The optimum performance is dominated by ${bf j} times {bf B}$ Lorentz body force acceleration, while it is increased with Joule heating and the ${bf u} times {bf B}$ term's contribution, which are detrimental to the propulsion performance. Moreover, friction forces resist the flow performance, particularly near the channel exit.      

Effects of hall current on MHD Couette flow in a rotating system with arbitrary magnetic field

Author has studied the MHD Couette flow in a rotating environment with non- conducting walls in the presence of an arbitrary magnetic field. The solution in dimensionless form contains four pertinent flow parameters, viz. the Hartmann number, the rotation parameter which is the reciprocal of the Ekman number, the Hall current parameter, and the angle of inclination of the magnetic field to the positive direction of the axis of rotation. An interplay of hydromagnetic force and Coriolis force with an inclusion of Hall current plays a significant role in determining the MHD flow behaviour. The velocity and induced magnetic field distributions are depicted graphically. Also, the numerical results of shear stresses and the rate of mass flows are presented graphically.     

The Modified Magnetohydrodynamical Equations

After finding the really self-consistent electromagnetic equations for a plasma, we proceed in a similar fashion to find how the magnetohydrodynamical equations have to be modified accordingly. Substantially this is done by replacing the “Lorentz“ force equation by the correct (in our case) force equation. Formally we have to use the vector potential instead of the magnetic field intensity. The appearance of the formulae presented is the one of classical vector analysis. We thus find a set of eight equations in eight unknowns, as previously known concerning the traditional MHD equations.     

Characteristics of high-power radiating imploding discharge with cold start

A qualitative model of the dynamics of a multiterawatt radiating Z-pinch with cold start and high rate of current rise is proposed. The model is used to analyze discharges with currents I ～ 2–5 MA (with dI/dt > 10¹³ A/s) through uniform or structured plasma-producing loads, including wire arrays. The most important consequence of cold start is that spatially nonuniform plasma production is prolonged to almost the entire current rise time. Under these conditions, the Ampére force begins to play a dominant role in the plasma dynamics before the plasma-producing load is completely transformed into an accelerated plasma. The results of computations of wire-array vaporization are presented. A formula is proposed for estimating the highest attainable velocity of plasma flow into a heterogeneous liner driven by the Ampére force. It is shown that local imbalance between radial motion of the produced plasma and supply of the plasma-producing substance to be ionized leads to axially nonuniform breakthrough of magnetic flux into the liner, which precedes plasma collapse. The magnetic-flux breakthrough gives rise to a chaotic azimuthal-axial plasma structure consisting of radial plasma jets of relatively small diameter, which is called a radial plasma rainstorm. The breaking-through azimuthal magnetic flux obstructs further current flow in the breakthrough region. Analyses of Z-pinch implosion based on the theory of Rayleigh-Taylor instability or the snowplow model are incorrect under the plasma-rainstorm conditions. The processes taking place in a stagnant Z-pinch include conversion of the energy carried by the current-generated magnetic field into turbulent MHD flow of the ion component of the plasma, its convective mixing with magnetic field, heating, energy transfer from ions to electrons, and emission from the plasma. Under typical experimental conditions, emission plays a key role in the energy balance in an imploding pinch. Z-pinch is modeled by an electric-circuit component that has a time-dependent nonlinear impedance and consumes the magnetic energy supplied by a generator through a magnetically insulated transmission line (MITL). The peak power reached in the circuit is comparable to the peak soft X-ray power output emitted by the pinch in terms of magnitude and timing. Optimum matching conditions are formulated for the generator-MITL-pinch circuit.     

Three-dimensional structures of MHD flow in a disk generator

Three-dimensional nonuniform plasmas and boundary layers have been studied numerically under an MHD interaction. The nonuniform plasma of “streamer” owing to weak ionization of seed material has a spiral structure in the r-&thetas; plane, and the plasma becomes almost uniform between the walls in the r-z plane. This structure is almost the same as that in our previous paper (1997), where the gas (heavy particle) properties are assumed to be invariant and steady. In addition to the streamer, the nonuniform plasma of “domain” owing to weak ionization of noble gas is revealed. The domain has the structure perpendicular to the streamer. In a strong MHD interaction case, the static pressure considerably increases in the upstream region of a generation channel, and the pseudo-shock waves appear in the generator, but the plasma is almost uniform along the &thetas; direction. The boundary layer in the strong MHD interaction is considerably thick, and in the broad region near the wall the Hall current flows reversely. In the weak MHD interaction case, the plasma forms a nonuniform structure along the &thetas; direction, and the Hall current does not always flow in the opposite direction even on the insulator wall since the azimuthal electric field is not zero