Next Article in Journal
Chaos in a Magnetized Modified Gravity Schwarzschild Spacetime
Next Article in Special Issue
Massive Neutron Stars and White Dwarfs as Noncommutative Fuzzy Spheres
Previous Article in Journal
Covariant Evolution of Gravitoelectromagnetism
Previous Article in Special Issue
Analogue Quantum Gravity in Hyperbolic Metamaterials
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Maximal Kinematical Invariance Group of Fluid Dynamics and Applications

Chennai Mathematical Institute, H1 SIPCOT IT Park, Kelambakkam, Chennai 603103, India
*
Author to whom correspondence should be addressed.
Universe 2022, 8(6), 319; https://doi.org/10.3390/universe8060319
Submission received: 31 March 2022 / Revised: 20 May 2022 / Accepted: 21 May 2022 / Published: 7 June 2022
(This article belongs to the Special Issue Quantum Gravity Phenomenology)

Abstract

:
The maximal kinematical invariance group of the Euler equations of fluid dynamics for the standard polytropic exponent is larger than the Galilei group. Specifically, the inversion transformation ( Σ : t 1 / t , x x / t ) leaves the Euler equation’s invariant. This duality has been used to explain the striking similarities observed in simulations of the supernova explosions and laboratory implosions induced in plasma by intense lasers. The inversion symmetry extends to discontinuous fluid flows as well. In this contribution, we provide a concise review of these ideas and discuss some applications. We also explicitly work out the implosion dual of the Sedov’s explosion solution.

1. Introduction

Surprises lurk in unexpected corners of physics. This review summarises a body of results that ensue from one such surprise, viz. the striking similarity between the earlier simulations of supernova explosions and the experimental evolution of an imploding plasma contained in a fusion capsule bombarded by high-intensity lasers [1]. It was hoped [1] that, over time, laser experiments would become more in line with actual supernovae behaviour. Hence, considerable efforts were devoted to simulating astrophysical systems in the laboratory. Modern supernova simulations have become much more complex; several new physical effects and numerical techniques are incorporated [2]. At the present stages of development, it is not clear how much one can learn about the astrophysical systems in the laboratory setting. However, the observations mentioned above have led to some intriguing theoretical developments. In this paper we concentrate on discussing the theoretical explanation [3,4,5] for the observed similarity [1], and some ramifications of the resulting analysis [6]. We limit our considerations mostly to references [3,4,5,6].
Earlier computational studies of the evolution of a supernova remnant (as cited in [1]) usually used initial conditions of dense pressure-free ejecta expanding ballistically outwards from the site of the explosion, taken for convenience to be the origin, and interacting with a stationary ambient medium of much lower density and negligible pressure. At early times, the bulk of the ejecta expands ballistically, except for a thin interaction region on the outside consisting of a forward shock running into the ambient medium, a zone of hot-shocked ambient medium, a contact discontinuity, a zone of shocked ejecta, and a reverse shock propagating into the ejecta. At later times, when the mass of the swept-up ambient medium becomes comparable with the ejecta mass, the reverse shock detaches itself from the contact discontinuity and implodes on the origin.
In the laboratory plasma, we have, initially, a stationary sphere of high density material surrounded by a low density converging flow. The inflowing gas has to decelerate at the shock front, building up pressure and driving a reverse shock which leads to an implosion. From an experimental point of view, a perfectly spherically symmetric explosion is not realistic. The ejecta emerging from a supernova explosion is also highly nonuniform on a wide range of scales making it computationally challenging to calculate its evolution.
From a theoretical point of view, it is interesting to study the underlying symmetry that enables one to map an explosion to an implosion. One can straightaway rule out time-reversal as an answer because the supernova explosions occur over astronomical time-scales, while the plasma implosions happen in a few nano-seconds.
Both an exploding star and an imploding plasma can be modelled by the equations of a perfect fluid, as we are taught in standard textbooks [7,8]. The explanation offered in Refs. [3,4] for the observed explosion–implosion duality stems from a hitherto unnoticed nonlinear symmetry of these equations, which we expand upon in the next section.
This analysis of ref. [4] highlights that the maximal kinematical invariance group of the Euler equations of fluid dynamics for the standard polytropic exponent is larger than the Galilei group. The techniques required to establish this find applications in other situations, viz. fluid flows in non-inertial frames. The Earth’s oceans and atmosphere are important examples of fluid flows in non-inertial reference frames, where the Earth’s rotation provides the underlying non-inertial frame. In order to describe oceanic and atmospheric fluid flows, it is natural to analyse fluid phenomena in the Earth’s frame. This requires adding Coriolis forces to the right hand side of the fluid equations. The Coriolis force terms lead to surprising phenomenon: weather storms and ocean’s currents. Going to the non-inertial reference frame allows us to separate out the rotational component of the fluid flows. We can then concentrate on the parts of flow patterns that matter the most.
Related situations arise when fluid flows are characterised by a large degree of expansion or contraction. Poludnenko and Khokhlov [6] considered the formulation of Euler equations of fluid dynamics in an expanding or contracting or possibly rotating reference frame. The motivation being that by going to an appropriately chosen frame we can discard the expanding or contracting or rotating nature of the fluid flows. The frame motion is adjusted to minimise the local fluid velocities. This method allows to accommodate efficiently large degrees of change in the flow extent, such as those encountered in astrophysical flows: supernovae, contracting stars, etc. Their work investigated numerical computations in such non-inertial reference frames.
As in the case of rigidly rotating reference frames, going to an expanding or contracting reference frame requires adjustments of the fluid flow equations. Ref. [6] argued that these adjustments do not come at any additional numerical cost: the new equations can be as easily implemented numerically using any of the standard numerical schemes. However, by separating out the global component of the fluid flow, it leads to significant improvement in the physical understanding of the fluid flows, which would be difficult to extract in inertial reference frame simulations. (More precisely, in numerical work it is important to work with smaller local fluid velocities. If a fluid flow is dominated by the global component associated with expansion or contraction or rotation, then it is inefficient to model such flows in inertial frames.) They carried out extensive numerical testing of the method for a variety of reference frames representative of realistic applications.
The rest of the article is organised as follows. In Section 2, we review the maximal kinematical invariance group of fluid dynamics, based on the original work of O’ Raifeartaigh and one of the authors [4]. In Section 3, we discuss the symmetries of discontinuous flows, based on our original work with Oliver Jahn [5]. In Section 4, we review the work of Poludnenko and Khokhlov [6], who considered the formulation of Euler equations of fluid dynamics in non-inertial reference frames. In Section 5, we present the conclusions. Appendix A explicitly works out the implosion dual of Sedov’s explosion solution. Throughout the review we will be concise, referring the reader to original references for further details.

