Next Article in Journal
Flow Dynamics in a Model of a Left Ventricle with Different Mitral Valve Orientations
Next Article in Special Issue
Extraction of Tangential Momentum and Normal Energy Accommodation Coefficients by Comparing Variational Solutions of the Boltzmann Equation with Experiments on Thermal Creep Gas Flow in Microchannels
Previous Article in Journal
Optimization of a Grid-Connected Microgrid Using Tidal and Wind Energy in Cook Strait
Previous Article in Special Issue
A Review on BGK Models for Gas Mixtures of Mono and Polyatomic Molecules
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Shock Structures Using the OBurnett Equations in Combination with the Holian Conjecture

Department of Mechanical Engineering, Indian Institute of Technology—Bombay Powai, Mumbai 400076, India
*
Author to whom correspondence should be addressed.
Fluids 2021, 6(12), 427; https://doi.org/10.3390/fluids6120427
Submission received: 26 September 2021 / Revised: 10 November 2021 / Accepted: 22 November 2021 / Published: 26 November 2021
(This article belongs to the Special Issue Rarefied Gas Dynamics)

Abstract

:
In the present work, we study the normal shock wave flow problem using a combination of the OBurnett equations and the Holian conjecture. The numerical results of the OBurnett equations for normal shocks established several fundamental aspects of the equations such as the thermodynamic consistency of the equations, and the existence of the heteroclinic trajectory and smooth shock structures at all Mach numbers. The shock profiles for the hydrodynamic field variables were found to be in quantitative agreement with the direct simulation Monte Carlo (DSMC) results in the upstream region, whereas further improvement was desirable in the downstream region of the shock. For the discrepancy in the downstream region, we conjecture that the viscosity–temperature relation ( μ T φ ) needs to be modified in order to achieve increased dissipation and thereby achieve better agreement with the benchmark results in the downstream region. In this respect, we examine the Holian conjecture (HC), wherein transport coefficients (absolute viscosity and thermal conductivity) are evaluated using the temperature in the direction of shock propagation rather than the average temperature. The results of the modified theory (OBurnett + HC) are compared against the benchmark results and we find that the modified theory improves upon the OBurnett results, especially in the case of the heat flux shock profile. We find that the accuracy gain is marginal at lower Mach numbers, while the shock profiles are described better using the modified theory for the case of strong shocks.

1. Introduction

In the continuum regime for vanishing Knudsen numbers ( K n ( = λ / L ) 0 ), the linear constitutive laws (Newton’s law of viscosity and Fourier’s law of heat conduction) employed in the Navier-Stokes equations are sufficient to discern the flow physics. However, as we start moving into the slip regime ( 0.001 < K n < 0.1 ) and the transition regime ( 0.1 < K n < 10 ) by virtue of the reduction in the characteristic length scale (L) or increase in the mean free path ( λ ), these linear constitutive laws are shown to be inaccurate [1,2,3,4,5]. Several non-equilibrium flow phenomena [6,7,8] start arising in the flow that cannot be described using the Navier-Stokes equations. For such flows, it becomes necessary to employ particle methods such as the direct simulation Monte Carlo (DSMC) technique [9,10] or formulate higher-order continuum theories starting from the Boltzmann kinetic equation, for example, conventional Burnett equations [11,12], Grad 13 moment equations [13,14], Onsager-Burnett equations [15], and Onsager 13 moment equations [16]. These higher-order continuum theories have shown promise in modeling simple flow problems such as force-driven plane Poiseuille flow [17,18], shear-driven Couette flow [19,20], and the recently proposed Grad’s second problem [21,22]. However, the fundamental limitations of some of these theories, such as conventional Burnett equations and Grad moment equations, have also been exposed for shock waves, which have somewhat restricted the mainstream adoption of these higher-order continuum theories.
For the shock wave flow problem, the Mach number and Knudsen number are two important dimensionless numbers that characterize the flow problem. The Mach number is greater than unity and is defined as the ratio of the speed of the shock to the adiabatic sound speed. The Knudsen number is of the order of unity and falls well within the transition regime ( 0.1 < K n < 10 ). With such a high degree of non-equilibrium across the shock, the Navier-Stokes equations fail to quantitatively agree with the benchmark results for Mach numbers ( M a ) greater than 1.8 [23,24,25,26,27,28]. The solution to the same problem within the higher-order continuum framework was also not successful. For example, the conventional Burnett equations become invalidated for M a > 2.69 owing to the unstable nature of the equations, possible violation of the second law of thermodynamics, and the non-existence of the heteroclinic trajectory [25,29,30,31]. In case of the hyperbolic Grad 13 moment equations, discontinuity (sub-shocks) is observed in the shock profiles for M a > 1.65 [4,32]. Furthermor, the Grad distribution function assumes unphysical negative values for M a > 4 , invalidating the moment equations at higher Mach numbers [33]. Only recently [34], we showed that the OBurnett equations based on the Onsager-consistent approach do not suffer from these limitations and, at the same time, provide significant improvement over the results of Navier-Stokes equations for all Mach numbers. Without tweaking the equations in any way, the equations accurately resolved the upstream region of the shock, whereas some discrepancy was observed in the downstream region. In the present work, we modified the OBurnett equations by incorporating the Holian conjecture (HC) [35,36,37] in order to achieve quantitative agreement in the downstream region. In the present work, the modified theory (OBurnett + HC) was tested for a wide Mach number range including the case of a very strong shock ( M a = 134 ), and the results were benchmarked against the DSMC/non-equilibrium molecular dynamics (NEMD) results for a dilute gas system composed of hard-sphere molecules.
The remainder of this paper is organized as follows: The OBurnett equations are introduced in Section 2, with a brief introduction to the Onsager-consistent approach. In Section 3, the significance and superiority of the OBurnett equations over other higher-order continuum theories are highlighted. The problem definition for the normal shock wave flow problem is introduced in Section 4, and the reduced form of the equations based on the modified theory (OBurnett + HC) is dscribed in Section 5. The results are presented in Section 6, followed by the important conclusions in Section 7.

