Next Article in Journal
Exploring the Impact of Technology 4.0 Driven Practice on Warehousing Performance: A Hybrid Approach
Next Article in Special Issue
Accurate Solutions to Non-Linear PDEs Underlying a Propulsion of Catalytic Microswimmers
Previous Article in Journal
Hamiltonian Cycles in Cayley Graphs of Gyrogroups
Previous Article in Special Issue
First-Order Approximate Mei Symmetries and Invariants of the Lagrangian
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Invariant Finite-Difference Schemes for Plane One-Dimensional MHD Flows That Preserve Conservation Laws

by
Vladimir Dorodnitsyn
1,* and
Evgeniy Kaptsov
2
1
Keldysh Institute of Applied Mathematics, Russian Academy of Science, Miusskaya Pl. 4, Moscow 125047, Russia
2
School of Mathematics, Institute of Science, Suranaree University of Technology, Nakhon Ratchasima 30000, Thailand
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(8), 1250; https://doi.org/10.3390/math10081250
Submission received: 1 March 2022 / Revised: 7 April 2022 / Accepted: 8 April 2022 / Published: 11 April 2022

Abstract

:
Invariant finite-difference schemes are considered for one-dimensional magnetohydrodynamics (MHD) equations in mass Lagrangian coordinates for the cases of finite and infinite conductivity. The construction of these schemes makes use of results of the group classification of MHD equations previously obtained by the authors. On the basis of the classical Samarskiy–Popov scheme, new schemes are constructed for the case of finite conductivity. These schemes admit all symmetries of the original differential model and have difference analogues of all of its local differential conservation laws. New, previously unknown, conservation laws are found using symmetries and direct calculations. In the case of infinite conductivity, conservative invariant schemes are constructed as well. For isentropic flows of a polytropic gas the proposed schemes possess the conservation law of energy and preserve entropy on two time layers. This is achieved by means of specially selected approximations for the equation of state of a polytropic gas. In addition, invariant difference schemes with additional conservation laws are proposed. A new scheme for the case of finite conductivity is tested numerically for various boundary conditions, which shows accurate preservation of difference conservation laws.

1. Introduction

Magnetic hydrodynamics equations describe the flows of electrically conductive fluids such as plasma, liquid metals, and electrolytes and are widely used in modeling processes in various fields from engineering to geophysics and astrophysics.
In the present publication, we restrict ourselves to considering plane one-dimensional MHD flows under the assumption that the medium is inviscid and thermally non-conducting. A group classification of the MHD equations under the above conditions was carried out recently in [1] (for some particular results see also [2,3,4,5]). The group classification splits into four essentially different cases according to whether the conductivity of the medium is finite or infinite, and the longitudinal component of the magnetic field vector is zero or a non-zero constant.
The MHD equations are nonlinear, so that even in the one-dimensional case only their particular solutions are known [6,7,8,9,10]. Therefore, numerical modeling in magnetohydrodynamics is of great practical interest. There are many approaches to numerically modeling MHD equations, including finite-difference, finite element, and finite volume methods (see, e.g., [11,12,13,14,15,16,17,18,19]). Further we consider finite-difference schemes taking as a starting point the classical Samarsky–Popov schemes [12,13] for the MHD equations for the case of finite conductivity. The main properties of the considered schemes are invariance, i.e., preservation of the symmetries of the original differential equations, and the presence of difference analogues of local differential conservation laws. It is known that there is a connection between the invariance of equations and the presence of conservation laws [20,21,22,23]. Recall that local conservation laws in the form of a divergence of some vector make it possible to calculate flows through the boundary of a region using the Gauss–Ostrogradsky theorem [24]. The latter has a clear physical meaning—the vector flux (electromagnetic field, etc.) through the computational domain is preserved if there are no sources inside the domain. The presence of a finite-difference analog of such a local conservation law also makes it possible to sum it over the entire computational domain and preserve the vector flow through the boundary of a domain. Thus, such a difference scheme retains the geometric and physical meaning inherent in the differential model.
Invariant schemes have been studied for a long time [25,26,27,28], and over the past decades, significant progress has been made in the development of methods for their construction and integration. For schemes for ordinary differential equations with Lagrangian or Hamiltonian functions, a number of methods [29] have been developed that make it possible to decrease the order or even integrate the schemes. A method based on the Lagrangian identity has also been developed for the case when the equations do not admit a variational formulation [30].
For partial difference schemes, the main methods used are the method of differential invariants [28,31,32] and the difference analogue of the direct method [28]. Using these methods, the authors have constructed invariant schemes for various shallow water models [33,34,35,36]. In addition, some previously known schemes have been investigated from a group analysis point of view. In particular, symmetries and conservation laws of the Samarskiy–Popov schemes for the one-dimensional gas dynamics equations of a polytropic gas have been investigated in [37,38,39]. Based on the results of the group classification [1] and Samarskiy–Popov schemes for the MHD equations, we further construct invariant finite-difference schemes possessing conservation laws. The set and number of conservation laws depend on the conductivity, the form of the magnetic field vector, and the equation of state of the medium.
This paper is organized as follows. In Section 2, the simplest version of the finite-conductivity MHD equations in mass Lagrangian coordinates in the case of one-dimensional plane flows is considered. Electric and magnetic fields are represented by one-component vectors, which greatly simplifies the form of the equations. This was the main case considered in Samarsky and Popov’s publications [12,13]. The section also provides basic notation and definitions. Then, symmetries and conservation laws of the Samarsky–Popov scheme for the MHD equations are investigated. In addition to the previously known conservation laws, the center-of-mass conservation law is given, as well as new conservation laws for the specific conductivity function, obtained on the basis of the group classification [1].
Section 3 is devoted to various generalizations of the scheme of Section 2. In Section 3.1, the scheme for arbitrary electric and magnetic fields is considered. Its symmetries are investigated and conservation laws are given. The case of infinite conductivity is considered in Section 3.2. It is shown that in this case the Samarsky–Popov scheme requires some additional modifications in order to possess the conservation law of angular momentum. In the case of a polytropic gas, it turns out to be possible to preserve not only energy, but also entropy along pathlines. This can be done using a specially selected equation of state for a polytropic gas. At the end of the section, an example of an invariant scheme is given that does not possess a conservation law of energy, but preserves entropy and has additional conservation laws in the case of isentropic flows. In Section 4, one of the invariant schemes for the case of finite conductivity is numerically implemented for the example of plasma bunch deceleration by crossed electromagnetic fields. The results are discussed in the Conclusions.

2. Conservative Schemes for MHD Equations with Finite Conductivity

Problems of continuum mechanics and plasma physics are often considered in mass Lagrangian coordinates [13,40] since for them the formulation of boundary conditions is greatly simplified. In particular, conservative Samarsky–Popov schemes for the equations of gas dynamics and magnetohydrodynamics have been constructed in mass Lagrangian coordinates.
In mass Lagrangian coordinates the MHD equations, describing the plane one-dimensional MHD flows, are [1,13]
1 ρ t = u s ,
u t = p + ( H y ) 2 + ( H z ) 2 2 s ,
v t = H 0 H s y , y t = v ,
w t = H 0 H s z , z t = w ,
H y ρ t = ( H 0 v + E z ) s ,
H z ρ t = ( H 0 w E y ) s ,
σ E y = ρ H s z , σ E z = ρ H s y ,
ε t = p u s + 1 ρ ( i · E ) ,
x t = u , x s = 1 / ρ ,
where t is time, s is the mass Lagrangian coordinate, x is the Eulerian coordinate, ρ is density, p is pressure, ε is internal energy, u = ( u , v , w ) is the velocity of a particle, E = ( E x , E y , E z ) is the electric field vector, H = ( H x , H y , H z ) is the magnetic field vector, and i = ( i x , i y , i z ) is the electric current. The conductivity σ is some function of p and ρ , i.e., σ = σ ( p , ρ ) .
Following [13] we first consider the simplest case of one-component electric and magnetic fields. Here, we also introduce the notation and some basic concepts. In the next section some generalizations are considered, including the case of infinite conductivity.
For simplicity, the longitudinal component of the magnetic field H is set to zero, and the coordinate system is chosen in such a way that H = ( 0 , H , 0 ) . Consequently, the electric current i and the electric field E are also one-component vectors, i.e., i = ( 0 , 0 , i ) , E = ( 0 , 0 , E ) . Electromagnetic force f = ( f , 0 , 0 ) acts in the x-direction, and the velocity is u = ( u , 0 , 0 ) (see Figure 1).
Given the above, the system of the one-dimensional MHD equations with a finite conductivity | σ | < in mass Lagrangian coordinates can be written as [13]
1 ρ t = u s ,
u t = p s + f , f = i H / ρ ,
H ρ t = E s ,
i = σ E = κ ρ H s ,
ε t = p u s + q , q = i E / ρ ,
x t = u , x s = 1 / ρ ,
where κ = 1 / ( 4 π ) and q is Joule heating per unit mass.
In particular, we consider a polytropic gas for which the following relation holds
ε = 1 γ 1 p ρ , γ = const > 1 .
Equation (2e) for the energy evolution can be rewritten in the semi-divergent form
ε + u 2 2 t = ( p u ) s + f u + q ,
or in the divergent form
ε + u 2 2 + κ H 2 2 ρ t = p + κ H 2 2 u s + κ E H s .
Note that the electromagnetic force f = i H / ρ can be represented in the divergent form f = κ ( H 2 / 2 ) s , and Equation (2b) can be rewritten as
u t = p + κ H 2 2 s .
Further, we assume κ = 1 since it can be discarded by means of the scaling transformation
s ˜ = κ s , p ˜ = κ p , ρ ˜ = κ ρ , σ ˜ = κ σ .

