Next Article in Journal
Weak Approximations of the Wright–Fisher Process
Next Article in Special Issue
First-Order Approximate Mei Symmetries and Invariants of the Lagrangian
Previous Article in Journal
The Influences of the Hyperbolic Two-Temperatures Theory on Waves Propagation in a Semiconductor Material Containing Spherical Cavity
Previous Article in Special Issue
Approximate Mei Symmetries and Invariants of the Hamiltonian
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

New Exact Solutions with a Linear Velocity Field for the Gas Dynamics Equations for Two Types of State Equations

by
Renata Nikonorova
,
Dilara Siraeva
* and
Yulia Yulmukhametova
Mavlyutov Institute of Mechanics, Subdivision of the Ufa Federal Research Centre of the Russian Academy of Sciences, 450054 Ufa, Russia
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(1), 123; https://doi.org/10.3390/math10010123
Submission received: 18 November 2021 / Revised: 14 December 2021 / Accepted: 27 December 2021 / Published: 1 January 2022

Abstract

:
In this paper, exact solutions with a linear velocity field are sought for the gas dynamics equations in the case of the special state equation and the state equation of a monatomic gas. These state equations extend the transformation group admitted by the system to 12 and 14 parameters, respectively. Invariant submodels of rank one are constructed from two three-dimensional subalgebras of the corresponding Lie algebras, and exact solutions with a linear velocity field with inhomogeneous deformation are obtained. On the one hand of the special state equation, the submodel describes an isochoric vortex motion of particles, isobaric along each world line and restricted by a moving plane. The motions of particles occur along parabolas and along rays in parallel planes. The spherical volume of particles turns into an ellipsoid at finite moments of time, and as time tends to infinity, the particles end up on an infinite strip of finite width. On the other hand of the state equation of a monatomic gas, the submodel describes vortex compaction to the origin and the subsequent expansion of gas particles in half-spaces. The motion of any allocated volume of gas retains a spherical shape. It is shown that for any positive moment of time, it is possible to choose the radius of a spherical volume such that the characteristic conoid beginning from its center never reaches particles outside this volume. As a result of the generalization of the solutions with a linear velocity field, exact solutions of a wider class are obtained without conditions of invariance of density and pressure with respect to the selected three-dimensional subalgebras.

1. Introduction

Group analysis of differential equations is a powerful tool for obtaining exact solutions to nonlinear differential equations [1,2]. The gas dynamics equations have been well studied from the point of view of group analysis within the framework of the program “Submodels” [3]. This became possible thanks to the outstanding scientist Ovsyannikov, who sparked interest in this field of research. As a result of the research, a group classification with respect to the state equation was carried out, and optimal systems of subalgebras [4,5,6,7,8,9] and submodels [10,11] were obtained for various state equations. In addition, many types of motions were investigated, for example, barochronous [12] and isothermal [13], with a linear velocity field [14,15,16]. In this paper, new exact solutions with a linear velocity field are sought.
The motion of a continuous medium with a linear velocity field had been studied since the XIX century by Dirichlet and Riemann [17,18]. In their works, they considered motions with the homogeneous deformation of an incompressible fluid. It was assumed that the liquid moves in a force field caused by the mutual attraction of particles according to Newton’s law of universal gravitation. Ovsyannikov and Dyson [19,20] independently showed that for a polytropic gas, the system of the gas dynamics equations is reduced to a system of nine ordinary differential equations of the second order. Several first integrals of the system were found. The development of the mathematical theory of these equations were obtained in the works of Andreev [21], Bogoyavlensky [22], Nemchinov [23], etc. Some of the last works were [24,25,26,27,28]. The general form of the solution to the gas dynamics equations with an arbitrary state equation with a linear velocity field with pressure and density depending on time is determined by an autonomous system of ordinary differential equations in [29]. In addition, there is the interesting paper where the author considers a solution possessing linearity with respect to a part of the space coordinates (one or two coordinates) [30].
In order to construct exact solutions, it is necessary to construct submodels, which are the gas dynamics equations written in terms of invariants. The rank of the submodel is the number of independent variables. In this paper, the gas dynamics equations are considered in the case of the special state equation (a pressure equals to the sum of two arbitrary functions of density and entropy) and in the case of the state equation of a monatomic gas. 12-dimensional Lie algebra L 12 and the 14-dimensional Lie algebra L 14 correspond to the transformation groups admitted by the gas dynamics equations with specified state equations. All subalgebras of these Lie algebras up to internal automorphisms were listed in optimal systems of nonsimilar subalgebras [6,7]. The constructed submodels and description of particle motion for small-dimensional subalgebras are given in [31,32,33,34,35,36,37]. In order to construct submodels, two three-dimensional subalgebras are selected from the Lie algebras L 12 and L 14 , and invariant submodels of rank one are constructed. Exact solutions of the gas dynamics equations with a linear velocity field are obtained from the submodels. Particle motion is investigated for the obtained exact solutions.
A different approach is also used to obtain exact solutions. If, initially, for the gas dynamics equations with an arbitrary state equation, a representation of the solution with a linear velocity field is specified, then a classification of submodels according to the state equations is obtained. In this paper, the submodel is selected which generalizes representations of velocities of obtained invariant submodels for the two types of state equations. Exact solutions are obtained and particle trajectories are plotted.

2. Symmetries of Gas Dynamics Equations and Solutions with a Linear Velocity Field

We consider the gas dynamics equations as follows [38]:
u t + ( u · ) u + ρ 1 p = 0 , ρ t + ( u · ) ρ + ρ · u = 0 , S t + ( u · ) S = 0 or p t + ( u · ) p + ρ a c 2 · u = 0 ,
where t and x are independent variables; u is a velocity; p is a pressure; ρ is a density; and S is an entropy, with an arbitrary state equation:
p = f ( ρ , S ) .
In Equation (1) in the Cartesian coordinate system, we have:
x = x i + y j + z k , = i x + j y + k z , u = u i + v j + w k ,
where i , j , and k is the orthonormal basis.
The sound velocity a c > 0 is defined by formula [38]:
a c 2 = d p d ρ | S = c o n s t = f ρ ( ρ , S ) .
The gas dynamics Equation (1) with an arbitrary state Equation (2) are invariant under the action Galilean group extended by uniform dilatation:
x = x + a ( space translations ) ; t = t + a 0 ( time translation ) ; x = O x ,   u = O u ,   O O T = E ,   det O = 1 ( rotations ) ; x = x + t b ,   u = u + b ( Galilean translations ) ; t = c t , x = c x ( uniform dilatation ) .
The transformation group (4) corresponds to an 11-dimensional Lie algebra L 11 with the basis generators in the Cartesian coordinate system [3,39]:
X 1 = x , X 2 = y , X 3 = z , X 4 = t x + u , X 5 = t y + v , X 6 = t z + w , X 7 = y z z y + v w w v , X 8 = z x x z + w u u w , X 9 = x y y x + u v v u , X 10 = t , X 11 = t t + x x + y y + z z .
In the case of special state equation [3]:
p = f ( ρ ) + h ( S )
the group (4) also includes the following transformation:
p = p + p 0 ( pressure translation ) .
The Lie algebra L 11 is extended to a 12-dimensional Lie algebra L 12 , and the generator
Y 1 = p
is added to the basis generators of L 11 (5). The Lie algebra L 12 decomposes into the direct sum of two ideals L 12 = L 11 Y 1 .
In the case of a monatomic gas [3], we have the state equation:
p = f ( S ) ρ 5 3 .
Remark 1.
The system (1), (2) admits the equivalence transformation S = K ( S ) , therefore, in what follows, we will assume in (6), (8) instead of h ( S ) , f ( S ) just S [40].
For the state Equation (8), the group (4) also includes the following transformations:
t = c 1 t , u = u c 1 , ρ = ρ c 1 3 , p = p c 1 5 ( dilatations ) ; ρ = g ρ , p = g p , S = S g 2 3 ( dilatations of thermodynamic parameters of the gas ) ; t = t 1 f t , x = x 1 f t , u = f x + ( 1 f t ) u , ρ = ( 1 f t ) 3 ρ , p = ( 1 f t ) 5 p ( projective transformation ) ;
and the Lie algebra L 11 is extended to a 14-dimensional Lie algebra L 14 . In this case, generators
X 12 = t 2 t + t x x + t y y + t z z + ( x t u ) u + ( y t v ) v + ( z t w ) w 3 t ρ ρ 5 t p p , X 13 = t t u u v v w w 3 ρ ρ 5 p p , X 14 = ρ ρ + p p 2 3 S S
are added to basis generators of L 11 (5).
A solution with a linear velocity field has the form:
u = A ( t ) x + u 0 ( t ) ,
where A ( t ) is matrix 3 × 3 , and u 0 ( t ) is a 3-dimensional vector. If u 0 = 0 , then (11) is a solution with uniform deformation. If u 0 0 , then (11) is a solution with an inhomogeneous deformation. It is assumed that the state equation is (2). Many of the solutions with a linear velocity field are obtained from invariants for 4-dimensional subalgebras. In [14], a complete classification of submodels with a solution in the form of a linear velocity field was carried out according to the types of state equations.

