Compressible Viscous Equations

Compressible viscous fluid flow is described by a set of coupled partial differential equations that represent the conservation laws of physics. These equations govern the behavior of fluids under conditions where both compressibility effects and viscous forces are significant. The complete set of equations includes the continuity equation (conservation of mass), momentum equations (Navier-Stokes), and the energy equation, along with constitutive relations.

Continuity Equation (Conservation of Mass)

The continuity equation expresses the principle of mass conservation for a compressible fluid:

Where:

• is the fluid density
• u is the velocity vector (u, u, u)
• represents the divergence operator
• t is time

This equation states that the rate of change of density within a control volume plus the net mass flux out of the volume equals zero.

Conservation of Momentum (Navier-Stokes Equations)

For compressible flows, the Navier-Stokes equations account for the variable density and temperature effects on viscosity. The momentum equation is:

Where:

• p is the pressure
• is the stress tensor
• g is the body force per unit mass (e.g., gravity)

The stress tensor for a Newtonian fluid is defined as:

Where:

• is the second coefficient of viscosity (often related to the bulk viscosity)
• is the dynamic viscosity
• I is the identity matrix
• (u)^T represents the transpose of the velocity gradient tensor

For compressible flows, the viscosity coefficients typically depend on temperature, following relationships such as Sutherland's law:

Where is the reference viscosity at temperature T, and S is Sutherland's constant (typically around 110K for air).

Conservation of Energy

The energy equation represents the first law of thermodynamics for a fluid element:

Or alternatively, in terms of total energy:

Where:

• c_v is the specific heat at constant volume
• T is the temperature
• e is the total energy per unit mass (internal + kinetic)
• k is the thermal conductivity
• is the viscous dissipation function, representing viscous heating
• S_rad represents radiation effects (if applicable)

The viscous dissipation function is given by:

Equation of State

To close the system, an equation of state relates pressure, density, and temperature. For a perfect gas:

Where R is the specific gas constant. For real gases, more complex equations of state like the Van der Waals equation or other cubic equations might be used:

Where a and b are Van der Waals constants specific to each gas.

Compressibility Effects

When the flow velocity approaches or exceeds the speed of sound, compressibility effects become significant. The Mach number (M) characterizes the importance of compressibility:

Where u is the flow velocity and c is the speed of sound, given by:

Where is the ratio of specific heats (c_p/c_v).

Boundary Conditions

Appropriate boundary conditions are essential for solving the compressible Navier-Stokes equations:

• No-slip condition: At a solid wall, u = 0 (no relative motion between fluid and wall)
• No-penetration condition: At an impermeable wall, un = 0 (normal component of velocity is zero)
• Temperature condition: Either a specified wall temperature (T = T_w) or specified wall heat flux (kT/n = q_w)
• Pressure condition: At an inlet/outlet, either pressure is specified or extrapolated based on interior values
• Far-field condition: As r , flow approaches freestream conditions

Non-dimensionalization

Non-dimensional forms of the equations highlight the important dimensionless parameters:

Where asterisk denotes dimensionless quantities, and:

• Re = UL/ (Reynolds number)
• M = U/(RT) (Freestream Mach number)
• Pr = c_p/k (Prandtl number)
• Fr = U/(gL) (Froude number)

Simplifications and Special Cases

Stokes Flow (Creeping Flow)

For very low Reynolds numbers (Re << 1), inertial terms can be neglected:

Boundary Layer Approximation (Prandtl)

For high Reynolds numbers, the Navier-Stokes equations can be simplified in thin regions near solid surfaces. For a 2D steady compressible boundary layer:

turbulence Modeling

For turbulent compressible flows, averaging techniques like Reynolds decomposition leads to the Reynolds-averaged Navier-Stokes (RANS) equations:

The Reynolds stress tensor (-u'_iu'_j) and turbulent heat flux (-h'_ju'_j) require turbulence modeling using approaches like k-, k-, or Reynolds stress models.

Numerical Solution Methods

Solving the compressible viscous equations numerically presents several challenges due to:

• Nonlinearity of the equations
• Coupling between different variables
• Shock waves and discontinuities in supersonic flows
• Multiscale nature of turbulent flows

Common numerical approaches include:

1. Finite Volume Method: Conservative formulation ideal for capturing shocks
2. Finite Element Method: Flexible for complex geometries
3. Spectral Method: High accuracy for smooth solutions
4. Discontinuous Galerkin Method: Combines advantages of finite volume and finite element methods

For compressible flows, specialized techniques are needed to handle shock waves, such as:

• Flux limiters and shock capturing schemes (e.g., TVD, MUSCL)
• Artificial viscosity methods
• Adaptive mesh refinement for resolving shock regions

Applications

The compressible viscous equations find applications in numerous engineering and scientific fields:

• Aerodynamics: Aircraft design, rocket propulsion, hypersonic vehicles
• Turbomachinery: Jet engines, gas turbines, compressors
• Internal combustion engines: Inlet and exhaust flow modeling
• Environmental flows: Atmospheric boundary layers, pollutant dispersion
• Biomedical flows: Airflow in respiratory systems
• Acoustics: Sound generation and propagation

Conclusion

The compressible viscous equations form the foundation for understanding and predicting fluid flow behavior where density variations and viscous effects are both important. These coupled nonlinear partial differential equations capture the essential physics of mass, momentum, and energy conservation in fluid systems. Despite their complexity, analytical solutions are possible only for highly simplified cases, and numerical methods are essential for solving practical problems. Advances in computing power and algorithm development continue to expand the range of problems that can be accurately modeled using these fundamental equations.