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.     

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稳定性分析、湍流的大涡模拟及直接数值模拟提供很好的选择。     

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

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

液态包层流动通道插件三维MHD效应数值模拟

用直接数值方法对高效液态锂铅包层内的金属流体三维MHD效应进行分析。用投影法对包含洛仑兹力源项的不可压Navier．Stokes方程求解,用相容守恒格式计算电磁力。研究了不同材料的流动通道插件（FCI）对金属磁流体流速、MHD压降和电流流线分布的影响。主要分析了以下三种情况：无FCI插件的通道内的流动状况;加入绝缘材料...     

Consistent projection methods for variable density incompressible Navier–Stokes equations with continuous surface forces on a rectangular collocated mesh

Two consistent projection methods of second-order temporal and spatial accuracy have been developed on a rectangular collocated mesh for variable density Navier–Stokes equations with a continuous surface force. Instead of the original projection methods (denoted as algorithms I and II in this paper), in which the updated cell center velocity from the intermediate velocity and the pressure gradient is not guaranteed solenoidal, the consistent projection methods (denoted as algorithms III and IV) obtain the cell center velocity based on an interpolation from a conservative fluxes with velocity unit on surrounding cell faces. Dependent on treatment of the continuous surface force, the pressure gradient in algorithm III or the sum of the pressure gradient and the surface force in algorithm IV at a cell center is then conducted from the difference between the updated velocity and the intermediate velocity in a consistent projection method. A non-viscous 3D static drop with serials of density ratios is numerically simulated. Using the consistent projection methods, the spurious currents can be greatly reduced and the pressure jump across the interface can be accurately captured without oscillations. The developed consistent projection method are also applied for simulation of interface evolution of an initial ellipse driven by the surface tension and of an initial sphere bubble driven by the buoyancy with good accuracy and good resolution.     

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

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

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.     

数值模拟磁场作用下气泡在导电流体中运动

在自适应网格上,采用VOF方法捕捉界面,相容守恒格式计算电流及电磁力,发展了金属流体自由界面MHD数值方法。通过数值模拟磁场作用下不同Hartmann数的气泡在导电溶液中的运动和变形,分析磁场对气泡以及流场的影响,同时给出诱导电场和电流的分布。为进一步深入研究冶金及热核聚变相关的金属流体在强磁场作用下的自由界面流打下基础。     

On the calculation of steady-state magnetohydrodynamic flows of liquid metals in circular ducts of a rectangular cross section

A numerical code has been designed to calculate two-dimensional steady-state magnetohydrodynamic (MHD) flows of incompressible conducting fluids (liquid metals) in linear and circular thin-wall ducts of a rectangular cross section. The flows are caused by the Lorentz force J × B that appears when an electric current passes through a fluid placed in a vertical uniform magnetic field. The code is the generalization of the well-known iteration Gauss-Seidel method to the case of a set of elliptical equations. The method proposed can be used to calculate steady-state flows over wide ranges of Hartmann (Ha = 1–10³) and Reynolds (Re = 1–10⁶) numbers.     

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.     

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

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

Approximation of the inductionless MHD problem using a stabilized finite element method

In this work, we present a stabilized formulation to solve the inductionless magnetohydrodynamic (MHD) problem using the finite element (FE) method. The MHD problem couples the Navier–Stokes equations and a Darcy-type system for the electric potential via Lorentz’s force in the momentum equation of the Navier–Stokes equations and the currents generated by the moving fluid in Ohm’s law. The key feature of the FE formulation resides in the design of the stabilization terms, which serve several purposes. First, the formulation is suitable for convection dominated flows. Second, there is no need to use interpolation spaces constrained to a compatibility condition in both sub-problems and therefore, equal-order interpolation spaces can be used for all the unknowns. Finally, this formulation leads to a coupled linear system; this monolithic approach is effective, since the coupling can be dealt by effective preconditioning and iterative solvers that allows to deal with high Hartmann numbers.     

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.     

Experimental and theoretical investigation of high-pressure arcs.II. The magnetically deflected arc (three-dimensional modeling)