3. Exact Solutions of the Gas Dynamics Equations with Special State Equation

We consider the basis generators of a 3-dimensional subalgebra 3.36 from the optimal system of nonsimilar subalgebras of Lie algebra L 12 [6]:
X 1 = x , X 3 + X 4 = z + t x + u , a X 2 + b X 3 + γ Y 1 = a y + b z + γ p , a 2 + b 2 = 1 .
The coefficient γ is equal to 1 in the case of Lie algebra L 12 and equal to 0 in the case of Lie algebra L 11 . For γ = 0 , we have the case of subalgebra 3.44 from L 11 [41].
The invariants of subalgebra (12) are as follows:
t , u + b a y z , v , w , ρ , p γ y a , a 0 .
The case a = 0 gives us a partially invariant submodel of rank 2 and defect 1.
Representation of the invariant solution from (13) has the form
u = u 1 ( t ) b a y + z , v = v ( t ) , w = w ( t ) , ρ = ρ ( t ) , p = p 1 ( t ) + γ y a , S = S 1 ( t ) + γ y a , γ = 0 1 , a 2 + b 2 = 1 .
Substituting (14) into (1) and (6), we obtain the invariant submodel:
u 1 t = b a v w , v t = ρ 1 γ a , w t = 0 , ρ t = 0 , S 1 t = γ a v .
The exact solution of the system of Equations (1) and (6) from (14) and (15) is:
u = b a y + z γ b 2 a 2 ρ 0 t 2 + b a v 0 w 0 t + u 0 , v = γ a ρ 0 t + v 0 , w = w 0 , ρ = ρ 0 , p = γ y a + γ 2 2 a 2 ρ 0 t 2 γ a v 0 t p 0 , S = γ y a + γ 2 2 a 2 ρ 0 t 2 γ a v 0 t + S 0 .
On applying Galilean translations (4) with b = ( u 0 , v 0 , w 0 ) and pressure translation (7) to the solution (16), we have that u 0 = v 0 = w 0 = p 0 = 0 .
The motion of the particles is given by the following equation [38]:
d x d t = u ( x , t ) .
Integral curves of the equation (17) are the world lines of particles in space R 4 ( t , x ) , the projection of which in R 3 ( x ) are particle trajectories. From Equations (16) and (17), world lines have the form:
x = ( a C 3 b C 2 ) t + a C 1 , y = γ 2 a ρ 0 t 2 + a C 2 , z = a C 3 .
From the generators of subalgebra (12) and Lie equations [1]
x ¯ a = ξ ( x ¯ ) , x ¯ | a = 0 = x
it follows that:
x ¯ = x + ξ ; x ¯ = x + η t , z ¯ = z + η , u ¯ = u + η ; y ¯ = y + a ζ , z ¯ = z + b ζ , p ¯ = p + ζ .
Applying the superposition of transformations (20) to Equation (18) with parameters
ξ = a C 1 , η = b C 2 a C 3 , ζ = C 2 ,
we obtain:
x = 0 , y = γ 2 a ρ 0 t 2 , z = 0 .
The trajectories (21) describe Oy axis with y 0 for a > 0 . From this, we can reproduce the following formulas using (20):
x = η t + ξ , y = γ 2 a ρ 0 t 2 + a ζ , z = b ζ + η .
Jacobian J = | x / x 0 | [41] of transformation (22) is equal to a 0 , where x 0 = ( ξ , ζ , η ) . This means that the world lines of particles do not intersect. The vortex is ω = ( w y v z , u z w x , v x u y ) = ( 0 , 1 , b / a ) [38]. The motion of particles is vortex and is restricted by moving plane y = γ 2 a ρ 0 t 2 for p 0 . The pressure is constant along the world line and has the value p = γ ζ .
For η = 0 , the trajectories (22) are rays. If the t variable is eliminated in (22) for η 0 , we obtain the following formulas of trajectories:
y = γ ( x ξ ) 2 2 a ρ 0 η 2 + a ζ , z = b ζ + η .
The trajectories (23) are parabolas in each plane z = c o n s t , in which the vertex of the parabola is the point ( ξ , a ζ ) . The parabolas (22) are illustrated in the Figure 1. For subalgebra 3.44 from L 11 , the world lines of particles are straight lines in a uniformly moving coordinate system [41].
Let the formulas (22) satisfy the initial conditions
x ( t 0 ) = x 0 , y ( t 0 ) = y 0 , z ( t 0 ) = z 0 ,
where x 0 , y 0 , z 0 are local Lagrangian coordinates. In this case, the formulas (22) have the form:
x = z 0 b a y 0 γ b t 0 2 2 a 2 ρ 0 t + x 0 t 0 z 0 + b a t 0 y 0 + γ b t 0 3 2 a 2 ρ 0 , y = γ t 2 2 a ρ 0 + y 0 + γ t 0 2 2 a ρ 0 , z = z 0 .
From now on, we make the following assumption: b = 0 , then a = 1 . Let the particles (24) be on the sphere with radius r and centre ( x 1 , y 1 , z 1 ) at time t = t 0 :
( x 0 x 1 ) 2 + ( y 0 y 1 ) 2 + ( z 0 z 1 ) 2 = r 2 .
We obtain from Equations (24) and (25) location of the particles
x 2 + y 2 + ( t t 0 ) 2 + 1 z 2 2 x z ( t t 0 ) 2 x 1 x + A y + B z = C , A = γ ( t 2 t 0 2 ) ρ 0 2 y 1 , B = 2 x 1 ( t t 0 ) 2 z 1 , C = r 2 x 1 2 y 1 2 z 1 2 + γ ( t 2 t 0 2 ) ρ 0 y 1 γ 2 ( t 2 t 0 2 ) 2 4 ρ 0 2 .
Let us rotate the coordinate axes by the α angle
x = x ˜ cos α z ˜ sin α , z = x ˜ sin α + z ˜ cos α .
After substituting (27) into (26), we find the angle α equating the coefficient at x ˜ z ˜ to zero
tan 2 α = 2 t t 0 α = 1 2 arctan 2 t t 0 + π 2 n , n Z .
If t , then α 0 . Let us choose α = 1 2 arctan 2 t t 0 from (28), then Equation (26) may be written in the following form:
A 1 x ˜ 2 + y 2 + B 1 z ˜ 2 + A 2 x ˜ + B 2 z ˜ + A y = C , A 1 = ( t t 0 ) 2 sin 2 α ( t t 0 ) sin 2 α + 1 , B 1 = ( t t 0 ) 2 cos 2 α + ( t t 0 ) sin 2 α + 1 , A 2 = B sin α 2 x 1 cos α , B 2 = B cos α + 2 x 1 sin α ; sin α = 1 2 1 t t 0 ( t t 0 ) 2 + 4 , cos α = 1 2 1 + t t 0 ( t t 0 ) 2 + 4 , sin 2 α = 2 ( t t 0 ) 2 + 4 , cos 2 α = t t 0 ( t t 0 ) 2 + 4 , sin 2 α = 1 2 t t 0 2 ( t t 0 ) 2 + 4 , cos 2 α = 1 2 + t t 0 2 ( t t 0 ) 2 + 4 .
From Equation (29), the canonical form of a second-order surface has the following form:
x ˜ + A 2 2 A 1 2 K A 1 + y + A 2 2 K + z ˜ + B 2 2 B 1 2 K B 1 = 1 , K = C + A 2 2 4 A 1 + B 2 2 4 B 1 + A 2 4 .
Equation (30) defines the ellipsoid at t t 0 (Figure 2). Calculations by the MAPLE computer mathematics system give us the following limits:
lim t K A 1 = + , lim t K = r , lim t K B 1 = 0 , lim t A 2 2 A 1 = sign ( z 1 ) , lim t B 2 2 B 1 = 0 , lim t A 2 = sign γ ρ 0 .
Thus, according to formulas (31), if t , the ellipsoid (30) turns into the plane strip. The volume of the ellipsoid (30) is V = 4 3 π r 3 at t .

