Next Article in Journal
Evaluation of SPH and FVM Models of Kinematically Prescribed Peristalsis-like Flow in a Tube
Next Article in Special Issue
Bulk Viscosity of Dilute Gases and Their Mixtures
Previous Article in Journal
Numerical Analysis Related to the ROCOM Pressurized Thermal Shock Benchmark
Previous Article in Special Issue
Six-Field Theory for a Polyatomic Gas Mixture: Extended Thermodynamics and Kinetic Models
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Navier–Stokes Equations and Bulk Viscosity for a Polyatomic Gas with Temperature-Dependent Specific Heats

1
Institute for Liberal Arts and Sciences, Kyoto University, Kyoto 606-8501, Japan
2
Department of Mathematics, National Cheng Kung University, Tainan 70101, Taiwan
*
Author to whom correspondence should be addressed.
Submission received: 11 November 2022 / Revised: 10 December 2022 / Accepted: 16 December 2022 / Published: 22 December 2022
(This article belongs to the Special Issue Bulk Viscosity and Relaxation Processes: Revisited)

Abstract

:
A system of Navier–Stokes-type equations with two temperatures is derived, for a polyatomic gas with temperature-dependent specific heats (thermally perfect gas), from the ellipsoidal statistical (ES) model of the Boltzmann equation extended to such a gas. Subsequently, the system is applied to the problem of shock-wave structure for a gas with large bulk viscosity (or, equivalently, with slow relaxation of the internal modes), and the numerical results are compared with those based on the ordinary Navier–Stokes equations. It is shown that the latter equations fail to describe the double-layer structure of shock profiles for a gas with large bulk viscosity.

1. Introduction

Nonequilibrium flows of polyatomic gases have attracted increasing attention in connection with practical problems in high-temperature gas dynamics [1,2,3,4]. The basis of the study of such flows is kinetic theory of gases based on the Boltzmann Equation [3,5,6,7,8]. However, since the Boltzmann equation has a complicated structure particularly for polyatomic gases, simplified Boltzmann equations for a polyatomic gas, in which the internal states are represented by a single variable associated with the energy of internal modes, have been proposed [9,10,11,12,13,14]. It should be mentioned that recent progress in the mathematical study of these Boltzmann equations is remarkable [12,14,15,16,17].
However, even the simplified versions of the Boltzmann equation are still not simple enough for practical applications. Therefore, alternative approaches using model Boltzmann equations [18,19,20,21,22,23,24,25,26,27] or macroscopic equations [28,29,30,31,32,33,34,35] have been proposed. The model Boltzmann equations bear much simpler structures while retaining certain basic properties of the original Boltzmann equation. The ellipsoidal-statistical (ES) model [21,22] is one of such models and is relevant herein. The macroscopic equations partially take account of kinetic effects of polyatomic gases, such as bulk viscosity, empirically or systematically on the basis of kinetic theory.
In a recent paper by the present authors, M. Bisi, and M. Groppi [36], a set of Navier–Stokes-type equations with two temperatures [two-temperature Navier–Stokes (NS) equations for short] was derived on the basis of kinetic theory and applied successfully to the structure of a plane shock wave for a polyatomic gas with large bulk viscosity. More specifically, the set was derived by the Chapman–Enskog procedure [37,38,39,40] using the ES model under the assumption that the relaxation of the internal modes is slow. It should be noted that, in this model, this assumption corresponds to large bulk viscosity.
The two-temperature NS equations in [36] are for a gas with constant specific heats (calorically perfect gas), since they were derived using the ES model for such a gas. In reality, however, the specific heats for most gases are temperature-dependent even when treated conceptually as ideal gases (thermally perfect gas). The effect of temperature-dependent specific heats becomes particularly important when there are significant variations in temperature in the flow field, such as gas flows containing strong shock waves.
For this reason, herein, we attempt to extend the two-temperature NS equations to a gas with temperature-dependent specific heats. Since the ES model has been extended to such a gas in [41], we will adopt this extended ES model as the starting point. Subsequently, an analysis can be conducted in parallel to [36], and given that the analysis details have already been published, we simply provide the result herein. We then shift our attention to a comparison between the ordinary Navier–Stokes (NS) equations with a single temperature, in which the bulk viscosity is contained, and the two-temperature NS equations, in which the bulk viscosity does not occur. The comparison is performed numerically for the structure of a plane shock wave.
Here, the following comments are in order. Navier–Stokes-type equations with two or multiple temperatures have widely been used in computational fluid dynamics for hypersonic flows of real gases in connection with reentry flight and contributed to the development of aerothermodynamics [2,3,42,43]. These flows comprise mixtures of many component gases with dissociation and chemical reactions. Because of the complexity of the phenomena, the equations are not strictly based on kinetic theory, and various modelings are introduced in many different processes. In this sense, the approach is different from ours that originates from kinetic theory in a simpler and idealized setting. Moreover, it is worth mentioning that the relation between multi-temperature fluid models and the bulk viscosity has been discussed on the basis of kinetic theory [44,45].
The rest of the paper is organized as follows. We summarize the extended ES model for a gas with temperature-dependent specific heats and its basic properties in Section 2. Next, the ordinary NS equations with a single temperature and the two-temperature NS equations, both derived from the extended ES model, are summarized in Section 3 and Section 4, respectively. Brief remarks on the relation between the bulk viscosity and the presence of the two temperatures, the translational and internal temperatures, are given in Section 5. In Section 6, the shock-structure problem is analyzed numerically for certain gases with large bulk viscosities using the two different NS equations, and the difference in the shock profiles is discussed. Section 7 provides concluding remarks.

2. ES Model for a Gas with Temperature-Dependent Specific Heats

Let us consider a polyatomic rarefied gas. Let t be the time variable, X (or X i ) the position vector in the physical space, ξ (or ξ i ) the molecular velocity, and E the energy associated with the internal modes per unit mass. We denote by
f ( t , X , ξ , E ) d X d ξ d E ,
the total mass of the gas molecules contained in an infinitesimal volume d X d ξ d E around a point ( X , ξ , E ) in the seven-dimensional space consisting of X , ξ , and E at time t. We may call f ( t , X , ξ , E ) the velocity–energy distribution function of the gas molecules.

2.1. ES Model

In the present paper, the gas is assumed to be thermally perfect, that is, the specific heat at constant volume c v and that at constant pressure c p are not constant but are functions of the temperature T. For such a gas, an ES model governing the velocity-energy distribution function f was proposed in [41]. It is a straightforward extension of the ES model for a gas with constant c v and c p [21]. In this subsection, we summarize the ES model proposed in [41].
We first consider an equilibrium state with a uniform and constant temperature T and assume for the moment that the specific heats are independent of T and that the classical equipartition law holds. Then, the internal energy E of the gas per unit mass can be expressed as E = ( 3 + δ 0 ) R T / 2 , where δ 0 is the (constant) number of internal degrees of freedom, and R is the gas constant per unit mass (note that R = k B / m with k B the Boltzmann constant and m the mass of a gas molecule). Now, we consider a more general gas for which the internal energy E is a given (monotonically increasing) function of the temperature T, i.e., E = E ( T ) . Then, we define a function D ( T ) as D ( T ) = 2 E ( T ) / R T 3 . This D ( T ) is a generalized number of internal degrees of freedom. Note that the relations c v ( T ) = d E ( T ) / d T and c p ( T ) = c v ( T ) + R hold, so that c v and c p depend on the temperature.
Next, we apply the functions E ( T ) , D ( T ) , c v ( T ) , and c p ( T ) , defined in the equilibrium, to the nonequilibrium situation and use them in the construction of the ES model for the thermally perfect gas. We omit the process of the construction (see [41]) and simply show the result. The ES model is given by the following equation:
f t + ξ i f X i = Q ( f ) ,
where
Q ( f ) = A c ( T ) ρ ( G f ) .
Here,
G = ρ E δ / 2 1 ( 2 π ) 3 / 2 [ det ( T ) ] 1 / 2 ( R T rel ) δ / 2 Γ ( δ / 2 ) × exp 1 2 ( T 1 ) i j ( ξ i v i ) ( ξ j v j ) E R T rel ,
( T ) i j = ( 1 θ ) [ ( 1 ν ) R T tr δ i j + ν p i j / ρ ] + θ R T δ i j ,
ρ = 0 f d E d ξ , v i = 1 ρ 0 ξ i f d E d ξ ,
p i j = 0 ( ξ i v i ) ( ξ j v j ) f d E d ξ ,
T = E 1 ( e ) , δ = D ( T ) = 2 E ( T ) / R T 3 = 2 e / R T 3 ,
T tr = 2 e tr / 3 R , T int = 2 e int / R δ , T rel = θ T + ( 1 θ ) T int ,
where e, e tr , and e int are defined by
e = e tr + e int , e tr = 1 2 ρ 0 | ξ v | 2 f d E d ξ , e int = 1 ρ 0 E f d E d ξ .
In (2)–(5), ρ is the density, v (or v i ) is the flow velocity, p i j is the stress tensor, e is the internal energy per unit mass, e tr is that associated with the translational motion, e int is that associated with the internal modes, T is the temperature, T tr is the temperature associated with the translational motion, T int is the temperature associated with the energy of the internal modes, d ξ = d ξ 1 d ξ 2 d ξ 3 , and the domain of integration with respect to ξ is the whole space of ξ . The symbol δ i j indicates the Kronecker delta, and ν [ 1 / 2 , 1 ) and θ [ 0 , 1 ] are parameters. In addition, A c ( T ) is a function of T such that A c ( T ) ρ is the collision frequency of the gas molecules, Γ ( z ) = 0 s z 1 e s d s is the gamma function, T is the 3 × 3 positive-definite symmetric matrix whose ( i , j ) component is defined by (4b), and det ( T ) and T 1 are, respectively, its determinant and inverse.
Note that all the macroscopic quantities contained in G are generated from f. To be more specific, (i) ρ , v , p i j , e tr , e int , and e are obtained by (4c), (4d) and (5); (ii) T and then δ are determined by (4e) using the inverse function E 1 of the function E; (iii) T tr , T int , and T rel are determined by (4f), and then T is established by (4b). Since e = e tr + e int = ( 3 T tr + δ T int ) R / 2 and also e = ( 3 + δ ) R T / 2 , we have the relation T = ( 3 T tr + δ T int ) / ( 3 + δ ) [note that δ depends on T: cf. (4e)].
The pressure p and the heat-flow vector q i are given by
p = R ρ T ,
q i = 0 ( ξ i v i ) 1 2 | ξ v | 2 + E f d E d ξ ,
where (6a) is the equation of state.

