Open access peer-reviewed chapter

Benchmarks for Non-Ideal Magnetohydrodynamics

By Sofronov Vasily, Zhmailo Vadim and Yanilkin Yury

Submitted: April 3rd 2017Reviewed: February 20th 2018Published: July 11th 2018

DOI: 10.5772/intechopen.75713

Downloaded: 201

Abstract

The paper presents an overview of benchmarks for non-ideal magnetohydrodynamics. These benchmarks include dissipative processes in the form of heat conduction, magnetic diffusion, and the Hall effect.

Keywords

  • ALE method
  • numerical simulation
  • magnetohydrodynamics
  • benchmark
  • verification

1. Introduction

Numerical modeling of magnetohydrodynamics (MHD) is an important and challenging problem addressed in numerous publications (e.g., see [1, 2]). This problem is further complicated in case of multi-flux models that account for the relative motion and interaction of particles of different nature (electrons, various species of ions, neutral atoms, and molecules) both with each other and with an external magnetic field.

This class of problems is generally solved using the fractional-step method, when complex operators are represented as a product of operators having a simpler structure. Thus, within the splitting method, the calculation of one-time step consists of a series of simpler procedures. It is obvious that difference schemes for each splitting stage should, where possible, preserve the properties of corresponding difference equations.

Note that the task of constructing reference solutions accounting for the whole range of physical processes is challenging (and often unfeasible). Existing benchmarks enable accuracy assessment of individual splitting stages rather than the simulation as a whole.

Magnetohydrodynamic problems are naturally divided into two groups: problems for an ideal infinitely conducting plasma and problems with dissipative processes in the form of heat conduction and magnetic viscosity.

Numerous publications on the construction of difference methods for ideal magnetohydrodynamics use a standard set of test problems. These include propagation of one-dimensional Alfven waves at various angles to grid lines [3, 4, 5], Riemann problem for MHD equations [6, 7, 8, 9], and various two-dimensional problems accounting for the presence of a uniform magnetic field [3, 5, 10]. In [11], a number of additional ideal MHD benchmarks are presented, which are basically shock-wave problems. A special class of tests includes problems with a weak magnetic field not affecting the medium motion. If there is an exact solution for a given hydrodynamic problem, the magnetic field “freezing-in” principle allows finding components of the field HHxHyHzat any time with the known medium displacements X=XX0t.

The representation in publications of the problem of testing the dissipative stage of MHD equations is much the worse. Possibly, this is owing to complexity problems that require accounting the interaction of the shock-wave processes, heat conduction, diffusion of magnetic field, and Joule heating.

Numerical simulations of some of the tests presented here have been done using the Lagrangian-Eulerian code EGIDA developed at VNIIEF [12, 13] for multi-material compressible flow simulations.

The magnetohydrodynamic equation system in one-temperature approximation modified by the Hall effect can be written in the following conservative form [2]:

ρt+divρu=0,ρudt+divρuu+P+PHI0.5HH=0,PH=0.5H2,∂Ht+divuHHu=rotνrotH+bH×rotH,ρnεt+divρnεu=0,∂Ξt+divΞ+P+PHuHuHκgradT=0,Ξ=ρe+0.5u2+PH,P=PρT,ε=ερT.E1

where ν=c2/4πσis the magnetic viscosity coefficient, κis the heat conduction factor, b=c/4πeneis a local exchange (Hall) parameter [2], and e and ne are charge and density of electrons. When writing Eq. (18), we assume that bias currents and electron inertia are negligibly small [2]. Equation system (1) differs from equation system for ideal MHD owing to diffusion terms present in the equations of energy and inductance of magnetic field.

2. A plane diffusion wave with regard to the Hall effect

Let the components of magnetic field depend on coordinate z only, i.e., H=HxtzHytzHz0. We neglect the medium motion. Then, the magnetic field equation (for components) is written in the following form:

dHxdt=ν2Hxz2+β2Hyz2,dHydt=ν2Hyz2β2Hxz2,dHzdt=0,β=bHz0.E2