4. Exact Solutions of the Gas Dynamics Equations with Monatomic Gas State Equation

We consider the basis generators of 3-dimensional subalgebra 3.4* from the optimal system of nonsimilar subalgebras of Lie algebra L 14 for Equation (8) [42]:
X 3 + X 5 = z + t y + v , a ( X 2 + X 6 ) + X 7 + X 10 + X 12 = = a ( y + t z + w ) + y z z y + v w w v + ( 1 + t 2 ) t + + t x x + t y y + t z z + ( x t u ) u + ( y t v ) v + ( z t w ) w 3 t ρ ρ 5 t p p , X 2 + X 6 + b X 14 = y + t z + w + b ( ρ ρ + p p 2 3 S S ) .
Invariants have the following form:
x 1 = x 1 + t 2 , u 1 + t 2 t x 1 + t 2 , v + t w 2 y t + z ( t 2 1 ) 1 + t 2 , w t v 2 z t y ( t 2 1 ) 1 + t 2 , ρ ( 1 + t 2 ) 3 2 e x p b ( a τ y + t z 1 + t 2 ) , p ( 1 + t 2 ) 5 2 e x p b ( a τ y + t z 1 + t 2 ) , τ = arctan t .
The representation of the invariant solution is given as follows:
u = u 1 ( x 1 ) 1 + t 2 + t x 1 + t 2 , v = v 1 ( x 1 ) w 1 ( x 1 ) t + y t z 1 + t 2 , w = w 1 ( x 1 ) + v 1 ( x 1 ) t + y + z t 1 + t 2 , ρ = ρ 1 ( x 1 ) ( 1 + t 2 ) 3 2 e x p b ( a τ + y + t z 1 + t 2 ) , p = p 1 ( x 1 ) ( 1 + t 2 ) 5 2 e x p b ( a τ + y + t z 1 + t 2 ) , S = S 1 ( x 1 ) e x p 2 3 b ( a τ + y + t z 1 + t 2 ) , S 1 ( x 1 ) = p 1 ( x 1 ) ρ 1 ( x 1 ) 5 3 .
Substituting (33) in (1) and (8), we obtain the invariant submodel:
u 1 u 1 x 1 + p 1 x 1 ρ 1 = x 1 , u 1 v 1 x 1 = 2 w 1 b p 1 ρ 1 , u 1 w 1 x 1 = 2 v 1 , u 1 ρ 1 x 1 + ρ 1 u 1 x 1 = b ρ 1 ( a v 1 ) , u 1 S 1 x 1 = 2 3 b S 1 ( a v 1 ) , S 1 = p 1 ρ 1 5 3 .
The system (34) has an integral ρ 1 u 1 S 1 3 2 = p 1 ρ 1 3 2 u 1 = K 0 and admits reflection a a , b b , v 1 v 1 , w 1 w 1 .
If the condition u 1 = 0 is satisfied, then there is a solution:
a b = 0 , u 1 = 0 , v 1 = 0 , w 1 = b p 1 2 ρ 1 , ρ 1 = p 1 x 1 x 1 ,
where p 1 ( x 1 ) is an arbitrary function.
Let us distinguish among them solutions with a linear velocity field w 1 = b x 1 p 1 2 p 1 x 1 = a 3 x 1 + b 3 :
1 )   b = 0 , u 1 = 0 , v 1 = 0 , w 1 = 0 , ρ 1 = p 1 x 1 x 1 ; p 1 ( x 1 ) is arbitrary function ,
2 )   a = 0 , b b 3 > 0 , u 1 = 0 , v 1 = 0 , w 1 = b 3 , ρ 1 = ρ 0 e b 4 b 3 x 1 2 , p 1 = 2 b 3 b ρ 0 e b 4 b 3 x 1 2 ;
3 )   a = 0 , a 3 0 , b ( a 3 x 1 + b 3 ) > 0 , u 1 = 0 , v 1 = 0 , w 1 = a 3 x 1 + b 3 , ρ 1 = ρ 0 | a 3 x 1 + b 3 | b b 3 2 a 3 2 1 e b 2 a 3 x 1 , p 1 = 2 | b | ρ 0 | a 3 x 1 + b 3 | b b 3 2 a 3 2 e b 2 a 3 x 1 .
If u 1 0 and we assume u 1 = ( u 1 , v 1 , w 1 ) = ( a 1 x 1 + b 1 , a 2 x 1 + b 2 , a 3 x 1 + b 3 ) with a linear velocity field, then there are no solutions.
Solution (35) produces solution to the gas dynamics Equations (1) and (8):
u = t x 1 + t 2 , v = y t z 1 + t 2 , w = y + z t 1 + t 2 , ρ = p 1 x 1 x 1 ( 1 + t 2 ) 3 2 , p = p 1 ( x 1 ) ( 1 + t 2 ) 5 2 , S = p ρ 5 3 ; p 1 ( x 1 ) is arbitrary function of x 1 = x 1 + t 2 ,
which was investigated in the work [43].

Analysis of the Exact Solution

