Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userShock capturing schemes based on nonuniform nonlinear weighted interpolation for conservation laws and their application as subcell limiters for FR/CPR
149
0
0.0
(
0
)
Added by
Huajun Zhu
Publication date
2021
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
A series of shock capturing schemes based on nonuniform nonlinear weighted interpolation on nonuniform points are developed for conservation laws. Smoothness indicator and discrete conservation laws are discussed. To make fair comparisons between different types of schemes, the properties of eigenvalues of spatial discretization matrices are proved. And the proposed schemes are compared with Weighted Compact Nonlinear Schemes (WCNS) and Flux Reconstruction or Correction Procedure via Reconstruction (FR/CPR) in dispersion, dissipation properties and numerical accuracy. Then, the proposed shock capturing schemes are used as subcell limiters for high-order FR/CPR and the hybrid scheme has superiority in data transformation and satisfying discrete conservation laws. Accuracy, discrete conservation laws and shock capturing properties are tested. Numerical results in one and two dimensions are provided to illustrate that the proposed schemes have good properties in shock capturing and can be applied as subcell limiters for FR/CPR.
rate research
Read More
A main disadvantage of many high-order methods for hyperbolic conservation laws lies in the famous Gibbs-Wilbraham phenomenon, once discontinuities appear in the solution. Due to the Gibbs-Wilbraham phenomenon, the numerical approximation will be polluted by spurious oscillations, which produce unphysical numerical solutions and might finally blow up the computation. In this work, we propose a new shock capturing procedure to stabilise high-order spectral element approximations. The procedure consists of going over from the original (polluted) approximation to a convex combination of the original approximation and its Bernstein reconstruction, yielding a stabilised approximation. The coefficient in the convex combination, and therefore the procedure, is steered by a discontinuity sensor and is only activated in troubled elements. Building up on classical Bernstein operators, we are thus able to prove that the resulting Bernstein procedure is total variation diminishing and preserves monotone (shock) profiles. Further, the procedure can be modified to not just preserve but also to enforce certain bounds for the solution, such as positivity. In contrast to other shock capturing methods, e.g. artificial viscosity methods, the new procedure does not reduce the time step or CFL condition and can be easily and efficiently implemented into any existing code. Numerical tests demonstrate that the proposed shock-capturing procedure is able to stabilise and enhance spectral element approximations in the presence of shocks.
A priori subcell limiting approach is developed for high-order flux reconstruction/correction procedure via reconstruction (FR/CPR) on two-dimensional unstructured quadrilateral meshes. Firstly, a modified indicator based on modal energy coefficients is proposed to detect troubled cells. Then, troubled cells are decomposed into nonuniform subcells and each subcell has one solution point. A second-order finite difference shock-capturing scheme based on nonuniform nonlinear weighted (NNW) interpolation is constructed to calculate troubled cells while smooth cells are calculated by the CPR method. Numerical investigations show that the subcell limiting strategy on unstructured quadrilateral meshes is robust in shock-capturing.
We propose an Exponential DG approach for numerically solving partial differential equations (PDEs). The idea is to decompose the governing PDE operators into linear (fast dynamics extracted by linearization) and nonlinear (the remaining after removing the former) parts, on which we apply the discontinuous Galerkin (DG) spatial discretization. The resulting semi-discrete system is then integrated using exponential time-integrators: exact for the former and approximate for the latter. By construction, our approach i) is stable with a large Courant number (Cr > 1); ii) supports high-order solutions both in time and space; iii) is computationally favorable compared to IMEX DG methods with no preconditioner; iv) requires comparable computational time compared to explicit RKDG methods, while having time stepsizes orders magnitude larger than maximal stable time stepsizes for explicit RKDG methods; v) is scalable in a modern massively parallel computing architecture by exploiting Krylov-subspace matrix-free exponential time integrators and compact communication stencil of DG methods. Various numerical results for both Burgers and Euler equations are presented to showcase these expected properties. For Burgers equation, we present detailed stability and convergence analyses for the exponential Euler DG scheme.
We formulate an oversampled radial basis function generated finite difference (RBF-FD) method to solve time-dependent nonlinear conservation laws. The analytic solutions of these problems are known to be discontinuous, which leads to occurrence of non-physical oscillations (Gibbs phenomenon) that pollute the numerical solutions and can make them unstable. We address these difficulties using a residual based artificial viscosity stabilization, where the residual of the conservation law indicates the approximate location of the shocks. The location is then used to locally apply an upwind viscosity term, which stabilizes the Gibbs phenomenon and does not smear the solution away from the shocks. The proposed method is numerically tested and proves to be robust and accurate when solving scalar conservation laws and systems of conservation laws, such as compressible Euler equations.
We build a multi-element variant of the smoothness increasing accuracy conserving (SIAC) shock capturing technique proposed for single element spectral methods by Wissink et al. (B.W. Wissink, G.B. Jacobs, J.K. Ryan, W.S. Don, and E.T.A. van der Weide. Shock regularization with smoothness-increasing accuracy-conserving Dirac-delta polynomial kernels. Journal of Scientific Computing, 77:579--596, 2018). In particular, the baseline scheme of our method is the nodal discontinuous Galerkin spectral element method (DGSEM) for approximating the solution of systems of conservation laws. It is well known that high-order methods generate spurious oscillations near discontinuities which can develop in the solution for nonlinear problems, even when the initial data is smooth. We propose a novel multi-element SIAC filtering technique applied to the DGSEM as a shock capturing method. We design the SIAC filtering such that the numerical scheme remains high-order accurate and that the shock capturing is applied adaptively throughout the domain. The shock capturing method is derived for general systems of conservation laws. We apply the novel SIAC filter to the two-dimensional Euler and ideal magnetohydrodynamics (MHD) equations to several standard test problems with a variety of boundary conditions.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
