Abstract:The generalized Riemann problem (GRP) scheme for the Euler equations and gas-kinetic scheme (GKS) for the Boltzmann equation are two high resolution shock capturing schemes for fluid simulations. The difference is that one is based on the characteristics of the inviscid Euler equations and their wave interactions, and the other is based on the particle transport and collisions. The similarity between them is that both methods can use identical MUSCL-type initial reconstructions around a cell interface, and the spatial slopes on both sides of a cell interface involve in the gas evolution process and the construction of a time-dependent flux function. Although both methods have been applied successfully to the inviscid compressible flow computations, their performances have never been compared. Since both methods use the same initial reconstruction, any difference is solely coming from different underlying mechanism in their flux evaluation. Therefore, such a comparison is important to help us to understand the correspondence between physical modeling and numerical performances. Since GRP is so faithfully solving the inviscid Euler equations, the comparison can be also used to show the validity of solving the Euler equations itself. The numerical comparison shows that the GRP exhibits a slightly better computational efficiency, and has comparable accuracy with GKS for the Euler solutions in 1D case, but the GKS is more robust than GRP. For the 2D high Mach number flow simulations, the GKS is absent from the shock instability and converges to the steady state solutions faster than the GRP. The GRP has carbuncle phenomena, likes a cloud hanging over exact Riemann solvers. The GRP and GKS use different physical processes to describe the flow motion starting from a discontinuity. One is based on the assumption of equilibrium state with infinite number of particle collisions, and the other starts from the non-equilibrium free transport process to evolve into an equilibrium one through particle collisions. The different mechanism in the flux evaluation deviates their numerical performance. Through this study, we may conclude scientifically that it may NOT be valid to use the Euler equations as governing equations to construct numerical fluxes in a discretized space with limited cell resolution. To adapt the Navier–Stokes (NS) equations is NOT valid either because the NS equations describe the flow behavior on the hydrodynamic scale and have no any corresponding physics starting from a discontinuity. This fact alludes to the consistency of the Euler and Navier–Stokes equations with the continuum assumption and the necessity of a direct modeling of the physical process in the discretized space in the construction of numerical scheme when modeling very high Mach number flows. The development of numerical algorithm is similar to the modeling process in deriving the governing equations, but the control volume here cannot be shrunk to zero.

Comparison of computational codes for direct numerical simulations of turbulent Rayleigh-Bénard convection

A Three-Dimensional Fictitious Domain Method for Direct Numerical Simulations of Particle-Laden Flows with Heat Transfer

Assessment of Numerical Accuracy and Parallel Performance of OpenFOAM and its Reacting Flow Extension EBIdnsFoam

3D DNS of laminar Rayleigh-Bénard convection in a cylinder for incompressible fluid flow

Rescaled Equations of Fluid Motion for Well-Conditioned Direct Numerical Simulations of Rapidly Rotating Convection

Direct Numerical Simulation of high Prandtl number fluids and supercritical carbon dioxide canonical flows using the spectral element method

A comparison of 3-D spherical shell thermal convection results at low to moderate Rayleigh number using ASPECT (version 2.2.0) and CitcomS (version 3.3.1)

High-resolution High-Order Upwind Compact Scheme-Based Numerical Computation of Natural Convection Flows in a Square Cavity

HIGHER-ORDER DIFFERENCING METHOD WITH A MULTIGRID APPROACH FOR THE SOLUTION OF THE INCOMPRESSIBLE FLOW EQUATIONS AT HIGH REYNOLDS NUMBERS

Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows

A generalized curvilinear solver for spherical shell Rayleigh-Bénard convection

Comparison of Pebble Bed Velocity Profiles Between High-Fidelity and Intermediate-Fidelity Codes

Data-driven spectral turbulence modeling for Rayleigh-Bénard convection

A-priori evaluation of data-driven models for large-eddy simulations in Rayleigh–Bénard convection

Lattice Boltzmann simulations of three-dimensional thermal convective flows at high Rayleigh number

Toward Direct Numerical Simulation of Turbulent and Transitional Flow in Hexagonal Subchannels for Helium Conditions

Comparison of the Generalized Riemann Solver and the Gas-Kinetic Scheme for Inviscid Compressible Flow Simulations

Compressible turbulent convection at very high Rayleigh numbers

Comparison of the performance of high-order schemes based on the gas-kinetic and HLLC fluxes

Flow organization in non-Oberbeck-Boussinesq Rayleigh-Benard convection in water

Steady Rayleigh–Bénard convection between no-slip boundaries