Invariant solution (36) produces exact solution to the gas dynamics Equations (1) and (8):
u = t x 1 + t 2 , v = y t z b 3 t 1 + t 2 , w = y + z t + b 3 1 + t 2 , ρ = ρ 0 ( 1 + t 2 ) 3 2 e x p b 1 + t 2 ( x 2 4 b 3 + y + t z ) , p = 2 b 3 ρ 0 b ( 1 + t 2 ) 5 2 e x p b 1 + t 2 ( x 2 4 b 3 + y + t z ) , S = 2 b 3 b ρ 0 2 3 e x p 2 b 3 ( 1 + t 2 ) ( x 2 4 b 3 + y + t z ) , b b 3 > 0 .
The vortex is ω = ( w y v z , u z w x , v x u y ) = 2 1 + t 2 , 0 , 0 . The motion of the particles is vortex.
Equations (17) and (39) give the world lines:
x ( t ) = x 0 1 + t 2 , y ( t ) = ( b 3 τ + z 0 ) t + y 0 , z ( t ) = y 0 t + b 3 τ + z 0 ,
which are space curves for x 0 0 .
Here, x 0 , y 0 , z 0 are the particle coordinates at the moment t = 0 (the global Lagrange coordinates). The Jacobian of transformation (40) is J = ( 1 + t 2 ) 3 2 0 . The motion of the particles occurs without blow-up. The particles lying in the plane x = 0 do not leave it (plane curves). The motion of particles occurs in half-spaces symmetrically relative to the plane x = 0 .
The velocity components (39) along world lines (40) are as follows:
u = t x 0 1 + t 2 , v = z 0 b 3 τ b 3 t 1 + t 2 , w = y 0 + b 3 1 + t 2 .
At t = 0 , we have u ( x 0 ; y 0 ; z 0 ) = ( 0 ; z 0 ; y 0 + b 3 ) ; at t ± , we have u ( x 0 ; y 0 ; z 0 ) = ( ± x 0 ; z 0 b 3 π 2 ; y 0 ) .
We introduce new axes in the plane O y z :
ξ = y + z t 1 + t 2 , η = z t y 1 + t 2 ,
which, at the time t = 0 , coincide with the axes O y and O z ( ξ 0 = y 0 and η 0 = z 0 , respectively), and at any other time t, they are rotated by an angle τ = arctan t counter-clockwise relative to the axes O y and O z . If t ± , then ξ ± z , η y , τ ± π 2 .
Thus, the formulas (40) are given by:
x ( t ) = x 0 1 + t 2 , ξ ( t ) = ξ 0 1 + t 2 , η ( t ) = ( η 0 + b 3 τ ) 1 + t 2 .
In the space ( x , ξ , η ) , trajectories (41) are plane curves. In polar coordinates x = r cos ϕ , ξ = r sin ϕ , Equations (41) have the following form:
r ( t ) = r 0 1 + t 2 , ϕ ( t ) = ϕ 0 , η ( t ) = ( η 0 + b 3 τ ) 1 + t 2 .
In the plane ( r , η ) , the trajectories (42) are given by the equation η = η 0 ± b 3 arccos r 0 r r r 0 , which define the curves with oblique asymptotes η a = η 0 + 0.5 sign ( t ) b 3 π r 0 r b 3 sign ( t ) at t ± . If we consider the curve as r = r ( η ) , then it has a derivative r η = r t η t = r 0 t t η 0 + b 3 ( 1 + t τ ) and an extremum at t = 0 in the point ( η 0 , r 0 ) . For any ϕ 0 , its graph is shown in the Figure 3.
Thus, the motion of particles can be represented as a complex motion consisting of the condensation of gas particles to the origin, followed by scattering and successive (uneven) rotation in the plane O y z by an angle π . Figure 4 represents the trajectories of the same particles in spaces ( x , ξ , η ) (or ( r , ϕ , η ) ) and ( x , y , z ) , showing how the rotation of the coordinate system around the axis O x changes the trajectories particles motion.
If the particle has coordinates ( x t 0 ; y t 0 ; z t 0 ) at the moment of time t = t 0 , then the Equation (40) may be rewritten as follows:
x ( t ) = x t 0 1 + t 0 2 1 + t 2 , y ( t ) = b 3 t ( τ τ 0 ) z t 0 y t 0 t 0 1 + t 0 2 t + y t 0 + z t 0 t 0 1 + t 0 2 , z ( t ) = b 3 ( τ τ 0 ) + y t 0 + z t 0 t 0 1 + t 0 2 t + z t 0 y t 0 t 0 1 + t 0 2 .
If at t = t 0 the particle with coordinates ( x t 0 , y t 0 , z t 0 ) is on a sphere with a center ( x 1 , y 1 , z 1 ) and radius R 0 > 0
( x t 0 x 1 ) 2 + ( y t 0 y 1 ) 2 + ( z t 0 z 1 ) 2 = R 0 2 ,
then at time moment t, it will be on the sphere
x x 1 1 + t 0 2 1 + t 2 2 + y + b 3 t ( τ τ 0 ) + z 1 y 1 t 0 1 + t 0 2 t y 1 + z 1 t 0 1 + t 0 2 2 + + z b 3 ( τ τ 0 ) y 1 + z 1 t 0 1 + t 0 2 t z 1 y 1 t 0 1 + t 0 2 2 = R 0 2 1 + t 0 2 ( 1 + t 2 ) ,
which has a radius R 0 1 + t 0 2 1 + t 2 and the center moves according to the law of motion (43).
The sound surface is given by the equation u 2 + v 2 + w 2 = a c 2 [38] a c 2 = 5 p 3 ρ = 10 b 3 3 b ( 1 + t 2 ) from ( 3 ) . It has the following form from (39):
t 2 x 2 1 + t 2 + y + b 3 1 t 2 1 + t 2 2 + z + b 3 2 t 1 + t 2 2 = 10 b 3 3 b .
For t 0 , Equation (45) defines an ellipsoid with a center 0 ; b 3 1 t 2 1 + t 2 ; b 3 2 t 1 + t 2 , one semiaxis is equal to 10 b 3 ( 1 + t 2 ) 3 b t 2 , and the other two semiaxes are equal to r c = 10 b 3 3 b .
For t = 0 , Equation (45) defines a circular cylinder ( y + b 3 ) 2 + z 2 = 10 b 3 3 b with radius r c and an axis parallel to the O x , passing through a point ( 0 ; b 3 ; 0 ) .
For t ± , Equation (45) defines a sphere x 2 + ( y b 3 ) 2 + z 2 = 10 b 3 3 b with a center ( 0 , b 3 , 0 ) and radius r c .
The motion of the volume of particles (44) has been considered. The surface of the volume of particles (sphere) is the contact characteristic of the gas dynamics equations. For a normal gas, the gas dynamics equations also have two sound characteristics C ± .
We consider the equations of sound characteristics [38]:
C ± : h t + u h x + v h y + w h z ± a h x 2 + h y 2 + h z 2 = 0 .
The equations of the bicharacteristics have the following form:
d x d t = u ± a h | h | ,
d h j d t = u j · h a j | h | , j = t , x , y , z .
Equations (39) and (47) arrive at the following formulas:
h x = h x 0 1 + t 2 , h y = h y 0 t h z 0 1 + t 2 , h z = h z 0 + t h y 0 1 + t 2 , | h | = | h 0 | 1 + t 2 , | h 0 | = h x 0 2 + h y 0 2 + h z 0 2 .
Equations (39), (46), and (48) give:
x = ± 10 b 3 3 b h x 0 | h 0 | ( τ τ 0 ) + x t 0 1 + t 0 2 1 + t 2 , y = y t 0 + t 0 z t 0 1 + t 0 2 + t 0 y t 0 z t 0 1 + t 0 2 t b 3 t ( τ τ 0 ) ± 10 b 3 3 b h y 0 t h z 0 | h 0 | ( τ τ 0 ) , z = z t 0 t 0 y t 0 1 + t 0 2 + y t 0 + t 0 z t 0 1 + t 0 2 t + b 3 ( τ τ 0 ) ± 10 b 3 3 b t h y 0 + h z 0 | h 0 | ( τ τ 0 ) .
Let us construct a characteristic conoid K ( P ) with a vertex P—a characteristic surface formed by all the bicharacteristics coming out of a given point P ( x t 0 , t 0 ) . In order to represent all such bicharacteristics, it is necessary to take into account the initial conditions to the system (46), (47) [38]:
x ( t 0 ) = x t 0 , h j ( t 0 ) = h j o ( j = t , x , y , z ) .
In order to satisfy (50), it is necessary to make a replacement h x 0 h x 0 1 + t 0 2 ; h y 0 h y 0 + t 0 h z 0 ; h z 0 h z 0 t 0 h y 0 ; | h 0 | | h 0 | 1 + t 0 2 in (48), (49), but still, we exclude h x 0 , h y 0 , h z 0 , | h 0 | from (49) and obtain the following equation:
x x t 0 1 + t 0 2 1 + t 2 2 + y y t 0 + t 0 z t 0 1 + t 0 2 t 0 y t 0 z t 0 1 + t 0 2 t + b 3 t ( τ τ 0 ) 2 + + z z t 0 t 0 y t 0 1 + t 0 2 y t 0 + t 0 z t 0 1 + t 0 2 t b 3 ( τ τ 0 ) 2 = 10 b 3 3 b ( τ τ 0 ) 2 ( 1 + t 2 ) .
The section of the conoid (51) by hyperplane t = c o n s t is a sphere in R 3 ( x ) , the center of which moves according to the law of motion (43), and the radius is equal to 10 b 3 3 b ( 1 + t 2 ) | τ τ 0 | .
At t , the sound velocity a c 0 and sound characteristic tends to the contact characteristic. Indeed, the conoid (51) tends to the sphere (44), which has a radius equal to R 0 = 10 b 3 3 b ( 1 + t 0 2 ) | π 2 τ 0 | . Thus, for any moment of time t > 0 , the characteristic conoid with a vertex P ( x t 0 , t 0 ) beginning from the center of a spherical volume with a radius R 0 = 10 b 3 3 b ( 1 + t 0 2 ) | π 2 τ 0 | never reaches particles outside this volume. Figure 5 shows the propagation of a characteristic conoid (51) with a vertex P ( x t o , t 0 ) and the motion of a spherical volume of gas (44) with the center at the same point at t = t 0 and radius R 0 .
Depending on where the point P ( x t 0 , t 0 ) is located (in the subsonic region, on the sound surface, or in the supersonic region), the sections of the characteristic conoid by planes t = c o n s t are nested into each other, touch at the point, or intersect each other, respectively (Figure 6).
At t = c o n s t , the pressure and density have the same values on parabolic cylinders.