2. OBurnett Equations

The Onsager-Burnett (OBurnett) equations belong to a class of higher-order continuum transport equations that form a super-set of Navier-Stokes equations. In the derivation of these equations, the Onsager-consistent approach forms the basis that is completely independent of the Chapman-Enskog approach or the Grad moment method. In the Onsager-consistent approach, we ignore the traditional phenomenological approach of deriving continuum theories; rather, we start with the fundamental Boltzmann equation from kinetic theory, which is given as:
f t + c k f x k + F k f c k = J ( f , f 1 )
where f ( x , c , t ) is the single particle distribution function, which represents the solution of the Boltzmann equation; x is th location of a particle in physical space; c is the molecular velocity; t is the time variable; F i is the external force per unit mass acting on the molecules; and J ( f , f 1 ) represents the collision integral.
The transition from microscopic to macroscopic dynamics can be achieved by taking the appropriate moments of the Boltzmann equation. The macroscopic properties of interest, such as density, bulk velocity, and temperature, are related to the distribution function as
ρ ( x , t ) = m f ( x , c , t ) d c
u i ( x , t ) = 1 n c i f ( x , c , t ) d c
T ( x , t ) = 1 3 m 2 ρ k B C 2 f ( x , c , t ) d c
where n is the number density of the molecules, m is the molecular mass, and k B is the Boltzmann constant.
Multiplying the Boltzmann equation by { m , m c i , m c 2 / 2 } and then integrating over the velocity space, the basic conservation laws of mass, momentum, and energy can be obtained as
ρ t + ρ u k x k = 0 ,
ρ u i t + ρ u k u i x k + p x i + σ i k x k = ρ F i ,
ρ ϵ t + ρ u k ϵ x k + q k x k + p u k x k + σ i j u i x j = 0 ,
where ρ is the mass density, u k is the bulk velocity vector, p is the thermodynamic pressure, q i is the heat flux vector, σ i j is the stress tensor, F i is the external body force per unit mass, ϵ ( = 3 R T / 2 ) is the internal energy, and T is the absolute temperature.
It is easily recognized that the above system of equations is not closed as we have 5 equations and 13 unknowns. To obtain a closed set of equations, we need to substitute the constitutive relations for the stress tensor and the heat flux vector in terms of the primary variables and their gradients. In the traditional phenomenological approach employed while deriving the Navier-Stokes equations, the necessary closure is obtained by using Newton’s law of viscosity and Fourier’s law of heat conduction. However, in the kinetic approach, we evaluate the stress tensor and the heat flux vector as
σ i j ( x , t ) = m C i C j f ( x , c , t ) d c
q i ( x , t ) = m 2 C i C 2 f ( x , c , t ) d c .
Note that C is peculiar velocity given as C i = c i u i and the angular bracket denotes the symmetric and trace-free part of the tensor, given as C i C j = C i C j 1 3 C 2 δ i j .
Expressions for the stress tensor (Equation (8)) and the heat flux vector (Equation (9)) can only be evaluated once we know the form of the distribution function, i.e., we need to solve the Boltzmann equation for the distribution function. The Chapman-Enskog approach, the Grad moment method, and the Onsager-consistent approach focus on obtaining the form of the distribution function, and all these approaches are independent of one another. In the Onsager-consistent approach [15], the distribution function is represented in terms of the thermodynamic forces and fluxes [16,38,39] as
f O = f ( 0 ) Y τ : X τ + Y q · X q I order correction + ( Y τ τ X τ ) : X τ + ( Y q q X q ) · X q II order correction
where f ( 0 ) represents the Maxwell–Boltzmann distribution function and the subscript O in f O signifies that the distribution function is consistent with Onsager’s reciprocity principle. When the series is truncated until the Maxwell–Boltzmann distribution, we can obtain the inviscid Euler equations. Taking into account the first-order correction terms, the Navier-Stokes equations can be obtained, which are first-order accurate for Knudsen number. As we extend one step further by considering the second-order correction terms, the OBurnett equations can be obtained, which form a super-set of the Navier-Stokes equations and are second-order accurate for Knudsen number.
In the Onsager-consistent approach, the formulation of the distribution function is performed carefully so that it is consistent with Onsager’s symmetry principle [40,41] and the H-theorem. As these important principles are satisfied at the level of the distribution function, we ultimately obtain thermodynamically consistent theory at the macroscopic level. Th Onsager-consistent distribution function also satisfies the linearized Boltzmann equation and the collision invariance property. This particular form of the distribution function is then utilized to derive the OBurnett [15] (Burnett-like) and O13 [16] (Grad-like) equations.
Substituting the expression for the Onsager-consistent distribution function (Equation (10)) in Equations (8) and (9), the OBurnett constitutive relations for the stress tensor and the heat flux vector can be obtained as
σ x x = σ x x N S + σ x x B = μ δ 1 u x + δ 2 v y + δ 2 w z + 4 μ 2 β ρ [ α 1 u x 2 + α 2 u y 2 + α 3 u z 2 + α 4 u y v x + α 5 u z w x + α 6 w x 2 + α 7 v x 2 + α 8 u x v y + α 9 v y 2 + α 10 w z 2 + α 11 v y w z + α 12 u x w z + α 13 v z w y + α 14 w y 2 + α 15 v z 2 ] ,
σ x y = σ x y N S + σ x y B = μ δ 3 u y + μ δ 3 v x + 4 μ 2 β ρ [ β 1 u x u y + β 2 v x v y + β 3 u z v z + β 4 u x v x + β 5 u y v y + β 6 w x w y + β 7 v z w x + β 8 u z w y + β 9 u y w z + β 10 v x w z ] ,
q x = q x N S + q x B = δ 4 k 1 2 R β 2 β x + 4 μ 2 β ρ [ γ 1 1 β g x u x + γ 2 1 β 2 β x v y + γ 3 1 β 2 β x w z + γ 4 1 β g y u y + γ 5 1 β g y v x + γ 6 1 β g z w x + γ 7 1 β 2 β x u x + γ 8 1 β 2 β y u y + γ 9 1 β 2 β z u z + γ 10 1 β 2 β y v x + γ 11 1 β 2 β z w x + γ 12 1 β g x v y + γ 13 1 β g x w z ] + 2 k ( γ 1 ) R γ 2 1 ρ β [ γ 14 β y v x + γ 15 β z w x + γ 16 β x u x + γ 17 β y u y + γ 18 β z u z + γ 19 β x v y + γ 20 β x w z ] ,
where u, v, and w are the x, y, and z components of the bulk velocity vector, respectively; μ is the absolute viscosity; k is the thermal conductivity of the gas; and γ is the specific heat ratio. The expressions for β and g are given as
β = 1 2 R T , g = l o g ρ β
The coefficients, α , β , γ , and δ , with numerical subscripts, are functions of the type of gas and the interaction potential between the molecules. The constants appearing in front of the derivatives of flow field variables are
δ 1 = 4 3 , δ 2 = 2 3 , α 1 = 125 γ 2 576 γ + φ ( 110 160 γ + 50 γ 2 ) + 643 40 , α 2 = 4 5 , α 3 = 4 5 , α 4 = 1 5 , α 5 = 1 5 , α 6 = 3 5 , α 7 = 3 5 , α 8 = 100 γ 2 484 γ + φ ( 90 140 γ + 50 γ 2 ) + 459 20 , α 9 = 125 γ 2 392 γ + φ ( 70 120 γ + 50 γ 2 ) + 291 40 , α 10 = α 9 , α 11 = α 9 + 18 40 , α 12 = α 11 , α 13 = 2 5 , α 14 = 1 5 , α 15 = α 14 ,
δ 3 = 1 , β 1 = 23 γ + 5 γ φ 5 φ 37 10 , β 2 = β 1 , β 3 = 1 , β 4 = β 1 14 10 , β 5 = β 4 , β 6 = 4 β 4 , β 7 = 3 10 , β 8 = β 7 , β 9 = β 1
δ 4 = 1 , γ 1 = 47 + 25 γ 8 , γ 2 = φ 49 35 γ 8 , γ 3 = γ 2 , γ 4 = 1 2 , γ 5 = γ 4 , γ 6 = γ 4 , γ 7 = + φ 77 35 γ 8 , γ 8 = 7 φ 4 , γ 9 = γ 8 , γ 10 = γ 9 , γ 11 = γ 8 , γ 12 = 39 + 25 γ 8 , γ 13 = γ 12 , γ 14 = 1 , γ 15 = γ 14 , γ 16 = 77 + 35 γ + 10 φ ( 1 + γ ) 8 , γ 17 = 7 4 , γ 18 = γ 17 , γ 19 = 59 + 35 γ + 10 φ ( 1 + γ ) 8 , γ 20 = γ 19 .
To evaluate the Burnett contribution for other components of the stress tensor and heat flux vector, we apply a suitable change in variables in an appropriate base equation as given in Table 1.