2.1. Conservative Samarskiy–Popov Schemes for System (2)

The family of Samarskiy–Popov conservative difference schemes for system (2) is
1 ρ t = u s ( 0.5 ) ,
u t = p s ¯ ( α ) + f , f = H H ^ 2 s ¯ = 1 2 [ i ^ H / ρ ^ + i H ^ / ρ ] ,
H ρ t = E s ( β ) ,
i = σ E = ρ H s ¯ ,
ε t = p ( α ) u s ( 0.5 ) + q , q = 1 2 [ ( i / ρ ) ( 0.5 ) E ( β ) + ( i + / ( ρ + ) ) ( 0.5 ) E + ( β ) ] ,
x t = u ( 0.5 ) , x s = 1 ρ ,
where 0 ( α , β ) 1 are free parameters. For arbitrary α and β , scheme (8) approximates system (2) up to O ( τ + h 2 ) , and for α = β = 0.5 , the scheme is of order O ( τ 2 + h 2 ) .
Here and further ϕ t , ϕ t ˇ and ϕ s , ϕ s ¯ denote finite-difference derivatives of some quantity ϕ = ϕ ( t n , s m , u m n , )
ϕ t = S + τ ( ϕ ) ϕ τ n , ϕ s = S + s ( ϕ ) ϕ h m , ϕ t ˇ = ϕ S τ ( ϕ ) τ n 1 , ϕ s ¯ = ϕ S s ( ϕ ) h m 1 ,
which are defined with the help of the finite-difference right and left shifts along the time and space axes correspondingly
S ± τ ( ϕ ( t n , s m , u m n , ) ) = ϕ ( t n ± 1 , s m , u m n ± 1 , ) , S ± s ( ϕ ( t n , s m , u m n , ) ) = ϕ ( t n , s m ± 1 , u m ± 1 n , ) .
The indices n and m are, respectively, changed along time and space axes t and s. The time and space steps are defined as follows
τ n = τ ^ = t n + 1 t n = t ^ t , τ n 1 = τ ˇ = t n t n 1 = t t ˇ , h m = h + = s m + 1 s m = s + s , h m 1 = h = s m s m 1 = s s .
Following the Samarskiy–Popov notation throughout the text we denote
S + s ( ϕ ) = ϕ + , S s ( ϕ ) = ϕ , S + τ ( ϕ ) = ϕ ^ , S τ ( ϕ ) = ϕ ˇ ,
ϕ ( α ) = α ϕ ^ + ( 1 α ) ϕ
and
ϕ = ( ϕ ) i j = h i ϕ i 1 / 2 j + h i 1 ϕ i + 1 / 2 j h i + h i 1 .
Note that on a uniform lattice h i = h = const in its integral nodes (13) becomes
ϕ = ϕ + ϕ 2 .
Remark 1.
The energy Equation (8e) can be reduced to one of the three following forms [13] using equivalent algebraic transformations:
ε t = p ( α ) u s ( 0.5 ) + q ,
ε + u 2 + u + 2 4 t = p ( α ) u ( 0.5 ) s + 1 2 f u ( 0.5 ) + f + u + ( 0.5 ) + q ,
ε + u 2 + u + 2 4 + H 2 2 ρ t + p ( α ) + ( H H ^ ) 2 u ( 0.5 ) E ( β ) H ( 0.5 ) s = 0 .
These different forms of equation reflect the balance of certain types of energy, i.e., they express the different physical aspects of energy conservation. To emphasize this property, such schemes are also called completely conservative.

2.2. Invariance of Samarskiy–Popov Schemes

System (2) can be rewritten in the following form that is more convenient for symmetry analysis
1 ρ t = u s ,
u t = p + H 2 2 s ,
H ρ t = E s ,
σ E = ρ H s ,
p t = γ ρ p u s + ( γ 1 ) σ E 2 ,
x t = u , x s = 1 / ρ .
Remark 2.
Note that for the polytropic gas with the state Equation (3), one can rewrite the energy evolution Equation (8e) of the Samarskiy–Popov scheme as
p t = ρ ^ ( p + ( γ 1 ) p ( α ) ) u ( 0.5 ) + ( γ 1 ) ρ ^ 2 σ E ρ ( 0.5 ) E ( β ) + ( i + ) E + ( ρ + ) ( 0.5 ) E + ( β ) .
In this form the energy evolution equation corresponds to Equation (18e).
Calculations show [1] that the Lie algebra admitted by the system for an arbitrary σ = σ ( p , ρ ) is as follows (here and further the notation f f is used):
X 1 = t , X 2 = s , X 3 = x , X 4 = t x + u .
The group generator
X = ξ t t + ξ s s + η x
is prolonged to the finite-difference space as follows [25,28]
X ˜ = k , l = S k τ S l s ( X ) ,
The scheme of the form
Φ ( t , t ˇ , t ^ , s , s + , s , u , u + , u , u ^ , u ˇ , u ^ + , u ˇ + , u ^ , u ˇ , ) = 0 ,
h + = h = h , τ ^ = τ ˇ = τ , ( τ , h ) = 0
defined on a uniform orthogonal mesh is invariant if the following criterion of invariance holds [28]
X ˜ Φ | ( 23 ) = 0 , X ˜ ( τ ^ τ ˇ ) | ( 23 ) = 0 , X ˜ ( h + h ) | ( 23 ) = 0 .
To preserve uniformness and orthogonality of the mesh it is also required that [25,28]
D + s D s ( ξ s ) = 0 , D + τ D τ ( ξ t ) = 0 ,
D ± s ( ξ t ) = D ± τ ( ξ s ) ,
where D ± τ and D ± s are finite-difference differentiation operators
D + τ = S + τ 1 τ n , D τ = 1 S τ τ n 1 , D + s = S + s 1 h m , D s = 1 S s h m 1 .
One can verify that scheme (8) is indeed invariant with respect to the generators X 1 , …, X 4 and all the generators (20) satisfy the mesh orthogonality and uniformness conditions (25) and (26). Hence, one can use an orthogonal uniform mesh (23b) that is an invariant one.

2.3. Conservation Laws Possessed by the Samarskiy–Popov Scheme