2.2. Basic Properties

In this subsection, we summarize the basic properties that follow from the model Equation (2) [41]. Some properties are different when the parameter θ is equal to zero.
  • Conservations: For an arbitrary function f ( t , X , ξ , E ) , the following relation holds:
    0 φ r Q ( f ) d E d ξ = 0 ,
    where φ r ( r = 0 , , 4 ) are φ 0 = 1 , φ i = ξ i ( i = 1 , 2 , 3 ), and φ 4 = 1 2 | ξ | 2 + E for θ 0 ; and φ r ( r = 0 , , 5 ) are φ 0 = 1 , φ i = ξ i ( i = 1 , 2 , 3 ), φ 4 = 1 2 | ξ | 2 , and φ 5 = E for θ = 0 .
  • Equilibrium: The vanishing of the collision term Q ( f ) = 0 is equivalent to the fact that f is the following local equilibrium distribution:
    f eq = ρ ¯ E δ ¯ / 2 1 ( 2 π R T ¯ ) 3 / 2 ( R T ¯ ) δ ¯ / 2 Γ ( δ ¯ / 2 ) exp | ξ v ¯ | 2 2 R T ¯ E R T ¯ ,
    for θ 0 , where ρ ¯ , v ¯ , and T ¯ are arbitrary functions of t and X , and δ ¯ = D ( T ¯ ) ; and
    f eq = ρ ¯ E δ ¯ / 2 1 ( 2 π R T ¯ tr ) 3 / 2 ( R T ¯ int ) δ ¯ / 2 Γ ( δ ¯ / 2 ) exp | ξ v ¯ | 2 2 R T ¯ tr E R T ¯ int ,
    for θ = 0 , where ρ ¯ , v ¯ , T ¯ tr , and T ¯ int are arbitrary functions of t and X , and δ ¯ and T ¯ are determined by the following coupled equations:
    δ ¯ = D ( T ¯ ) , T ¯ = E 1 ( 3 R T ¯ tr / 2 + δ ¯ R T ¯ int / 2 ) .
    The solution ( δ ¯ , T ¯ ) of (10) exists. In particular, it is unique when T ¯ int T ¯ tr .
  • Entropy inequality: For an arbitrary function f ( t , X , ξ , E ) , the following inequality holds:
    0 ln f E δ / 2 1 Q ( f ) d E d ξ 0 ,
    where δ is defined by (4e) and (5). The equality sign holds if and only if f = f eq in (8) ( θ 0 ) or (9) ( θ = 0 ).
  • Mean free path: The mean free path l 0 of the gas molecules in the equilibrium state at rest, at density ρ 0 , and temperature T 0 is given by
    l 0 = 2 π ( 2 R T 0 ) 1 / 2 A c ( T 0 ) ρ 0 ,
    for the ES model (2), since A c ( T 0 ) ρ 0 is the collision frequency at this equilibrium state.
  • Relaxation: In the space homogeneous case, where f = f ( t , ξ , E ) , the macroscopic quantities ρ , v , and T are constant, and the temperatures T tr and T int relax to T in the following manner:
    T tr = T + ( T tr * T ) e θ A c ( T ) ρ t ,
    T int = T + ( T int * T ) e θ A c ( T ) ρ t .
    where T tr * and T int * are, respectively, the values of T tr and T int at the initial time t = 0 satisfying the relation ( 3 T tr * + δ T int * ) / ( 3 + δ ) = T with δ = D ( T ) . Since 1 / A c ( T ) ρ is the mean free time, the time scale of relaxation of the temperatures is the mean free time devided by θ . Therefore, we may regard θ as the parameter controlling the speed of relaxation of the internal modes. This situation is the same as in the case of the ES model with constant δ [36].
It should be remarked that the inequality (11) enables us to show the H-theorem only for the spatially homogeneous case [41]. It is a drawback of this model, which should be the consequence of the simple and convenient definition of the temperature T and the resulting δ using the given function E [cf. (4e)].

3. Ordinary Navier–Stokes Equations with a Single Temperature

When the mean free path l 0 of the gas molecules in the reference equilibrium state at rest with density ρ 0 and temperature T 0 is much shorter than the characteristic length L of the system, i.e., when the Knudsen number Kn = l 0 / L is small, one can formally derive fluid-dynamic equations by a classical procedure, such as the Chapman–Enskog expansion, from the ES model (2). In this section, restricting ourselves to the case where θ is of the order of unity, we summarize the Navier–Stokes equations derived in [41] by the standard Chapman–Enskog procedure.
The macroscopic quantities ρ , v i , and T are governed by the following Navier–Stokes equations for a compressible fluid:
ρ t + ρ v j X j = 0 , ρ v i t + ρ v i v j X j = p X i + X j μ ( T ) v i X j + v j X i 2 3 v k X k δ i j
+ X i μ b ( T ) v k X k , t ρ 3 + δ 2 R T + 1 2 | v | 2 + X j ρ v j 5 + δ 2 R T + 1 2 | v | 2 = X j λ ( T ) T X j + X j μ ( T ) v i v i X j + v j X i 2 3 v k X k δ i j
+ X j μ b ( T ) v j v k X k ,
with p = R ρ T [cf. (6a)] and δ = D ( T ) . Here, the viscosity μ ( T ) , the bulk viscosity μ b ( T ) , and the thermal conductivity λ ( T ) are given by
μ ( T ) = Pr R T A c ( T ) , μ b ( T ) = 1 θ 5 3 γ ( T ) 1 Pr μ ( T ) , λ ( T ) = γ ( T ) γ ( T ) 1 R R T A c ( T ) , Pr = 1 1 ν + θ ν ,
where γ ( T ) = c p ( T ) / c v ( T ) = [ c v ( T ) + R ] / c v ( T ) is the ratio of the specific heats, and Pr = c p ( T ) μ ( T ) / λ ( T ) is the Prandtl number. In addition, the translational and internal temperatures T tr and T int are, respectively, given as
T tr = T 1 R μ b ( T ) 1 ρ v k X k , T int = T + 3 R δ μ b ( T ) 1 ρ v k X k .
One of the advantages of the present ES model is that the resulting transport coefficients (15) are simple and explicit, so that the Navier–Stokes eqautions (14) have a wide applicability for practical flows of a polyatomic gas with temperature-dependent specific heats.

4. Navier–Stokes Equations with Two Temperatures

In this section, we consider the case where the parameter θ contained in the ES model (2) is small. As was seen in Section 2.2, this corresponds to the case with slow relaxation of the internal modes. On the other hand, Equation (15) indicates that small values of θ correspond to large values of the ratio μ b / μ of the bulk viscosity to the viscosity.
In [36], the case where θ is small and of the same order of magnitude as the Knudsen number Kn ( 1 ) was considered for the ES model for a gas with constant specific heats, and fluid-dynamic equations were derived by an adapted Chapman–Enskog expansion. The fluid-dynamic equations are of Navier–Stokes type containing two temperatures, the temperature associated with the energy of the internal modes and that associated with the translational energy, and relaxation terms of the temperatures.
If we let
θ = α ϵ , ϵ = ( π / 2 ) Kn 1 ,
where α is a positive constant of the order of unity, in the ES model (2) and carry out almost the same analysis using the Chapman–Enskog expansion for small ϵ as in [36], we obtain the equations of Navier–Stokes type with two temperatures. In this paper, the process of analysis is omitted, and only the result is shown. The dimensional form of the resulting equations is as follows:
ρ t + ρ v j X j = 0 ,
ρ v i t + ρ v i v j X j = ρ R T tr X i + X j μ ( T , T tr ) v i X j + v j X i 2 3 v k X k δ i j , t ρ 3 2 R T tr + 1 2 | v | 2 + X j ρ v j 5 2 R T tr + 1 2 | v | 2 + θ A c ( T ) ρ 2 · 3 2 R T tr T
= X j λ ( T , T tr ) T tr X j + X j μ ( T , T tr ) v i v i X j + v j X i 2 3 v k X k δ i j , t ρ δ 2 R T int + X j ρ v j δ 2 R T int + θ A c ( T ) ρ 2 · δ 2 R T int T
= 1 5 X j λ ( T , T tr ) δ T int X j ,
with
T = 3 T tr + δ T int 3 + δ , δ = D ( T ) ,
where the viscosity μ ( T , T tr ) and the thermal conductivity λ ( T , T tr ) are given, respectively, by
μ ( T , T tr ) = 1 1 ν R T tr A c ( T ) , λ ( T , T tr ) = 5 2 R R T tr A c ( T ) .
With the known functions A c ( T ) and D ( T ) , the system (18)–(20) form the equations for ρ , v i , T, T tr , and T int . It is shown in Appendix A in [41] that a solution ( T , δ ) of (19) exists for given T tr and T int under the assumption that E ( T ) > 3 R T / 2 . When T int < T tr , T is such that T int < T < T tr , and the solution ( T , δ ) is unique; when T tr < T int , T is such that T tr < T < T int , and the possibility that the solution is not unique is not excluded depending on the form of E ( T ) (see Appendix A in [41]). If the solution is not unique, it is not clear if there is a criterion to choose one of the solutions. In the numerical examples of shock-wave structure in Section 6.3, T int < T tr holds in all cases. Finally, it should be recalled that θ in (18) is small as indicated by (17).