3. Significance of the OBurnett Equations

The conventional Burnett equations, as is now well-known [31,39,42], are unstable, known to violate the second law of thermodynamics, and require additional boundary conditions for their complete solution, which are generally unknown. As such, the theory is far from complete and cannot be applied to complex boundary value problems. All these limitations associated with the conventional Burnett equations do not arise in the OBurnett equations. With respect to the boundary conditions part, the OBurnett equations require the same number of boundary conditions as the Navier-Stokes equations. This is evident when we notice the absence of second- and higher-order derivatives of velocity and temperature in the OBurnett constitutive relations. Stability is another important aspect wherein the OBurnett equations score above the conventional Burnett equations. The OBurnett equations are unconditionally stable and, hence, accurate solutions can be obtained on finer grids, which is not possible in the case of conventional Burnett equations. Hence, it is safe to say that the OBurnett equations form a complete theory and can serve as a monolithic solver for simulating flows starting from the continuum regime to the transition regime without any restrictions.
The accuracy of the OBurnett equations in describing non-equilibrium flow phenomena was successfully demonstrated in the force-driven plane compressible flow problem [18]. The equations were shown to capture the non-constant pressure profile, the tangential heat flux, and the presence of the normal stresses, in agreement with the DSMC results. After this initial success, the equations were applied to the shock wave flow problem, which is regarded as the most stringent test case for higher-order continuum theories [34,43]. This flow problem helped to establish several fundamental aspects of the OBurnett equations by providing strong evidence for the following:
1.
Smooth shock structures at all Mach numbers;
2.
Existence of the heteroclinic trajectory (a curve connecting the upstream and downstream equilibrium states);
3.
Positive entropy generation at all Mach numbers, implying the thermodynamic consistency of the equations.
It is important to highlight that neither the conventional Burnett equations nor the Grad 13 moment equations satisfy all these claims. For example, the heteroclinic trajectory does not exist for the conventional Burnett equations for M a > 2.69 and the equations predict negative entropy generation in the upstream region of the shock [25,30,43]. In the case of moment equations, since the equations are hyperbolic, the issue of subshocks arise and we cannot obtain smooth shock structures for M a > 1.65 [2,4,32,33]. This properties help to establish the superiority of the OBurnett equations over the conventional Burnett equations and the Grad moment equations.
While describing the shock structures for the hydrodynamic field variables using the OBurnett equations [34,43], the equations provide significant improvement over the results of the Navier-Stokes equations for a wide range of Mach numbers. Further comparison with the benchmark data (DSMC/NEMD results) revealed that the agreement is quantitative in the upstream region, whereas further improvement is desirable in the downstream region. This result was a bit surprising since the more rarefied upstream region is supposedly more difficult to resolve than the downstream region. In order to explain this behavior, we conjecture that the relation between the viscosity and the absolute temperature ( μ T φ ) is insufficient to capture the flow physics in the high-temperature, downstream region. To be more specific, the dissipation needs to be increased in the downstream region either by enhancing the viscosity index, as in [27,28], or by employing the Holian conjecture [35,36,37]. In the present work, we apply the Holian conjecture in the OBurnett equations and aim to verify whether the modified theory reasonably agrees with the benchmark data in the downstream region of the shock.

