A brief introduction to bulk viscosity of fluids

Bhanuday Sharma,Rakesh Kumar
2023-03-15
Abstract:Fluid flows are typically studied by solving the Navier--Stokes equation. One of the fundamental assumptions of this equation is Stokes' hypothesis. This hypothesis assumes bulk viscosity, to be identically zero. The Stokes' hypothesis is a reasonable approximation for commonly observed fluid flows; therefore, Navier--Stokes equation gives satisfactory results in these situations. However, there are circumstances where this hypothesis does not hold good, and therefore, the classical Navier--Stokes equation becomes inapt. These situations include the absorption of sound waves, hypersonic flows, turbulent flows, and flow of Martian air etc. Reliable analytical and computational studies of these flows requires account of bulk viscosity in the governing equations. In this article, we provide a brief review of the subject of bulk viscosity. We start with a brief background of this topic. Then we discuss the underlying microscopic mechanisms that give rise to bulk viscosity effects. It was followed by a review of methods available in the literature for estimation of this parameter. Finally, a review of the studies that analyze the effects of bulk viscosity in various fluid flows is provided.
Fluid Dynamics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is about the importance and influence of fluid bulk viscosity under specific conditions. Specifically, the paper explores the problem that the Stokes assumption (i.e., considering the bulk viscosity to be zero) is no longer applicable in some special cases. These special cases include acoustic wave absorption, hypersonic flow, turbulence, and flow in the Martian atmosphere. In these cases, the influence of bulk viscosity on the flow behavior cannot be ignored. Therefore, the role of bulk viscosity needs to be considered in the control equations in order to obtain more accurate analysis results. ### Summary of the main problems in the paper: 1. **Limitations of the Stokes assumption**: The traditional Navier - Stokes equations are based on the Stokes assumption, which considers the bulk viscosity to be zero. However, under certain specific conditions, this assumption no longer holds, resulting in the inapplicability of the classical Navier - Stokes equations. 2. **Importance of bulk viscosity**: In cases such as acoustic wave absorption, hypersonic flow, turbulence, and Martian atmospheric flow, bulk viscosity has a significant impact on the flow behavior and needs to be considered in theoretical and computational research. 3. **Microscopic mechanisms of bulk viscosity**: The paper discusses the microscopic mechanisms that generate bulk viscosity, including energy exchange of internal degrees of freedom of molecules and molecular structure rearrangement. 4. **Methods for determining bulk viscosity**: The paper reviews existing theoretical, experimental, and numerical methods for estimating and measuring bulk viscosity, including theoretical formulas, acoustic wave absorption and dispersion experiments, Rayleigh - Brillouin scattering, shock wave thickness measurement, and molecular dynamics simulation. ### Formula summary: - **Stress - strain rate relationship**: \[ \sigma_{ik} = -P_{\text{thermo}} \delta_{ik} + \mu \left( \frac{\partial u_i}{\partial x_k} + \frac{\partial u_k}{\partial x_i} \right) + \lambda \frac{\partial u_j}{\partial x_j} \delta_{ik} \] where \(\sigma_{ik}\) is the Cauchy stress tensor, \(\delta_{ik}\) is the Kronecker delta, \(u_i\) is the fluid velocity, \(x_i\) is the spatial coordinate, \(P_{\text{thermo}}\) is the thermodynamic pressure, \(\mu\) is the shear viscosity, and \(\lambda\) is the longitudinal viscosity. - **Mechanical pressure**: \[ P_{\text{mech}} = -\frac{1}{3} (\sigma_{11} + \sigma_{22} + \sigma_{33}) = P_{\text{thermo}} - \mu_b \nabla \cdot \vec{u} \] - **Tisza formula**: \[ \mu_b = \rho_{\text{eq}} a^2 \frac{\gamma - 1}{\gamma} \sum_i \frac{c_v,i}{c_v} \tau_i \] where \(\rho_{\text{eq}}\) is the gas density in the equilibrium state, \(a\) is the sound speed in the non - viscous case, \(\gamma\) is the specific heat ratio in the equilibrium state, \(c_v,i\) is the heat capacity of the \(i\) - th internal mode, \(c_v\) is the total heat capacity of the gas, and \(\tau_i\) is the relaxation time of the \(i\) - th internal mode. - **Kustova formula**: \[ \mu_b = \frac{p R (c_{\text{int}} / c_v)^2 (c_{\text{rot}} / \tau_{\text{rot}} + c_{\text{vib}} / \tau_{\text{vib}})^{-1}} \]