Abstract:Generalized hydrodynamic equations were originally proposed to describe the rarefied non-equilibrium flows beyond the Navier–Stokes–Fourier (NSF) equations by constructing a non-equilibrium canonical distribution function for the mesoscopic Boltzmann equation. Subsequently, nonlinear coupled constitutive relations (NCCR) were developed under the adiabatic assumption, the Eu’s closure, and Myong’s simplification. NCCR+ was also proposed to include the omitted terms in Myong’s simplification. The goal of this paper is to assess the improvements in the accuracy due to NCCR+ and the influence of the bulk viscosity in one-dimensional steady shock wave structure for monatomic and diatomic gases. In order to solve NCCR+ equations, a coupled solution process based on the time-independent method for non-conserved variables is employed, which is different from the previous uncoupled solution process used for the NCCR equations. Shock structures in argon and nitrogen are calculated up to Mach 50, where the shock profile, inverse shock thickness, asymmetry parameter, and temperature–density separation distance are validated by DSMC and available experimental measurements. The results show that NCCR+ could not provide much improvement in accuracy compared to NCCR but adds to the computational cost, suggesting that Myong’s simplification used in NCCR is satisfactory. Both NCCR and NCCR+ perform better than NSF in computing the one-dimensional shock wave structure at high Mach numbers. It is also shown that the bulk viscosity has significant influence on the accuracy of prediction of the shock wave structure in a diatomic gas using both NSF and NCCR. Stokes’ hypothesis in conventional NSF is valid only for flows in a monatomic gas and for very low Mach number flows in a diatomic gas. Additionally, it is also found that the viscosity exponent s = 0.81 in variable-hard-sphere model provides a good fit with the experimental data for shock wave in argon using NCCR.

Molecular dynamics of shock fronts and their transitions in Lennard-Jonesium and Tin

Anisotropic Responses and Initial Decomposition of Condensed-Phase Β-Hmx under Shock Loadings Via Molecular Dynamics Simulations in Conjunction with Multiscale Shock Technique.

Computation of Shock Wave Structure Using a Simpler Set of Generalized Hydrodynamic Equations Based on Nonlinear Coupled Constitutive Relations

Molecular Dynamics Simulation on Stability of Converging Shocks

Quantum molecular dynamics simulation of shock-wave experiments in aluminum

Investigating shock wave propagation, evolution, and anisotropy using a moving window concurrent atomistic-continuum framework

Molecular Dynamics Simulation of Strong Shock Waves Propagating in Dense Deuterium, Taking into Consideration Effects of Excited Electrons

Molecular dynamics simulation of the shock response of materials: A tutorial

Molecular Dynamics Simulation of Strong Shock Waves Propagating in Dense Deuterium With the Effect of Excited Electrons

Ab initio calculation of the shock Hugoniot of bulk silicon

Multiscale Modeling of Shock Wave Localization in Porous Energetic Material

Discontinuous Molecular Dynamics for Rigid Bodies: Applications

Dynamics of Bond Breaking and Formation in Polyethylene Near Shock Front.

Numerical study of dynamic phase transitions in shock tube

Molecular dynamics simulation of shock-wave loading of copper and titanium

A Coupled Phase-Field and Crystal Plasticity Model for Understanding Shock-Induced Phase Transition of Iron

Shock-induced plasticity of semi-coherent 111 Cu-Ni multilayers

Shock Hugoniot calculations using on-the-fly machine learned force fields with ab initio accuracy

Shock interactions with heterogeneous energetic materials

Shear relaxation behind the shock front in $\langle$110$\rangle$ Molybdenum - From the Atomic Scale to Continuous Dislocation Fields

Analysis of simulation technique for steady shock waves in materials with analytical equations of state