4. Problem Definition

A typical normal shock wave consists of two equilibrium states, namely the upstream state before the shock and the downstream state after the shock, as shown in Figure 1. The shock wave propagates with a constant velocity; hence, modeling the problem in the shock wave reference frame removes the time dependency in the governing equations. This simplification reduces the complexity of the problem, and an important aspect of the time of shock formation does not apply. Accordingly, all the hydrodynamic variables are now functions of the x coordinate only; hence, the problem can be modeled in a one-dimensional framework. As such, the velocity vector ( u i ), the stress tensor ( σ i j ), and the heat flux vector ( q i ) have only their x components, which are u, σ x x , and q x , respectively.
The upstream state is characterized by supersonic velocity ( u 0 ), whereas the velocity at the downstream state is subsonic ( u 1 ). The thermodynamic variables, density and temperature, rise sharply from the upstream state ( ρ 0 , T 0 ) to the downstream state ( ρ 1 , T 1 ) across the narrow region of the shock. These strong gradients of the field variables induces dissipative effects inside the shock region leading to the formation of a strong non-equilibrium region across the narrow width of the shock. Accordingly, the mass, momentum, and energy conservation equations (Equations (5)–(7)) can be readily obtained as
d ρ d x = 0 ,
ρ u d u d x = d p d x d σ d x ,
u d ρ ϵ d x = d q d x p d u d x σ d u d x
where the normal stress σ x x and x component of the heat flux vector q x are represented by σ and q, respectively. The thermodynamic pressure p appearing in the momentum and energy equation can be replaced using ideal gas equation as
p = ρ R T
Furthermore, for a dilute, monatomic gas system, the internal energy of a monatomic gas without internal degrees of freedom is given as
ϵ = 3 2 R T
Substituting for p and ϵ , the conservation Equations (15)–(17) can be written in a more convenient form as
d d x ρ u = 0 ,
d d x ρ u 2 + ρ R T + σ = 0 ,
d d x ρ u 3 + 5 ρ R T u + 2 u σ + 2 q = 0
Noting that the stresses and heat fluxes are absent before and after the shock ( σ = 0 and q = 0 ) owing to the zero gradients of the field variables, the Rankine-Hugoniot conditions can be obtained by performing integration between the upstream and downstream equilibrium states as
ρ 0 u 0 = ρ 1 u 1
ρ 0 u 0 2 + ρ 0 R T 0 = ρ 1 u 1 2 + ρ 1 R T 1
u 0 2 + 5 R T 0 = u 1 2 + 5 R T 1
where we use Equation (23) to obtain Equation (25).
To obtain the variation in the hydrodynamic field variables across the shock, we need to carry out the integration between the upstream state and any point inside the shock as
ρ 0 u 0 = ρ u ,
ρ 0 u 0 2 + ρ 0 R T 0 = ρ u 2 + ρ R T + σ ,
ρ 0 u 0 3 + 5 ρ 0 R T 0 u 0 = ρ u 3 + 5 ρ R T u + 2 u σ + 2 q
Note that the above equations are not in a closed form, and we need to substitute the constitutive relations for the normal stress and the heat flux. The OBurnett constitutive relations for the normal shock wave flow problem can be obtained as
σ = 4 3 μ d u d x + 2 α 1 μ 2 ρ R T d u d x 2
q = k d T d x + 4 γ 1 μ ρ 2 d ρ d x d u d x + 4 γ 1 γ 7 γ 16 P r 2 μ 2 ρ T d u d x d T d x
The space coordinate and other field variables are non-dimensionalized, as in [29,36]
s x l , l = 5 m 12 ρ 0 σ 2 π , τ k T m u 0 2 , u * u u 0 , ρ * ρ ρ 0 , τ 0 p 0 ρ 0 u 0 2 , σ * σ ρ 0 u 0 2 , q * q ρ 0 u 0 3
The quantities l and λ are different and related to each other as l = 1.04 λ . Note that λ is defined at the upstream point for the case of hard-sphere molecules, and the transport coefficients are proportional to square root of temperature as [12]
μ = 5 c μ 16 σ 2 m k T ( x ) π 1 / 2 , κ = 75 c λ 64 σ 2 k 3 T ( x ) π m 1 / 2
where σ is the particle diameter and the coefficients c μ and c λ are taken as 1.
By substituting the OBurnett constitutive relations for normal stress (Equation (29)) and heat flux (Equation (30)), the closed form of the momentum (27) and energy (28) equations can be obtained. The final non-dimensionalized form of the reduced OBurnett equations was derived in our earlier work (see Equations (28) and (29) in [34]). Without repeating the mathematical details, the final equations read as
τ d u * d s ̲ 9 8 α 1 u * d u * d s 2 = τ u * τ 0 1 + u *
45 16 τ d τ d s ̲ + 9 4 γ 1 τ d u * d s 2 9 4 Ψ u * d τ d s d u * d s = 3 2 ( τ τ 0 ) 1 2 1 u * 2 τ 0 1 u *
where Ψ = γ 1 γ 7 γ 16 / Pr 2 , where Pr is the Prandtl number, which is 2 / 3 for monatomic gases. The relevant OBurnett coefficients for the case of hard-sphere molecules ( φ = 0.5 ) are obtained as
α 1 = 8 15 , γ 1 = 2 3 , Ψ = 119 48 .
The Rankine-Hugoniot conditions in the non-dimensionalized form can be obtained as
u 1 * = 5 4 τ 0 + 1 4 , τ 1 = 7 8 τ 0 + 3 16 5 16 τ 0 2 , u 0 * = 1