While Part I deals with cylindrical arcs, Part II studies the influence of transverse magnetic fields on the arc column for ambient pressures of 0.1-5.0 MPa. If exposed to a magnetic induction of several millitesla, the column of an arc is deflected by the Lorentz forces. In this paper, heat transfer and fluid flow with coupled electromagnetic forces are modeled for the magnetically deflected arc. To verify the predictions, the three-dimensional temperature distributions of the arc column are determined by line and continuum radiation measurements using tomographic methods. These temperature maps are compared with the results of the numerical simulations. To gain insight into the physical professes of the discharge and to make the arc properties available which are not readily measured, a self-consistent numerical model of the arc column is applied to the time-dependent and three-dimensional case. The temperature, velocity, pressure, and current densities are predicted by solving the conservation equations for mass, momentum, and energy, and Ohm's and Biot-Savart's law using material functions of the plasma. A control volume approach facilitates a numerically conservative scheme for solving the coupled partial differential equations. The predictions are in fair agreement with experimental results. A time-dependent fully implicit simulation of the arc was used to investigate the arc instabilities for large magnetic inductions     

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.     

Division of the energy and of the momentum of electromagnetic waves in linear media into electromagnetic and material parts

We defend a natural division of the energy density, energy flux and momentum density of electromagnetic waves in linear media in electromagnetic and material parts. In this division, the electromagnetic part of these quantities have the same form as in vacuum when written in terms of the macroscopic electric and magnetic fields, the material momentum is calculated directly from the Lorentz force that acts on the charges of the medium, the material energy is the sum of the kinetic and potential energies of the charges of the medium and the material energy flux results from the interaction of the electric field with the magnetized medium. We present reasonable models for linear dispersive non-absorptive dielectric and magnetic media that agree with this division. We also argue that the electromagnetic momentum of our division can be associated with the electromagnetic relativistic momentum, inspired on the recent work of Barnett [Phys. Rev. Lett. 104 (2010) 070401] that showed that the Abraham momentum is associated with the kinetic momentum and the Minkowski momentum is associated with the canonical momentum.     

Circuit-field coupled finite element analysis method for an electromagnetic acoustic transducer under pulsed voltage excitation

This paper presents an analytical method for electromagnetic acoustic transducers (EMATs) under voltage excitation and considers the non-uniform distribution of the biased magnetic field. A complete model of EMATs including the non-uniform biased magnetic field, a pulsed eddy current field and the acoustic field is built up. The pulsed voltage excitation is transformed to the frequency domain by fast Fourier transformation (FFT). In terms of the time harmonic field equations of the EMAT system, the impedances of the coils under different frequencies are calculated according to the circuit-field coupling method and Poynting’s theorem. Then the currents under different frequencies are calculated according to Ohm’s law and the pulsed current excitation is obtained by inverse fast Fourier transformation (IFFT). Lastly, the sequentially coupled finite element method (FEM) is used to calculate the Lorentz force in the EMATs under the current excitation. An actual EMAT with a two-layer two-bundle printed circuit board (PCB) coil, a rectangular permanent magnet and an aluminium specimen is analysed. The coil impedances and the pulsed current are calculated and compared with the experimental results. Their agreement verified the validity of the proposed method. Further- more, the influences of lift-off distances and the non-uniform static magnetic field on the Lorentz force under pulsed voltage excitation are studied.     

Force law in material media,hidden momentum and quantum phases

We address to the force law in classical electrodynamics of material media, paying attention on the force term due to time variation of hidden momentum of magnetic dipoles. We highlight that the emergence of this force component is required by the general theorem, deriving zero total momentum for any static configuration of charges/currents. At the same time, we disclose the impossibility to add this force term covariantly to the Lorentz force law in material media. We further show that the adoption of the Einstein–Laub force law does not resolve the issue, because for a small electric/magnetic dipole, the density of Einstein–Laub force integrates exactly to the same equation, like the Lorentz force with the inclusion of hidden momentum contribution. Thus, none of the available expressions for the force on a moving dipole is compatible with the relativistic transformation of force, and we support this statement with a number of particular examples. In this respect, we suggest applying the Lagrangian approach to the derivation of the force law in a magnetized/polarized medium. In the framework of this approach we obtain the novel expression for the force on a small electric/magnetic dipole, with the novel expression for its generalized momentum. The latter expression implies two novel quantum effects with non-topological phases, when an electric dipole is moving in an electric field, and when a magnetic dipole is moving in a magnetic field. These phases, in general, are not related to dynamical effects, because they are not equal to zero, when the classical force on a dipole is vanishing. The implications of the obtained results are discussed.     

Can Torsion Play a Role in Angular Momentum Conservation Law?

In Einstein-Cartan theory, by the use of thegeneral Noether theorem, the general covariantangular-momentum conservation law is obtained withrespect to the local Lorentz transformations. Thecorresponding conservative Noether current is interpreted asthe angular momentum tensor of the gravity-matter systemincluding the spin density. It is pointed out that,assuming the tetrad transformation given by eq. (15), then torsion does not play a role in theconservation law of angular momentum.