All the conservation laws of system (18) have their finite-difference counterparts for the Samarskiy–Popov scheme. They are given in Table 1. For convenience, the conservation laws are numbered (see column “#”).
An additional conservation law
s H ρ t + ( H s E ) s = 0 .
only occurs in case σ = ρ [1]. In this case, system (18) admits two more symmetries, namely
X 5 = s s x x u u + 2 ρ ρ E E , X 6 = 2 t t + 2 s s 2 u u + 2 p p + 2 ρ ρ 3 E E H H .
Conservation law (27) has its finite-difference counterpart, which can be found by direct calculations.
Note that neither conservation law (27) nor the center-of-mass law
t u x t + t p + H 2 2 s = 0
was mentioned in [13]. Perhaps the authors of [13] knew the finite-difference analogue of (29).

3. Generalizations of the Samarskiy–Popov Schemes for MHD Equations

3.1. The Case of Finite Conductivity

We consider a more general case H = ( H 0 , H y , H z ) , E = ( 0 , E y , E z ) , i = ( 0 , i y , i z ) , u = ( u , v , w ) , x = ( x , y , z ) , and H 0 = const . Here we used the fact that the coordinate system can always be chosen in such a way that the first component of the vector field E is equal to zero.
Further, we consider Equation (1), where, by analogy with (18), the energy evolution Equation (1h) is written as
p t = γ p ρ u s + ( γ 1 ) σ ( ( E y ) 2 + ( E z ) 2 ) ,
A generalization of scheme (8) for E = ( 0 , E y , 0 ) and H = ( H 0 , 0 , H z ) is given in [13]. Since the MHD equations are almost symmetric in terms of the components E y , E z and H y , H z , one can extend the scheme proposed in [13] as follows
1 ρ t = u s ( 0.5 ) ,
u t = p ( α ) + H y H ^ y + H z H ^ z 2 s ¯ , v t = H 0 ( H y ) s ¯ ( 0.5 ) , w t = H 0 ( H z ) s ¯ ( 0.5 ) ,
H y ρ t = H 0 v s ( 0.5 ) + ( E z ) s ( β 1 ) , H z ρ t = H 0 w s ( 0.5 ) ( E y ) s ( β 2 ) ,
i y = σ E y = ρ H s ¯ z , i z = σ E z = ρ H s ¯ y ,
ε t = p ( α ) u s ( 0.5 ) + q y + q z ,
x t = u ( 0.5 ) , y t = v ( 0.5 ) , z t = w ( 0.5 ) , x s = 1 ρ ,
where 0 ( α , β 1 , β 2 ) 1 and
q y = 1 2 [ ( i y / ρ ) ( 0.5 ) ( E y ) ( β 2 ) + ( i + y / ( ρ + ) ) ( 0.5 ) ( E y ) + ( β 2 ) ] ,
q z = 1 2 [ ( i z / ρ ) ( 0.5 ) ( E z ) ( β 1 ) + ( i + z / ( ρ + ) ) ( 0.5 ) ( E z ) + ( β 1 ) ] .
Note that this generalization of the scheme was discussed in [13] but it was not given explicitly.
Remark 3.
One can generalize (19) for scheme (31), (3) as follows
p t + ρ ^ u s ( 0.5 ) ( p + ( γ 1 ) p ( α ) ) + ( 1 γ ) ρ ^ ( q y + q z ) = 0 .
According to [1], the symmetries admitted by system (1) are the following.
  • If H 0 = 0 and σ is arbitrary, then the admitted Lie algebra is
    X 1 = t , X 2 = s , X 3 = x , X 4 = t x + u , X 5 = E z E y E y E z + H z H y H y H z .
    In case σ = ρ , there are two additional generators that are admitted, namely
    X 1 a = s s x x u u + 2 ρ ρ E y E y E z E z , X 2 a = 2 t t + 2 x x 2 p p 2 ρ ρ E y E y E z E z H y H y H z H z .
    There are also two more conservation laws in the latter case (see Table 2). Additional conservation laws do not occur for any other forms of the function σ .
  • If H 0 0 and σ is arbitrary then the admitted Lie algebra is
    X 1 = t , X 2 = s , X 3 = x , X 4 = t x + u , X 5 = E z E y E y E z + H z H y H y H z + w v v w + z y y z , X 6 = t y + v , X 7 = t z + w , X 8 = h 1 ( s ) y , X 9 = h 2 ( s ) z ,
    where h 1 and h 2 are arbitrary functions of s.
    Additional conservation laws do not occur for any specific σ .
Table 2. Differential and difference conservation laws for the extended scheme.
Table 2. Differential and difference conservation laws for the extended scheme.
#Conservation Laws of the SystemConservation Laws of the SchemePhysics Interpretation
σ = σ ( p , ρ )
1 1 ρ t u s = 0 1 ρ t u ( 0.5 ) s = 0 Mass conservation
2 H y ρ t ( E z + H 0 v ) s = 0 H y ρ t ( E z ) ( β 1 ) + H 0 v ( 0.5 ) s = 0 Magnetic flux conservation
3 H z ρ t + ( E y H 0 w ) s = 0 H z ρ t + ( E y ) ( β 2 ) H 0 w ( 0.5 ) s = 0 Magnetic flux conservation
4 u t + p + ( H y ) 2 + ( H z ) 2 2 s = 0 u t + p ( α ) + H y H ^ y + H z H ^ z 2 s = 0 Momentum conservation
5 v t ( H 0 H y ) s = 0 v t H 0 ( H y ) ( 0.5 ) s = 0 Momentum conservation
6 w t ( H 0 H z ) s = 0 w t H 0 ( H z ) ( 0.5 ) s = 0 Momentum conservation
7 t u x t + t p + ( H y ) 2 + ( H z ) 2 2 s = 0 t + t ˇ 2 u x t + t p ( α ) + H y H ^ y + H z H ^ z 2 s = 0 Center of mass law
8 ( t v y ) t ( t H 0 H y ) s = 0 t + t ˇ 2 v y t t H 0 ( H y ) ( 0.5 ) s = 0 Center of mass law
9 ( t w z ) t ( t H 0 H z ) s = 0 t + t ˇ 2 w z t t H 0 ( H z ) ( 0.5 ) s = 0 Center of mass law
10 ( ε + u 2 + v 2 + w 2 2 + ( H y ) 2 + ( H z ) 2 2 ρ ) t + [ ( p + ( H y ) 2 + ( H z ) 2 2 ) u + E y H z E z H y H 0 ( v H y + w H z ) ] s = 0 ( ε + u 2 + u + 2 + v 2 + v + 2 + w 2 + w + 2 4 + ( H y ) 2 + ( H z ) 2 2 ρ ) t + [ ( p ( α ) + ( H y H ^ y + H z H ^ z ) 2 ) u ( 0.5 ) + ( E y ) ( β 2 ) ( H z ) ( 0.5 ) ( E z ) ( β 1 ) ( H y ) ( 0.5 ) H 0 ( v ( 0.5 ) ( H y ) ( 0.5 ) + w ( 0.5 ) ( H z ) ( 0.5 ) ) ] s = 0 Energy conservation
H 0 = 0
11 z v y w t = 0 ( z v ( 0.5 ) y w ( 0.5 ) ) t ˇ = 0 Angular momentum conservation
H 0 = 0 , σ = ρ
12 s H y ρ t + ( H y s E z ) s = 0 s H y ρ t + ( ( H y ) ( β 1 ) s ( E z ) ( β 1 ) ) s = 0 Unknown
13 s H z ρ t + ( H z + s E y ) s = 0 s H z ρ t + ( ( H z ) ( β 2 ) + s ( E y ) ( β 2 ) ) s = 0 Unknown
In both the cases above, scheme (31) is invariant. The rotation generator X 5 is only admitted for β 1 = β 2 . The remaining generators are admitted by the scheme for any set of parameters α , β 1 , β 2 and γ > 1 .
The conservation laws possessed by system (31) and their finite-difference counterparts are given in Table 2. Here and further, conservation laws whose fluxes vanish for H 0 = 0 are marked with †. In case H 0 = 0 , their densities preserve along the pathlines.

3.2. The Case of Infinite Conductivity σ

In this case, system (1) reduces to
1 ρ t = u s ,
u t = p + ( H y ) 2 + ( H z ) 2 2 s , x t = u ,
v t = H 0 H s y , y t = v ,
w t = H 0 H s z , z t = w ,
H y ρ t = ( H 0 v ) s ,
H z ρ t = ( H 0 w ) s ,
p t = γ p ρ u s ,
where the internal energy is given by (3).
In addition to the analogues of conservation laws presented in the previous section, system (36) possesses the conservation law of angular momentum, namely
( z v y w ) t + H 0 ( y H z z H y ) s = 0 .
As σ , scheme (31) becomes
1 ρ t = u s ( 0.5 ) ,
u t = p ( α ) + H y H ^ y + H z H ^ z 2 s ¯ , v t = H 0 ( H y ) s ¯ ( 0.5 ) , w t = H 0 ( H z ) s ¯ ( 0.5 ) ,
H y ρ t = H 0 v s ( 0.5 ) , H z ρ t = H 0 w s ( 0.5 ) ,
ε t = p ( α ) u s ( 0.5 ) ,
x t = u ( 0.5 ) , y t = v ( 0.5 ) , z t = w ( 0.5 ) , x s = 1 ρ .
One can verify that scheme (38) is an invariant one. As the symmetries of (36) and the corresponding difference schemes are reviewed in Section 3.2.3, we defer our discussion until then.

3.2.1. Conservation of Angular Momentum and Energy

Apparently, the latter scheme does not preserve angular momentum, i.e., it does not possess a difference analogue of the conservation law (37). One can verify it by algebraic manipulations with the scheme or with the help of the finite-difference analogue of the direct method [33]. We overcome this issue by modifying the latter scheme as follows
1 ρ t = u s ( 0.5 ) ,
u t = p ( α ) + H y H ^ y + H z H ^ z 2 s ¯ , v t = H 0 H ^ s ¯ y , w t = H 0 H ^ s ¯ z ,
H y ρ t = H 0 v s , H z ρ t = H 0 w s ,
ε t = p ( α ) u s ( 0.5 ) ,
x t = u ( 0.5 ) , y t = v , z t = w , x s = 1 ρ .
The latter modification allows one to obtain the whole set of finite-difference analogues of the conservation laws of Equation (36) excluding the conservation of the entropy along the pathlines. The conservation laws are presented in Table 3. Note that the three-layer conservation law of energy given in the table can be rewritten in the following two-layer form by means of (39c)
ε + u 2 + u + 2 + v 2 + v + 2 + w 2 + w + 2 4 + ( H y ) 2 + ( H z ) 2 2 ρ + τ 2 H 0 ( H y v s + H z w s ) t + p ( α ) + ( H y H ^ y + H z H ^ z ) 2 u ( 0.5 ) H 0 ( v ( 0.5 ) H ^ y + w ( 0.5 ) H ^ z ) s = 0
In addition, in order to verify the conservation law (37), one has to consider the following equations, which can be obtained by integration of (39c)
y s = H y H 0 ρ , z s = H z H 0 ρ .
We also note that the modified scheme (39) is still invariant and a completely conservative one.

3.2.2. Conservation of the Entropy along the Pathlines

From the latter system (36) it follows that
p ρ γ t = S t = 0 .
This represents the conservation of the entropy S along pathlines, which is a crucial difference between the finite and infinite conductivity cases.
It is known [37] that the Samarskiy–Popov scheme for polytropic gas does not preserve the entropy S for arbitrary γ . However, the following relation holds on solutions of the system
p ^ p p ( α ) = γ ρ ^ ρ ρ ( α ) ,
which approximates the differential relation
d p p = γ d ρ ρ .
The latter relation holds along trajectories of the particles up to O ( τ ) for α 0.5 or up to O ( τ 2 ) for α = 0.5 .
In [37], an entropy-preserving invariant scheme for gas dynamics equations in the case of a polytropic gas with γ = 3 was proposed. This scheme conserves the entropy along the pathlines but has only one conservation law, namely the conservation law of mass. It seems that the conservation of entropy by the difference scheme usually leads to the “loss” of some other conservation laws.
Here, we propose a way of preserving the entropy along the pathlines for polytropic gas with integer values of adiabatic exponent γ 2 for scheme (39). We show that this can be done by choosing appropriate approximations of the state Equation (3).
We note that by means of (3) and (36a), Equation (36g) can be represented as the identity
1 γ 1 ρ γ 1 t = ρ γ 2 ρ t .
In the finite-difference case, the rules of differentiation are different. As a result, not every approximation of the latter identity is a finite-difference identity. For a proper discrete analogue of (45) the right hand side of the identity should also be expressed in the divergent form as well as the left hand side. Choosing the difference approximation for the scheme in the case of a polytropic gas, one has an additional “degree of freedom”: the choice of approximation for the state Equation (3). This should be done so that both the left and right hand sides of the resulting approximation for (45) are divergent expressions. Note that this does not affect the conservativeness of the total energy conservation law equation since it does not depend on any specific form of the equation of state.
Further, we consider the shifted version of Equation (38d)
ε ˇ t = p ˇ ( α ) u ˇ s ( 0.5 ) = p ˇ ( α ) 1 ρ ˇ t .
First, we choose the following approximation of the state Equation (3) for γ = 2 ,
ε = p ( α ) ρ ^ .
Substituting (47) into (46), one derives
p ( α ) ρ ^ p ˇ ( α ) ρ = τ p ˇ ( α ) 1 ρ ˇ t .
Solving with resect to p ˇ ( α ) , one obtains
p ˇ ( α ) = p ( α ) ρ ˇ ρ ^ .
The latter equation can be rewritten as
p ( α ) p ˇ ( α ) = ρ ρ ^ ρ ˇ ρ .
Equation (50) can be integrated, i.e.,
p ˇ ( α ) ρ ρ ˇ t = S t = 0 .
This means conservation of entropy S along pathlines for γ = 2 on two time layers. We have achieved the integrability of the difference analogue of Equation (36g) by choosing a suitable approximation for the state equation.
In a similar way one can arrive at the conservation of entropy for γ = 3 , namely
ε = ρ p ( α ) ρ ^ ( ρ ^ + ρ ) , 2 p ˇ ( α ) ρ ρ ˇ ( ρ + ρ ˇ ) t = 0 .
Similarly, for γ = 4
ε = ρ 2 p ( α ) ρ ^ ( ρ ^ 2 + ρ ρ ^ + ρ 2 ) , 3 p ˇ ( α ) ρ ρ ˇ ( ρ 2 + ρ ρ ˇ + ρ ˇ 2 ) t = 0 ,
etc.
Thus, by induction, one establishes the following general formula for an arbitrary natural γ 2
ε = p ( α ) k = 0 γ 2 ρ ^ γ k 1 ρ k γ + 2 , ( γ 1 ) p ˇ ( α ) k = 0 γ 2 ρ γ k 1 ρ ˇ k + 1 t = 0 .
Entropy-preservation Formulas (54) are presented in Table 3 among the other conservation laws.
Remark 4.
From the preservation of entropy S t = 0 in the differential case it follows that (for simplicity, we consider the specific case γ = 2 )
0 T p ρ 2 t d t = p ( T , s ) ρ 2 ( T , s ) p ( 0 , s ) ρ 2 ( 0 , s ) = const .
Since the constant can be omitted, this means
p ( 0 , s ) ρ 2 ( 0 , s ) = p ( T , s ) ρ 2 ( T , s ) .
In the finite-difference case, by means of (51), one derives the following analogue of (55)
k = 0 N α p m n + k + ( 1 α ) p m n 1 + k ρ m n + k ρ m n 1 + k t τ = α p m n + N + ( 1 α ) p m n 1 + N ρ m n + N ρ m n 1 + N α p m n + ( 1 α ) p m n 1 ρ m n ρ m n 1 = 0 ,
where N = T / τ and we recall that p ( α ) = α p ^ + ( 1 α ) p . Similar to the differential case, the latter gives
α p m n + ( 1 α ) p m n 1 ρ m n ρ m n 1 = α p m n + N + ( 1 α ) p m n 1 + N ρ m n + N ρ m n 1 + N
which means entropy preservation for a given liquid particle.
Remark 5.
The approach described above can also lead to entropy conservation for rational values of γ. Without proof of the existence of a general formula we present the result for γ = 5 / 3 , which occurs for one-atomic ideal gas. One can verify that for γ = 5 / 3 the approximation
ε = p ( α ) ρ ^ 2 / 3 + ( ρ ρ ^ ) 1 / 3 + ρ 2 / 3 ρ 1 / 3 ρ ^ ( ρ ^ 1 / 3 + ρ 1 / 3 ) = 3 2 p ρ + O ( τ )
for the internal energy ε leads to the following preservation of entropy
2 3 p ˇ ( α ) ρ ˇ 2 / 3 + ( ρ ρ ˇ ) 1 / 3 + ρ 2 / 3 ρ ρ ˇ ( ρ ˇ 1 / 3 + ρ 1 / 3 ) t = 0 .
Remark 6.
Note that in the case H 0 = 0 , according to Table 3, scheme (39) possesses an infinite set of conservation laws for the following form
Φ ( γ 1 ) p ( α ) k = 0 γ 2 ρ ^ γ k 1 ρ k + 1 , H y ρ , H z ρ , v , w , y t v , z t w t = 0
where γ N \ { 1 } and Φ is an arbitrary function of its arguments.
Remark 7.
From (51) it follows that
p ˇ ( α ) ρ ρ ˇ S = 0 .
The Taylor series expansion of the latter equation is
p ρ 2 S + p ρ t ρ 3 + ( α 1 ) p t ρ 2 τ + O ( τ 2 ) = 0 .
Equation (43) for γ = 2 can be represented as
p ( α ) ( ρ ( α ) ) 2 ρ t 2 p t = 0 .
The corresponding expansion is
p ρ 2 ρ t 2 p t + ρ p t 2 ρ t p ρ 3 α + ( ρ t p t t 2 p t ρ t t ) ρ t 2 p t 2 τ + O ( τ 2 ) = 0 .
Equations (61) and (59) approximate the conservation of entropy with the same order O ( τ ) . In contrast to (61), approximation (59) can be written in a divergent form. Thus, it represents a conservation law of the scheme, while (61) does not. This gives an advantage in the case of isentropic flows when additional conservation laws include entropy. Then, the expression for the entropy given by Equation (59) can be considered as a constant and included in conserved quantities. Invariant schemes and their conservation laws in the case of isentropic flows are discussed in the next section.

3.2.3. On Specific Symmetries and Conservation Laws in the Case of Isentropic Flows ( S = const )

According to [1], in case S = p / ρ γ = const , Equation (36) admit the following symmetries.
  • If H 0 0 , the admitted Lie algebra is
    X 1 = t , X 2 = s , X 3 = x , X 4 = t x + u , X 5 = z y y z + w v v w + E z E y E y E z + H z H y H y H z , X 6 = t t + 2 s s u u v v w w + 2 ρ ρ , X 7 = s s + x x + y y + z z + u u + v v + w w 2 ρ ρ , X 8 = q 1 ( s ) y , X 9 = q 2 ( s ) z , X 10 = t y + v , X 11 = t z + w ,
    where q 1 , q 2 are arbitrary functions of s.
    Scheme (39) supplemented by the state Equation (54) admits all the generators (63).
  • In case H 0 = 0 , the admitted Lie algebra is
    X 1 = t , X 2 = s , X 3 = x , X 4 = t x + u , X 5 = q 3 ( s ) ( H z H y H y H z ) , X 6 = t t + 2 s s u u + 2 ρ ρ , X 7 = s s + x x + u u 2 ρ ρ , X 8 = 2 s s + 2 ρ ρ + p p + H y H y + H z H z .
    In case γ = 2 , there are two additional generators, namely
    X 9 = q 4 ( s ) ρ ( H y H y p ) , X 10 = q 5 ( s ) ρ ( H z H z p ) .
    Here, q 3 , q 4 , and q 5 are arbitrary functions of s.
    Scheme (39), (54) admits all the generators (64). However, the scheme does not admit the generators X 9 and X 10 .
There are the following additional conservation laws for system (36).
(a)
In case H 0 0 , there is an additional conservation law that corresponds to the generator s
u ρ + v H y + w H z H 0 ρ t + γ S γ 1 ρ γ 1 u 2 + v 2 + w 2 2 s = 0 .
(b)
Case H 0 = 0 .
  • The conservation law corresponding to the generator s is
    u ρ t + γ S 0 γ 1 ρ γ 1 u 2 2 + ρ B 0 s = 0
    provided
    S 0 = S = const and B 0 = ( H y ) 2 + ( H z ) 2 ρ 2 = const .
    The latter follows from system (36). When conductivity of the medium tends to infinity, the phenomenon of frozen-in magnetic field is observed (see, e.g., [41]). In this case, in the absence of the longitudinal component H 0 of the magnetic field, the quantity B 0 , which is proportional to the magnetic pressure, turns out to be preserved along the pathlines.
  • In case γ = 2 , the admitted generator
    5 X 6 + 3 X 7 4 X 8
    corresponds to the conservation law
    5 t ρ A 0 + 5 t u 2 2 s u ρ 3 x u t + ( ( 5 t u 3 x ) ρ 2 2 s ρ ) A 0 + s u 2 2 s = 0
    provided
    A 0 = S γ 1 + ( H y ) 2 + ( H z ) 2 2 ρ 2 γ = 2 = p ρ 2 + ( H y ) 2 + ( H z ) 2 2 ρ 2 = const
    which follows from system (36).
Remark 8.
Conservation law (70) is a basis one. Its partial derivative with respect to s is equivalent to (67), and its partial derivative with respect to t is
A 0 ( ρ t + ρ 2 u s ) + u ( u t + ( A 0 ρ 2 ) s ) = 0
which is a combination of (36a) and (36b) provided (71).
By virtue of the content of Remark 6, one can verify that the finite-difference analogues of (68) and (71) hold along the pathlines for scheme (39), namely
( H y ) 2 + ( H z ) 2 ρ 2 t = 0 or H y H ˇ y + H z H ˇ z ρ ρ ˇ t = 0 if H 0 = 0 ,
and
p ( α ) ρ ˇ ρ + ( H y ) 2 + ( H z ) 2 2 ρ 2 t = 0 or p ( α ) ρ ˇ ρ + H y H ˇ y + H z H ˇ z 2 ρ ρ ˇ t = 0 if H 0 = 0 , γ = 2 .
Scheme (39) also admits the generators s and (69) under the same conditions as for the differential case.
Analyzing scheme (39), one can conclude that for the additional conservation laws (66), (67), and (70) there are no approximations in terms of rational expressions. This means that construction of finite difference analogues of the mentioned conservation laws is extremely hard.
Further, we restrict ourselves to the case γ = 2 and S = S 1 = const , and consider another invariant scheme on an extended finite-difference stencil.
We introduce the pressure for the polytropic gas as
p = S 1 ρ ^ ρ ^ + .
Then, the conservation law of entropy
p ˇ ρ ρ + t = ( S 1 ) t = 0
is defined by the following invariant expression
p ρ ^ ρ ^ + = p ˇ ρ ρ + = p ˇ ρ ρ = p ρ ^ ρ ^ = S 1 .
The scheme under consideration is based on scheme (39) and it has the following form
1 ρ t = u s ,
u t = p + H ^ + y H ^ y + H ^ + z H ^ z 2 s ¯ , v t = H 0 H ^ s y , w t = H 0 H ^ s z ,
H y ρ t = H 0 v s ¯ , H z ρ t = H 0 w s ¯ ,
p ρ ^ ρ ^ + = p ˇ ρ ρ + = p ˇ ρ ρ = p ρ ^ ρ ^ = S 1 ,
x t = u , y t + = v , z t + = w , x s = 1 ρ .
One can verify that the latter scheme is invariant. It admits the same symmetries as scheme (39) and (54). The following quantities hold for (78)
H + y H y + H + z H z ρ ρ + t = ( B 1 ) t = 0 if H 0 = 0 ,
p ˇ ρ ρ + + H + y H y + H + z H z 2 ρ ρ + t = 1 2 ( 2 S 1 + B 1 ) t = 0 if H 0 = 0 , γ = 2 .
Scheme (78) possesses the difference analogues of (66) and (67), namely
u ρ + v H ^ y + w H ^ z H 0 ρ ^ t + 2 ρ ^ S 1 u u + v ^ 2 + w ^ 2 2 s = 0 ,
u ρ t + ρ ^ ( 2 S 1 + B 1 ) u u 2 s = 0 .
To construct the latter conservation law one should use the following relation
H ^ + y H ^ y + H ^ + z H ^ z 2 = H ^ + y H ^ y + H ^ + z H ^ z 2 ρ ^ ρ ^ + = 1 2 ρ ^ ρ ^ + B 1 .
Remark 9.
The angular momentum and center of mass conservation laws are
( z v y w ) t + H 0 ( y ^ H ^ z z ^ H ^ y ) s = 0 ,
( t u x ) t + t ^ p + H ^ y H ^ y + H ^ z H ^ z 2 s = 0 ,
( t v y + ) t t ^ H 0 H ^ y s = 0 ,
( t u z + ) t t ^ H 0 H ^ z s = 0 .
The remaining conservation laws of mass, momentum, magnetic flux, and entropy follow directly from the scheme as it is written in a divergent form.

4. Numerical Experiments

In this section, we consider the problem of deceleration of a plasma bunch in a crossed electromagnetic field under the presence and absence of a longitudinal component of magnetic field. We use scheme (31), and consider how the conservation laws hold on the solutions of this scheme. In addition to the transverse component H y of the magnetic field, we also consider the case of the presence of a longitudinal magnetic field H 0 0 .
A plasma bunch is considered, which moves from left to right in a railgun channel. The channel is filled with a relatively cold weakly conducting gas. With the help of an external electric circuit, a strong transverse magnetic field is generated in the channel, which causes the bunch to decelerate. During its motion, the plasma bunch closes the electric circuit and moves along the background gas; therefore, the magnetic field and pressure at the left boundary of the computational domain are considered equal to zero. The differential boundary conditions are as follows
p ( 0 , t ) = 0 , H y ( 0 , t ) = 0 , H z ( 0 , t ) = 0 ,
u ( S , t ) = 0 , E y ( S , t ) = 0 , H y ( S , t ) = 4 π J ( t ) ,
L 0 d J d t + R 0 J V ( t ) + E z ( S , t ) = 0 ,
d V d t = J / C 0 , V ( 0 ) = V 0 , J ( 0 ) = 0 ,
where 0 s S and 0 t t max , S is the total mass of the gas, J and V are current and voltage, and C 0 , L 0 , R 0 are the external circuit parameters. The boundary conditions (88c) and (88d) are approximated in the same way as in [13], namely
L 0 J t + R 0 J ( 0.5 ) V ( 0.5 ) + E M z ( 0.5 ) = 0 , V t = J ( 0.5 ) / C 0 ,
where M = S / h .
All calculations were carried out using the dimensionless version of scheme (31) with the value of the coefficient κ = 4 π . For the dimensionless form of the scheme, the initial conditions are: ρ 0 = 1.0 , p 0 = 0.0056 , R 0 = 1.17 , C 0 = 1.64 , L 0 = 0.0035 , S = 4.0 , the temperature of the plasma T 0 = 3.0 , and the initial speed of the plasma bunch u 0 = 0.75 . The gas is considered polytropic with γ = 5 / 3 . The uniform mesh steps are h = 0.067 and τ = 0.003 , and t max = 0.7 . The initial voltage V 0 is varied between 1.67 and 2.6 which approximately correspond to the voltage 650 and 1000 V. In experiments where the longitudinal magnetic field H 0 is present, a value close to 1 is taken for H 0 . In the calculations, a linear artificial viscosity is used, with a viscosity coefficient ν = 2 h .
The problem under consideration is close to the problem described in [42] (see also [43]) in which, however, tabulated real plasma parameters, including electrical conductivity, were used. In our problem, we used the ideal gas equation and an exponential conductivity function.
Scheme (31) is implemented using the iterative methods described in [13]. In this case, the scheme equations are divided into two parts, dynamic and magnetic. The dynamic part is preliminarily linearized using the Newton method, and for the magnetic part a flow version of the sweep method is used [44], which is well suited for the case of finite conductivity, especially when its values are small. The bunch motion is modeled by a shock wave. Conductivity σ of the plasma bunch is proportional to T 3 / 2 , and the conductivity function σ = σ ( ρ , T ) is very sensitive to the density ρ in such a way that in the rarefied background gas region it has values close to zero.
Three essentially different cases are considered:
  • The bunch is decelerated using a transverse magnetic field H y at a relatively low voltage in the circuit.
  • The bunch is decelerated using a transverse magnetic field at a high voltage in the circuit.
  • A rather strong longitudinal magnetic field H 0 is added to the previous case. (Calculations show that a weak longitudinal magnetic field has little effect on the experimental results.)
In all cases, at the initial moment of time, the gas particles are given a small constant transverse velocity v > 0 . This is necessary in order to track the influence of the longitudinal magnetic field on the transverse component of the particle velocity, which should be observed only in the third numerical experiment.
Figure 2 shows the evolution of the magnetic field and plasma temperature in the first experiment. The magnetic field is not strong enough to stop the bunch. If the bunch reaches the right boundary of the computational domain, the reflection of the wave can be observed due to the boundary condition u ( S , t ) = 0 . Figure 3 shows the second case where the transverse magnetic field is strong enough. The plasma bunch is decelerated by the magnetic field and after a short period of time begins to move backward. Adding a sufficiently strong longitudinal magnetic field H 0 to the previous experiment leads to an intermediate picture: the magnetic field is “smeared” over the computational domain, the plasma deceleration process is not as intense as in the previous case, and is inhomogeneous along the mass coordinate, which leads to a kind of fragmentation of the temperature profile (see Figure 4).
In Figure 5, Figure 6 and Figure 7 the evolution of the work A = p u s performed on the gas and the electromagnetic force f x in the direction of the x axis are shown. At the beginning of the process at the shock wave front, A is positive, which corresponds to gas heating due to compression. At the left boundary of the computational domain, A is negative, and the rarefaction wave cools the gas. The electromagnetic force f is mainly localized in the front part of the bunch. With an increase in the total current, the electromagnetic force increases, which leads to the deceleration of the wave and further to its stop and reverse motion. Note that the component f y = H 0 H s ¯ y appears only in the third experiment ( H 0 0 ). Its evolution is depicted in Figure 8. In addition, recall that, as follows from (31),
f x = 1 2 σ ^ ρ ^ ( E ^ z H y E ^ y H z ) + σ ρ ( E z H ^ y E y H ^ z ) .
In Figure 9 the trajectories of particles under the action of magnetic fields are shown. The left part (Figure 9a–c) shows x-trajectories of particles for three experiments. The right side of the figure shows y-trajectories associated with the transverse velocity component v. Figure 9d corresponds to the first and second experiments where v has a constant value and H 0 = 0 . Figure 9e,f correspond to the third experiment at H 0 > 0 and H 0 < 0 where under the action of the longitudinal magnetic field the transverse velocity component increases or slows down accordingly. Note that the choice of sign of the value H 0 otherwise does not affect the results of the third experiment.
In Figure 10, Figure 11 and Figure 12 the finite-difference conservation laws of energy, magnetic flux (along the y axis), momentum, and center-of-mass motion (along the x axis) are given for the selected moment in time, when the interaction of magnetic fields and the plasma bunch is already quite intense. The results are provided only for the third experiment, since in other cases the control of conservation laws gives similar results. The accurate enough preservation of the conservation laws on solutions is due to the conservativeness of scheme (31).

5. Conclusions

Finite-difference schemes for MHD equations in the case of plane one-dimensional flows were considered. The Samarsky–Popov classical scheme for the case of finite conductivity was taken as a starting point. Symmetries and conservation laws of this scheme were investigated. It was shown that the scheme admits the same symmetries as the original differential model. It also has difference analogues of the conservation laws of the original model. In addition to the conservation laws previously known for the scheme, new conservation laws were given, which were obtained on the basis of the group classification recently carried out in [1].
The classical Samarskiy–Popov scheme was generalized to the case of arbitrary vectors of electric and magnetic fields, as well as to the case of infinite conductivity. In the case of finite conductivity the scheme possesses difference analogues of all differential local conservation laws obtained in [1], some of which were not previously known. In the case of infinite conductivity, straightforward generalization of the scheme leads to a scheme that does not preserve angular momentum. The proposed modification makes it possible to obtain an invariant scheme that also possesses the conservation law of angular momentum. In addition, it was shown how to approximate the equation of state for a polytropic gas to preserve the entropy along the pathlines on the extended stencil for two time layers.
A numerical implementation of the generalized Samarsky–Popov scheme for the case of finite conductivity was performed for the problem of deceleration of a plasma bunch by crossed electromagnetic fields. Various cases of the action of fields on a plasma were considered. Calculations showed that the finite-difference conservation laws are preserved on the solutions of the scheme quite accurately.

Author Contributions

Conceptualization, V.D.; methodology, V.D. and E.K.; software, E.K.; validation, E.K.; investigation, V.D. and E.K.; data curation, E.K.; writing—original draft preparation, V.D. and E.K.; writing—review and editing, V.D. and E.K.; visualization, E.K.; supervision, V.D.; project administration, V.D.; funding acquisition, V.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by Russian Science Foundation Grant No. 18-11-00238 “Hydrodynamics-type equations: symmetries, conservation laws, invariant difference schemes”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors thank E. Schulz and S.V. Meleshko for valuable discussions. E.K. sincerely appreciates the hospitality of the Suranaree University of Technology.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dorodnitsyn, V.A.; Kaptsov, E.I.; Kozlov, R.V.; Meleshko, S.V.; Mukdasanit, P. Plane one-dimensional MHD flows: Symmetries and conservation laws. Int. J. Non-Linear Mech. 2022, 140, 103899. [Google Scholar] [CrossRef]
  2. Gridnev, N. Study of the magnetohydrodynamics equations’ group properties and invariant solutions. J. Appl. Mech. Tech. Phys. 1968, 6, 103–107. (In Russian) [Google Scholar]
  3. Dorodnitsyn, V.A. On invariant solutions of one-dimensional nonstationary magnetohydrodynamics with finite conductivity. Keldysh Inst. Prepr. 1976, 143. Available online: https://library.keldysh.ru/preprint.asp?id=1976-143 (accessed on 28 February 2022). (In Russian).
  4. Rogers, C. Invariant Transformations in Non-Steady Gasdynamics and Magneto-Gasdynamics. Zeit Angew. Math. Phys. 1969, 20, 370–382. [Google Scholar] [CrossRef]
  5. Ibragimov, N.H. (Ed.) CRC Handbook of Lie Group Analysis of Differential Equations; CRC Press: Boca Raton, FL, USA, 1995; Volume 2. [Google Scholar]
  6. Oliveri, F.; Speciale, M.P. Exact solutions to the ideal magneto-gas-dynamics equations through Lie group analysis and substitution principles. J. Phys. A Math. Gen. 2005, 38, 8803–8820. [Google Scholar] [CrossRef]
  7. Picard, P. Some exact solutions of the ideal MHD equations through symmetry reduction. J. Math. Anal. Appl. 2008, 337, 360–385. [Google Scholar] [CrossRef] [Green Version]
  8. Golovin, S. Regular partially invariant solutions of defect 1 of the equations of ideal magnetohydrodynamics. J. Appl. Mech. Tech. Phys. 2019, 50, 171–180. [Google Scholar] [CrossRef]
  9. Golovin, S. Natural curvilinear coordinates for ideal MHD equations. Non-stationary flows with constant total pressure. Phys. Lett. Sect. A Gen. At. Solid State Phys. 2011, 375, 283–290. [Google Scholar] [CrossRef] [Green Version]
  10. Golovin, S.; Sesma, L. Exact Solutions of Stationary Equations of Ideal Magnetohydrodynamics in the Natural Coordinate System. J. Appl. Mech. Tech. Phys. 2019, 60, 234–247. [Google Scholar] [CrossRef]
  11. Toro, E.F. Riemann Solvers and Numerical Methods for Fluid Dynamics; Springer: Berlin/Heidelberg, Germany, 1997. [Google Scholar]
  12. Samarskiy, A.A.; Popov, Y.P. Completely conservative difference schemes for the equations of magneto-hydrodynamics. USSR Comput. Math. Math. Phys. 1970, 10, 233–243. [Google Scholar]
  13. Samarskii, A.A.; Popov, Y.P. Difference Methods for Solving Problems of Gas Dynamics; Nauka: Moscow, Russian, 1980. (In Russian) [Google Scholar]
  14. Falle, S.A.E.G.; Komissarov, S.S.; Joarder, P. A multidimensional upwind scheme for magnetohydrodynamics. Mon. Not. R. Astron. Soc. 1998, 297, 265–277. [Google Scholar] [CrossRef] [Green Version]
  15. Powell, K.G.; Roe, P.L.; Linde, T.J.; Gombosi, T.I.; De Zeeuw, D.L. A Solution-Adaptive Upwind Scheme for Ideal Magnetohydrodynamics. J. Comput. Phys. 1999, 154, 284–309. [Google Scholar] [CrossRef]
  16. Yakovlev, S.; Xu, L.; Li, F. Locally divergence-free central discontinuous Galerkin methods for ideal MHD equations. J. Comput. Sci. 2013, 4, 80–91. [Google Scholar] [CrossRef]
  17. Yang, Y.; Feng, X.S.; Jiang, C.W. A high-order CESE scheme with a new divergence-free method for MHD numerical simulation. J. Comput. Phys. 2017, 349, 561–581. [Google Scholar] [CrossRef]
  18. Hirabayashi, K.; Hoshino, M.; Amano, T. A new framework for magnetohydrodynamic simulations with anisotropic pressure. J. Comput. Phys. 2016, 327, 851–872. [Google Scholar] [CrossRef] [Green Version]
  19. Ryu, D.; Jones, T.W. Numerical Magnetohydrodynamics in Astrophysics: Algorithm and Tests for One-dimensional Flow. ApJL 1995, 442, 228. [Google Scholar] [CrossRef]
  20. Noether, E. Invariante Variationsprobleme. Konigliche Gesellschaft der Wissenschaften zu Gottingen, Nachrichten, Mathematisch-Physikalische Klasse Heft 2. Transp. Theory Stat. Phys. 1971, 1, 183–207. [Google Scholar]
  21. Ibragimov, N.H. Transformation Groups Applied to Mathematical Physics; Reidel: Boston, MA, USA, 1985. [Google Scholar]
  22. Olver, P.J. Applications of Lie Groups to Differential Equations; Springer: New York, NY, USA, 1986. [Google Scholar]
  23. Bluman, G.W.; Anco, S.C. Symmetry and Integration Methods for Differential Equations; Springer: New York, NY, USA, 2002; p. 422. [Google Scholar]
  24. Kreyszig, E. Advanced Engineering Mathematics; Wiley: New York, NY, USA, 1983. [Google Scholar]
  25. Dorodnitsyn, V.A. Transformation groups in net spaces. J. Sov. Math. 1991, 55, 1490–1517. [Google Scholar] [CrossRef]
  26. Maeda, S. Extension of discrete Noether theorem. Math. Jpn. 1985, 26, 85–90. [Google Scholar]
  27. Maeda, S. The similarity method for difference equations. J. Inst. Math. Appl. 1987, 38, 129–134. [Google Scholar] [CrossRef]
  28. Dorodnitsyn, V.A. Applications of Lie Groups to Difference Equations; CRC Press: Boca Raton, FL, USA, 2011. [Google Scholar]
  29. Dorodnitsyn, V.A.; Kozlov, R.V.; Winternitz, P. Continuous symmetries of Lagrangians and exact solutions of discrete equations. J. Math. Phys. 2004, 45, 336–359. [Google Scholar] [CrossRef] [Green Version]
  30. Dorodnitsyn, V.A.; Kaptsov, E.I.; Kozlov, R.V.; Winternitz, P. The adjoint equation method for constructing first integrals of difference equations. J. Phys. A Math. Theor. 2015, 48, 055202. [Google Scholar] [CrossRef]
  31. Bourlioux, A.; Cyr-Gagnon, C.; Winternitz, P. Difference schemes with point symmetries and their numerical tests. J. Phys. A Math. Gen. 2006, 39, 6877–6896. [Google Scholar] [CrossRef]
  32. Bourlioux, A.; Rebelo, R.; Winternitz, P. Symmetry Preserving Discretization of SL(2, ℝ) Invariant Equations. J. Nonlinear Math. Phys. 2007, 15, 362–372. [Google Scholar] [CrossRef] [Green Version]
  33. Dorodnitsyn, V.A.; Kaptsov, E.I. Shallow water equations in Lagrangian coordinates: Symmetries, conservation laws and its preservation in difference models. Commun. Nonlinear. Sci. Numer. Simulat. 2020, 89, 105343. [Google Scholar] [CrossRef]
  34. Dorodnitsyn, V.A.; Kaptsov, E.I.; Meleshko, S.V. Symmetries, conservation laws, invariant solutions and difference schemes of the one-dimensional Green-Naghdi equations. J. Nonlinear Math. Phys. 2020, Submitted. [Google Scholar] [CrossRef]
  35. Dorodnitsyn, V.A.; Kaptsov, E.I. Discrete shallow water equations preserving symmetries and conservation laws. J. Math. Phys. 2021, 62, 083508. [Google Scholar] [CrossRef]
  36. Kaptsov, E.I.; Dorodnitsyn, V.A.; Meleshko, S.V. Conservative invariant finite-difference schemes for the modified shallow water equations in Lagrangian coordinates. Stud. Appl. Math. 2021, submitted. [Google Scholar]
  37. Dorodnitsyn, V.A.; Kozlov, R.; Meleshko, S.V. One-dimensional gas dynamics equations of a polytropic gas in Lagrangian coordinates: Symmetry classification, conservation laws, difference schemes. Commun. Nonlinear. Sci. Numer. Simulat. 2019, 74, 201–218. [Google Scholar] [CrossRef] [Green Version]
  38. Kozlov, R. Conservative difference schemes for one-dimensional flows of polytropic gas. Commun. Nonlinear. Sci. Numer. Simulat. 2019, 78, 104864. [Google Scholar] [CrossRef] [Green Version]
  39. Dorodnitsyn, V.A.; Kozlov, R.V.; Meleshko, S.V. One-dimensional flows of a polytropic gas: Lie group classification, conservation laws, invariant and conservative difference schemes. In Symmetries and Applications of Differential Equations: In Memory of Nail H. Ibragimov (1939–2018); Albert, C.J.L., Gazizov, R.K., Eds.; Nonlinear Physical Science; Springer: Singapore, 2021. [Google Scholar]
  40. Rojdestvenskiy, B.L.; Yanenko, N.N. Systems of Quasilinear Equations and Their Applications to Gas Dynamics; Nauka: Moscow, Russian, 1968. (In Russian) [Google Scholar]
  41. Kulikovskiy, A.G.; Lyubimov, G.A. Magnetohydrodynamics; International Series in Physics, Reading, Mass; Addison-Wesley: Boston, MA, USA, 1965. [Google Scholar]
  42. Danilova, G.V. Dorodnitsyn, V.A.; Kurdyumov, S.P.; Popov, Y.P.; Smarskiy, A.A.; Tsaryova, L.S. Interaction of a plasma clot with a magnetic field in a railgun channel. Keldysh Inst. Prepr. 1973, 62. Available online: https://library.keldysh.ru/preprint.asp?id=1973-63 (accessed on 28 February 2022). (In Russian).
  43. Samarskiy, A.A.; Dorodnitsyn, V.A.; Kurdyumov, S.P.; Popov, Y.P. Formation of T-layers in the process of plasma deceleration by a magnetic field. Dokl. Akad. Nauk SSSR 1974, 216, 1254–1257. (In Russian) [Google Scholar]
  44. Degtyarev, L.M.; Favorskii, A.P. A flow variant of the sweep method. USSR Comput. Math. Math. Phys. 1968, 8, 252–261. [Google Scholar] [CrossRef]
Figure 1. A plane one-dimensional flow for the chosen coordinates.
Figure 1. A plane one-dimensional flow for the chosen coordinates.
Mathematics 10 01250 g001
Figure 2. Evolution of the magnetic field H y and the temperature T for the first experiment.
Figure 2. Evolution of the magnetic field H y and the temperature T for the first experiment.
Mathematics 10 01250 g002
Figure 3. Evolution of the magnetic field H y and the temperature T for the second experiment.
Figure 3. Evolution of the magnetic field H y and the temperature T for the second experiment.
Mathematics 10 01250 g003
Figure 4. Evolution of the magnetic field H y and the temperature T for the third experiment.
Figure 4. Evolution of the magnetic field H y and the temperature T for the third experiment.
Mathematics 10 01250 g004
Figure 5. Evolution of the electromagnetic force f x and the work on gas A for the first experiment.
Figure 5. Evolution of the electromagnetic force f x and the work on gas A for the first experiment.
Mathematics 10 01250 g005
Figure 6. Evolution of the electromagnetic force f x and the work on gas A for the second experiment.
Figure 6. Evolution of the electromagnetic force f x and the work on gas A for the second experiment.
Mathematics 10 01250 g006
Figure 7. Evolution of the electromagnetic force f x and the work on gas A for the third experiment.
Figure 7. Evolution of the electromagnetic force f x and the work on gas A for the third experiment.
Mathematics 10 01250 g007
Figure 8. Evolution of the electromagnetic force f y for the third experiment ( H 0 > 0 ).
Figure 8. Evolution of the electromagnetic force f y for the third experiment ( H 0 > 0 ).
Mathematics 10 01250 g008
Figure 9. The trajectories of particles under the action of magnetic fields: (ac) show x-trajectories for cases 1, 2, and 3; (d) y-trajectories for H 0 = 0 ; (e) y-trajectories for H 0 > 0 ; (f) y-trajectories for H 0 < 0 .
Figure 9. The trajectories of particles under the action of magnetic fields: (ac) show x-trajectories for cases 1, 2, and 3; (d) y-trajectories for H 0 = 0 ; (e) y-trajectories for H 0 > 0 ; (f) y-trajectories for H 0 < 0 .
Mathematics 10 01250 g009
Figure 10. Conservation laws of energy (solid line) and y-flux (dashed line) at t = 0.64 .
Figure 10. Conservation laws of energy (solid line) and y-flux (dashed line) at t = 0.64 .
Mathematics 10 01250 g010
Figure 11. Conservation law of x-momentum at t = 0.64 .
Figure 11. Conservation law of x-momentum at t = 0.64 .
Mathematics 10 01250 g011
Figure 12. Conservation law of center-of-mass along axis x at t = 0.64 .
Figure 12. Conservation law of center-of-mass along axis x at t = 0.64 .
Mathematics 10 01250 g012
Table 1. Differential and difference conservation laws for system (18) and scheme (8).
Table 1. Differential and difference conservation laws for system (18) and scheme (8).
#Conservation Laws of the SystemConservation Laws of the SchemePhysics Interpretation
σ = σ ( p , ρ )
1 1 ρ t u s = 0 1 ρ t u ( 0.5 ) s = 0 Mass conservation
2 H ρ t E s = 0 H ρ t ( E ( β ) ) s = 0 Magnetic flux conservation
3 u t + p + H 2 2 s = 0 u t + p ( α ) + H H ^ 2 s = 0 Momentum conservation
4 t u x t + t p + H 2 2 s = 0 t + t ˇ 2 u x t + t p ( α ) + H H ^ 2 s = 0 Center of mass law
5 ε + u 2 2 + H 2 2 ρ t + + p + H 2 2 u E H s = 0 ε + u 2 + u + 2 4 + H 2 2 ρ t + + p ( α ) + ( H H ^ ) 2 u ( 0.5 ) E ( β ) H ( 0.5 ) s = 0 Energy conservation
σ = ρ
6 s H ρ t + ( H s E ) s = 0 s H ρ t + ( H ( β ) s E ( β ) ) s = 0 Unknown
Table 3. Differential and difference conservation laws for the modified scheme (39) (and the corresponding system (36) for an arbitrary entropy S ( s ) in case σ .
Table 3. Differential and difference conservation laws for the modified scheme (39) (and the corresponding system (36) for an arbitrary entropy S ( s ) in case σ .
#Conservation Laws of the SystemConservation Laws of the SchemePhysics Interpretation
1 1 ρ t u s = 0 1 ρ t u ( 0.5 ) s = 0 Mass conservation
2 H y ρ t ( H 0 v ) s = 0 H y ρ t H 0 v s = 0 Magnetic flux conservation
3 H z ρ t ( H 0 w ) s = 0 H z ρ t H 0 w s = 0 Magnetic flux conservation
4 u t + p + ( H y ) 2 + ( H z ) 2 2 s = 0 u t + p ( α ) + H y H ^ y + H z H ^ z 2 s = 0 Momentum conservation
5 v t ( H 0 H y ) s = 0 v t H 0 H y s = 0 Momentum conservation
6 w t ( H 0 H z ) s = 0 w t H 0 H z s = 0 Momentum conservation
7 t u x t + t p + ( H y ) 2 + ( H z ) 2 2 s = 0 t + t ˇ 2 u x t + t p ( α ) + H y H ^ y + H z H ^ z 2 s = 0 Center of mass law
8 ( t v y ) t ( t H 0 H y ) s = 0 t v y t t ^ H 0 H ^ y s = 0 Center of mass law
9 ( t w z ) t ( t H 0 H z ) s = 0 t w z t t ^ H 0 H ^ z s = 0 Center of mass law
10 ( ε + u 2 + v 2 + w 2 2 + ( H y ) 2 + ( H z ) 2 2 ρ ) t + [ ( p + ( H y ) 2 + ( H z ) 2 2 ) u + E y H z E z H y H 0 ( v H y + w H z ) ] s = 0 ( ε ˇ + u ˇ 2 + u ˇ + 2 + v ˇ 2 + v ˇ + 2 + w ˇ 2 + w ˇ + 2 4 + H y H ˇ y + H z H ˇ z 2 ρ ) t + [ ( p ˇ ( α ) + ( H y H ˇ y + H z H ˇ z ) 2 ) u ˇ ( 0.5 ) H 0 ( v ˇ ( 0.5 ) H y + w ˇ ( 0.5 ) H z ) ] s = 0 Energy conservation
11 ( z v y w ) t + H 0 ( y H z z H y ) s = 0 . ( z v y w ) t + H 0 ( y ^ H ^ z z ^ H ^ y ) s = 0 Angular momentum conservation
ε = p ( α ) k = 0 γ 2 ρ ^ γ k 1 ρ k γ + 2 , γ N \ { 1 }
12 p ρ γ t = 0 ( γ 1 ) p ˇ ( α ) k = 0 γ 2 ρ γ k 1 ρ ˇ k + 1 t = 0 Entropy conservation
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dorodnitsyn, V.; Kaptsov, E. Invariant Finite-Difference Schemes for Plane One-Dimensional MHD Flows That Preserve Conservation Laws. Mathematics 2022, 10, 1250. https://doi.org/10.3390/math10081250

AMA Style

Dorodnitsyn V, Kaptsov E. Invariant Finite-Difference Schemes for Plane One-Dimensional MHD Flows That Preserve Conservation Laws. Mathematics. 2022; 10(8):1250. https://doi.org/10.3390/math10081250

Chicago/Turabian Style

Dorodnitsyn, Vladimir, and Evgeniy Kaptsov. 2022. "Invariant Finite-Difference Schemes for Plane One-Dimensional MHD Flows That Preserve Conservation Laws" Mathematics 10, no. 8: 1250. https://doi.org/10.3390/math10081250

APA Style

Dorodnitsyn, V., & Kaptsov, E. (2022). Invariant Finite-Difference Schemes for Plane One-Dimensional MHD Flows That Preserve Conservation Laws. Mathematics, 10(8), 1250. https://doi.org/10.3390/math10081250

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop