A Legendre spectral viscosity (LSV) method applied to shock capturing for high-order flux reconstruction schemes

• Published in: Journal of computational physics 2022-07, Vol.460, p.111157, Article 111157
• Main Authors: C. B. Sousa, Victor, Scalo, Carlo
• Format: Article
• Language: English
• Subjects: Basis functions; Burgers equation; Computational physics; Energy dissipation; Flux-reconstruction; High-order methods; Large eddy simulation; Mach reflection; Mathematical analysis; Polynomials; Reconstruction; Robustness (mathematics); Shock capturing; Spectral viscosity; Viscosity; Wave interaction
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2022.111157

•A new perspective on shock capturing for high-order flux reconstruction is introduced.•A dynamic model, active only in cells where nonlinear action is strong, is developed.•Spectral modulation of dissipation is performed with concentration at small scales.•Polynomial orders higher than 120 have been used for one-dimensional shocks. A novel approach to shock capturing for high-order flux reconstruction schemes is derived based on the mathematical formalism of the filtered governing equations. While the latter perspective is only typically used for turbulence modeling in the context of Large-Eddy Simulations (LES), the novel Legendre Spectral Viscosity (LSV) subfilter scale (SFS) closure model is capable of performing simulations in the presence of shock-discontinuities. The LSV method exploits the set of hierarchical basis functions formed by the Legendre polynomials to extract the information on the energy content near the resolution limit and estimate the overall magnitude of the required SFS dissipative terms, resulting in a scheme that dynamically activates only in cells where nonlinear behavior is important. Additionally, the modulation of such terms in the Legendre spectral space allows for the concentration of the dissipative action at small scales. The proposed method is tested in canonical shock-dominated flow setups in both one and two dimensions. These include the 1D Burgers' problem, a 1D shock tube, a 1D shock-entropy wave interaction, a 2D inviscid shock-vortex interaction and a 2D double Mach reflection. Results showcase a high-degree of resolution power, achieving accurate results with a small number of degrees of freedom, and robustness, being able to capture shocks associated with the Burgers' equation and the 1D shock tube within a single cell with orders 120 and higher.