Consider the problem of a diffusion wave propagating in an unbounded medium with the given boundary and initial conditions:

Htz=0=H1,Htz=H0,Ht=0z=H0,H0=00Hz0,H1=Hx0Hy0Hz0.E3

Let γ=ν2+β2. A general solution to Eq. (2) for the self-similar variable ξ=z/4γtlooks like

Hx=Hx0+C1Φξ+C2Ψξ,Hy=Hy0+C1ΨξC2Φξ,Hz=Hz0

where Φξ=0ξexpνx2/γsinβx2/γdx, Ψξ=0ξexpνx2/γcosβx2/γdx.

Since Φ=Γ12γν2γ=πγν2γ, Ψ=Γ12γ+ν2γ=πγ+ν2γconstants C1,C2with regard to boundary conditions can be found from equations

C1=Hx0ΦHy0ΨΦ+Ψ=2πHx0γν2γ+Hy0γ+ν2γ,
C2=Hx0Ψ+Hy0ΦΦ+Ψ=2πHx0γ+ν2γHy0γν2γ.

Simulation setup: the initial data is described by Eq. (3). A bounded computational domain 0 < z < L, L = 1 is considered. For this reason, the magnetic field value taken from the analytical solution Hxtz=L=1+C1Φ+C2ΨL4γt, Hytz=L=1+C1ΨC2ΦL4γt, Hztz=L=1is imposed on the right boundary z = L. On the left boundary z = 0, the field components take constant values according to Eq. (3). In simulations with 2D and 3D codes, boundary conditions H/n=0(nis a normal vector to a face) are imposed on lateral faces. By varying parameters ν and β, we can study the effect of the diffusion and Hall terms in Eq. (2) on the diffusion wave parameters. Consider an option with the Hall effect dominating over the diffusion effect: ν = 0, β = 1, Hx0 = Hy0 = Hz0 = 1. Profiles of magnetic field’s components Hx, Hy at time t = 0.01, 0.1 are shown in Figures 1 and 2. With the grid refinement, convergence to the reference solution takes place.

Figure 1.

Profiles of field components at time t = 0.01: (a) Hx and (b) Hy.

Figure 2.

Profiles of field components at time t = 0.1: (a) Hx and (b) Hy.

3. Diffusion of magnetic field in an immovable plane layer of plasma with regard to joule heating and its effect on the diffusion and heat conduction coefficients

The problem of magnetic diffusion in a plane layer of material has many applications in practice [14]. In its detailed formulation, the problem was considered in paper [14] for mega gauss fields. Hydrodynamic motion, magnetic diffusion, heat conduction by electrons, and radiant heat exchange in the “back and forth” approximation were taken into account. Since finding an exact solution to such a problem causes difficulties, the original formulation needs to be simplified. Self-similar solutions to the problem obtained with simplifying assumptions were also presented in [15].

A model problem is considered with the following assumptions:

  • plasma is immovable, it has a constant heat capacity,

  • plasma has Coulomb conductivity,

  • heat conduction is absent.

With such assumptions, the problem is reduced to solving equations

dHdt=rotνrotH, ρdTdt=γ1νrotHrotH,E4

where ν=c2/4πσ,σ=σ0T/Ryα,α=3/2,γ1=R/CV,R=1,CV=1.5

Here, ρ is the density of plasma, γ is the heat capacity ratio (adiabatic index), σ0is the conductivity at T = Ry (it is expressed via atomic constants), and Ry is the Rydberg constant. Energy units have been chosen for temperature.

At initial time t = 0, all quantities depend on one space coordinate. It is assumed that the magnetic field has only one component, H=00Hz. The solution is considered for the problem with initial conditions Eq. (5) and boundary conditions Eq. (6):

Hzxt=0=0ifx<0H0ifx>0,Txt=0=T0,ρxt=0=ρ0E5
Hzxt=0,Hzxt=H0,Tx±t=T0.E6

For dimensionless variables, hz=Hz/H0, τ=T/T0Eq. (4) are reduced to the form:

dhzdt=xντhzx,dt=ηντhzx2,ντ=ν0τα,ν0=c24πσ0T0Ryα,η=γ1H02ρT0.E7