5. Some Remarks

The ordinary NS equations (14) are based on the Navier–Stokes constitutive relations, i.e., the following expressions of the stress tensor p i j and the heat-flow vector q i :
p i j = p δ i j μ ( T ) v i X j + v j X i 2 3 v k X k δ i j μ b ( T ) v k X k δ i j ,
q i = λ ( T ) T X i .
Here, we should recall that Equation (21) has been derived under the assumption that θ is of the order of unity. On the other hand, the viscosity μ and the bulk viscosity μ b make sense only for the stress tensor p i j of the form of (21a). Suppose that we are able to measure, experimentally, the tangential and normal stresses as well as the flow velocity, pressure, and temperature. Then, the viscosity and the bulk viscosity are determined empirically with the help of (21a). It is known that there are many gases that have large values of the ratio μ b / μ . However, it follows from (15) that large μ b / μ means small θ . Therefore, for the gases with large μ b / μ , the appropriate fluid-dynamic equations are not (14), derived under the assumption θ = O ( 1 ) , but (18), derived under (17). This fact was already pointed out in [36] for a gas with constant specific heats. In the next section, we will compare the shock-wave structure based on (18) and that based on (14) for some gases with large μ b / μ numerically and find that they are significantly different.
The Chapman–Enskog expansion has been carried out for the dimensionless form of (2) in terms of the small parameter ϵ in (17). Equations (14) and (18) retain the terms of O ( 1 ) and O ( ϵ ) though it is not visible from the dimensional Equations (14) and (18). It is noted that the terms containing μ ( T ) , μ b ( T ) , and λ ( T ) in (14) and those containing μ ( T , T tr ) , λ ( T , T tr ) , and θ in (18) correspond to the terms of O ( ϵ ) . Equation (14) with p = R ρ T and δ = D ( T ) determines ρ , v , T, and p; then, Equation (16) gives T tr and T int as small perturbations of O ( ϵ ) from T. In this way, T tr and T int play a subsidiary role.
When μ b / μ is large, as can be guessed from (16), T tr and T int are not small perturbations from T and play a principal role. Therefore, we need the principal equations governing T tr and T int , which are given by the system (18). The transition from (18) to (14) for θ such that O ( ϵ ) θ O ( 1 ) can be shown by following the procedure in [44], as in the case of a gas with constant specific heats (see Appendix B in [46]).

6. Shock-Wave Structure

The problem of shock-wave structure is one of the most fundamental problems in kinetic theory and has been investigated for long time [6,39,40,47]. In this section, we investigate the problem of shock-wave structure for some gases with large values of μ b / μ using the two-temperature NS equations (18) as well as the ordinary NS equations (14) and compare the results.
A typical example of a polyatomic gas with large μ b / μ is CO 2 , and the study of shock-wave structure for CO 2 gas, which has a long history [48,49,50,51], has been revived in recent years [31,33,35,52,53,54,55,56,57]. For instance, an interesting classification of shock profiles depending on the upstream Mach number was provided in [52,53] using macroscopic equations based on extended thermodynamics. The results in [31,52,53] were subsequently confirmed by kinetic theory, more precisely, by the numerical analysis using the ES model in [41,58].
The first aim of this section is to confirm that the two-temperature NS equations (18) are able to reproduce the shock profiles obtained by the ES model (2) for a gas with temperature dependent specific heats [41]. The second aim is to assess the validity of the ordinary NS equations (14) with large bulk viscosity by comparing the obtained shock profiles with those based on the two-temperature NS equations (18).

6.1. Problem

We consider a stationary plane shock wave in a uni-directional flow of an ideal polyatomic gas with large values of μ b / μ , i.e., with slow relaxation of the internal modes, perpendicular to the shock wave. The X 1 axis of the coordinate system is taken in the flow direction. The gas at upstream infinity ( X 1 ) is in an equilibrium state with density ρ , flow velocity v = ( v , 0 , 0 ) , and temperature T , while the gas at downstream infinity ( X 1 + ) is in another equilibrium state with density ρ + , flow velocity v + = ( v + , 0 , 0 ) , and temperature T + . We investigate the steady behavior of the gas under the following assumptions: (i) the specific heat at constant volume c v is a known function of the temperature T; (ii) the behavior of the gas is described by the Navier–Stokes equations with a single temperature (14) or those with two temperatures (18); (iii) the problem is spatially one dimensional, so that the physical quantities are independent of X 2 and X 3 .
Let us denote by M the Mach number at upstream infinity and by γ the ratio of the specific heats at the temperature T , i.e.,
M = v γ R T , γ γ ( T ) = c v ( T ) + R c v ( T ) .
When θ 0 , the downstream parameters ρ + , v + , T + and the upstream ones ρ , v , T are related by the following Rankine–Hugoniot relations [41]:
ρ + ρ = v + v 1 , v + v = 1 + γ M 2 2 γ M 2 d ^ E ( τ ) + 1 γ M 2 , T + T = τ ,
where
d ^ E ( x ) = 1 R 1 x c v ( T s ) d s ,
and τ is a solution such that τ > 1 of the following equation:
τ + 2 d ^ E ( τ ) + 1 γ M 2 1 γ M + γ M 2 d ^ E ( τ ) + 1 γ M 2 = 0 .
If c v ( T ) is the increasing function of T, the solution satisfying τ > 1 is unique [41].
We try to obtain the numerical solutions of the steady and one-dimensional versions of (14) and (18), where / t = / X 2 = / X 3 = 0 and v 2 = v 3 = 0 are assumed, under the following boundary conditions:
( ρ , v 1 , T ) ( ρ , v , T ) , as X 1 , [ for ( 14 ) ] ,
( ρ , v 1 , T tr , T int ) ( ρ , v , T , T ) , as X 1 , [ for ( 18 ) ] .
In the actual computation, however, we solve corresponding time-dependent problems and obtain the steady shock profiles as the long-time limit of the unsteady solutions. Since the numerical method is essentially the same as that in [36], we omit its description here.

6.2. Properties of Some Gases with Large Bulk Viscosities

In Table IV in [46], thermophysical properties of some gases with large ratios of the bulk viscosity to the viscosity are listed. In this paper, we choose SF 6 , CH 4 , and CO 2 from the list and investigate the shock-wave structure for these gases.
We first consider the specific heat at constant volume c v ( T ) according to the data of that at constant pressure c p ( T ) in [59,60]. We use the following fitting formula for c v ( T ) :
c v ( T ) / R = c 0 + c 1 T + c 2 T 2 + c 3 T 3 + c 4 T 4 ,
where c 0 , , c 4 are fitting coefficients having units such that each term on the right-hand side is dimensionless. The values of c 0 , , c 4 listed in Table 1 are those determined by using the data in [59,60] for SF 6 and CH 4 and those taken from [31] for CO 2 . In the table, the values of c v ( T ) / R , γ [ = γ ( T ) ] , and δ D ( T ) when T = 300 K are also shown. Figure 1 shows the comparison between the fitting formula (27) with the values in Table 1 and the data in [59,60]. It is seen that the formula (27) reproduces the data well for the temperature range between 220 K and 1100 K. If T = 300 K and M = 5 , then T + = 1.938 T 581 K for SF 6 , T + = 3.408 T 1022 K for CH 4 , and T + = 3.708 T 1112 K for CO 2 .
Other thermophysical properties for SF 6 , CH 4 , and CO 2 are included in Table IV in [46]. The data for these three gases are shown in Table 2. The data for the viscosity μ , the ratio of the specific heats γ , and the Prandtl number Pr, which are at T = 300 K and p = 101 kPa, are taken from [59], whereas those for the ratio μ b / μ , which are at T = 300 K, are taken from [61]. Note that for the ES model, μ depends only on T. From these data, one obtains the values of θ and ν in Table 2 by using (15). As expected, the values of θ are small.
In [46], the viscosity μ ( T ) was assumed to be in the form
μ ( T ) = μ ( 300 K ) × ( T / 300 K ) s ,
and the exponent s was determined by using the method of least squares on the basis of the data provided in [59]. The obtained value of s for each gas is shown in Table 2. Equation (28) with these values of s reproduces the data of μ given in [62] quite well in the temperature range 250–1300 K for CO 2 and 250–600 K for SF 6 and CH 4 . The function A c ( T ) is then determined by (15) as
A c ( T ) = A c ( 300 K ) × T 300 K 1 s , A c ( 300 K ) = ( R Pr ) × 300 K μ ( 300 K ) .
With this A c ( T ) , the reference mean free path l 0 is obtained by (12).