2. The Fluid Equations

This section is based on the original work of O’Raifeartaigh and one of the authors [4]. The general fluid dynamic equations in n-dimensional space are (see, e.g., textbooks [7,8])
D ρ = ρ · u ,
ρ D u = p + V ,
D ϵ = ( ϵ + p ) · u ,
where the convective derivative D and the viscosity terms V are defined by
D = t + u · ,
V i = j η ( j u i + i u j 2 n δ i j k u k ) + i ( ζ k u k ) ,
respectively. In the above equations ρ , u , p , ϵ stand for the density, the velocity vector field, the pressure, and the energy density respectively, and η and ζ are the bulk and shear viscosity fields. These partial differential equations are augmented by an algebraic condition called the equation of state that relates the pressure and energy density. According to the polytropic equation of state
p = ( γ 0 1 ) ϵ p + ϵ = γ 0 ϵ ,
where the constant γ 0 is called the polytropic exponent. This equation can be used to eliminate p from the fluid equations. Further, by making the substitution
ϵ = χ ρ γ 0
the equations can be rewritten in the form
D ρ = ρ · u ,
ρ D u = ( γ 0 1 ) ( χ ρ γ 0 ) + V ,
D χ = 0 .

2.1. Action Formulation

The fluid equations may be given an action formulation by switching off the viscosity fields, i.e., η = ζ = 0 . Such fluids are called inviscid, or perfect fluids, and the equations are called Euler’s equations. We next set χ = 1 without loss of consistency, representing the isentropicity condition. Further, the Clebsch parametrisation [9,10,11]
u = ϕ ν θ ,
allows us to isolate the irrotational parts by setting ν = θ = 0 . The resulting action for inviscid, isentropic, irrotational flows in three dimensions is given by,
S = d 3 x d t ρ ϕ ˙ 1 2 ( ϕ ) 2 ρ γ 0 .
It may be mentioned in passing that the terms contained in the parentheses represent the Hamilton–Jacobi function for a free particle.
The symmetries of the aforementioned special flows, represented by the action, follow from the requirement of its form-invariance. The transformation properties of the fields in the general fluid equations may be extracted from these transformation properties once again by requiring the equations to transform covariantly.

2.2. Symmetries

We begin a priori with the most general transformations involving the coordinates and fields. The transformed coordinates ξ and τ are defined by
ξ = ξ ( x , t ) , τ = τ ( x , t ) , ϕ ˜ = ϕ ˜ ( ξ , τ , ϕ ) , ρ ˜ = ρ ˜ ( ξ , τ , ρ ) .
Substituting the transformations and demanding the form-invariance of the action produces a set of equations which can be solved exactly to yield [4],
ξ = f ( t ) ( R x + a + v t ) ,
ρ ˜ = f 3 ( t ) ρ ,
ϕ ˜ = ϕ + λ ( ξ , τ ) ,
where
τ = α t + β γ t + δ , f ( t ) = 1 γ t + δ , α δ β γ = 1 ,
λ τ 1 2 λ ξ 2 = 0 .
In the above, R represents the usual rotation matrix, a the translations, v the boosts, and f ( t ) , a time-dependent scale parameter. The α , β , γ , δ represent parameters of the S L ( 2 , R ) group, a non-relativistic remnant of the special conformal group, to be discussed in Section 2.4. For details on λ ( ξ , τ ) , we refer the reader to the original reference [4]. We note that the following discrete symmetries are permitted:
( α , β , γ , δ ) ( α , β , γ , δ ) ( α , β , γ , δ ) ( α , β , γ , δ ) .

2.3. Transformation Functions for General Flows

The transformation functions for general flows may be obtained by requiring the general fluid equations to transform covariantly. It is straightforward to see that for general non-isentropic flows, the equations transform covariantly if χ is a scalar under the coordinate transformations.
The velocity vector transforms inhomogeneously as,
u ˜ = ( γ t + δ ) u γ x
We note that these transformations do not preserve the condition · u = 0 , implying that unlike Galilean symmetry, the above symmetry is valid only when the fluid is compressible.
The viscosity fields transform similar to scalar densities [4]. This implies that the symmetry we are discussing is broken in the case of Navier–Stokes equations which approximate the viscosity to be constant.

2.4. The Maximal Invariance Group