In an infinite region (<x<), the problem has a self-similar solution depending on the variable ξ=x/ν0t. The solution can be obtained by integrating the system of ordinary differential equations:

ξ2dhz+dταdhz=0,ξ2+ηταdhz2=0.E8

with boundary conditions

hzξ=0,hzξ=1,τξ±=1.E9

Note that for a linear case, α = 0, the solution of Eqs. (8), (9) can be found in quadratures

hzξ=0.51+signξerfξ/2,τξ=1η4πEiξ2/4.E10

Since Eix=C+lnx+i1ixiii!, temperature in the vicinity of interface ξ2~0 for the linear case α = 0 has the logarithmic profile τξηlnξ2/4π.

In general, if α > 0, one does not manage to establish the asymptotic law for temperature in the vicinity of ξ2~0, because of no integral curves of Eq. (8) satisfying the boundary conditions at infinity Eq. (9).

Now, let us build the reference solution to the problem with regard to heat conduction. In this case temperature near the interface takes a finite value. The diffusion equations and energy equation of magnetic field with regard to Joule heating and heat conduction are considered. As it was assumed earlier, all quantities depend on one space coordinate, and the magnetic field has only one component, H=00Hz. For dimensionless variables, hz = Hz/H0 and τ = T/T0 equations are reduced to the forms

dhzdt=xν0ταhzx,dt=ην0ταhzx2+xκ0τβτx.E11

A self-similar solution depending on the variable ξ=x/ν0tcan be obtained by integrating the system of ordinary differential equations:

ξ2dhz+dταdhz=0,ξ2+ηταdhz2+adτβ=0,a=κ0D0,E12

with boundary conditions:

hzξ=0,hzξ=1,τξ±=1.E13

To find the reference solution, it is convenient to use the first-order system with an increased number of unknowns instead of the second-order system Eq. (12). The first-order system relative to variables hz,τ,Ψ=ταdhz/,w=aτβ/looks like

dhz=Ψτα,=wτβa,=ξΨτα2,dw=ξwτβ2a+ηΨ2τα.E14

Consider the numerical solution of Eq. (14) for the right half plane (0 < ξ < ∞). The solution in the left half plane <ξ<0follows from the symmetry conditions:

hzξ=1hzξ,τξ=τξ,Ψξ=Ψξ,wξ=wξ.

Consider the numerical solution of Eq. (14) in a bounded domain 0<ξ<ξ1. To formulate boundary conditions for this bounded domain, it is required to find the asymptotic behavior of functions with ξ → ∞. Asymptotic laws can be formulated, if we assume a = 0.5. In this case boundary conditions have the form:

hzξ=0.51+erfξ/2,τξ=1,Ψξ=c1expξ2/4,wξ=c2+ηξc12expξ2/2.E15

Constants C1, C2 are taken so that the following conditions are satisfied on the left boundary of the computational domain:

hzξ=0=0.5,wξ=0=0.

Confine oneself to the consideration of case a = 0.5. Introduce new variables Ψξ=Ψξexpξ2/4, Wξ=wξexpξ2/2with regard to boundary conditions Eq. (15). The replacement of variables gives us the equation system:

dhz=Ψταexpξ24,=Wτβaexpξ22,dΨ=ξΨ1τα2,dW=ξW1τβ2a+ηΨ2ταE16

with boundary conditions:

hzξ=ξ1=0.51+erfξ1/2,τξ=ξ1=1,Ψξ=ξ1=c1,Wξ=ξ1=c2+ηξ1c12,hzξ=0=0.5,wξ=0=0.E17

The set of Eqs. (16), (17) was solved numerically with the methods of automatically selecting an integration step. The following values of parameters were used in simulations: ξ1 = 10, η = 20/3, α = 3/2, D0 = 1, and k0 = aD0 = 1/2. The values of constants satisfying the boundary conditions Eq. (17) were obtained: C1 = 0.10231 and C2 = 1.79474. Since the behavior of functions near the right boundary corresponds to asymptotic laws Eq. (17), the reduction of parameter ξ1 from ξ1 = 10 to ξ1 = 1 does not affect simulation results.