6.3. Numerical Results

In this section, we show some numerical results for the three gases mentioned in Section 6.2, i.e., CO 2 , SF 6 , and CH 4 . The shock-wave structure based on the two-temperature NS equations (18) is compared with that based on the ordinary NS equations (14). For CO 2 gas, it is also compared with the result based on the original ES model (2) taken from [41].
As in [41,58], we show the profiles of the density ρ , the flow velocity v 1 (the X 1 component), and the temperatures T, T tr , and T int normalized in the conventional way, i.e.,
ρ ˇ = ρ ρ ρ + ρ , v ˇ = v 1 v + v v + , T ˇ = T T T + T , T ˇ tr = T tr T T + T , T ˇ int = T int T T + T .
In this normalization, ρ ˇ , T ˇ , T ˇ tr , and T ˇ int varies from 0 (upstream infinity) to 1 (downstream infinity), whereas v ˇ from 1 (upstream infinity) to 0 (downstream infinity). In the figures in this section, the X 1 coordinate is normalized by the mean free path l at the equilibrium state at rest with density ρ and temperature T , i.e., l = ( 2 / π ) ( 2 R T ) 1 / 2 / A c ( T ) ρ [cf. (12)], and the origin X 1 = 0 is set in such a way that ρ ˇ obtained by (18) and that by (14) coincide at X 1 = 0 . Otherwise, the absolute location of the origin is unimportant in the present problem. Here, it is put around the position where the density exhibits the steepest change, so that we can use the grid points concentrating on the origin in the numerical analysis.

6.3.1. CO 2 Gas

First, it should be noted that the calculation for CO 2 gas is carried out for the parameters slightly different from those shown in Table 1 and Table 2 because one of the purposes here is to compare the result with that based on the original ES model (2) reported in [41]. In [41], T and p are set to be T = 295 K and p = 69 mmHg, respectively, in accordance with [31,53], and this T gives c v ( T ) / R = 3.456 , γ = 1.289 , and δ = 3.913 . In addition, although the fitting formula (27) with the data in Table 1 is used for c v ( T ) / R , the exponent s in (28) and (29) and Pr are set as s = 0.935 and Pr = 0.73 , respectively.
One of the interests in this study is to observe the change of the shock-wave structure as the ratio μ b / μ increases. For this purpose, we do not fix the value of μ b / μ according to the thermophysical data, such as 3849 in Table 2, but specify it arbitrary; we assume μ b / μ = 10 , 100 , and 1000 in the following. In other words, our CO 2 gas is a fictitious one that we called the pseudo-CO 2 gas in [58]. Correspondingly, the parameter θ becomes θ = 5.169 × 10 2 , 5.169 × 10 3 , and 5.169 × 10 4 , respectively.
Figure 2, Figure 3 and Figure 4 show the profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int inside the shock wave at M = 1.3 in the case of the pseudo-CO 2 gas; Figure 2 is for μ b / μ = 10 , Figure 3 for μ b / μ = 100 , and Figure 4 for μ b / μ = 1000 . Note that T + / T = 1.141 and ρ + / ρ = ( p + / p ) ( T / T + ) = 1.566 at M = 1.3 . In each of Figure 2, Figure 3 and Figure 4, panel (a) shows the profiles of ρ ˇ , v ˇ , and T ˇ , and panel (b) those of T ˇ , T ˇ tr , and T ˇ int ; in Figure 4, panels (c) and (d) are, respectively, the magnified figures of panels (a) and (b). The solid line indicates the result based on the two-temperature NS equations (18), and the dot-dashed line that based on the ordinary NS equation (14). The respective macroscopic quantities are shown in different colors, which can be recognized from the colors of the labels on the vertical axis. To be more specific, the red, green, and blue lines indicate, respectively, ρ ˇ , v ˇ , and T ˇ in panels (a) and (c) and T ˇ tr , T ˇ int , and T ˇ in panels (b) and (d). In Figure 3 and Figure 4, ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int obtained by using the ES model (2) [41] is also shown by the black dashed line, and in Figure 4c,d, the result based on the ordinary NS equations is not included. In addition, in Figure 4, the result for θ = 0 (or μ b / μ = ), obtained by the two-temperature NS equations, is shown by dotted lines. It should be noted that the Rankine–Hugoniot relations are different from (23) when θ = 0 (see Appendix B.2 in [41]).
Figure 5, Figure 6 and Figure 7 show the profiles at M = 3 for μ b / μ = 10 (Figure 5), 100 (Figure 6), and 1000 (Figure 7). As for the profiles at M = 5 , only the result for μ b / μ = 1000 is shown in Figure 8. The explanation of Figure 5, Figure 6 and Figure 7 and that of Figure 8 are basically the same as the explanation of Figure 2, Figure 3 and Figure 4 and that of Figure 4, respectively, so that they are omitted. It is noted that T + / T = 2.078 and ρ + / ρ = 4.932 at M = 3 , and T + / T = 3.723 and ρ + / ρ = 7.819 at M = 5 .
Figure 2, Figure 3 and Figure 4 correspond to the case where the shock wave is rather weak ( M = 1.3 ). Even for μ b / μ = 10 (Figure 2), the profiles obtained by the two-temperature NS equations deviates significantly from those by the ordinary NS equations. For μ b / μ = 100 , the result based on the two-temperature NS equations, which agrees perfectly with that based on the ES model (2), shows separation of the profiles into the thin front layer with a steep change, which may be called a subshock, and the thick rear layer with a slow relaxation. This separation is not observed in the result based on the ordinary NS equations. When μ b / μ increases up to 1000 (Figure 4), the separation becomes more prominent in the result based on the two-temperature NS equations as well as that based on the ES model (2), and the thickness of the shock wave is more than 5000 mean free paths. For a monatomic gas, it is known that the Navier–Stokes equations can describe the structure of weak shock waves relatively well. In contrast, for a polyatomic gas with large bulk viscosity, the ordinary NS equations fail to describe the structure correctly even for weak shock waves. Contrarily, the two-temperature NS equations correctly describe the shock-wave structure for such a gas. As seen from Figure 4c,d, the subshock for large μ b / μ coincides with the shock wave with θ = 0 (or μ b / μ = ).
Figure 5, Figure 6 and Figure 7 are for a higher upstream Mach number M = 3 and correspond to Figure 2, Figure 3 and Figure 4, respectively. As in Figure 2, Figure 3 and Figure 4, the profiles based on the two-temperature NS equations show the double-layer structure more eminently as μ b / μ increases. However, the thickness of the rear layer for μ b / μ = 1000 , which is about 1000 mean free paths, is much smaller than that at M = 1.3 (cf. Figure 4). As seen from the magnified figures, Figure 7c,d, the two-temperature NS equations give a thinner and sharper subshock than the ES model (2) though the profiles based on both equations agree very well in the rear layer with slow relaxation. As Figure 8 shows, the deviation between the profiles based on the two-temperature NS equations and those based on the ES model in the subshock becomes larger for a relatively strong shock wave at M = 5 , but their agreement in the thick rear layer is still very good. Figure 7c,d and Figure 8c,d show that the profiles of the subshock for large μ b / μ coincide with those of the shock wave with θ = 0 (or μ b / μ = ).
The double-layer structure of the shock profiles for CO 2 gas had long been known and was revisited using the macroscopic equations based on extended thermodynamics in [31,52,53], which motivated subsequent studies using kinetic theory [36,41,58]. The profiles with the double-layer structure was classified as type C in [31,52,53] (see also [41,58]). It should be mentioned that a similar double-layer structure is also observed for a shock wave in a ternary mixture of Euler gases with disparate masses [63].
We summarize the conclusion for the pseudo-CO 2 gas with large μ b / μ obtained in the present computation as follows:
(i)
The two-temperature NS equations are able to reproduce correct shock profiles for weak shock waves ( M = 1.3 ).
(ii)
For moderate and strong shock waves ( M = 3 and 5), although the profiles obtained by the two-temperature NS equations deviate from those based on the ES model in the thin front layer (subshock), they agree very well in the thick rear layer. The two-temperature NS equations tend to give a thinner and sharper subshock.
(iii)
The profiles of the subshock coincide with those of the shock wave with μ b / μ = .
(iv)
The ordinary NS equations (with bulk viscosity) cannot be used to correctly describe the shock-wave structure even for weak shock waves ( M = 1.3 ).
The properties (i)–(iii) have already been observed in cases where the specific heats are constant [36]. Therefore, focusing our attention on the property (iv), we investigate the shock-wave structure for other gases on the basis of the two-temperature and ordinary NS equations.

6.3.2. SF 6 Gas

For SF 6 gas, we set T = 300 K and use the values of c m ( m = 0 , , 4 ), c v ( T ) / R , γ , and δ in Table 1 and those of s, Pr, and μ b / μ in Table 2. Then, ν and θ are obtained as ν = 0.2412 and θ = 2.224 × 10 3 (note that these values are slightly different from those in Table 2). The pressure p or the density ρ at upstream infinity is not specified here. Once it is given, l , which appears in the figures as X 1 / l , is determined using (12) (with l 0 = l , T 0 = T , and ρ 0 = ρ ) and (29) [with μ ( 300 K ) in Table 2].
Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13 show the profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int inside the shock wave at M = 1.05 , 1.234 , 1.3 , 3, and 5; in Figure 13, panels (c) and (d) are, respectively, the magnified figures of panels (a) and (b) (without the result based on the ordinary NS equations). The types and colors of the lines in the figures are the same as in Figure 2, so that the explanation is omitted. It is remarked that ( T + / T , ρ + / ρ ) = ( 1.009 , 1.098 ) ( M = 1.05 ); ( 1.039 , 1.492 ) ( M = 1.234 ); ( 1.049 , 1.643 ) ( M = 1.3 ); ( 1.368 , 6.875 ) ( M = 3 ); and ( 1.938 , 13.57 ) ( M = 5 ).
When M = 1.05 (Figure 9), which corresponds to a very weak shock, the profiles are smooth, and the profiles based on the ordinary NS equations (14) do not differ much from those based on the two-temperature NS equations (18). These smooth profiles are classified as type A in [31,52,53] (see also [41,58]). However, when M = 1.234 (Figure 10), the profiles start suddenly except for T ˇ int and are quite different from those in Figure 9 though M has been increased slightly. These profiles with a sudden start are classified as type B in [31,52,53] (see also [41,58]). Although the shock wave is still weak, the ordinary NS equations cannot describe the sudden start of the profiles. The value M = 1.234 corresponds to the sonic point M ˜ = 1 of a different upstream Mach number M ˜ = v / 5 R T / 3 , that is, the upstream Mach number when the gas is regarded as a monatomic gas, or when the internal modes are frozen (see [58]). At M = 1.3 (Figure 11), one can observe the appearance of a thin front layer with a steep change (subshock). As in the case of pseudo-CO 2 gas (see Figure 4), the ordinary NS equations cannot describe the profiles with a subshock. As M is increased to 3 (Figure 12) and further to 5 (Figure 13), the amplitude of the subshock increases, and the change of physical quantities becomes steeper there. We recall that in the case of the pseudo-CO 2 gas (Figure 7 and Figure 8), the subshock based on the two-temperature NS equations (18) is much thinner than that based on the original ES model (2) for moderately strong shock waves. Therefore, we cannot expect that the profiles inside the subshock in Figure 13c,d are accurate. Nevertheless, the jump over the subshock and the profiles in the thick rear layer are expected to be correct in Figure 13. With the exception of Figure 9 corresponding to the profiles of type A, the ordinary NS equations (with the bulk viscosity) provide completely incorrect profiles.

6.3.3. CH 4 Gas

Finally, we present some results for CH 4 gas. We also set T = 300 K here and use the values of c m ( m = 0 , , 4 ), c v ( T ) / R , γ , and δ in Table 1 and those of s, Pr, and μ b / μ in Table 2. Then, ν and θ are obtained as ν = 0.3112 and θ = 1.997 × 10 3 (note that these values are slightly different from those in Table 2). The pressure p or the density ρ at upstream infinity is not specified here (see the first paragraph in Section 6.3.2). Figure 14, Figure 15, Figure 16, Figure 17 and Figure 18 show the profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int inside the shock wave at M = 1.05 , 1.131 , 1.3 , 3, and 5; in Figure 18, panels (c) and (d) are, respectively, the magnified figures of panels (a) and (b) (without the result based on the ordinary NS equations). The different types and colors of the lines in the figures are the same as in Figure 2, thus the explanation is omitted. We remark that ( T + / T , ρ + / ρ ) = ( 1.026 , 1.089 ) ( M = 1.05 ); ( 1.067 , 1.240 ) ( M = 1.131 ); ( 1.147 , 1.563 ) ( M = 1.3 ); ( 2.060 , 5.042 ) ( M = 3 ); and ( 3.408 , 8.747 ) ( M = 5 ).
The value M = 1.131 in Figure 15 corresponds to M ˜ = 1 , as in Figure 10. Although the local values of the profile are different from those for SF 6 gas shown in Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13, the qualitative behavior at respective values of M is similar to SF 6 gas, and we therefore do not repeat similar explanations here.

7. Concluding Remarks

In our recent paper [36], we derived two-temperature Navier–Stokes equations for a polyatomic gas with slow relaxation of the internal modes from the ES model [21,22]. In this model, slow relaxation of the internal modes is equivalent to large values of the ratio of bulk viscosity to viscosity μ b / μ . The resulting Navier–Stokes equations were then successfully applied to the numerical analysis of the shock-structure problem. All these results pertain to a gas with constant specific heats (calorically perfect gas).
In the present paper, we have extended the two-temperature Navier–Stokes equations to a polyatomic gas whose specific heats are not constant but depend on temperature (thermally perfect gas) using the ES model for such a gas proposed in [41]. Then, using the extended two-temperature Navier–Stokes equations, we numerically analyzed the structure of a plane shock wave for certain gases bearing large values of μ b / μ , i.e., CO 2 , SF 6 , and CH 4 . The same problem has also been analyzed by using the ordinary Navier–Stokes equations with a single temperature and bulk viscosity derived from the same extended ES model.
We first confirmed for CO 2 gas (more precisely, pseudo-CO 2 gas) that, as in the case of a calorically perfect gas [36], the two-temperature Navier–Stokes equations recover the shock profiles obtained by the numerical analysis of the ES model very well for weak shock waves and sufficiently well even for moderately strong shock waves. Subsequently, we showed that the shock profiles obtained by the ordinary Navier–Stokes equations differ significantly from those based on the two-temperature Navier–Stokes equations even for weak shock waves. Notably, the double-layer structure with a thin front layer (subshock) and a thick rear layer, as obtained by the two-temperature Navier–Stokes equations, cannot be described by the ordinary Navier–Stokes equations even qualitatively. This drawback of the ordinary Navier–Stokes equations was also demonstrated for SF 6 and CH 4 gases.
Therefore, we may conclude that the ordinary Navier–Stokes equations with bulk viscosity are inappropriate for polyatomic gases with large values of μ b / μ , even at the so-called Navier–Stokes level. For such gases, one should use more appropriate macroscopic equations, such as the two-temperature Navier–Stokes equations derived in this paper.
The present two-temperature Navier–Stokes equations are derived from the ES model (2) with a single energy variable E for all the internal modes. Therefore, if, for instance, one of the vibrational modes relaxes slowly, we suppose that the internal modes relax slowly and that the two-temperature Navier–Stokes equations are applicable. In this sense, the latter equations are rather primitive. Recently, more sophisticated ES models using two separate energy variables for rotational and vibrational modes have been proposed [26,27]. With the help of these new models and of an appropriate parameter setting, one may derive multi-temperature Navier–Stokes equations with the translational, rotational, and vibrational temperatures. The derivation of such equations would be an interesting and important subject as the next step.

Author Contributions

Conceptualization, S.K. and K.A.; methodology, S.K. and K.A.; formal analysis, S.K. and K.A; numerical analysis, S.K.; writing—original draft preparation, K.A.; writing—review and editing, S.K.; All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Anderson, J.D., Jr. Hypersonic and High Temperature Gas Dynamics; McGraw-Hill: New York, NY, USA, 1989. [Google Scholar]
  2. Park, C. Nonequilibrium Hypersonic Aerothermodynamics; John Wiley & Sons: New York, NY, USA, 1990. [Google Scholar]
  3. Nagnibeda, E.; Kustova, E. Non-Equilibrium Reacting Gas Flows: Kinetic Theory of Transport and Relaxation Processes; Springer: Berlin, Germany, 2009. [Google Scholar]
  4. Boyd, I.D.; Schwartzentruber, T.E. Nonequilibrium Gas Dynamics and Molecular Simulation; Cambridge University Press: Cambridge, UK, 2017. [Google Scholar]
  5. Wang Chang, C.S.; Uhlenbeck, G.E. Transport Phenomena in Polyatomic Gases; Engineering Research Institute Report CM-681; University of Michigan: Ann Arbor, MI, USA, 1951. [Google Scholar]
  6. Ferziger, J.H.; Kaper, H.G. Mathematical Theory of Transport Processes in Gases; North Holland: Amsterdam, The Netherlands, 1972. [Google Scholar]
  7. McCourt, F.R.W.; Beenakker, J.J.M.; Köhler, W.E.; Kuščer, I. Nonequilibrium Phenomena in Polyatomic Gases, Volume 1: Dilute Gases; Clarendon: Oxford, UK, 1990. [Google Scholar]
  8. Giovangigli, V. Multicomponent Flow Modeling; Birkhäuser: Boston, MA, USA, 1999. [Google Scholar]
  9. Borgnakke, C.; Larsen, P.S. Statistical collision model for Monte Carlo simulation of polyatomic gas mixture. J. Comp. Phys. 1975, 18, 405–420. [Google Scholar] [CrossRef]
  10. Bourgat, J.-F.; Desvillettes, L.; Le Tallec, P.; Perthame, B. Microreversible collisions for polyatomic gases and Boltzmann’s theorem. Eur. J. Mech. B Fluids 1994, 13, 237–254. [Google Scholar]
  11. Desvillettes, L.; Monaco, R.; Salvarani, F. A kinetic model allowing to obtain the energy law of polytropic gases in the presence of chemical reactions. Eur. J. Mech. B Fluids 2005, 24, 219–236. [Google Scholar] [CrossRef]
  12. Borsoni, T.; Bisi, M.; Groppi, M. A general framework for the kinetic modelling of polyatomic gases. Commun. Math. Phys. 2022, 393, 215–266. [Google Scholar] [CrossRef]
  13. Pavić-Čolić, M.; Simić, S. Kinetic description of polyatomic gases with temperature-dependent specific heats. Phys. Rev. Fluids 2022, 7, 083401. [Google Scholar] [CrossRef]
  14. Boudin, L.; Rossi, A.; Salvarani, F. A kinetic model of polyatomic gas with resonant collisions. Ricerche Mat. 2022. [Google Scholar] [CrossRef]
  15. Gamba, I.M.; Pavić-Čolić, M. On the Cauchy problem for Boltzmann equation modelling a polyatomic gas. arXiv Prepr. 2020, arXiv:2005.01017. [Google Scholar]
  16. Brull, S.; Shahine, M.; Thieullen, P. Compactness property of the linearized Boltzmann operator for a diatomic single gas model. Netw. Heterog. Media 2022, 17, 847–861. [Google Scholar] [CrossRef]
  17. Bernhoff, N. Linearized Boltzmann collision operator: II. Polyatomic molecules modeled by a continuous internal energy variable. arXiv Prepr. 2022, arXiv:2201.01377. [Google Scholar]
  18. Morse, T.F. Kinetic model for gases with internal degrees of freedom. Phys. Fluids 1964, 7, 159–169. [Google Scholar] [CrossRef]
  19. Holway, L.H., Jr. New statistical models for kinetic theory: Methods of construction. Phys. Fluids 1966, 9, 1658–1673. [Google Scholar] [CrossRef]
  20. Rykov, V.A. A model kinetic equation for a gas with rotational degrees of freedom. Fluid Dyn. 1975, 10, 959–966. [Google Scholar] [CrossRef]
  21. Andries, P.; Le Tallec, P.; Perlat, J.-P.; Perthame, B. The Gaussian-BGK model of Boltzmann equation with small Prandtl number. Eur. J. Mech. B Fluids 2000, 19, 813–830. [Google Scholar] [CrossRef]
  22. Brull, S.; Schneider, J. On the ellipsoidal statistical model for polyatomic gases. Contin. Mech. Thermodyn. 2009, 20, 489–508. [Google Scholar] [CrossRef] [Green Version]
  23. Rahimi, B.; Struchtrup, H. Capturing non-equilibrium phenomena in rarefied polyatomic gases: A high-order macroscopic model. Phys. Fluids 2014, 26, 052001. [Google Scholar] [CrossRef]
  24. Bisi, M.; Spiga, G. On kinetic models for polyatomic gases and their hydrodynamic limits. Ric. Mat. 2017, 66, 113–124. [Google Scholar] [CrossRef]
  25. Baranger, C.; Dauvois, Y.; Marois, G.; Mathé, J.; Mathiaud, J.; Mieussens, L. A BGK model for high temperature rarefied gas flows. Eur. J. Mech. B Fluids 2020, 80, 1–12. [Google Scholar] [CrossRef] [Green Version]
  26. Dauvois, Y.; Mathiaud, J.; Mieussens, L. An ES-BGK model for polyatomic gases in rotational and vibrational nonequilibrium. Eur. J. Mech. B Fluids 2021, 88, 1–16. [Google Scholar] [CrossRef]
  27. Mathiaud, J.; Mieussens, L.; Pfeiffer, M. An ES-BGK model for diatomic gases with correct relaxation rates for internal energies. Eur. J. Mech. B Fluids 2022, 96, 65–77. [Google Scholar] [CrossRef]
  28. Myong, R.S. Coupled nonlinear constitutive models for rarefied and microscale gas flows: Subtle interplay of kinematics and dissipation effects. Contin. Mech. Thermodyn. 2009, 21, 389–399. [Google Scholar] [CrossRef]
  29. Arima, T.; Taniguchi, S.; Ruggeri, T.; Sugiyama, M. Extended thermodynamics of real gases with dynamic pressure: An extension of Meixner’s theory. Phys. Lett. A 2012, 376, 2799–2803. [Google Scholar] [CrossRef]
  30. Pavić, M.; Ruggeri, T.; Simić, S. Maximum entropy principle for rarefied polyatomic gases. Phys. A 2013, 392, 1302–1317. [Google Scholar] [CrossRef] [Green Version]
  31. Taniguchi, S.; Arima, T.; Ruggeri, T.; Sugiyama, M. Overshoot of the non-equilibrium temperature in the shock wave structure of a rarefied polyatomic gas subject to the dynamic pressure. Int. J. Non-Linear Mech. 2016, 79, 66–75. [Google Scholar] [CrossRef]
  32. Arima, T.; Ruggeri, T.; Sugiyama, M. Rational extended thermodynamics of a rarefied polyatomic gas with molecular relaxation processes. Phys. Rev. E 2017, 96, 042143. [Google Scholar] [CrossRef] [Green Version]
  33. Pavić-Čolić, M.; Madjarević, D.; Simić, S. Polyatomic gases with dynamic pressure: Kinetic non-linear closure and the shock structure. Int. J. Non-Linear Mech. 2017, 92, 160–175. [Google Scholar] [CrossRef]
  34. Bisi, M.; Ruggeri, T.; Spiga, G. Dynamical pressure in a polyatomic gas: Interplay between kinetic theory and extended thermodynamics. Kin. Rel. Models 2018, 11, 71–95. [Google Scholar] [CrossRef] [Green Version]
  35. Ruggeri, T.; Sugiyama, M. Classical and Relativistic Rational Extended Thermodynamics of Gases; Springer: Cham, Switzerland, 2021. [Google Scholar]
  36. Aoki, K.; Bisi, M.; Groppi, M.; Kosuge, S. Two-temperature Navier–Stokes equations for a polyatomic gas derived from kinetic theory. Phys. Rev. E 2020, 102, 023104. [Google Scholar] [CrossRef]
  37. Chapman, S.; Cowling, T.G. The Mathematical Theory of Non-Uniform Gases, 3rd ed.; Cambridge University Press: Cambridge, UK, 1991. [Google Scholar]
  38. Grad, H. Principles of the kinetic theory of gases. In Handbuch der Physik; Band XII; Flügge, S., Ed.; Springer: Berlin, Germany, 1958; pp. 205–294. [Google Scholar]
  39. Cercignani, C. The Boltzmann Equation and Its Applications; Springer: Berlin, Germany, 1988. [Google Scholar]
  40. Sone, Y. Molecular Gas Dynamics: Theory, Techniques, and Applications; Birkhäuser: Boston, MA, USA, 2007; Supplementary Notes and Errata; Available online: http://hdl.handle.net/2433/66098 (accessed on 21 July 2019).
  41. Kosuge, S.; Kuo, H.-W.; Aoki, K. A kinetic model for a polyatomic gas with temperature-dependent specific heats and its application to shock-wave structure. J. Stat. Phys. 2019, 177, 209–251. [Google Scholar] [CrossRef]
  42. Park, C.; Yoon, S. Calculation of real-gas effects on blunt-body trim angles. AIAA J. 1992, 30, 999–1007. [Google Scholar] [CrossRef]
  43. Park, C.; Lee, S.-H. Validation of multitemperature nozzle flow code. J. Thermophys. Heat Trans. 1995, 9, 9–16. [Google Scholar] [CrossRef]
  44. Bruno, D.; Giovangigli, V. Relaxation of internal temperature and volume viscosity. Phys. Fluids 2011, 23, 093104, Erratum in Phys. Fluids 2013, 25, 039902. [Google Scholar] [CrossRef] [Green Version]
  45. Bruno, D.; Giovangigli, V. Internal energy relaxation processes and bulk viscosities in fluids. Fluids 2022, 7, 356. [Google Scholar] [CrossRef]
  46. Kosuge, S.; Aoki, K.; Bisi, M.; Groppi, M.; Martalò, G. Boundary conditions for two-temperature Navier–Stokes equations for a polyatomic gas. Phys. Rev. Fluids 2021, 6, 083401. [Google Scholar] [CrossRef]
  47. Bird, G.A. Molecular Gas Dynamics and the Direct Simulation of Gas Flows; Oxford University Press: Oxford, UK, 1994. [Google Scholar]
  48. Smiley, E.F.; Winkler, E.H.; Slawsky, Z.I. Measurement of the vibrational relaxation effect in CO2 by means of shock tube interferograms. J. Chem. Phys. 1952, 20, 923–924. [Google Scholar] [CrossRef]
  49. Smiley, E.F.; Winkler, E.H. Shock-tube measurements of vibrational relaxation. J. Chem. Phys. 1954, 22, 2018–2022. [Google Scholar] [CrossRef]
  50. Griffith, W.D.; Kenny, A. On fully-dispersed shock waves in carbon dioxide. J. Fluid Mech. 1957, 3, 286–289. [Google Scholar] [CrossRef]
  51. Johannesen, N.H.; Zienkiewicz, H.K.; Blythe, P.A.; Gerrard, J.H. Experimental and theoretical analysis of vibrational relaxation regions in carbon dioxide. J. Fluid Mech. 1962, 13, 213–225. [Google Scholar] [CrossRef]
  52. Taniguchi, S.; Arima, T.; Ruggeri, T.; Sugiyama, M. Effect of the dynamic pressure on the shock wave structure in a rarefied polyatomic gas. Phys. Fluids 2014, 26, 016103. [Google Scholar] [CrossRef]
  53. Taniguchi, S.; Arima, T.; Ruggeri, T.; Sugiyama, M. Thermodynamic theory of the shock wave structure in a rarefied polyatomic gas: Beyond the Bethe-Teller theory. Phys. Rev. E 2014, 89, 013025. [Google Scholar] [CrossRef]
  54. Alekseev, I.V.; Kosareva, A.A.; Kustova, E.V.; Nagnibeda, E.A. Various continuum approaches for studying shock wave structure in carbon dioxide. In Proceedings of the The Eighth Polyakhov’s Reading: Proceedings of the International Scientific Conference on Mechanics, Saint Petersburg, Russia, 29 January–2 February 2018; Kustova, E., Leonov, G., Morosov, N., Yushkov, M., Mekhonoshina, M., Eds.; AIP: Melville, NY, USA, 2018; p. 060001. [Google Scholar]
  55. Alekseev, I.V.; Kosareva, A.A.; Kustova, E.V.; Nagnibeda, E.A. Shock waves in carbon dioxide: Simulations using different kinetic-theory models. In Proceedings of the 31st International Symposium on Rarefied Gas Dynamics, Glasgow, UK, 23–27 July 2018; Zhang, Y., Emerson, D.R., Lockerby, D., Wu, L., Eds.; AIP: Melville, NY, USA, 2019; p. 060005. [Google Scholar]
  56. Kustova, E.; Mekhonoshina, M.; Kosareva, A. Relaxation processes in carbon dioxide. Phys. Fluids 2019, 31, 046104. [Google Scholar] [CrossRef]
  57. Alekseev, I.; Kustova, E. Extended continuum models for shock waves in CO2. Phys. Fluids 2021, 33, 096101. [Google Scholar] [CrossRef]
  58. Kosuge, S.; Aoki, K. Shock-wave structure for a polyatomic gas with large bulk viscosity. Phys. Rev. Fluids 2018, 3, 023401. [Google Scholar] [CrossRef]
  59. Japan Society of Thermophysical Properties (Ed.) Thermophysical Properties Handbook; Yokendo: Tokyo, Japan, 1990. (In Japanese) [Google Scholar]
  60. Uribe, F.J.; Mason, E.A.; Kestin, J. Thermal conductivity of nine polyatomic gases at low density. J. Phys. Chem. Ref. Data 1990, 19, 1123–1136. [Google Scholar] [CrossRef]
  61. Cramer, M.S. Numerical estimates for the bulk viscosity of ideal gases. Phys. Fluids 2012, 24, 066102. [Google Scholar] [CrossRef] [Green Version]
  62. Boushehri, A.; Bzowski, J.; Kestin, J.; Mason, E.A. Equilibrium and transport properties of eleven polyatomic gases at low density. J. Phys. Chem. Ref. Data 1987, 16, 445–466, Erratum in J. Phys. Chem. Ref. Data 1988, 17, 255. [Google Scholar] [CrossRef]
  63. Madjarević, D.; Pavić-Čolić, M.; Simić, S. Shock structure and relaxation in the multi-component mixture of Euler fluids. Symmetry 2021, 13, 955. [Google Scholar] [CrossRef]
Figure 1. c v ( T ) / R vs. T. (a) SF 6 , (b) CH 4 , and (c) CO 2 . The solid line indicates the fitting Formula (27) with the data in Table 1, the red circles represent the data from [60], and the blue circles denote those from [59].
Figure 1. c v ( T ) / R vs. T. (a) SF 6 , (b) CH 4 , and (c) CO 2 . The solid line indicates the fitting Formula (27) with the data in Table 1, the red circles represent the data from [60], and the blue circles denote those from [59].
Fluids 08 00005 g001
Figure 2. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of the pseudo-CO 2 gas at M = 1.3 for μ b / μ = 10 . (a) ρ ˇ , v ˇ , and T ˇ ; (b) T ˇ tr , T ˇ int , and T ˇ . The solid line indicates the result based on the two-temperature NS equations (18), and the dot-dashed line that based on the ordinary NS equations (14). The red, green, and blue lines indicate, respectively, ρ ˇ , v ˇ , and T ˇ in panel (a) and T ˇ tr , T ˇ int , and T ˇ in panel (b).
Figure 2. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of the pseudo-CO 2 gas at M = 1.3 for μ b / μ = 10 . (a) ρ ˇ , v ˇ , and T ˇ ; (b) T ˇ tr , T ˇ int , and T ˇ . The solid line indicates the result based on the two-temperature NS equations (18), and the dot-dashed line that based on the ordinary NS equations (14). The red, green, and blue lines indicate, respectively, ρ ˇ , v ˇ , and T ˇ in panel (a) and T ˇ tr , T ˇ int , and T ˇ in panel (b).
Fluids 08 00005 g002
Figure 3. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of the pseudo-CO 2 gas at M = 1.3 for μ b / μ = 100 . (a) ρ ˇ , v ˇ , and T ˇ ; (b) T ˇ tr , T ˇ int , and T ˇ . The solid line indicates the result based on the two-temperature NS equations (18), and the dot-dashed line that based on the ordinary NS equations (14). The red, green, and blue lines indicate, respectively, ρ ˇ , v ˇ , and T ˇ in panel (a) and T ˇ tr , T ˇ int , and T ˇ in panel (b). The profiles based on the ES model (2) are shown by the black dashed line for all quantities.
Figure 3. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of the pseudo-CO 2 gas at M = 1.3 for μ b / μ = 100 . (a) ρ ˇ , v ˇ , and T ˇ ; (b) T ˇ tr , T ˇ int , and T ˇ . The solid line indicates the result based on the two-temperature NS equations (18), and the dot-dashed line that based on the ordinary NS equations (14). The red, green, and blue lines indicate, respectively, ρ ˇ , v ˇ , and T ˇ in panel (a) and T ˇ tr , T ˇ int , and T ˇ in panel (b). The profiles based on the ES model (2) are shown by the black dashed line for all quantities.
Fluids 08 00005 g003
Figure 4. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of the pseudo-CO 2 gas at M = 1.3 for μ b / μ = 1000 . (a,c) ρ ˇ , v ˇ , and T ˇ ; (b,d) T ˇ tr , T ˇ int , and T ˇ . Panels (c,d) are, respectively, the magnified figures of panels (a,b) [except that the dot-dashed line is not included in panels (c,d)]. The solid line indicates the result based on the two-temperature NS equations (18), the dotted line that for θ = 0 (or μ b / μ = ) based on the same equations, and the dot-dashed line that based on the ordinary NS equations (14). The red, green, and blue lines indicate, respectively, ρ ˇ , v ˇ , and T ˇ in panels (a,c) and T ˇ tr , T ˇ int , and T ˇ in panels (b,d). The profiles based on the ES model (2) are shown by the black dashed line for all quantities.
Figure 4. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of the pseudo-CO 2 gas at M = 1.3 for μ b / μ = 1000 . (a,c) ρ ˇ , v ˇ , and T ˇ ; (b,d) T ˇ tr , T ˇ int , and T ˇ . Panels (c,d) are, respectively, the magnified figures of panels (a,b) [except that the dot-dashed line is not included in panels (c,d)]. The solid line indicates the result based on the two-temperature NS equations (18), the dotted line that for θ = 0 (or μ b / μ = ) based on the same equations, and the dot-dashed line that based on the ordinary NS equations (14). The red, green, and blue lines indicate, respectively, ρ ˇ , v ˇ , and T ˇ in panels (a,c) and T ˇ tr , T ˇ int , and T ˇ in panels (b,d). The profiles based on the ES model (2) are shown by the black dashed line for all quantities.
Fluids 08 00005 g004
Figure 5. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of the pseudo-CO 2 gas at M = 3 for μ b / μ = 10 . (a) ρ ˇ , v ˇ , and T ˇ ; (b) T ˇ tr , T ˇ int , and T ˇ . See the caption of Figure 2.
Figure 5. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of the pseudo-CO 2 gas at M = 3 for μ b / μ = 10 . (a) ρ ˇ , v ˇ , and T ˇ ; (b) T ˇ tr , T ˇ int , and T ˇ . See the caption of Figure 2.
Fluids 08 00005 g005
Figure 6. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of the pseudo-CO 2 gas at M = 3 for μ b / μ = 100 . (a) ρ ˇ , v ˇ , and T ˇ ; (b) T ˇ tr , T ˇ int , and T ˇ . See the caption of Figure 3.
Figure 6. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of the pseudo-CO 2 gas at M = 3 for μ b / μ = 100 . (a) ρ ˇ , v ˇ , and T ˇ ; (b) T ˇ tr , T ˇ int , and T ˇ . See the caption of Figure 3.
Fluids 08 00005 g006
Figure 7. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of the pseudo-CO 2 gas at M = 3 for μ b / μ = 1000 . (a,c) ρ ˇ , v ˇ , and T ˇ ; (b,d) T ˇ tr , T ˇ int , and T ˇ . See the caption of Figure 4.
Figure 7. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of the pseudo-CO 2 gas at M = 3 for μ b / μ = 1000 . (a,c) ρ ˇ , v ˇ , and T ˇ ; (b,d) T ˇ tr , T ˇ int , and T ˇ . See the caption of Figure 4.
Fluids 08 00005 g007
Figure 8. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of the pseudo-CO 2 gas at M = 5 for μ b / μ = 1000 . (a,c) ρ ˇ , v ˇ , and T ˇ ; (b,d) T ˇ tr , T ˇ int , and T ˇ . See the caption of Figure 4.
Figure 8. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of the pseudo-CO 2 gas at M = 5 for μ b / μ = 1000 . (a,c) ρ ˇ , v ˇ , and T ˇ ; (b,d) T ˇ tr , T ˇ int , and T ˇ . See the caption of Figure 4.
Fluids 08 00005 g008
Figure 9. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of SF 6 gas at M = 1.05 . (a) ρ ˇ , v ˇ , and T ˇ ; (b) T ˇ tr , T ˇ int , and T ˇ . See the caption of Figure 2.
Figure 9. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of SF 6 gas at M = 1.05 . (a) ρ ˇ , v ˇ , and T ˇ ; (b) T ˇ tr , T ˇ int , and T ˇ . See the caption of Figure 2.
Fluids 08 00005 g009
Figure 10. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of SF 6 gas at M = 1.234 . (a) ρ ˇ , v ˇ , and T ˇ ; (b) T ˇ tr , T ˇ int , and T ˇ . See the caption of Figure 2.
Figure 10. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of SF 6 gas at M = 1.234 . (a) ρ ˇ , v ˇ , and T ˇ ; (b) T ˇ tr , T ˇ int , and T ˇ . See the caption of Figure 2.
Fluids 08 00005 g010
Figure 11. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of SF 6 gas at M = 1.3 . (a) ρ ˇ , v ˇ , and T ˇ ; (b) T ˇ tr , T ˇ int , and T ˇ . See the caption of Figure 2.
Figure 11. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of SF 6 gas at M = 1.3 . (a) ρ ˇ , v ˇ , and T ˇ ; (b) T ˇ tr , T ˇ int , and T ˇ . See the caption of Figure 2.
Fluids 08 00005 g011
Figure 12. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of SF 6 gas at M = 3 . (a) ρ ˇ , v ˇ , and T ˇ ; (b) T ˇ tr , T ˇ int , and T ˇ . See the caption of Figure 2.
Figure 12. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of SF 6 gas at M = 3 . (a) ρ ˇ , v ˇ , and T ˇ ; (b) T ˇ tr , T ˇ int , and T ˇ . See the caption of Figure 2.
Fluids 08 00005 g012
Figure 13. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of SF 6 gas at M = 5 . (a,c) ρ ˇ , v ˇ , and T ˇ ; (b,d) T ˇ tr , T ˇ int , and T ˇ . Panels (c,d) are, respectively, the magnified figures of panels (a,b). See the caption of Figure 2.
Figure 13. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of SF 6 gas at M = 5 . (a,c) ρ ˇ , v ˇ , and T ˇ ; (b,d) T ˇ tr , T ˇ int , and T ˇ . Panels (c,d) are, respectively, the magnified figures of panels (a,b). See the caption of Figure 2.
Fluids 08 00005 g013
Figure 14. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of CH 4 gas at M = 1.05 . (a) ρ ˇ , v ˇ , and T ˇ ; (b) T ˇ tr , T ˇ int , and T ˇ . See the caption of Figure 2.
Figure 14. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of CH 4 gas at M = 1.05 . (a) ρ ˇ , v ˇ , and T ˇ ; (b) T ˇ tr , T ˇ int , and T ˇ . See the caption of Figure 2.
Fluids 08 00005 g014
Figure 15. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of CH 4 gas at M = 1.131 . (a) ρ ˇ , v ˇ , and T ˇ ; (b) T ˇ tr , T ˇ int , and T ˇ . See the caption of Figure 2.
Figure 15. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of CH 4 gas at M = 1.131 . (a) ρ ˇ , v ˇ , and T ˇ ; (b) T ˇ tr , T ˇ int , and T ˇ . See the caption of Figure 2.
Fluids 08 00005 g015
Figure 16. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of CH 4 gas at M = 1.3 . (a) ρ ˇ , v ˇ , and T ˇ ; (b) T ˇ tr , T ˇ int , and T ˇ . See the caption of Figure 2.
Figure 16. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of CH 4 gas at M = 1.3 . (a) ρ ˇ , v ˇ , and T ˇ ; (b) T ˇ tr , T ˇ int , and T ˇ . See the caption of Figure 2.
Fluids 08 00005 g016
Figure 17. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of CH 4 gas at M = 3 . (a) ρ ˇ , v ˇ , and T ˇ ; (b) T ˇ tr , T ˇ int , and T ˇ . See the caption of Figure 2.
Figure 17. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of CH 4 gas at M = 3 . (a) ρ ˇ , v ˇ , and T ˇ ; (b) T ˇ tr , T ˇ int , and T ˇ . See the caption of Figure 2.
Fluids 08 00005 g017
Figure 18. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of CH 4 gas at M = 5 . (a,c) ρ ˇ , v ˇ , and T ˇ ; (b,d) T ˇ tr , T ˇ int , and T ˇ . Panels (c,d) are, respectively, the magnified figures of panels (a,b). See the caption of Figure 2.
Figure 18. Profiles of ρ ˇ , v ˇ , T ˇ , T ˇ tr , and T ˇ int in the case of CH 4 gas at M = 5 . (a,c) ρ ˇ , v ˇ , and T ˇ ; (b,d) T ˇ tr , T ˇ int , and T ˇ . Panels (c,d) are, respectively, the magnified figures of panels (a,b). See the caption of Figure 2.
Fluids 08 00005 g018
Table 1. The coefficients in the fitting formula (27) and c v ( T ) , γ , and δ for T = 300 K for SF 6 , CH 4 , and CO 2 .
Table 1. The coefficients in the fitting formula (27) and c v ( T ) , γ , and δ for T = 300 K for SF 6 , CH 4 , and CO 2 .
c 0 c 1 c 2 c 3 c 4 c v ( T ) / R γ δ
SF 6 4.017 7.670 ( 2 ) a 1.130 ( 4 ) 7.745 ( 8 ) 2.014 ( 11 ) 10.751.09318.50
CH 4 3.505 7.952 ( 3 ) 3.221 ( 5 ) 2.829 ( 8 ) 8.280 ( 12 ) 3.3221.3013.643
CO 2 1.412 b 8.697 ( 3 ) b 6.575 ( 6 ) b 1.987 ( 9 ) b03.4831.2873.966
a Read as 7.670 × 10−2; b Value in [31].
Table 2. Thermophysical properties of some gases with large bulk viscosities.
Table 2. Thermophysical properties of some gases with large bulk viscosities.
μ [ 10 6 Pa · s ] s γ Pr μ b / μ θ × 10 3 ν
SF 6 15.24 a0.901.097 a0.806 a320 b 2.21 0.241
CH 4 11.43 a0.831.305 a0.7630 a240 b 1.98 0.311
CO 2 14.91 a0.831.293 a0.767 a3849 b 0.127 0.304
a Value at 300 K, 101 kPa in [59]; b Value at 300 K in [61].
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kosuge, S.; Aoki, K. Navier–Stokes Equations and Bulk Viscosity for a Polyatomic Gas with Temperature-Dependent Specific Heats. Fluids 2023, 8, 5. https://doi.org/10.3390/fluids8010005

AMA Style

Kosuge S, Aoki K. Navier–Stokes Equations and Bulk Viscosity for a Polyatomic Gas with Temperature-Dependent Specific Heats. Fluids. 2023; 8(1):5. https://doi.org/10.3390/fluids8010005

Chicago/Turabian Style

Kosuge, Shingo, and Kazuo Aoki. 2023. "Navier–Stokes Equations and Bulk Viscosity for a Polyatomic Gas with Temperature-Dependent Specific Heats" Fluids 8, no. 1: 5. https://doi.org/10.3390/fluids8010005

APA Style

Kosuge, S., & Aoki, K. (2023). Navier–Stokes Equations and Bulk Viscosity for a Polyatomic Gas with Temperature-Dependent Specific Heats. Fluids, 8(1), 5. https://doi.org/10.3390/fluids8010005

Article Metrics

Back to TopTop