5. Modified Theory: OBurnett + HC

In this section, we combine OBurnett theory and the Holian conjecture and present the reduced form of the equations for the modified theory. The main idea of the Holian conjecture [35,36] is to modify the viscosity–temperature relation ( μ T φ ) by replacing the average temperature (T) with the component of temperature in the direction of the shock propagation ( T x x ). The average temperature is related to the directional temperatures as
T = T x x + T y y + T z z 3 .
For a monatomic gas, the temperature component T x x can be extracted from the normal component of pressure tensor P x x as
P x x = ρ R T x x .
Note that the pressure tensor is related to the stress tensor as P i j = σ i j + p δ i j , so that the relation P x x = σ x x + p can be obtained.
The temperature in the direction of the shock propagation ( T x x ) is always greater than the average temperature (T). By replacing T with T x x in the viscosity-temperature relation, we are essentially increasing the viscosity across the shock region. This results in increased dissipation, thereby broadening the shock profiles for the case of strong shocks. This modification when applied to the Navier-Stokes equations provides reasonable agreement with the NEMD results for strong shocks ( M a = 134 ), as shown in [36].
Here, we apply the Holian conjecture in OBurnett theory by correcting only the Navier-Stokes-order terms of the stress tensor and the heat flux vector. We do not attempt to correct the Burnett-order terms since the magnitude of these terms will be small compared to that of the Navier-Stokes terms and most of the Burnett-order terms are actually independent of temperature, as also noted in [37]. Following this procedure, it can be shown that the term τ appearing in the underlined terms of Equations (32) and (33) needs to be replaced with τ x x , and the expression for τ x x is given as [35,36,44],
τ x x = u * ( τ 0 + 1 u * ) 1 / 2 .
The final form of the system of the ordinary differential equations for the modified theory (OBurnett + HC) then assumes the form
u * ( τ 0 + 1 u * ) 1 / 2 d u * d s 9 8 α 1 u * d u * d s 2 = τ u * τ 0 1 + u *
45 16 u * ( τ 0 + 1 u * ) 1 / 2 d τ d s + 9 4 γ 1 τ d u * d s 2 9 4 Ψ u * d τ d s d u * d s = 3 2 ( τ τ 0 ) 1 2 1 u * 2 τ 0 1 u *
The Rankine-Hugoniot boundary conditions given in Equation (35) remain same for the modified theory. The system of Equations (39) and (40) is tackled as an initial value problem as performed in the literature [25,29,36,45]. The adopted numerical procedure is the same as followed in our earlier works [34,43]. The MATLAB solver ode15i of order 5 was selected to solve the implicit system of ordinary differential equations with the relative and absolute errors set to 10 13 and 10 14 , respectively. Once we obtain the solution for velocity and temperature, other shock profiles for density, normal stress, and heat flux can be readily obtained as
ρ * = 1 u *
σ * = u * ( τ 0 + 1 u * ) 1 / 2 d u * d s + 9 8 α 1 u * d u * d s 2
q * = 45 16 u * ( τ 0 + 1 u * ) 1 / 2 d τ d s 9 4 γ 1 τ d u * d s 2 + 9 4 Ψ u * d τ d s d u * d s
The results of the shock profiles for three Mach numbers ( M a = 3, 9, and 134) are presented in the ensuing sections.

6. Results