Results of simulations are illustrated in Figures 3 and 4. With the use of such regularity method (with regard to heat conduction), temperature at the central point of the computational domain takes its finite value. Note that with t = 1/ν0 the space coordinate coincides with the self-similar coordinate, x = ξ.

Figure 3.

Profiles of self-similar functions: (a) W and (b) Ψ.

Figure 4.

Profiles: (a) magnetic field and (b) temperature.

4. A point explosion in a perfectly non-conducting atmosphere

Let us consider the problem of a point blast in the presence of a uniform magnetic field (for definiteness) along the z axis (Hz = Hz0 = 0.01) in a perfectly non-conducting atmosphere. Initial data are chosen in such a way that the field has no effect on the motion of matter ε0r2/r13>>Hz02/ρ0. It is assumed that behind the shock front, the medium becomes perfectly conducting. If a self-similar solution to the problem of a point explosion is known, then one can calculate the magnetic field components at some time t>0. In an external domain r>rFt, the magnetic field’s vector potential is the solution to stationary equation rotrotΨ=0. With regard to conditions at infinity, this solution takes the form [16]:

Ψrrϑt=0,Ψφrϑt=0.5Hz0r1Ct/r3sinϑ,Ψϑrϑt=0.E18

Here, unknown constant Ct=brF3tcan be found from the condition of coupling with the solution in an internal domain r<rFt. Write components of magnetic field H=rotΨ:

Hrrϑt=Hz01Ct/r3sinϑ,Hφrϑt=0,Hϑrϑt=Hz01+0.5Ct/r3cosϑ.

The solution in the internal domain (r < rF(t)) is found from the freezing-in condition of the magnetic field:

ddtHrρ=Hrρurr,ddtHϑρ=Hϑρurr,ddtHφρ=Hφρurr.

The integration of these equations with regard to the solution in external domain Eq. (18) gives us

Hrrϑt=Hz0hrrtcosϑ,Hϑrϑt=Hz0hϑrtsinϑ,Hφrϑt=0,E19

where

hrrt=1βrFtr3,r>rFt1βr0rtr2,rrFt,hϑrt=1β2rFtr3,r>rFt1+β2ρrtrγ1ρ0r0rtγ+1,rrFt.

Here, the functions r0rt, ρrtare defined from the self-similar solution to the point blast problem [17].

Unknown constant βcan be found from the condition of the solenoidal distribution of magnetic field in internal domain.

divH=Hz0cosϑ1r2r2hrr+2hϑr=2Hz0r0cosϑ1βr02r0r2r1+β2ρρ0γ1γ+1=0.

Since ρρ0=r02r0r2r, then 1β=1+β2γ1γ+1. That is why β=43γ+1.

Note that in external domain this condition is satisfied automatically. In Cartesian coordinates, the solution of Eq. (19) looks like

Hxxyzt=Hz0xzr2hrrt+hϑrt,Hyxyzt=Hz0yzr2hrrt+hϑrt,Hzxyzt=Hz0z2r2hrrt+hϑrthϑrt.E20

It is convenient to compare the numerical and exact solutions using the field components depending on one space coordinate:

hrrt=HrHz0cosϑ=1Hz0Hzxyzt+xHxxyzt+yHyxyztz,hϑrt=HϑHz0sinϑ=1Hz0HzxyztzxHxxyzt+yHyxyztr2z2hφrt=HφHz0sinϑ=1Hz0yHxxyzt+xHyxyztr=0.,E21

The magnetic field lines can be obtained by integrating equations

dxdz=xzhrrt+hϑrtr2hϑrt+z2hrrt+hϑrt,dydz=yzhrrt+hϑrtr2hϑrt+z2hrrt+hϑrt.E22

We consider the process stage, at which the numerical simulation becomes self-similar. In this case, the shock wave is considerably far (compared to the energy release region) from the blast center. For example, at the final time t = 3, the wave front is located at a distance of RF = 13.467.