Let g be a general element of the symmetry group G obtained by setting β = γ = 0 , α = 1 . It follows
ξ = R x + v t + a , τ = t
In this case,
ρ ˜ = ρ and u ˜ = u + v
We notice that these correspond to the static Galilei transformations.
Let σ denote an element of the S L ( 2 , R ) group obtained by setting a = v = 0 . In this case,
ξ = f ( t ) x , τ = α t + β γ t + δ ,
ρ ˜ = ( γ t + δ ) 3 ρ and u ˜ = ( γ t + δ ) u γ x .
This represents a combination of dilatations and inversions, which are a nonrelativistic limit of the special conformal group.
It is a straightforward exercise to check that G is an invariant subgroup:
σ 1 · g ( R , a , v ) · σ = g ( R , a σ , v σ ) ,
where
a σ v σ = δ β γ α a v .
It therefore follows that the full group under which the fluid equations are invariant under the specified transformation properties for the coordinates and the fields is a semi-direct product
G = S L ( 2 , R ) G .
A special element of the group viz. Σ : ( α , β , γ , δ ) = ( 0 , 1 , 1 , 0 ) corresponds to a composition of an inversion and reversal of time, and plays an important role in the explosion–implosion duality discovered by Drury and Mendonça in [3]. A curious result follows immediately: Σ 2 = P , where P is the parity operator.
Cosets defined using these elements, g Σ ( R , a , v ) = Σ · g ( R , a , v ) have the interesting property that they are fourth roots of Galilean transformations,
g Σ 4 ( R , a , v ) = g R 4 , ( R 2 1 ) ( R a v ) , ( R 2 1 ) ( R v + a )
Since ( R a v ) and ( R v a ) are linearly independent, it follows that every Galilean transformation is a fourth power of a coset transformation [4].

2.5. Conservation Laws

Euler’s equations for a perfect fluid can be expressed in the form of conservation laws for mass, momentum, and energy, as follows:
ρ t = x j ( ρ u j ) ,
t ( ρ u i ) = x j ( ρ u i u j + δ i j p ) ,
t 1 2 ρ u 2 + ϵ = x j 1 2 ρ u 2 + ϵ + p u j .
These equations can be expressed succinctly as follows:
μ J ρ μ = μ J P μ = μ J H μ = 0 .
The zeroth components of the above currents namely, ρ , ρ u , 1 2 ρ u 2 + ϵ , give the charge densities which, when integrated over all space, give the conserved charges. The corresponding current densities are
J ρ j = ρ u j ,
J P i j = ( ρ u i u j + δ i j p ) ,
J H j = 1 2 ρ u 2 + ϵ + p u j .
The conservation laws corresponding to rotations ( δ x i = ω i j x j ), boosts ( δ x i = v i t ), dilatations ( δ x i = λ x i , δ t = 2 λ t ), and expansions ( δ x i = μ t x i , δ t = μ t 2 ) can be stated similarly,
μ J L μ = μ J K μ = μ J D μ = μ J A μ = 0 ,
where the appropriate charge densities are [5],
L = P × x , K = P t ρ x , D = 2 t H + x · P , A = t 2 H t x · P + ρ 2 2 x 2 .
The corresponding current densities can also be written down in a straightforward manner
J L i = ϵ i k l x k J P l ,
J K i = t J P i x i J ρ ,
J D = x i J P i 2 t J H ,
J A = 1 2 x 2 J ρ t x i J P i + t 2 J H .
These laws follow as a direct consequence of Noether’s theorem. They will be useful in studying flows with discontinuities, a topic to which we now pass.

3. Discontinuous Flows

This section is based on our original work with Oliver Jahn [5]. As long as the flows are smooth, i.e., the functions ρ , u , p , ϵ are smooth functions of x , t , Equations (1)–(3) and (31) are equivalent. Real flows, however, may develop discontinuities as they evolve. Such flows are described by weak solutions of differential Equations [12]. A weak solution is generally piecewise smooth. The smooth parts satisfy the differential equations in the usual strong form. The entire solution required to specify the course of motion of the initial conditions is obtained by supplementing the strong solutions by jump conditions. The jump conditions are derived from the conservation laws. We briefly review these concepts in the next two subsections. For pedagogical discussions on these topics we refer the reader to [7,8,13].

3.1. Weak Solutions and Jump Conditions

By definition, any, possibly non-smooth, function J μ ( x , t ) that satisfies
μ ω ( x , t ) J μ ( x , t ) d 3 x d t = 0 ,
for all test functions ω ( x , t ) is called a weak solution of the differential equation μ J μ = 0 .
Suppose J μ ( x , t ) has a jump discontinuity across a hypersurface S in x , t space while otherwise being continuously differentiable in some neighbourhood N of S ; see Figure 1. Let ω ( x , t ) be a test function with support in region N . Let R be the part of the region N that lies on one side of S , say to the right. We take ω ( x , t ) = 0 on the boundary of R , except on S . Then by Gauss’s theorem,
R μ ω J μ d 3 x d t + R ω μ J μ d 3 x d t = R μ ( ω J μ ) d 3 x d t = S ω n μ J μ d S ,
since ω ( x , t ) = 0 on the boundary of R , except on S . Here n μ is the outward normal to the hypersurface S . The second integral on the left hand side is zero because the conservation law holds in the strong sense in the interior of R . If we integrate similarly over the left part of the support ω ( x , t ) , and add the two results we obtain for a weak solution
0 = S ω n μ Δ J μ d S ,
where Δ f denotes the difference of the two limiting values of the the function f on the two sides of the hypersurface S , i.e., the jump of the function. This result follows as the normal which points outwards by convention, flips its sign when we move from the right to the left side of the hypersurface. Since ω ( x , t ) is an arbitrary function, it follows that
n μ Δ J μ = 0 on S .
This equation is the jump condition we mentioned earlier.

3.2. Rankine–Hugoniot Conditions

The general expression for the jump condition derived above can be applied to the conservation laws derived earlier. The conservation laws associated with mass, momentum, and energy yield
n μ Δ J ( ρ ) μ = 0 ,
n μ Δ J ( P ) μ = 0 ,
n μ Δ J ( H ) μ = 0 ,
and are called the Rankine–Hugoniot conditions in the fluid dynamics literature [7,8,13,14].
Similar equations can be derived for the other conservation laws, viz.
n μ Δ J ( L ) μ = 0 ,
n μ Δ J ( K ) μ = 0 ,
n μ Δ J ( D ) μ = 0 ,
n μ Δ J ( A ) μ = 0 .
These new jump conditions are associated with angular momentum, boosts, dilatations, and expansions.

3.3. Dual Rankine–Hugoniot Conditions

The new set of jump conditions are identically true because the current densities associated with angular momentum, boosts, dilatations, and expansions are linear combinations of the current densities associated with mass, momentum and energy conservation as shown in Equation (40). This suggests that the Rankine–Hugoniot conditions are invariant under the full kinematical invariance group of smooth flows including the S L ( 2 , R ) part. We examine this point in what follows.
If an (abstract) symmetry generator T r transforms under the S L ( 2 , R ) transformation σ as
T r = σ 1 T r σ = s M r s ( σ ) T s ,
then the corresponding currents transform as [5]
J r μ ( x ) = det x x x μ x ν s M r s ( σ ) J s ν ( x ) .
The fact that the currents transform similar to vector densities can be appreciated by looking at the temporal components, which pick up the multiplicative factor ( γ t + δ ) 3 .
Arranging the currents in a column J μ = J ( ρ ) μ , J ( K ) μ , J ( P ) μ , J ( A ) μ , J ( D ) μ , J ( H ) μ T , one has the following transformation matrix
M = 1 0 0 0 0 0 0 α β 0 0 0 0 γ δ 0 0 0 0 0 0 α α β β 2 0 0 0 2 α γ ( β γ + α δ ) 2 β γ 0 0 0 γ 2 γ δ δ 2 .
Using α δ β γ = 1 , it is easy to check that the matrix M has unit determinant.
The temporal components transform according to the transformations,
ρ = ( γ t + δ ) 3 ρ ,
P = ( γ t + δ ) 3 ( δ P + γ K ) ,
H = ( γ t + δ ) 3 ( γ 2 A δ γ D + δ 2 H ) .
Thus, ρ transforms under the singlet representation of S L ( 2 , R ) as a scalar density. The translations and boosts constitute the doublet representation and transform similar to vector densities. Likewise the Hamiltonian, and the generators of dilatations and expansions transform similar to densities under the triplet (adjoint) representation of S L ( 2 , R ) . The transformation properties of the spatial components of the current can similarly be read off from the above matrix.
The dual Rankine-Hugoniot conditions are now easily obtained. The normal vector n μ appearing in the jump condition (44) transforms like a covector
n μ = x ν x μ n ν .
Thus, the transformed jump conditions for the currents are
n μ Δ J r μ det x x s M r s ( σ ) n μ Δ J s μ ( x ) = 0 on S .
Since the determinant is smooth across the hypersurface S , the factor in front can be omitted, and the transformed jump condition is a linear combination of the original jump conditions.
In particular, the transformed conditions for J ( ρ ) μ , J ( P ) μ , J ( H ) μ (the Rankine–Hugoniot conditions), become linear combinations of the jump conditions of J ( ρ ) μ , J ( P ) μ , J ( H ) μ , J ( K ) μ , J ( D ) μ and J ( A ) μ . Specifically,
n μ Δ J ( ρ ) μ n μ Δ J ( ρ ) μ
n μ Δ J ( P ) μ n μ ( γ Δ J ( K ) μ + δ Δ J ( P ) μ )
n μ Δ J ( H ) μ n μ ( δ 2 Δ J ( H ) μ γ δ Δ J ( D ) μ + γ 2 Δ J ( A ) μ )
The original Rankine-Hugoniot conditions, in conjunction with the new conditions (48)–(51), imply that the right hand sides of the above equations are identically zero, i.e., the Rankine-Hugoniot conditions are form-invariant [5].
In particular, this holds for the Drury–Mendonça transformation [3] t 1 / t , x x / t used to relate the explosion and implosion problems. This corresponds to the choice ( α , β , γ , δ ) = ( 0 , 1 , 1 , 0 ) . We conclude that, if an explosion is described by the standard Rankine–Hugoniot conditions, the corresponding implosion is described by the dual Rankine–Hugoniot conditions.

3.4. Physical Conditions

As already mentioned, for a polytropic gas, ϵ = χ ρ γ 0 , which enables us to write the third of Euler’s Equation (3) as D χ = 0 . χ transforms similar to a scalar. For a polytropic gas, it is well known [14] that χ is related to the specific entropy (entropy per unit mass) as follows:
S S 0 = C V log [ χ ( ρ , V ) γ 0 ] ,
where C V = R / ( γ 0 1 ) , R being the universal gas constant divided by the molecular weight, V, the volume, and S 0 , an appropriate constant. Since χ transforms similar to a scalar, it follows that the specific entropy of a moving particle remains constant under an S L ( 2 , R ) transformation. Hence a physical shock is mapped to a physical shock under the transformation.
The requirement that the transformation preserves the physicality of a shock puts a condition on the viscosity viz. that its positivity is preserved. As already pointed out, the viscosity fields transform as scalar densities, similar to ρ . It follows that the total viscosity, such as mass, is an invariant under the transformation.

4. Fluid Equations in Non-Inertial Frames

In this section, we review the work of Poludnenko and Khokhlov [6], who considered the formulation of Euler equations of fluid dynamics in non-inertial reference frames. We start with the inertial reference frame Euler Equations (1)–(3)
t ρ + i ( ρ u i ) = 0 ,
t ( ρ u i ) + j ( ρ u i u j ) + i p = 0 ,
t ϵ + i ( ( ϵ + p ) u i ) = 0 ,
where ρ is the density, u i the fluid velocity, p the pressure, and ϵ the total energy density. The total energy ϵ is related to the internal energy per unit mass e as,
ϵ = e ρ + 1 2 ρ u 2 .
Consider a non-inertial reference frame { x ˜ i , τ } that expands or contracts with respect to the inertial frame { t , x i } :
x ˜ i = x i a ( t ) , τ = 0 t d t a ( t ) β + 1 ,
where β is a constant and the scale factor a ( t ) is a smooth non-vanishing function of time. (Poludnenko and Khokhlov [6] also considered additional rotational terms in transformation (68).) For simplicity we do not consider such terms; essential ideas are all captured by the simplified transformation). We use quantities with tilde signs to refer to quantities in the non-inertial reference frame { x ˜ i , τ } . Time t is the physical time, and τ is the computational time.
For the density, pressure, and energy density fields we introduce the scaling,
ρ ˜ ( x ˜ , τ ) = a α ρ ( x , t ) ,
p ˜ ( x ˜ , τ ) = a α + 2 β p ( x , t ) ,
e ˜ ( x ˜ , τ ) = a 2 β e ( x , t ) ,
where α and β are constant scaling exponents. A short calculation shows that
u i = d d t x i ( t ) = a β d ln a d τ x ˜ i + a β u ˜ i ,
and
a ¨ = 1 a 2 β + 1 d 2 ln a d τ 2 β d ln a d τ 2 .
As a result, the mass conservation Equation (64) in the non-intertial frame (68) become,
τ ρ ˜ + ˜ i ( ρ ˜ u ˜ i ) = ( α n ) d ln a d τ ρ ˜ ,
where ˜ i are partial derivatives with respect to x ˜ i and n is the dimension of space. The momentum conservation Equation (65) becomes
τ ( ρ ˜ u ˜ i ) + ˜ j ( ρ ˜ u ˜ i u ˜ j ) + ˜ i p ˜ = ( α n + β 1 ) d ln a d τ ρ ˜ u ˜ i a 2 β + 1 a ¨ ρ ˜ x ˜ i .
The transformation of the energy Equation (66) to the non-intertial frame (68) is quite tedious. When the dust settles one finds,
τ ϵ ˜ + ˜ i ( ( ϵ ˜ + p ˜ ) u ˜ i ) = d ln a d τ ( α n + 2 β ) ϵ ˜ n p ˜ ρ ˜ u ˜ i u ˜ i a 2 β + 1 a ¨ ρ ˜ u ˜ i x ˜ i .

4.1. Conditions for Invariance of the Fluid Equations

There is subclass of transformations (68) that preserve the form of the Euler’s equation. Let us look at this subclass in relation to the discussion of the previous sections. For the form invariance of mass conservation Equation (64) we require α = n from Equation (74). For the invariance of momentum and energy conservation (65) and (66), we require from Equations (75) and (76) , a ¨ = 0 , together with
α = n , β = 1 ,
and
ϵ = n 2 p + 1 2 ρ u 2 .
Recalling that the total energy ϵ is related to the internal energy e via (67), we have
e = n p 2 ρ
This is a restriction on the equation of state. For a polytropic gas p = χ ρ γ 0 we have the general relation,
e = p ( γ 0 1 ) ρ .
Thus, we conclude that the form invariance of Equations (74)–(76) singles out a special value of the polytropic index,
γ 0 = 1 + 2 n .
For n = 3 , γ 0 = 5 / 3 . Low mass white-dwarf stars are well approximated by this polytropic index. These results are perfectly consistent with [4] reviewed in the previous sections. a ¨ = 0 implies,
a ( t ) = c t + d .
Thus, transformation (68) becomes
x i x i c t + d , t 1 c ( γ t + d ) .
Comparing this with general S L ( 2 , R ) transformations [4]
x i x i c t + d , t a t + b c t + d , a d b c = 1 ,
we have a = 0 , b = 1 / c . The scaling of density, pressure, and energy density (69)–(71) are also compatible with the scalings in [4], and so is the value of the polytropic index.

4.2. Non-Invariant Terms as Sources

Poludnenko and Khokhlov argue that the above formulation based on general scaling of the fluid variables provides a degree of flexibility, provided we treat the non-invariant terms as sources. They consider values of exponents other than in Equation (77) that do not leave the form of the equations invariant. For example, a set of exponents can be obtained by demanding the invariance of the first law of thermodynamics. For isentropic flows, the first law of thermodynamics in inertial frames take the form
d s = 0 d e = p d 1 ρ ,
which for fluid flows implies,
t e + u i i e = p ρ 2 ( t ρ + u i i ρ ) .
In non-inertial frames (68), Equation (86) becomes
τ e ˜ + u ˜ i ˜ i e ˜ = p ˜ ρ ˜ 2 ( ˜ τ ρ ˜ + u ˜ i ˜ i ρ ˜ ) ( α p ˜ 2 β e ˜ ρ ˜ ) 1 ρ ˜ d ln a d τ .
Using the polytropic equation of state (80), it simplifies to
τ e ˜ + u ˜ i ˜ i e ˜ = p ˜ ρ ˜ 2 ( ˜ τ ρ ˜ + u ˜ i ˜ i ρ ˜ ) + ( 2 β α ( γ 0 1 ) ) d ln a d τ e ˜ .
The choice of exponents
α = n , β = n ( γ 0 1 ) 2 ,
ensures the invariance of the first law of thermodynamics together with the mass conservation for all values of the polytropic index. The momentum conservation equation no longer takes the conservative form, however. With this choice of exponents (89), the first of the source terms of the momentum conservation Equation (76),
( α n + 2 β ) ϵ ˜ n p ˜ ρ ˜ u ˜ i u ˜ i
simplifies to
( β 1 ) ρ ˜ u ˜ i u ˜ i .
This source term is proportional to the kinetic energy. This new set of exponents may be a preferred choice in numerical simulations if the thermal energy dominates the local kinetic energy in the non-inertial frame.

4.3. Primitive Fields as Simulation Variables

The key drawback in working with exponents (77) or (89) is the fact that they modify the primitive fields (69)–(71). From the general transformed equations (74)–(76), we immediately note that the homogeneous part of the equations is form-invariant. This allows for straightforward numerical implementations of the transformed equations for any value of the scaling exponents α and β treating the right hand side terms in Equations (74)–(76) as sources. The equations no longer take the form of conservation laws, but this is not a problem. In most practical situations this is a necessity. For example, if systems are governed by a non-polytropic equation of state, then we must work with source terms irrespective of the choice of the scaling exponents. We may as well work with the primitive fields as simulation variables, that is, we choose
α = 0 , β = 0 .
The use of primitive fields as simulation variables has the advantage of direct interpretation.

4.4. Numerical Results

In numerical work, source terms are frequently treated. Depending on the problem under consideration, source terms representing gravitational forces, electromagnetic forces, energy release due to radiation, etc are routinely added. Therefore, numerical computation in a moving frame can be performed at virtually no extra technical complication and at no extra computational cost. Poludnenko and Khokhlov mostly focus on tests of moving frame formulation of the fluid flow Equations (74)–(76) with zero exponents (92). They only briefly discuss other choices of scaling parameters. They perform their numerical simulations in a variety of frames for diverse physical problems. The details can be found in their paper. The key points are summarised as follows:
  • They consider several types of non-inertial reference frames: accelerating, expanding, contracting, oscillating (sinusoidal) reference frames, etc.
  • They treat in detail simulations of blast solutions (e.g., Sedov solution), converging shock solutions (e.g., Guderley blast wave solution), problems involving expansion of a gas sphere in vacuum, etc. They work in different reference frames best suited for the problem at hand.
  • They note that the computation in moving frames does not introduce systematic errors. Numerical solutions properly converge to the exact ones when they are known, e.g., the Sedov solution.
  • The method accuracy is valid even when solving the fluid equations for non-zero values of the exponents α and β .
There are some problems for which numerical simulations in non-inertial frames are not ideally suited. Such problems typically involve stationary regions of fluids in an inertial frame. The canonical example being the strong explosion in an otherwise stationary environment. (Expanding or collapsing environments where the ambient conditions are vacuous or dynamically unimportant can be optimally treated in non-inertial frames. In such problems, the ambient fluid can be set to be stationary in the computational–non-inertial–frame.) However, the demonstration in [6] that the numerical solution converges to the correct analytic Sedov solution is a crucial test of the method. The success of the simulation clearly shows that the non-conservative nature of the method does not introduce systematic errors and the Rankine–Hugoniot conditions are valid in the transformed reference frame too. The Rankine–Hugoniot conditions were shown to be valid in a subclass of transformed reference frames with scaling exponents Equation (77) in [5], as reviewed above. Further generalisation for different scaling exponents has not yet been carried out; however, given the numerical results it is likely that a useful formulation exists in more general situations.

5. Conclusions

In this article, we reviewed that the maximal kinematical invariance group of the Euler equations of fluid dynamics is larger than the Galilei group. Specifically, the inversion transformation ( Σ : t 1 / t , x x / t ) leaves the Euler equations invariant. This duality has been used to explain the striking similarities observed in simulations of the supernova explosions and laboratory implosions induced in plasma by intense lasers. It is quite remarkable that the inversion symmetry extends to discontinuous fluid flows as well. We also reviewed how this comes about.
We summarised the work of Poludnenko and Khokhlov [6]. They presented methods for computation of fluid flows characterised by large degree of expansion or contraction. The key idea is the transformation to a non-inertial reference frame. The scaling transformation of the primitive fluid variables provides a degree of flexibility. The use of non-inertial frames often leads to non-conservative formulation of the fluid equations, however, this does not affect the accuracy of the numerical work. For many problems of astrophysical and geophysical interests, going to an appropriate non-inertial frame allows for a cleaner extraction of relevant physics. We focused on [6] because of its close connection to the invariance properties of fluid equations.
There are several other papers addressing these issues, see, for example [15,16] and references therein. Similar ideas are frequently used in simulations of galaxies and the large-scale structure in an expanding universe and in atmospheric simulations.

Author Contributions

All authors have contributed equally. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Implosion Dual of the Sedov Explosion Solution

In this appendix we write down the explicit implosion solution dual of the Sedov explosion solution. To the best of our knowledge this has not been carried out before. (Drury and Mendonça in [3] have made some elementary remarks. They comment that Dwarkadas and Drury would publish details on the implosion solution dual to the Sedov solution in a separate paper. However, we are unable to locate the relevant references. Perhaps the results were not communicated to a journal. We will be glad if someone can point out relevant references to us.) We start with a brief review of the Sedov solution. We then apply the duality transformation to write the dual implosion solution. Some properties of the dual solution are presented. For the Sedov solution our presentation closely follows ([7], Section 99).

Appendix A.1. Sedov Solution