In this section, the variation in the conserved variables ( u * , ρ * , and τ ) and the non-conserved variable ( q * ) across the narrow region of the shock, as obtained by combination of the OBurnett equations and the Holian conjecture, is compared with the benchmark DSMC results. To generate the DSMC results, Bird’s DSMC code [9] for the normal shock wave flow problem was implemented; details are provided in our earlier work [43]. By employing the combination of the OBurnett equations and the Holian conjecture, we first verified the thermodynamic consistency of the modified theory and the existence of the heteroclinic trajectory. The results indicate that the modified theory satisfies all these fundamental aspects, similar to the OBurnett equations [34].
The shock profiles for the velocity, density, temperature, and the heat flux for two Mach numbers, M a = 3 and 9, are shown in Figure 2 and Figure 3. The effect of the increased dissipation that is achieved by employing the Holian conjecture is not that significant for the conserved variables, as is evident from the velocity (Figure 2a and Figure 3a), density (Figure 2b and Figure 3b), and temperature (Figure 2c and Figure 3c) shock profiles. The modified theory marginally improves upon the results of the OBurnett equations in the downstream region due to increased dissipation. However, the same cannot be said for the upstream region, where the results of the modified theory and the OBurnett equations are almost indistinguishable. For the heat flux (Figure 2d and Figure 3d), the modified theory accurately captures the peak of the heat flux while improving upon the results of the OBurnett equations. An interesting phenomenon of temperature overshoot is observed in the downstream region of the shock where the temperature of the gas molecules exceeds that downstream. This behavior is confirmed by the experiments as well as the DSMC method and is probably due to the delayed relaxation of the transversal kinetic energy of the gas molecules. Although higher-order terms are added in the constitutive relation of heat flux vector in OBurnett theory, we observe a monotonous rise in the temperature across the shock as per the modified theory.
The case of a very strong shock with M a = 134 is also studied using the modified theory (Figure 4) and the results are compared against the non-equilibrium molecular dynamics (NEMD) results, as obtained by Salamons and Mareschal [26]. For the velocity (Figure 4a) and the density (Figure 4b) profiles, we observe that the modified theory provides significant improvement in the downstream region compared with the results of the OBurnett equations. However, the temperature profile is described much better by the OBurnett equations than the modified theory, both in the upstream and the downstream regions. With respect to the heat flux profile (Figure 4d), the results of the modified theory are in good agreement with those of the NEMD with substantial improvement achieved compared to the OBurnett results.
For all the three cases studied here ( M a = 3, 9, and 134), we observe that the shock thickness increases only slightly with the Mach number. This observation is consistent with the experimental and the DSMC results in the literature [46,47], where it is shown that the for higher Mach numbers, the shock thickness based on density profiles increases only slightly with Mach number.
In summary, the Holian conjecture when applied to OBurnett theory, a modified theory, performs well for the case of a very strong shock. However, for weak shocks, the modified theory provides marginal improvement over the results of the OBurnett equations. Hence, our hypothesis that the modified theory will be able to resolve the downstream region well only for the case of strong shocks. This also implies that the applicability of the Holian conjecture is only for strong shocks, as also noted in the works of [37,48].

7. Conclusions

In the present work, we proposed a modified theory by combining the OBurnett equations and the Holian conjecture (OBurnett + HC) to model the normal shock waves for a dilute gas system of hard-sphere molecules. Rather than using the average temperature, the longitudinal temperature in the direction of the shock propagation was used in the evaluation of the transport coefficients. The sole aim of this modification was to achieve increased dissipation in the hot region of the shock and thereby obtain better quantitative agreement with the benchmark results than the previously obtained OBurnett results. The modified theory still formed an entropic consistent system, a clear heteroclinic trajectory existed, and smooth shock structures were obtained at all Mach numbers, similar to OBurnett theory. For two Mach numbers tested ( M a = 3 and 9), the modified theory improved upon the OBurnett results for the heat flux shock profiles. However, the accuracy gain for the other hydrodynamic fields ( u * , ρ * , and τ ) was only marginal. The results for the case of a very strong shock ( M a = 134 ) were significantly improved using the modified theory, with excellent agreement observed for the heat flux profile. These results suggest that accounting for the Holian conjecture in the OBurnett equations works better for the case of strong shocks, whereas the original OBurnett theory seems to be a better choice for modeling weak shocks.

Author Contributions

Conceptualization and methodology, R.S.J. and A.A.; formal analysis, investigation, and writing—original draft preparation, R.S.J.; writing—review and editing, A.A.; supervision, A.A.; project administration and funding acquisition, A.A. All authors have read and agreed to the published version of the manuscript.

Funding

This work was performed as part of the Department of Atomic Energy Science Research Council Outstanding Investigator Award (21/01/2015-BRNS/35038).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
OBurnettOnsager-Burnett
HCHolian Conjecture
DSMCDirect Simulation Monte Carlo
NEMDNon-equilibrium Molecular Dynamics

