Abstract:In this work, we propose a novel selective discontinuity sensor approach for numerical simulations of the compressible Navier-Stokes equations. Since transformation to characteristic space is already a common approach to reduce high-frequency oscillations during interpolation to cell interfaces, we exploit the characteristic wave structure of the Euler equations to selectively treat the various waves that the equations comprise. The approach uses the Ducros shock sensing criterion to detect and limit oscillations due to shocks while applying a different criterion to detect and limit oscillations due to contact discontinuities. Furthermore, the method is general in the sense that it can be applied to any method that employs characteristic transformation and shock sensors. However, in the present work, we focus on the Gradient-Based Reconstruction family of schemes. A series of inviscid and viscous test cases containing various types of discontinuities are carried out. The proposed method is shown to markedly reduce high-frequency oscillations that arise due to improper treatment of the various discontinuities; i.e., applying the Ducros shock sensor in a flow where a strong contact discontinuity is present. Moreover, the proposed method is shown to predict similar volume-averaged kinetic energy and enstrophy profiles for the Taylor-Green vortex simulation compared to the base Ducros sensor, indicating that it does not introduce unnecessary numerical dissipation when there are no contact discontinuities in the flow.

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to develop a new method for selective discontinuity sensors for the numerical simulation of the Navier - Stokes equations in compressible fluids. Specifically, the author proposes a selective discontinuity detector based on the characteristic wave structure of the Euler equations, aiming to handle various types of discontinuities more accurately, especially shock waves and contact discontinuities. ### Main Problems 1. **Accurate Handling of Shock Waves and Contact Discontinuities**: - Shock waves: They are formed when the local flow velocity exceeds the local speed of sound. - Contact discontinuities: They are formed when fluids with different velocities meet or at the interface between two fluids with different densities. 2. **Deficiencies of Existing Methods**: - Existing discontinuity detectors (such as the Ducros sensor) are mainly aimed at shock waves, but have limited detection ability for contact discontinuities. - Excessive dissipation: Unnecessary numerical dissipation is introduced in regions without discontinuities, affecting the accuracy of the solution. ### Solutions 1. **Characteristic - Wave Selective Discontinuity Sensor**: - Utilize the characteristic wave structure of the Euler equations, including acoustic waves, entropy waves, and shear/vortex waves. - Handle shock waves and contact discontinuities separately through different detection criteria. - Use the Ducros sensor to detect shock waves and use the contact discontinuity detector proposed by Frahan et al. 2. **Extension of the Gradient - Based Reconstruction (GBR) Method**: - Combine with the Ducros sensor to reduce high - frequency oscillations caused by shock waves. - Propose a new selective discontinuity detector to handle contact discontinuities more accurately. ### Method Overview 1. **Characteristic Variable Transformation**: - Calculate the Roe - averaged variables and construct the left and right eigenvectors of the characteristic space. - Transform the conservative variables and their derivatives into the characteristic space. 2. **Discontinuity Detection**: - Use the Ducros sensor to detect shock waves. - Use the Frahan detector to detect contact discontinuities. - Select an appropriate limiting strategy according to the detection results. 3. **Numerical Flux Reconstruction**: - Perform unrestricted interpolation in the characteristic space. - Apply a non - linear limiting strategy according to the detection results. - Project the characteristic variables back to the physical field. ### Numerical Results - The effectiveness of the new method has been verified through multiple test cases, especially in handling shock waves and contact discontinuities in multi - dimensional flows. - Compared with the existing TENO scheme, the new method can better preserve small - scale vortex structures and reduce numerical dissipation. ### Conclusions - The new method performs well in handling shock waves and contact discontinuities and can accurately capture various types of discontinuities while maintaining low numerical dissipation. - This method has broad application prospects, especially in high - precision numerical simulations.

A Wave Appropriate Discontinuity Sensor Approach for Compressible Flows

Stable, entropy-consistent, and localized artificial-diffusivity method for capturing discontinuities

Novel, simple and robust contact-discontinuity capturing schemes for high speed compressible flows

Valid Physical Processes from Numerical Discontinuities in Computational Fluid Dynamics

Consistent Interface Capturing Adaptive Reconstruction Approach for Viscous Compressible Multicomponent Flows

Discontinuity-Capturing Finite Element Computation of Unsteady Flow with Adaptive Unstructured Mesh

Discontinuity-resolving shock-capturing schemes on unstructured grids

A generalized adaptive central-upwind scheme for compressible flow simulations and preventing spurious vortices

Assessment of a high-order shock-capturing central-difference scheme for hypersonic turbulent flow simulations

A Time-Accurate Dynamic Mesh Approach for High-Order Shock Computation

A high order discontinuous Galerkin method for the symmetric form of the anisotropic viscoelastic wave equation

High-Resolution Viscous Terms Discretization and ILW Solid Wall Boundary Treatment for the Navier–Stokes Equations

An Immersed Boundary Formulation for Simulating High-Speed Compressible Viscous Flows with Moving Solids

A Dissipation‐free Numerical Method to Solve One‐dimensional Hyperbolic Flow Equations

A Moving Embedded Boundary Approach For The Compressible Navier-Stokes Equations In A Block-Structured Adaptive Refinement Framework

An Adaptive Discontinuity Fitting Technique on Unstructured Dynamic Grids

Centralized Gradient-Based Reconstruction for Wall Modelled Large Eddy Simulations of Hypersonic Boundary Layer Transition

An advanced algorithm for solving incompressible fluid dynamics: from Navier–Stokes to Poisson equations

A displacement-based material point method for weakly compressible free-surface flows

Application of Central-Weighted Essentially Non-Oscillatory Finite-Volume Interface-Capturing Schemes for Modeling Cavitation Induced by an Underwater Explosion

A note on reducing spurious pressure oscillations in fully conservative discontinuous Galerkin simulations of multicomponent flows