Sedov solution refers to an exact solution of compressible fluid dynamics equations in which spherical shock of great intensity propagates radially outwards as a result of a strong explosion. Strong explosion is characterised by the instantaneous release of energy E at the center. The equation of fluid dynamics in spherical symmetric situations take the form
ρ t + ( ρ u ) r + 2 ρ u r = 0 , u t + u u r + 1 ρ p r = 0 , s t + u s r = 0 ,
where s = ln ( p ρ γ ) . The last equation is the conservation of entropy. For the Sedov solution, the pressure discontinuity is very large: the pressure behind the shock is p 1 is much larger than the pressure in front of the shock p 0 . Sedov solution neglects p 0 everywhere (more on this later). The flow pattern is completely determined by the energy released E and the ambient density ρ 0 . The ratio of densities just behind and front of the shock is obtained by the Rankine–Hugoniot condition,
ρ 1 = γ + 1 γ 1 ρ 0 ,
assuming p 1 p 0 . The shock front is defined by
R ( t ) = ξ 0 E ρ 0 1 / 5 t 2 / 5 .
The propagation velocity of the shock is
D = d R d t = 2 5 R t .
Now, using the other two Rankine–Hugoniot conditions, which determine the gas velocity u 1 and pressure p 1 immediately behind the shock front, we obtain
p 1 = 2 γ + 1 ρ 0 D 2 , u 1 = 2 γ + 1 D .
Note that as the shock expands p 1 and u 1 also change as a function of time. We define
ξ = r R ( t )
and define
p ( r , t ) = 8 ρ 0 25 ( γ + 1 ) · r 2 t 2 · p ˇ ( ξ ) ,
u ( r , t ) = 4 5 ( γ + 1 ) · r t · u ˇ ( ξ ) ,
ρ ( r , t ) = ρ 0 γ + 1 γ 1 · ρ ˇ ( ξ ) .
Variables p ˇ , u ˇ , ρ ˇ are dimensionless pressure, velocity, and density. These are functions of the dimensionless variable ξ . The Rankine–Hugoniot jumps conditions in terms of these functions become
p ˇ = u ˇ = ρ ˇ = 1 a t ξ = ξ 0 .
These are the new boundary conditions. We warn the reader that there is a huge variation in the literature on the use of dimensionless variables for pressure, density, and velocity.
Using self-similarity of the solution one can argue that energy contained in a sphere of constant ξ remains constant in time. This gives an integral of motion. The argument proceeds as follows. For more details see [7]. Consider a spherical volume size r at constant ξ . It expands at the rate 2 r 5 t . The energy exiting the sphere in time d t due to the motion of the fluid is
4 π r 2 · ρ v s . h + 1 2 v 2 · d t .
This energy must be equal to the increase in the internal energy of the sphere in time d t due to its expansion
4 π r 2 · ρ ϵ + 1 2 v 2 · 2 r 5 t · d t .
Equating the two expressions give the first integral,
p ˇ ( ξ ) ρ ˇ ( ξ ) = γ + 1 2 u ˇ ( ξ ) 2 γ u ˇ ( ξ ) γ 1 u ˇ 2 ( ξ ) .
We obtain the remaining equations from the mass and entropy conservation equations. These equations simplify in the form
d u ˇ d ln ξ + u ˇ γ + 1 2 d ln ρ ˇ d ln ξ = 3 u ˇ ,
and
d d ln ξ ln p ˇ ρ ˇ γ = 5 ( γ + 1 ) 4 u ˇ 2 u ˇ ( γ + 1 ) .
These equations can be integrated to give implicitly the functions p ˇ ( ξ ) , u ˇ ( ξ ) , ρ ˇ ( ξ ) . They take the form
ξ 0 ξ 5 = u ˇ 2 5 ( γ + 1 ) 2 ( 3 γ 1 ) u ˇ 7 γ ν 1 2 γ u ˇ γ 1 γ 1 ν 2 ,
ρ ˇ = 2 γ u ˇ γ 1 γ 1 ν 3 5 ( γ + 1 ) 2 ( 3 γ 1 ) u ˇ 7 γ ν 4 γ + 1 2 u ˇ γ 1 ν 5 ,
where
ν 1 = 13 γ 2 7 γ + 12 ( 3 γ 1 ) ( 2 γ + 1 ) ,
ν 2 = 5 ( γ 1 ) 2 γ + 1 ,
ν 3 = 3 2 γ + 1 ,
ν 4 = 13 γ 2 7 γ + 12 ( 2 γ ) ( 3 γ 1 ) ( 2 γ + 1 ) ,
ν 5 = 2 γ 2 .
Here we have corrected a few typos from [7] (the ν 5 there has a typo). The parameter ξ 0 is determined by the requirement that the total energy of the gas up to radius R ( t ) is E. Other details of the solution can be found in [7,13]. Although reference [13] does not discuss the explicit solution, the discussion on the physical properties of the solution is very thorough and lucid.

Appendix A.2. Duality in Spherical Coordinates

In order to work out the implosion dual it is instructive to first work out the invariance of the simplified fluid Equation (A1) in spherical coordinates. To this end, we define the new time (now we use the notation t ˜ for the new time as opposed to τ ) and the new radial variable r ˜
r ˜ = r a ( t ) , t ˜ = 0 t d t a ( t ) β + 1 ,
and rescale pressure and density as,
ρ ˜ ( r ˜ , t ˜ ) = a α ρ ( r , t ) ,
p ˜ ( r ˜ , t ˜ ) = a α + 2 β p ( r , t ) .
These transformations give
u = a β u ˜ + r ˜ d ln a d t ˜ .
As for the derivatives we need to use
t X = X t ˜ · t ˜ t + X r ˜ · r ˜ t = a β 1 X t ˜ r ˜ a β 1 d ln a d t ˜ X r ˜ ,
and
r X = X r ˜ · r ˜ r = a 1 X r ˜ .
Now it is not difficult to verify that:
(1)
The continuity equation is form invariant for α = 3 .
(2)
The momentum equation is form invariant provided a ¨ = 0 and β = 1 .
(3)
The entropy equation is form invariant for γ = 5 / 3 .

Appendix A.3. Implosion Dual