References

  1. Gad-el Hak, M. The fluid mechanics of microdevices—The Freeman scholar lecture. J. Fluids Eng. 1999, 121, 5–33. [Google Scholar] [CrossRef]
  2. Torrilhon, M. Modeling Nonequilibrium Gas Flow Based on Moment Equations. Annu. Rev. Fluid Mech. 2016, 48, 429–458. [Google Scholar] [CrossRef]
  3. Uribe, F.; Garcia, A. Burnett description for plane Poiseuille flow. Phys. Rev. E 1999, 60, 4063–4078. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Grad, H. The profile of a steady plane shock wave. Commun. Pure Appl. Math. 1952, 5, 257–300. [Google Scholar] [CrossRef]
  5. Myong, R. Theoretical description of the gaseous Knudsen layer in Couette flow based on the second-order constitutive and slip-jump models. Phys. Fluids 2016, 28, 012002. [Google Scholar] [CrossRef]
  6. Balaj, M.; Roohi, E.; Mohammadzadeh, A. Regulation of anti-Fourier heat transfer for non-equilibrium gas flows through micro/nanochannels. Int. J. Therm. Sci. 2017, 118, 24–39. [Google Scholar] [CrossRef]
  7. Mohammadzadeh, A.; Rana, A.S.; Struchtrup, H. Thermal stress vs. thermal transpiration: A competition in thermally driven cavity flows. Phys. Fluids 2015, 27, 112001. [Google Scholar] [CrossRef] [Green Version]
  8. Hemadri, V.; Biradar, G.S.; Shah, N.; Garg, R.; Bhandarkar, U.V.; Agrawal, A. Experimental study of heat transfer in rarefied gas flow in a circular tube with constant wall temperature. Exp. Therm. Fluid Sci. 2018, 93, 326–333. [Google Scholar] [CrossRef]
  9. Bird, G.A.; Brady, J. Molecular Gas Dynamics and the Direct Simulation of Gas Flows; Clarendon press Oxford: Oxford, UK, 1994; Volume 42. [Google Scholar]
  10. Bird, G. The DSMC Method; CreateSpace Independent Publishing Platform: Scotts Valley, CA, USA, 2013. [Google Scholar]
  11. Burnett, D. The distribution of molecular velocities and the mean motion in a non-uniform gas. Proc. Lond. Math. Soc. 1936, 2, 382–435. [Google Scholar] [CrossRef]
  12. Chapman, S.; Cowling, T.G. The Mathematical Theory of Non-Uniform Gases: An Account of the Kinetic Theory of Viscosity, Thermal Conduction and Diffusion in Gases; Cambridge University Press: Cambridge, UK, 1970. [Google Scholar]
  13. Grad, H. On the kinetic theory of rarefied gases. Commun. Pure Appl. Math. 1949, 2, 331–407. [Google Scholar] [CrossRef]
  14. Grad, H. Principles of the Kinetic Theory of Gases. In Thermodynamik der Gase/Thermodynamics of Gases; Springer: Berlin/Heidelberg, Germany, 1958; pp. 205–294. [Google Scholar]
  15. Singh, N.; Jadhav, R.S.; Agrawal, A. Derivation of stable Burnett equations for rarefied gas flows. Phys. Rev. E 2017, 96, 013106. [Google Scholar] [CrossRef]
  16. Singh, N.; Agrawal, A. Onsager’s-principle-consistent 13-moment transport equations. Phys. Rev. E 2016, 93, 063111. [Google Scholar] [CrossRef] [PubMed]
  17. Tij, M.; Santos, A. Perturbation analysis of a stationary nonequilibrium flow generated by an external force. J. Stat. Phys. 1994, 76, 1399–1414. [Google Scholar] [CrossRef]
  18. Jadhav, R.S.; Singh, N.; Agrawal, A. Force-driven compressible plane Poiseuille flow by Onsager-Burnett equations. Phys. Fluids 2017, 29, 102002. [Google Scholar] [CrossRef]
  19. Lockerby, D.A.; Reese, J.M. High-resolution Burnett simulations of micro Couette flow and heat transfer. J. Comput. Phys. 2003, 188, 333–347. [Google Scholar] [CrossRef]
  20. Singh, N.; Gavasane, A.; Agrawal, A. Analytical solution of plane Couette flow in the transition regime and comparison with Direct Simulation Monte Carlo data. Comput. Fluids 2014, 97, 177–187. [Google Scholar] [CrossRef]
  21. Jadhav, R.S.; Agrawal, A. Grad’s second problem and its solution within the framework of Burnett hydrodynamics. J. Heat Transf. 2020, 142, 102105. [Google Scholar] [CrossRef]
  22. Jadhav, R.S.; Agrawal, A. Evaluation of Grad’s Second Problem Using Different Higher Order Continuum Theories. J. Heat Transf. 2021, 143, 012102. [Google Scholar] [CrossRef]
  23. Becker, R. Impact Waves and Detonation; Zeitschrift für Physik Transtation, N.A.C.A.-T.M. No. 505 (1929); Springer: Berlin/Heidelberg, Germany, 1929; Volume 8, pp. 321–362. [Google Scholar]
  24. Thomas, L. Note on Becker’s theory of the shock front. J. Chem. Phys. 1944, 12, 449–453. [Google Scholar] [CrossRef]
  25. Uribe, F.; Velasco, R.; Garcia-Colin, L.; Diaz-Herrera, E. Shock wave profiles in the Burnett approximation. Phys. Rev. E 2000, 62, 6648. [Google Scholar] [CrossRef] [PubMed]
  26. Salomons, E.; Mareschal, M. Usefulness of the Burnett description of strong shock waves. Phys. Rev. Lett. 1992, 69, 269. [Google Scholar] [CrossRef] [PubMed]
  27. Greenshields, C.J.; Reese, J.M. The structure of shock waves as a test of Brenner’s modifications to the Navier–Stokes equations. J. Fluid Mech. 2007, 580, 407–429. [Google Scholar] [CrossRef] [Green Version]
  28. Uribe, F.; Velasco, R. Shock-wave structure based on the Navier-Stokes-Fourier equations. Phys. Rev. E 2018, 97, 043117. [Google Scholar] [CrossRef] [PubMed]
  29. Uribe, F.J.; Velasco, R.M.; García-Colín, L.S. Burnett Description of Strong Shock Waves. Phys. Rev. Lett. 1998, 81, 2044–2047. [Google Scholar] [CrossRef]
  30. Comeaux, K.A.; Chapman, D.R.; MacCormack, R.W. An analysis of the Burnett equations based on the second law of thermodynamics. In Proceedings of the 33rd Aerospace Sciences Meeting and Exhibit, Reno, NV, USA, 9–12 January 1995; p. 415. [Google Scholar]
  31. Garcia-Colin, L.; Velasco, R.; Uribe, F. Beyond the Navier-Stokes equations: Burnett hydrodynamics. Phys. Rep. 2008, 465, 149–189. [Google Scholar] [CrossRef]
  32. Weiss, W. Continuous shock structure in extended thermodynamics. Phys. Rev. E 1995, 52, R5760–R5763. [Google Scholar] [CrossRef] [PubMed]
  33. Torrilhon, M.; Struchtrup, H. Regularized 13-moment equations: Shock structure calculations and comparison to Burnett models. J. Fluid Mech. 2004, 513, 171–198. [Google Scholar] [CrossRef] [Green Version]
  34. Jadhav, R.S.; Agrawal, A. Strong shock as a stringent test for Onsager-Burnett equations. Phys. Rev. E 2020, 102, 063111. [Google Scholar] [CrossRef]
  35. Holian, B.L. Modeling shock-wave deformation via molecular dynamics. Phys. Rev. A 1988, 37, 2562. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Holian, B.L.; Patterson, C.; Mareschal, M.; Salomons, E. Modeling shock waves in an ideal gas: Going beyond the Navier-Stokes level. Phys. Rev. E 1993, 47, R24. [Google Scholar] [CrossRef] [PubMed]
  37. He, Y.; Tang, X.; Pu, Y. Modeling shock waves in an ideal gas: Combining the Burnett approximation and Holian’s conjecture. Phys. Rev. E 2008, 78, 017301. [Google Scholar] [CrossRef] [PubMed]
  38. Mahendra, A.K.; Singh, R.K. Onsager reciprocity principle for kinetic models and kinetic schemes. arXiv 2013, arXiv:1308.4119. [Google Scholar]
  39. Agrawal, A.; Kushwaha, H.M.; Jadhav, R.S. Microscale Flow and Heat Transfer: Mathematical Modelling and Flow Physics; Springer: Cham, Switzerland, 2020. [Google Scholar]
  40. Onsager, L. Reciprocal relations in irreversible processes. I. Phys. Rev. 1931, 37, 405. [Google Scholar] [CrossRef]
  41. Onsager, L. Reciprocal relations in irreversible processes. II. Phys. Rev. 1931, 38, 2265. [Google Scholar] [CrossRef] [Green Version]
  42. Struchtrup, H. Macroscopic Transport Equations for Rarefied Gas Flows; Springer: Berlin/Heidelberg, Germany, 2005. [Google Scholar]
  43. Jadhav, R.S.; Gavasane, A.; Agrawal, A. Improved theory for shock waves using the OBurnett equations. J. Fluid Mech. 2021, 929, A37. [Google Scholar] [CrossRef]
  44. Holian, B. A History of constitutive modeling via molecular dynamics: Shock waves in fluids and gases. In Proceedings of the EPJ Web of Conferences. EDP Sciences, Poitiers, France, 4–9 July 2010; Volume 10, p. 00002. [Google Scholar]
  45. Gilbarg, D.; Paolucci, D. The Structure of Shock Waves in the Continuum Theory of Fluids. J. Ration. Mech. Anal. 1953, 2, 617–642. [Google Scholar]
  46. Alsmeyer, H. Density profiles in argon and nitrogen shock waves measured by the absorption of an electron beam. J. Fluid Mech. 1976, 74, 497–513. [Google Scholar] [CrossRef]
  47. Solovchuk, M.A.; Sheu, T.W. Prediction of strong-shock structure using the bimodal distribution function. Phys. Rev. E 2011, 83, 026301. [Google Scholar] [CrossRef] [Green Version]
  48. Velasco, R.; Uribe, F. A study on the Holian conjecture and Linear Irreversible Thermodynamics for shock-wave structure. Wave Motion 2021, 100, 102684. [Google Scholar] [CrossRef]