5. Generalization of Exact Solutions of the Gas Dynamics Equations with a Linear Velocity Field

We consider the solution of the gas dynamics equations with a linear velocity field (11). Let us substitute the velocity representation (11) into the gas dynamics Equations (1) and (2) and express the time derivative and the gradient of the pressure function. On equating the shifted derivatives, we obtain a classifying relation for the matrix B = A + A 2 . Matrix B is decomposed into the sum of symmetric S and antisymmetric parts. We obtain one or another submodel, depending on the rank of the matrix S . It is shown in [14] that there are 11 such submodels. Let us choose one of them so that the state Equations (6) and (8) satisfy the state equation of the selected submodel:
S + 2 S A = ( 1 γ ) S tr A , A + A 2 = S , S = S T , v + A T v + S u 0 = ( 1 γ ) v tr A , u 0 + A u 0 = v , ρ = 2 e tr A d t R ( I ) , I = ( x · S x + 2 v · x ) e ( 1 γ ) tr A d t 2 u 0 · v e ( 1 γ ) tr A d t d t , p = ρ γ h 0 ( S ) + P 0 1 ρ γ γ ,
where γ , P 0 are arbitrary constants; R, h 0 are arbitrary function; tr A is matrix trace A.

5.1. Generalization of Exact Solution of the Gas Dynamics Equations with Special State Equation

We consider an invariant solution (14) for the velocity components from subalgebra 3.36 (12), then matrix A and vector u 0 from (11) have the form:
A = 0 b a 1 0 0 0 0 0 0 , u 0 = u 1 ( t ) v 1 ( t ) w 1 ( t ) .
The state equation from (52) at γ 0 tends to:
p = h 0 ( S ) P 0 ln ρ .
Since S = A + A 2 0 , substituting the matrix A into the differential equations of the submodel (52), the differential equation for the matrix S is fulfilled. After substituting the matrices A and S into the differential equation for vector u 0 ( t ) from (52), we obtain a system of six-order differential equations for unknown functions u 1 ( t ) , v 1 ( t ) , w 1 ( t ) :
u 1 b a v 1 + w 1 = 0 ,
v 1 b a v 1 + b 2 a 2 v 1 b a w 1 = 0 ,
w 1 + u 1 b a v 1 + w 1 = 0 .
The solution of the last system has the form:
u 1 = 1 a 3 u ˜ 03 t 3 6 + u ¯ 04 a b u ¯ 01 t 2 2 + ( u ¯ 05 u ¯ 03 a 2 ) t + u ¯ 06 a ,
v 1 = b a 2 u ˜ 03 t 2 2 + u ¯ 01 a b a 2 u ¯ 04 t + u ¯ 02 b a 2 u ¯ 05 , w 1 = 1 a u ¯ 04 t + u ˜ 03 t 2 2 + u ¯ 05 ,
where u ¯ 01 , u ¯ 02 , u ¯ 03 , u ¯ 04 , u ¯ 05 , u ¯ 06 are arbitrary constants; u ˜ 03 = b u ¯ 02 u ¯ 03 .
Let us write the components of the velocity vector using the last formulas for u 1 , v 1 , and w 1 (14) and apply the Galilean translations with b = ( u ¯ 06 a 1 , u ¯ 05 b a 2 , u ¯ 05 a 1 ) and dilatation from (4) with c = a :
u = a z b y u ˜ 03 t 3 6 u ¯ 04 a b u ¯ 01 t 2 2 + u ¯ 03 t , v = b u ˜ 03 t 2 2 + u ¯ 01 b a u ¯ 04 t + u ¯ 02 , w = u ¯ 04 t + a u ˜ 03 t 2 2 .
The density function from the formulas (52) has the form:
ρ = 2 R ( I ) ,
where variable I has the form:
I = 2 u ˜ 03 x 2 b a ( a u ˜ 03 t + u ¯ 04 ) u ¯ 01 y + 2 ( a u ˜ 03 t + u ¯ 04 ) z u ˜ 03 2 3 t 4 + + 4 3 u ˜ 03 a ( a b u ¯ 01 u ¯ 04 ) t 3 u ¯ 03 2 b 2 u ¯ 02 2 + u ¯ 04 2 + ( u ¯ 01 b a u ¯ 04 ) 2 t 2 2 u ¯ 02 u ¯ 01 b a u ¯ 04 t .
We need to solve the system of differential equations for the function p from the gas dynamics equations in order to define the function of pressure:
p = ρ ( S x + v ) , p t = ρ ( A x + u 0 ) · ( S x + v ) ρ a c 2 ( p , ρ ) tr A .
Substituting the known function ρ , matrix S, and vector v into the first equation in (57), we obtain:
p = 2 R ( I ) v
or
p = R ( I ) I .
The integration one of the two last differential equations gives us the form of the pressure:
p = R ( I ) + p ¯ 0 ( t ) ,
where p ¯ 0 is an arbitrary function. Without loss of generality, we can assume in (58) that p ¯ 0 ( t ) is equal to zero because after substituting (58) in the second differential equation from (57), we find that p ¯ 0 ( t ) is constant p ¯ 0 . Then, the pressure (58) has the form:
p = R ( I ) .
Thus, the solution of the gas dynamics Equations (1) with the state Equation (53) is given by the formulas (54), (55), (56), and (59).
Remark 2.
The solution (54), (55), (56), and (59) coincides with the solution (16) after using uniform dilatation from (4) with c = a and if u ¯ 01 = γ ρ 0 1 , u ¯ 02 = u ¯ 03 = u ¯ 04 = u ˜ 03 = 0 , R ( I ) = 2 1 ρ 0 I , and P 0 = 0 .
The world lines of gas particles (17) are defined from solution (54):
x = u ˜ 03 t 2 2 + ( a z 0 b y 0 ) t + x 0 , y = b u ˜ 03 6 t 3 + ( u ¯ 01 b a u ¯ 04 ) t 2 2 + u ¯ 02 t + y 0 , z = a u ˜ 03 6 t 3 + u ¯ 04 2 t 2 + z 0 ,
where x 0 , y 0 , and z 0 are global Lagrangian coordinates. The Jacobian J = | x / x 0 | [41] of transformations (60) is equal to 1, where x 0 = ( x 0 , y 0 , z 0 ) . This fact means that the motion of particles does not have singularities. Trajectories of the gas particles (60) are built in the Figure 7.

5.2. Generalization of Exact Solution of the Gas Dynamics Equations with State Equation of the Monatomic Gas

We consider the representation of invariant solution (33) for subalgebra 3.4* (32) when the functions u 1 , v 1 , and w 1 are linear
u 1 = a 1 x 1 + t 2 + b 1 , v 1 = a 2 x 1 + t 2 + b 2 , w 1 = a 3 x 1 + t 2 + b 3 ,
where a 1 , a 2 , a 3 , b 1 , b 2 , and b 3 are arbitrary constants. Then, matrix A and u 0 from submodel (52) have the form:
A = a 1 + t 1 + t 2 0 0 a 2 a 3 t ( 1 + t 2 ) 3 / 2 t 1 + t 2 1 1 + t 2 a 3 + a 2 t ( 1 + t 2 ) 3 / 2 1 1 + t 2 t 1 + t 2 , u 0 = b 1 1 + t 2 b 2 t b 3 1 + t 2 b 3 + b 2 t 1 + t 2 .
The state equation from (52) with γ = 5 3 , P 0 = 0 is equal to p = ρ 5 / 3 h 0 ( S ) . Matrix A + A 2 has the form:
A + A 2 = a 1 2 + 1 ( 1 + t 2 ) 2 0 0 2 a 2 t + a 2 a 1 2 a 3 a 3 a 1 t ( 1 + t 2 ) 5 / 2 0 0 2 a 2 2 a 3 t + a 3 a 1 + a 2 a 1 t ( 1 + t 2 ) 5 / 2 0 0
From the symmetry condition of the matrix A + A 2 , we obtain that a 2 = a 3 = 0 . After substituting matrices S and A into the submodel Equation (52), we obtain a relation from which it follows that a 1 = 0 .
Next, we substitute the known matrices into a vector differential equation for vector u 0 from the submodel (52). The vector form u 0 is taken from (61). From the three obtained relations, we find that b 1 = 0 and b 2 and b 3 are arbitrary constants. The components of the velocity vector have the form:
u = t x 1 + t 2 , v = t ( y b 3 ) z + b 2 1 + t 2 , w = y + b 3 + t ( z + b 2 ) 1 + t 2 .
Let us define the type of density function from (52). We receive:
ρ = 2 R ( I ) ( 1 + t 2 ) 3 / 2 ,
where
I = x 2 4 y ( b 3 + t b 2 ) 4 z ( b 3 t b 2 ) 1 + t 2 .
From (57), we obtain the pressure
p = R ( I ) ( 1 + t 2 ) 5 / 2 + p ¯ 0 ( t ) .
Function p ¯ 0 ( t ) is defined from the second Equation (57), where a c 2 ( p , ρ ) = 5 3 p ρ 1 from (3). As a result, the pressure function is:
p = R ( I ) ( 1 + t 2 ) 5 / 2 .
Thus, the solution of the gas dynamics Equations (1) with the state Equation (8) is given by the formulas (62), (63), (64), and (65).
Remark 3.
The solution (62), (63), (64), and (65) coincides with the solution (38) if b 2 = b 3 = 0 , R ( I ) = p ( x 1 ) , I = x 1 2 .
Remark 4.
The solution (62), (63), (64), and (65) coincides with the solution (39) if b 2 = 0 , R ( I ) = 2 b 3 ρ 0 b e b 4 b 3 I .
Remark 5.
The solution (62), (63), (64), and (65) does not coincide with the solution (33), (37) because the matrix B = A + A 2 has an antisymmetrical part which contradicts the formulas (52).
The world lines of gas particles (17) are defined from solution (62):
x = x 0 1 + t 2 , y = ( b 2 t b 3 ) τ z 0 t + y 0 , z = ( t b 2 + b 3 ) τ + y 0 t + z 0 ,
where x 0 , y 0 , and z 0 are global Lagrangian coordinates.
The Jacobian J = | x / x 0 | [41] of transformations (66) is equal to ( 1 + t 2 ) 3 2 0 . At the initial moment of time t = 0 , the gas particle is located at the point with coordinates ( x 0 , y 0 , z 0 ) . The trajectories of the gas particles (66) are built in the Figure 8.

6. Conclusions

In this work, the gas dynamics equations have been considered in the case of the state equation of a monatomic gas and the special state equation (a pressure equals to the sum of two arbitrary functions of density and entropy). One three-dimensional subalgebra has been selected for each specified state equation. Invariant submodels of rank one have been constructed for them, and exact solutions with a linear velocity field with inhomogeneous deformation have been obtained.
In the case of the special state equation, the submodel describes the isochoric motion of particles with constant pressure along each world line. The motion is vortex, limited by a moving plane. The motion of particles occurs along parabolas and along rays in parallel planes. The motion of the particles volume in the case of the subalgebra without one parameter is considered. The spherical volume of particles turns into an ellipsoid at finite moments of time, and as time tends to infinity, the particles end up on an infinite strip with a width equals to the diameter of the original sphere.
In the case of the state equation of a monatomic gas, the submodel describes a complex vortex condensation (compaction) to the origin and subsequent expansion of gas particles in half-spaces relative to the plane x = 0 . The motion of any allocated spherical volume of gas retains a spherical shape and expands (narrows) for t > 0 (for t < 0 ). An equation defining a characteristic conoid with an arbitrary vertex in R 4 is obtained. Its projection into R 3 at each moment of time gives the equation of the sphere. It is shown that for any point R 3 at time t = t 0 > 0 , it is possible to obtain a radius of the sphere depending on t 0 . The perturbation coming from the center does not reach the particles outside this sphere.
The obtained exact solutions with a linear velocity field are generalized if the representations for the velocities are chosen according to the considered invariant solutions. From the previously performed classification of submodels, a submodel with a state equation including the above is selected. As a result, exact solutions of a wider class have been obtained, since density and pressure may not be invariant with respect to the generators of the selected three-dimensional subalgebras.

Author Contributions

Conceptualization, R.N., D.S. and Y.Y.; investigation, R.N., D.S. and Y.Y.; writing—original draft preparation, R.N., D.S. and Y.Y.; writing—review and editing, R.N., D.S. and Y.Y. All authors have read and agreed to the published version of the manuscript.

Funding