We choose a ( t ) = t . Then,
r ˜ = r t , t ˜ = 1 t .
Since the radial variable does not undergo an inversion, the interior Sedov solution is mapped to an interior solution, and the exterior solution is mapped to an exterior solution. Due to this, the dual solution does not satisfy any physically interesting boundary conditions, i.e., it cannot be compared with standard implosion solutions of the sort discussed in say, chapter XII of [13]. For a simple physically interesting laboratory realisable implosion solution, one would require the interior to be stationary at fixed density and negligible pressure. This is certainly not the case for the dual solution. On the contrary the exterior solution is at zero pressure (hence there the speed of sound is zero) and the fluid is moving.
We ask in what sense is the solution an implosion solution. Does it satisfy expected properties, specifically the Rankine–Hugoniot jumps conditions? We take 0 < r ˜ < and < t ˜ < 0 . For the initially stationary exterior region the velocity transformation gives
u ˜ ( r ˜ , t ˜ ) = r ˜ t ˜ .
Since t ˜ is negative, the velocity of the exterior fluid is directed inwards; an implosion. The location of the shock is
R ˜ ( t ˜ ) = R ( t ) t = t ˜ R ( t ) = ξ 0 E ρ 0 1 / 5 ( t ˜ ) 3 / 5 .
As t ˜ increases from negative value towards zero, R ˜ ( t ˜ ) decreases, i.e., it represents an implosion. The velocity of the shock surface is
d R ( t ˜ ) d t ˜ = 3 5 R ˜ ( t ˜ ) t ˜ .
Since the inward velocity of the shock surface is smaller than the inward velocity of ambient fluid just outside, the fluid is injected into the interior region through the shock surface. Thus, in this set-up the exterior fluid of low (zero) pressure is getting compressed at the shock surface into the interior region.
Let us now calculate the velocity of the fluid just behind the shock surface. Recall ξ = r R ( t ) . It follows that,
ξ = r ˜ R ˜ ( t ˜ ) ,
so the interpretation of ξ as a dimensionless variable remains the same. We have
u ˜ ( r ˜ , t ˜ ) = 4 5 ( γ + 1 ) · r ˜ t ˜ · u ˇ ( ξ ) + r ˜ t ˜ ,
for the interior region. Similarly, other variable can be constructed. At the shock surface,
u ˜ = 4 5 ( γ + 1 ) R ˜ t ˜ + R ˜ t ˜ = 1 + 5 γ 5 ( 1 + γ ) R ˜ t ˜ .
the frame of the shock, we can confirm that these velocities satisfy the Rankine–Hugoniot conditions. The other Rankine–Hugoniot conditions can also be checked similarly. The post-shock pressure increases as ( t ˜ ) 19 / 5 . These results are all consistent with the comments in Drury and Mendonça in [3].

References

  1. Remington, B. Supernova Hydrodynamics Upclose. Science and Technology Review. Lawrence Livermore Library. 2000. Available online: https://str.llnl.gov/content/pages/past-issues-pdfs/2000.01.pdf (accessed on 31 March 2022).
  2. Available online: https://www.mpa-garching.mpg.de/84411/Core-collapse-supernovae (accessed on 31 March 2022).
  3. Drury, L.O.; Mendonça, J.T. Explosion implosion duality and the laboratory simulation of astrophysical systems. Phys. Plasmas 2000, 7, 5148. [Google Scholar] [CrossRef] [Green Version]
  4. O’Raifeartaigh, L.; Sreedhar, V.V. The Maximal kinematical invariance group of fluid dynamics and explosion–implosion duality. Ann. Phys. 2001, 293, 215–227. [Google Scholar] [CrossRef] [Green Version]
  5. Jahn, O.; Sreedhar, V.V.; Virmani, A. Symmetries of discontinuous flows and the dual Rankine–Hugoniot conditions in fluid dynamics. Ann. Phys. 2005, 316, 30–43. [Google Scholar] [CrossRef] [Green Version]
  6. Poludnenko, A.Y.; Khokhlov, A.M. Computation of fluid flows in non-inertial contracting, expanding, and rotating reference frames. J. Comput. Phys. 2007, 220, 678–711. [Google Scholar] [CrossRef] [Green Version]
  7. Landau, L.; Lifshitz, E. Fluid Mechanics. In Course of Theoretical Physics; Pergamon Press: Oxford, UK, 2013; Volume 6. [Google Scholar]
  8. Thorne, K.S.; Blandford, R.D. Modern Classical Physics; Princeton University Press: Princeton, NJ, USA, 2017. [Google Scholar]
  9. Lamb, H. Hydrodynamics; Cambridge University Press: Cambridge, UK, 1942. [Google Scholar]
  10. Deser, S.; Jackiw, R.; Polychronakos, A.P. Clebsch (string) decomposition in d=3 field theory. Phys. Lett. A 2001, 279, 151–153. [Google Scholar] [CrossRef] [Green Version]
  11. Jackiw, R.; Nair, V.P.; Pi, S.Y. Chern-Simons reduction and nonAbelian fluid mechanics. Phys. Rev. D 2000, 62, 085018. [Google Scholar] [CrossRef] [Green Version]
  12. Richtmeyer, R.D. Principles of Advanced Mathematical Physics; Springer: New York, NY, USA, 1978; Volume 1. [Google Scholar]
  13. Zel’dovich, Y.B.; Raizer, Y.P. Physics of Shock Waves and High-Temperature Hydrodynamic Phenomena; Hayes, W.D., Probstein, R.F., Eds.; Dover Publication: Mineola, NY, USA, 2002. [Google Scholar]
  14. Courant, R.; Friedrichs, K.O. Supersonic Flow and Shock Waves; Interscience Publishers: New York, NY, USA, 1948. [Google Scholar]
  15. Trac, H.; Pen, U.L. A moving frame algorithm for high Mach number hydrodynamics. New Astron. 2004, 9, 443–465. [Google Scholar] [CrossRef]
  16. Gledhill, I.M.A.; Roohani, H.; Forsberg, K.; Eliasson, P.; Skews, B.W.; Nordström, J. Theoretical treatment of fluid flow for accelerating bodies. Theor. Comput. Fluid Dyn. 2016, 30, 449–467. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Diagram for the jump condition.
Figure 1. Diagram for the jump condition.
Universe 08 00319 g001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sreedhar, V.V.; Virmani, A. Maximal Kinematical Invariance Group of Fluid Dynamics and Applications. Universe 2022, 8, 319. https://doi.org/10.3390/universe8060319

AMA Style

Sreedhar VV, Virmani A. Maximal Kinematical Invariance Group of Fluid Dynamics and Applications. Universe. 2022; 8(6):319. https://doi.org/10.3390/universe8060319

Chicago/Turabian Style

Sreedhar, V. V., and Amitabh Virmani. 2022. "Maximal Kinematical Invariance Group of Fluid Dynamics and Applications" Universe 8, no. 6: 319. https://doi.org/10.3390/universe8060319

APA Style

Sreedhar, V. V., & Virmani, A. (2022). Maximal Kinematical Invariance Group of Fluid Dynamics and Applications. Universe, 8(6), 319. https://doi.org/10.3390/universe8060319

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