Figure 1. Schematic of the internal structure of a one-dimensional normal shock wave.
Figure 1. Schematic of the internal structure of a one-dimensional normal shock wave.
Fluids 06 00427 g001
Figure 2. Variation in hydrodynamics field (a) velocity, (b) density, (c) temperature, and (d) heat flux across the shock for M a = 3 .
Figure 2. Variation in hydrodynamics field (a) velocity, (b) density, (c) temperature, and (d) heat flux across the shock for M a = 3 .
Fluids 06 00427 g002
Figure 3. Variation in hydrodynamics field (a) velocity, (b) density, (c) temperature, and (d) heat flux across the shock for M a = 9 .
Figure 3. Variation in hydrodynamics field (a) velocity, (b) density, (c) temperature, and (d) heat flux across the shock for M a = 9 .
Fluids 06 00427 g003aFluids 06 00427 g003b
Figure 4. Variation in hydrodynamics field (a) velocity, (b) density, (c) temperature, and (d) heat flux across the shock for M a = 134 . Non-equilibrium molecular dynamics (NEMD) results were obtained from [26].
Figure 4. Variation in hydrodynamics field (a) velocity, (b) density, (c) temperature, and (d) heat flux across the shock for M a = 134 . Non-equilibrium molecular dynamics (NEMD) results were obtained from [26].
Fluids 06 00427 g004aFluids 06 00427 g004b
Table 1. Base equation and change in variables to be followed while evaluating the Burnett-order contribution for other components of the stress tensor and heat flux vector.
Table 1. Base equation and change in variables to be followed while evaluating the Burnett-order contribution for other components of the stress tensor and heat flux vector.
VariableBase EquationChange in Variables
σ y y σ x x u v , x y , v u and y x
σ z z σ x x u w , x z , w u and z x
σ y z σ x y u v , x y , v w , y z , w u and z x
σ z x σ x y u w , x z , v u , y x , w v and z y
q y q x u v , x y , v u and y x
q z q x u w , x z , w u and z x
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jadhav, R.S.; Agrawal, A. Shock Structures Using the OBurnett Equations in Combination with the Holian Conjecture. Fluids 2021, 6, 427. https://doi.org/10.3390/fluids6120427

AMA Style

Jadhav RS, Agrawal A. Shock Structures Using the OBurnett Equations in Combination with the Holian Conjecture. Fluids. 2021; 6(12):427. https://doi.org/10.3390/fluids6120427

Chicago/Turabian Style

Jadhav, Ravi Sudam, and Amit Agrawal. 2021. "Shock Structures Using the OBurnett Equations in Combination with the Holian Conjecture" Fluids 6, no. 12: 427. https://doi.org/10.3390/fluids6120427

APA Style

Jadhav, R. S., & Agrawal, A. (2021). Shock Structures Using the OBurnett Equations in Combination with the Holian Conjecture. Fluids, 6(12), 427. https://doi.org/10.3390/fluids6120427

Article Metrics

Back to TopTop