The authors were supported by the Russian Foundation for Basic Research (project No. 18-29-10071) and partially from the Federal Budget by the State Target (project No. 0246-2019-0052).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ovsyannikov, L.V. Group Analysis of Differential Equations; Academic Press: New York, NY, USA, 1982. [Google Scholar]
  2. Bluman, G.W.; Kumei, S. Symmetries and Differential Equations; Springer: New York, NY, USA, 1989. [Google Scholar] [CrossRef]
  3. Ovsyannikov, L.V. The “podmodeli” program. Gas dynamics. J. Appl. Math. Mech. 1994, 58, 601–627. [Google Scholar] [CrossRef]
  4. Khabirov, S.V. Optimal system for sum of two ideals admitted by hydrodynamic type equations. Ufa Math. J. 2014, 6, 97–101. [Google Scholar] [CrossRef] [Green Version]
  5. Makarevich, E.V. Optimal system of subalgebras admitted by the gas dynamics equations in case of state equation with separated density. Sib. Electron. Math. Rep. 2011, 8, 19–38. (In Russian) [Google Scholar]
  6. Siraeva, D.T. Optimal system of non-similar subalgebras of sum of two ideals. Ufa Math. J. 2014, 6, 90–103. [Google Scholar] [CrossRef] [Green Version]
  7. Cherevko, A.A. The Optimal System of Subalgebras for Lie Algebra of Generators Admitted by System of Gas Dynamics Equations with the State Equation p = f(S)ρ5/3, Preprint 4–96; Institute of Hydrodynamics: Novosibirsk, Russia, 1996; 39p. (In Russian) [Google Scholar]
  8. Golovin, S.V. The Optimal System of Subalgebras for the Lie Algebra of Generators Admitted by the Gas Dynamics Equations in the Case of a Polytropic Gas, Preprint 5–96; Institute of Hydrodynamics: Novosibirsk, Russia, 1996; 33p. (In Russian) [Google Scholar]
  9. Khabirov, S.V. Optimal Systems of Subalgebras Admitted by Equations of Gas Dynamics Preprint; Institute of Mechanics of the UNC RAS: Ufa, Russia, 1998; 33p. (In Russian) [Google Scholar]
  10. Mamontov, E.V. Invariant submodels of rank two of the equations of gas dynamics. J. Appl. Mech. Tech. Phys. 1999, 40, 232–237. [Google Scholar] [CrossRef]
  11. Khabirov, S.V. Group Analysis of the Plane Steady Vortex Submodel of Ideal Gas with Varying Entropy. Mathematics 2021, 9, 2006. [Google Scholar] [CrossRef]
  12. Chupakhin, A.P. Barokhronic Gas Motions. General Properties and Submodels Type (1.2) and (1.1), Preprint; Institute of Hydrodynamics: Novosibirsk, Russia, 1998; 40p. (In Russian) [Google Scholar]
  13. Khabirov, S.V. The plane isothermal motions of an ideal gas without expansions. J. Appl. Math. Mech. 2014, 78, 287–297. [Google Scholar] [CrossRef]
  14. Yulmukhametova, Y.V. Submodels in gas dynamics with linear field of velocity. Sib. Electron. Math. Rep. 2012, 9, 208–226. (In Russian) [Google Scholar]
  15. Tarasova, Y.V. Classification of submodels with a linear velocity field in gas dynamics. J. Appl. Ind. Math. 2010, 4, 570–577. [Google Scholar] [CrossRef]
  16. Urazbakhtina, L.Z. Integrable hydrodynamic submodels with a linear velocity field. J. Appl. Ind. Math. 2013, 7, 117–126. [Google Scholar] [CrossRef]
  17. Dirichlet, G.L. Untersuchunger uber eih Problem der Gydrodynamik. J. Reine Angev. Math. 1860, 58, 181. [Google Scholar]
  18. Riemann, B. Collected Works; GITTL: Moscow, Russia; Leningrad, Russia, 1948; pp. 339–366. (In Russian) [Google Scholar]
  19. Ovsiannikov, L.V. New solution of hydrodynamic equations. Dokl. Akad. Nauk SSSR 1956, 3, 47–49. [Google Scholar]
  20. Dyson, J.F. Dynamics of a spinning gas cloud. J. Math. Mech. 1968, 18, 91–101. [Google Scholar] [CrossRef]
  21. Andreev, V.K. Unsteady motion of a gas jet with a linear field of velocity. Sib. J. Ind. Math. 2002, 5, 23–35. [Google Scholar]
  22. Bogoyavlensky, O.I. Methods in the Qualitative Theory of Dynamical Systems in Astrophysics and Gas Dynamics; Springer Series in Soviet Mathematics; Springer: Berlin/Heidelberg, Germany, 1985. [Google Scholar]
  23. Nemchinov, I. Expansion of a tri-axial gas ellipsoid in a regular behavior. J. Appl. Math. Mech. 1965, 29, 143–150. [Google Scholar] [CrossRef]
  24. Rosensteel, G. Rapidly rotating nuclei as Riemann ellipsoids. Ann. Phys. 1988, 186, 230–291. [Google Scholar] [CrossRef]
  25. Borisov, A.V.; Kilin, A.; Mamaev, I. The Hamiltonian dynamics of self-gravitating liquid and gas ellipsoids. Regul. Chaotic Dyn. 2009, 14, 179–217. [Google Scholar] [CrossRef]
  26. Giron, J.F.; Ramsey, S.D.; Baty, R.S. Nemchinov-Dyson Solutions of the Two-Dimensional Axisymmetric Inviscid Compressible Flow Equations. Phys. Fluids 2020, 32, 127116. [Google Scholar] [CrossRef]
  27. Guo, Y.; Hadžić, M.; Jang, J. Continued Gravitational Collapse for Newtonian Stars. Arch. Rational Mech. Anal. 2021, 239, 431–552. [Google Scholar] [CrossRef]
  28. Turzynsky, M. Properties of solutions to the gas dynamics equations on a rotating plane, corresponding to motions with homogeneous deformation. Mosc. Univ. Mech. Bull. 2020, 75, 37–43. [Google Scholar] [CrossRef]
  29. Polyanin, A.D.; Zaitsev, V.F. Handbook of Nonlinear Partial Differential Equations, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2012; pp. 1219–1220. [Google Scholar]
  30. Sidorov, A.F. Two classes of solutions of the fluid and gas mechanics equations and their connection to traveling wave theory. J. Appl. Math. Mech. 1989, 30, 197–203. [Google Scholar] [CrossRef]
  31. Shayakhmetova, R.F. Vortex scattering of monatomic gas along plane curves. J. Appl. Mech. Tech. Phys. 2018, 59, 241–250. [Google Scholar] [CrossRef]
  32. Khabirov, S.V.; Shayakhmetova, R.F. Simple solutions of rank 2 invariant submodel for monatomic gas. Chelyabinsk Phys. Math. J. 2018, 3, 353–373. (In Russian) [Google Scholar] [CrossRef]
  33. Siraeva, D.T. Classification of rank 2 stationary submodels of ideal hydrodynamics. Chelyabinsk Phys. Math. J. 2019, 4, 18–32. (In Russian) [Google Scholar] [CrossRef]
  34. Siraeva, D.T. The canonical form of the rank 2 invariant submodels of evolutionary type in ideal hydrodynamics. J. Appl. Ind. Math. 2019, 13, 340–349. [Google Scholar] [CrossRef]
  35. Siraeva, D.T.; Khabirov, S.V. Invariant submodel of rank 2 on subalgebra of translations linear combinations for a hydrodynamic type model. Chelyabinsk Phys. Math. J. 2018, 3, 38–57. (In Russian) [Google Scholar]
  36. Siraeva, D.T. Two invariant submodels of rank 1 of the hydrodynamic type equations and exact solutions. J. Phys. Conf. Ser. 2020, 1666, 012049. [Google Scholar] [CrossRef]
  37. Siraeva, D.T. Reduction of partially invariant submodels of rank 3 defect 1 to invariant submodels. Multiph. Syst. 2018, 13, 59–63. (In Russian) [Google Scholar] [CrossRef]
  38. Ovsyannikov, L.V. Lecture on Foundation of Gas Dynamics; Nauka: Moscow, Russia, 1981; 368p. (In Russian) [Google Scholar]
  39. Siraeva, D.T.; Yulmukhametova, Y.V. Transformations of gas dynamics equations and basis operators of a admitted 11-dimensional Lie algebra. Multiph. Syst. 2020, 3–4, 217–222. (In Russian) [Google Scholar] [CrossRef]
  40. Khabirov, S.V. Nonisomorphic Lie algebras admitted by gasdynamic models. Ufa Math. J. 2011, 3, 85–88. Available online: https://matem.anrb.ru/sites/default/files/files/vupe10/Khabirov.pdf (accessed on 26 December 2021).
  41. Khabirov, S.V. Lectures Analytical Methods in Gas Dynamics; BSU: Ufa, Russia, 2013; 224p. (In Russian) [Google Scholar]
  42. Nikonorova, R.F. The lowest-rank monatomic gas submodels constructed on the basis of three-dimensional symmetry subalgebras. Sib. Electron. Math. Rep. 2018, 15, 1216–1226. (In Russian) [Google Scholar] [CrossRef]
  43. Shayakhmetova, R.F. Vortex expansion of a monatomic gas. Proc. Mavlyutov Inst. Mech. 2014, 10, 110–113. (In Russian) [Google Scholar] [CrossRef]
Figure 1. The trajectories (22) with γ = 1 , a = 1 , b = 0 , ρ 0 = 1 2 , ξ = 1 ; 2 ; 3 , η = 1 , ζ = 1 . For t = 1..0, the trajectories are indicated by a dotted line; for t = 0 ..1, the trajectories are indicated by a solid line.
Figure 1. The trajectories (22) with γ = 1 , a = 1 , b = 0 , ρ 0 = 1 2 , ξ = 1 ; 2 ; 3 , η = 1 , ζ = 1 . For t = 1..0, the trajectories are indicated by a dotted line; for t = 0 ..1, the trajectories are indicated by a solid line.
Mathematics 10 00123 g001
Figure 2. (left) The motion of the particles volume (26) with x 1 = 0.1 , y 1 = 0.1 , z 1 = 0.1 , r = 1 , t 0 = 2 , ρ 0 = 1 , γ = 1 ; t = 2 ; 1 (grey sphere and ellipsoid); t = 0 (black ellipsoid); and t = 1 ; 2 ; 3 (red ellipsoids). (right) The motion of the particles volume (26) with x 1 = 0.1 , y 1 = 0.1 , z 1 = 4 , r = 1 , t 0 = 2 , ρ 0 = 1 , γ = 1 ; t = 2 ; 1 (grey sphere and ellipsoid); t = 0 (black ellipsoid); and t = 1 ; 2 ; 3 ; 4 (red ellipsoids).
Figure 2. (left) The motion of the particles volume (26) with x 1 = 0.1 , y 1 = 0.1 , z 1 = 0.1 , r = 1 , t 0 = 2 , ρ 0 = 1 , γ = 1 ; t = 2 ; 1 (grey sphere and ellipsoid); t = 0 (black ellipsoid); and t = 1 ; 2 ; 3 (red ellipsoids). (right) The motion of the particles volume (26) with x 1 = 0.1 , y 1 = 0.1 , z 1 = 4 , r = 1 , t 0 = 2 , ρ 0 = 1 , γ = 1 ; t = 2 ; 1 (grey sphere and ellipsoid); t = 0 (black ellipsoid); and t = 1 ; 2 ; 3 ; 4 (red ellipsoids).
Mathematics 10 00123 g002
Figure 3. The trajectory (42) in the plane ( r , η ) at r 0 = 2 , η 0 = 1 , b 3 = 1 .
Figure 3. The trajectory (42) in the plane ( r , η ) at r 0 = 2 , η 0 = 1 , b 3 = 1 .
Mathematics 10 00123 g003
Figure 4. (left) The trajectories (41), (42) in the space ( x , ξ , η ) at b 3 = 2 , t = 0 . . 4 with coordinates ( x 0 ; y 0 ; z 0 ) = ( 1 ; 1 ; 1 ) ; ( 3 ; 3 ; 1 ) ; ( 5 ; 5 ; 1 ) . (right) The trajectories (40) in the space ( x , y , z ) with the same parameters.
Figure 4. (left) The trajectories (41), (42) in the space ( x , ξ , η ) at b 3 = 2 , t = 0 . . 4 with coordinates ( x 0 ; y 0 ; z 0 ) = ( 1 ; 1 ; 1 ) ; ( 3 ; 3 ; 1 ) ; ( 5 ; 5 ; 1 ) . (right) The trajectories (40) in the space ( x , y , z ) with the same parameters.
Mathematics 10 00123 g004
Figure 5. The characteristic conoid (51) with a vertex P ( 2 ; 2 ; 1 ; 1 ) is represented in green color, the motion of the sphere (44) with center ( x 1 , y 1 , z 1 ) = ( 2 , 2 , 1 ) and with the radius R 0 = 15 π 6 at t 0 = 1 is red color; b 3 = 1 , b = 1 ; t = 1 ; 4 ; 15 .
Figure 5. The characteristic conoid (51) with a vertex P ( 2 ; 2 ; 1 ; 1 ) is represented in green color, the motion of the sphere (44) with center ( x 1 , y 1 , z 1 ) = ( 2 , 2 , 1 ) and with the radius R 0 = 15 π 6 at t 0 = 1 is red color; b 3 = 1 , b = 1 ; t = 1 ; 4 ; 15 .
Mathematics 10 00123 g005
Figure 6. (left) The characteristic conoid with a vertex P ( 0 ; 1 ; 1 ; 1 ) in the subsonic region is represented at t = 1.2 ; 1.6 ; 2.0 ; 2.5 ; 3.0 ; b 3 = 3 , b = 1 . (right) The characteristic conoid with a vertex P ( 2 ; 2 ; 1 ; 1 ) in the supersonic region is represented at t = 1.2 ; 1.6 ; 2.0 ; 2.5 ; 3.0 ; b 3 = 3 , b = 1 .
Figure 6. (left) The characteristic conoid with a vertex P ( 0 ; 1 ; 1 ; 1 ) in the subsonic region is represented at t = 1.2 ; 1.6 ; 2.0 ; 2.5 ; 3.0 ; b 3 = 3 , b = 1 . (right) The characteristic conoid with a vertex P ( 2 ; 2 ; 1 ; 1 ) in the supersonic region is represented at t = 1.2 ; 1.6 ; 2.0 ; 2.5 ; 3.0 ; b 3 = 3 , b = 1 .
Mathematics 10 00123 g006
Figure 7. The trajectories (60) with parameters a = 1 , b = 0 , u ¯ 0 i = 1 , i = 1..4, x 0 = 1 , 2 , 3 ; y 0 = z 0 = 1 , and t = 3..3.
Figure 7. The trajectories (60) with parameters a = 1 , b = 0 , u ¯ 0 i = 1 , i = 1..4, x 0 = 1 , 2 , 3 ; y 0 = z 0 = 1 , and t = 3..3.
Mathematics 10 00123 g007
Figure 8. The trajectories (66) are represented with b 2 = 5 , b 3 = 2 , ( x 0 , y 0 , z 0 ) = ( 1 , 1 , 1 ) ; ( 3 , 3 , 1 ) ; ( 5 , 5 , 1 ) , and t = 0..4.
Figure 8. The trajectories (66) are represented with b 2 = 5 , b 3 = 2 , ( x 0 , y 0 , z 0 ) = ( 1 , 1 , 1 ) ; ( 3 , 3 , 1 ) ; ( 5 , 5 , 1 ) , and t = 0..4.
Mathematics 10 00123 g008
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nikonorova, R.; Siraeva, D.; Yulmukhametova, Y. New Exact Solutions with a Linear Velocity Field for the Gas Dynamics Equations for Two Types of State Equations. Mathematics 2022, 10, 123. https://doi.org/10.3390/math10010123

AMA Style

Nikonorova R, Siraeva D, Yulmukhametova Y. New Exact Solutions with a Linear Velocity Field for the Gas Dynamics Equations for Two Types of State Equations. Mathematics. 2022; 10(1):123. https://doi.org/10.3390/math10010123

Chicago/Turabian Style

Nikonorova, Renata, Dilara Siraeva, and Yulia Yulmukhametova. 2022. "New Exact Solutions with a Linear Velocity Field for the Gas Dynamics Equations for Two Types of State Equations" Mathematics 10, no. 1: 123. https://doi.org/10.3390/math10010123

APA Style

Nikonorova, R., Siraeva, D., & Yulmukhametova, Y. (2022). New Exact Solutions with a Linear Velocity Field for the Gas Dynamics Equations for Two Types of State Equations. Mathematics, 10(1), 123. https://doi.org/10.3390/math10010123

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