The flow parameters in this problem depend on a single spatial variable, r, and the field components, on two variables, r and θ. One can restrict the consideration to any plane passing through the z axis. Magnetic field lines in the plane y = 0 at t = 3 are shown in Figure 5. It follows from this figure that these lines of force stretch along axis z and, thereby, prevent the spread of gas in the direction orthogonal to this axis. This effect is small in the given problem because of the field smallness. With an increased strength of the field, the pressure zone gets out of its spherical shape due to occurrence of the singled out direction.

Figure 5.

The magnetic field lines at time t = 3 in plane y = 0. A dashed line shows the shock front. Strong blast in a perfectly nonconducting atmosphere (a) and in a uniform conducting atmosphere (b) [11].

Profiles of the field components Hy and Hz along the line x = z in this plane are presented in Figure 6. The self-similar functions of fluid parameters and radial and angular field components hr and hθ depend on one spatial variable r. These functions are shown in Figure 7. In testing numerical methods, the values of grid functions for all cells in the domain can be mapped onto the figures. Such a comparison of reference versus numerical solution indicates whether the spherical symmetry is preserved during numerical simulations.

Figure 6.

Profiles of field components along line x = z and y = 0 at time t = 3: (a) Hx and (b) Hz.

Figure 7.

Profiles of the non-dimensionalized field components at time t = 3: (a) radial hr and (b) angular hθ.

Simulation setup: The energy release region is a sphere of radius r1 = 0.1, in which the initial specific energy per unit mass is set to ε = ε0 = 107 and the initial density is set to ρ = ρ0 = 1. In the spherical layer r1 < R < r2 = 15, the specific energy and the density are equal to ε = 0, ρ = ρ0 = 1, respectively. The EOS is P = ρ(γ−1)ε, γ = 1.4. The problem domain in the three-dimensional setup is a cube L × L× L. All boundary faces of the domain are rigid walls.

This problem requires taking into account the magnetic field diffusion in external domain (outside the shock front). It is assumed that behind the shock front, the medium becomes perfectly conducting due to ionization effects. The magnetic viscosity is approximated by the following dependence:

νε=c24πσ=ν1=105,ε0,ν11ε/ε1,0<ε<ε1=1,0,ε1ε.

Parameter ε1is chosen to provide that the magnetic viscosity behind the shock front is always zero, i.e., the condition ε1<εFtε0r1/rFt3,ttk=3is satisfied. Accounting of diffusion in external domain leads to the necessity of increasing the size of computational domain (rF(tk)/L))3 < <1) in comparison with the ideal MHD problems to be able to set boundary conditions corresponding to the initial undisturbed state. Thus, the size of computational domain must be almost five times larger, L75.

5. Diffusion of magnetic field into a spherical plasma cloud

The problem formulation and its analytical solution have been taken from [16]. In contrast to this paper, consider the diffusion problem (the plasma cloud motion is neglected):

∂Ht=rotνrotH+bH×rotH.E23

It is assumed that the magnetic field at infinity is uniform and directed along axis z: H00H0=2(see Figure 8). The magnetic viscosity coefficient is constant inside and outside the cloud: νr=ν1,r<r0=1ν2,r>r0.

Figure 8.

The problem of diffusion into a plasma cloud [16].

5.1. Diffusion of magnetic field in the absence of the Hall effect

Assume that the Hall effect contribution is small, bH0/ν1. Write the equation of diffusion relative to vector potential H=rotΨ:

∂Ψt=νrotrotΨ.E24

This is an axially symmetric problem, and, therefore, it is convenient to use the polar coordinate system rϑφ. With no Hall effect and with regard to the conditions at infinity, the initial data for vector potential takes the form [16]

Ψrrϑt=0=0,Ψϑrϑt=0=0,Ψφrϑt=0=rfrt=0sinϑ,E25
frt=0=H020r<r01r0/r3r>r0.E26

The solution of Eq. (24) with initial data Eq. (25) is reduced to solving Eq. (27) with initial data Eq. (26) and boundary conditions Eq. (28):

dfdt=νrr4rr4fr,E27
frr=0t=0,frrt=0.E28

The solution has the form

Ψrrϑt=0=0,Ψϑrϑt=0,Ψφrϑt=rfrtsinϑ.E29

Paper [16] considers the plasma cloud interaction with the magnetic field of vacuum, and, therefore, ν2is assumed. For this special case, the solution to Eq. (27) in quadratures has been obtained, and it has the forms:

frt=H0216ζ2n=11nπn2Tntcosπnζsinπnζπnζ,0<ζ=rr0<116ζ2n=11πn2Tnt,1<ζ,Tnt=expπn2ν1tr02.E30

For finite values of conductivity in external domain ν2>0, the limit numerical solution of Eq. (27) has been taken for the reference solution.

Components of magnetic field H=rotΨare found by differentiating vector potential Eq. (29). If function frtis known, these components are calculated using formulas.

Hrrϑt=hrrtcosϑ,Hϑrϑt=hϑrtsinϑ,Hφrϑt=hφsinϑ,hrrt=2frt,hϑrt=r2f/rr,hφrt=0E31

In Cartesian coordinates the field components have the forms

Hzxyzt=hϑrt+z2r2hrrthϑrt,Hyxyzt=zyr2hrrthϑrt,Hxxyzt=xzr2hrrthϑrt,E32

Since the problem is axially symmetric, any plane coming across axis z can be taken to calculate the magnetic field lines. For example, for plane y=0, the differential equation describing the slope of the magnetic field lines looks like

dxdz=xzhrhϑr2hϑ+z2hrhϑ.

The magnetic field lines for the reference solution at time t = 0.01 are shown in Figure 9.

Figure 9.

The magnetic field lines in the exact solution with parameters ν1 = 1 and ν2 → ∞. A dashed line shows the plasma cloud position.

Results of Eq. (32) are the formulas for the radial and angular components of the field depending on a single space coordinate:

hrrt=Hzxyzt+xHxxyzt+yHyxyztz,hϑrt=HzxyztzxHxxyzt+yHyxyztr2z2,hφrt=xHyxyztyHxxyztr=0.

Figures 10 and 11 show profiles of field components for different magnetic viscosity values in external domain (r > r0). Note that there is a small difference between the profiles obtained with ν2 = 50 and ν2 → ∞ corresponding to the simulation of the plasma cloud interaction with the magnetic field of vacuum ν2 → ∞. For the problem with ν2 → ∞, profiles of the non-dimensionalized field components for the initial phase of diffusion, t < r03/ν1, are given (see Figure 12).

Figure 10.

Profiles of the non-dimensionalized field components at time t = 0.1: (a) radial hr and (b) angular hθ.

Figure 11.

Profiles of the field components along line x = z and y = 0 at time t = 0.1: (a) Hx and (b) Hz.

Figure 12.

Evolution of non-dimensionalized field components in the problem with parameters ν1 = 1 and ν2 → ∞: (a) radial hr and (b) angular hθ.

Simulation setup: A computational domain x0.5Ly0.5Lz0.5Lis a cube with edges L = 10. Boundary conditions corresponding to the initial undisturbed state are imposed on its lateral faces for the components of field HrtrΓ=Hrt=0. The initial data can be set either for the magnetic field components Eq. (33) or the components of vector potential Eq. (34).

Hzxyzt=0=2frt=0+z21+rfrt=0/r2,Hyxyzt=0=xyfrt=0/r,Hxxyzt=0=xzfrt=0/r,E33
Ψzxyzt=0=0,Ψyxyzt=0=zfrt=0,Ψxxyzt=0=yfrt=0.E34

The EGIDA code uses a scheme preserving the field divergence at one step because difference operators DIV and ROT (div and rot) [18] determined at nodes and in cells of a grid, respectively, satisfy the vector analysis identities: DIVrot = 0 and ROTdiv = 0.

It has been found that for the first set of initial data the divergence norm depends on errors induced by the initial distribution of the Hfield components in the vicinity of sphere r = r0. Though the difference scheme does not change the magnetic field divergence, initial errors lead to a significantly distorted numerical solution (see Figure 13).

