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.
Existence and asymptotic stability of relaxation discrete shock profiles