Figure 13.

Calculation of the magnetic field diffusion into a spherical plasma cloud for the first set of initial data (33). Profiles of the field components at time t = 0.1: (a) radial hr and (b) angular hθ. Distribution of transverse field component Hy in section z = 0 (c).

In the second case, the magnetic field components are determined using the operator numerically differentiating the vector potential, and, hence, the magnetic field divergence norm equals zero at initial time and at all later times. For this case, a good agreement between the calculated results and the exact solution has been achieved even on the coarsest grid (see Figure 14).

Figure 14.

Calculation of the magnetic field diffusion into a spherical plasma cloud for the second set of initial data (34). Profiles of the field components at time t = 0.1: (a) radial hr and (b) angular hθ. Distribution of transverse field component Hy in section z = 0 (c).

5.2. Diffusion of magnetic field with a low Hall effect

Assume that the Hall effect contribution is small, bH0/ν1, but finite. For this reason, it is required to take into account the Hall term in the diffusion equation:

∂Ψt=νrotrotΨbrotΨ×rotrotΨ.

The Hall effect leads to the occurrence of the azimuthal component, Hφ, of magnetic field in plasma [16]:

Hrrϑt=hrrtcosϑ,Hϑrϑt=hϑrtsinϑ,Hφrϑt=hφsinϑ,hrrt=2frt,hϑrt=r2f/rr,hφrt=Ψrt/r.

In view of the smallness of parameter bH0/ν1, we have Eqs. (27), (28) to calculate function frt, while the calculation of small additive Ψis reduced to solving the following boundary value problem:

dt=rνΨr6νΨr22fr2rbr4rr4fr,Ψrt=0=0,Ψr=0t=0,Ψrt=0.

Figure 15 shows profiles of the non-dimensionalized azimuthal field component hφat early times, t<r02/ν1. With small values of parameter bH0/ν1, the rest two components—hr,hϑ—remain unchanged, and they are shown in Figure 12. Note that, as it has been shown in [16], if the motion of plasma is accounted, the Hall effect may lead to the occurrence of azimuthal velocity, i.e., to the plasma cloud rotation.

The changeover to Cartesian coordinates is performed using formulas

Hzxyzt=hϑrt+z2hrrthϑrt/r2,Hyxyzt=zyhrrthϑrt/r2+xhφrt/r,Hxxyzt=xzhrrthϑrt/r2yhφrt/r.

Simulation setup: The initial data is given in the previous section. This problem requires accounting the Hall effect. A local exchange parameter is b = 0.01. Figure 15 shows profiles of the azimuthal component of non-dimensionalized field.

Figure 15.

Profiles of the azimuthal component of non-dimensionalized field in the problem with parameters ν1=1,ν2→∞,b=0.01.

6. Conclusion

An important property of difference schemes in multidimensional flow simulations is that they keep the magnetic field divergence-free in difference solutions. An adverse aspect of this defect is the unphysical transport of matter orthogonal to the field H[2].

Note that the zero-divergence requirements to difference schemes get more stringent as applied to the solution of diffusion problems. A violation of this requirement results in the accumulation of errors and loss of solution structure, especially in problems with high conductivity gradients.

Acknowledgments

Some EGIDA calculations have been performed under contract no. 1239349 between Sandia National Laboratories and RFNC-VNIIEF. The authors thank sincerely M. Pokoleva for her assistance in simulations and J. Kamm and A. Robinson for their assistance in formulating some benchmarks and references.

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Sofronov Vasily, Zhmailo Vadim and Yanilkin Yury (July 11th 2018). Benchmarks for Non-Ideal Magnetohydrodynamics, Numerical Simulations in Engineering and Science, Srinivas P. Rao, IntechOpen, DOI: 10.5772/intechopen.75713. Available from:

chapter statistics

201total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

A Numerical Simulation of the Shallow Water Flow on a Complex Topography

By Alexander Khoperskov and Sergey Khrapov

Related Book

First chapter

Monte Carlo Simulations in NDT

By Frank Sukowski and Norman Uhlmann

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us