1. |
- Gassner, Gregor J, et al.
(författare)
-
The BR1 scheme is stable for the compressible Navier–Stokes equations
- 2018
-
Ingår i: Journal of Scientific Computing. - : Springer. - 0885-7474 .- 1573-7691. ; 77:1, s. 154-200
-
Tidskriftsartikel (refereegranskat)abstract
- We show how to modify the original Bassi and Rebay scheme (BR1) [F. Bassi and S. Rebay, A High Order Accurate Discontinuous Finite Element Method for the Numerical Solution of the Compressible Navier-Stokes Equations, Journal of Computational Physics, 131:267–279, 1997] to get a provably stable discontinuous Galerkin collocation spectral element method (DGSEM) with Gauss-Lobatto (GL) nodes for the compressible Navier-Stokes equations (NSE) on three dimensional curvilinear meshes.Specifically, we show that the BR1 scheme can be provably stable if the metric identities are discretely satisfied, a two-point average for the metric terms is used for the contravariant fluxes in the volume, an entropy conserving split form is used for the advective volume integrals, the auxiliary gradients for the viscous terms are computed from gradients of entropy variables, and the BR1 scheme is used for the interface fluxes.Our analysis shows that even with three dimensional curvilinear grids, the BR1 fluxes do not add artificial dissipation at the interior element faces. Thus, the BR1 interface fluxes preserve the stability of the discretization of the advection terms and we get either energy stability or entropy-stability for the linear or nonlinear compressible NSE, respectively.
|
|
2. |
- Winters, Andrew Ross, et al.
(författare)
-
A comparative study on polynomial dealiasing and split form discontinuous Galerkin schemes for under-resolved turbulence computations
- 2018
-
Ingår i: Journal of Computational Physics. - : Elsevier. - 0021-9991 .- 1090-2716. ; 372, s. 1-21
-
Tidskriftsartikel (refereegranskat)abstract
- This work focuses on the accuracy and stability of high-order nodal discontinuous Galerkin (DG) methods for under-resolved turbulence computations. In particular we consider the inviscid Taylor-Green vortex (TGV) flow to analyse the implicit large eddy simulation (iLES) capabilities of DG methods at very high Reynolds numbers. The governing equations are discretised in two ways in order to suppress aliasing errors introduced into the discrete variational forms due to the under-integration of non-linear terms. The first, more straightforward way relies on consistent/over-integration, where quadrature accuracy is improved by using a larger number of integration points, consistent with the degree of the non-linearities. The second strategy, originally applied in the high-order finite difference community, relies on a split (or skew-symmetric) form of the governing equations. Different split forms are available depending on how the variables in the non-linear terms are grouped. The desired split form is then built by averaging conservative and non-conservative forms of the governing equations, although conservativity of the DG scheme is fully preserved. A preliminary analysis based on Burgers’ turbulence in one spatial dimension is conducted and shows the potential of split forms in keeping the energy of higher-order polynomial modes close to the expected levels. This indicates that the favourable dealiasing properties observed from split-form approaches in more classical schemes seem to hold for DG. The remainder of the study considers a comprehensive set of (under-resolved) computations of the inviscid TGV flow and compares the accuracy and robustness of consistent/over-integration and split form discretisations based on the local Lax-Friedrichs and Roe-type Riemann solvers. Recent works showed that relevant split forms can stabilize higher-order inviscid TGV test cases otherwise unstable even with consistent integration. Here we show that stable high-order cases achievable with both strategies have comparable accuracy, further supporting the good dealiasing properties of split form DG. The higher-order cases achieved only with split form schemes also displayed all the main features expected from consistent/over-integration. Among test cases with the same number of degrees of freedom, best solution quality is obtained with Roe-type fluxes at moderately high orders (around sixth order). Solutions obtained with very high polynomial orders displayed spurious features attributed to a sharper dissipation in wavenumber space. Accuracy differences between the two dealiasing strategies considered were, however, observed for the low-order cases, which also yielded reduced solution quality compared to high-order results.
|
|
3. |
- Bohm, Marvin, et al.
(författare)
-
An entropy stable nodal discontinuous Galerkin method for the resistive MHD equations. Part I : Theory and numerical verification
- 2018
-
Ingår i: Journal of Computational Physics. - : Elsevier. - 0021-9991 .- 1090-2716.
-
Tidskriftsartikel (refereegranskat)abstract
- The first paper of this series presents a discretely entropy stable discontinuous Galerkin (DG) method for the resistive magnetohydrodynamics (MHD) equations on three-dimensional curvilinear unstructured hexahedral meshes. Compared to other fluid dynamics systems such as the shallow water equations or the compressible Navier-Stokes equations, the resistive MHD equations need special considerations because of the divergence-free constraint on the magnetic field. For instance, it is well known that for the symmetrization of the ideal MHD system as well as the continuous entropy analysis a non-conservative term proportional to the divergence of the magnetic field, typically referred to as the Powell term, must be included. As a consequence, the mimicry of the continuous entropy analysis in the discrete sense demands a suitable DG approximation of the non-conservative terms in addition to the ideal MHD terms.This paper focuses on the resistive MHD equations: Our first contribution is a proof that the resistive terms are symmetric and positive-definite when formulated in entropy space as gradients of the entropy variables, which enables us to show that the entropy inequality holds for the resistive MHD equations. This continuous analysis is the key for our DG discretization and guides the path for the construction of an approximation that discretely mimics the entropy inequality, typically termed entropy stability. Our second contribution is a detailed derivation and analysis of the discretization on three-dimensional curvilinear meshes. The discrete analysis relies on the summation-by-parts property, which is satisfied by the DG spectral element method (DGSEM) with Legendre-Gauss-Lobatto (LGL) nodes. Although the divergence-free constraint is included in the non-conservative terms, the resulting method has no particular treatment of the magnetic field divergence errors, which might pollute the solution quality. Our final contribution is the extension of the standard resistive MHD equations and our DG approximation with a divergence cleaning mechanism that is based on a generalized Lagrange multiplier (GLM).As a conclusion to the first part of this series, we provide detailed numerical validations of our DGSEM method that underline our theoretical derivations. In addition, we show a numerical example where the entropy stable DGSEM demonstrates increased robustness compared to the standard DGSEM.
|
|
4. |
- Bohm, Marvin, et al.
(författare)
-
Multi-element SIAC Filter for Shock Capturing Applied to High-Order Discontinuous Galerkin Spectral Element Methods
- 2019
-
Ingår i: Journal of Scientific Computing. - : Springer-Verlag New York. - 0885-7474 .- 1573-7691. ; 81:2, s. 820-844
-
Tidskriftsartikel (refereegranskat)abstract
- 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. (J Sci Comput 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 equations to several standard test problems with a variety of boundary conditions.
|
|
5. |
- Clausen, Frederik Banch, et al.
(författare)
-
Recommendation for validation and quality assurance of non-invasive prenatal testing for foetal blood groups and implications for IVD risk classification according to EU regulations
- 2022
-
Ingår i: Vox Sanguinis. - : Wiley. - 0042-9007 .- 1423-0410. ; 117:2, s. 157-165
-
Forskningsöversikt (refereegranskat)abstract
- Background and Objectives: Non-invasive assays for predicting foetal blood group status in pregnancy serve as valuable clinical tools in the management of pregnancies at risk of detrimental consequences due to blood group antigen incompatibility. To secure clinical applicability, assays for non-invasive prenatal testing of foetal blood groups need to follow strict rules for validation and quality assurance. Here, we present a multi-national position paper with specific recommendations for validation and quality assurance for such assays and discuss their risk classification according to EU regulations. Materials and Methods: We reviewed the literature covering validation for in-vitro diagnostic (IVD) assays in general and for non-invasive foetal RHD genotyping in particular. Recommendations were based on the result of discussions between co-authors. Results: In relation to Annex VIII of the In-Vitro-Diagnostic Medical Device Regulation 2017/746 of the European Parliament and the Council, assays for non-invasive prenatal testing of foetal blood groups are risk class D devices. In our opinion, screening for targeted anti-D prophylaxis for non-immunized RhD negative women should be placed under risk class C. To ensure high quality of non-invasive foetal blood group assays within and beyond the European Union, we present specific recommendations for validation and quality assurance in terms of analytical detection limit, range and linearity, precision, robustness, pre-analytics and use of controls in routine testing. With respect to immunized women, different requirements for validation and IVD risk classification are discussed. Conclusion: These recommendations should be followed to ensure appropriate assay performance and applicability for clinical use of both commercial and in-house assays.
|
|
6. |
- Derigs, Dominik, et al.
(författare)
-
A novel averaging technique for discrete entropy-stable dissipation operators for ideal MHD
- 2017
-
Ingår i: Journal of Computational Physics. - : Elsevier. - 0021-9991 .- 1090-2716. ; 330, s. 624-632
-
Tidskriftsartikel (refereegranskat)abstract
- Entropy stable schemes can be constructed with a specific choice of the numerical flux function. First, an entropy conserving flux is constructed. Secondly, an entropy stable dissipation term is added to this flux to guarantee dissipation of the discrete entropy. Present works in the field of entropy stable numerical schemes are concerned with thorough derivations of entropy conservative fluxes for ideal MHD. However, as we show in this work, if the dissipation operator is not constructed in a very specific way, it cannot lead to a generally stable numerical scheme. The two main findings presented in this paper are that the entropy conserving flux of Ismail & Roe can easily break down for certain initial conditions commonly found in astrophysical simulations, and that special care must be taken in the derivation of a discrete dissipation matrix for an entropy stable numerical scheme to be robust. We present a convenient novel averaging procedure to evaluate the entropy Jacobians of the ideal MHD and the compressible Euler equations that yields a discretization with favorable robustness properties.
|
|
7. |
- Derigs, Dominik, et al.
(författare)
-
A novel high-order, entropy stable, 3D AMR MHD solver with guaranteed positive pressure
- 2016
-
Ingår i: Journal of Computational Physics. - : Elsevier. - 0021-9991 .- 1090-2716. ; 317, s. 223-256
-
Tidskriftsartikel (refereegranskat)abstract
- We describe a high-order numerical magnetohydrodynamics (MHD) solver built upon a novel non-linear entropy stable numerical flux function that supports eight travelling wave solutions. By construction the solver conserves mass, momentum, and energy and is entropy stable. The method is designed to treat the divergence-free constraint on the magnetic field in a similar fashion to a hyperbolic divergence cleaning technique. The solver described herein is especially well-suited for flows involving strong discontinuities. Furthermore, we present a new formulation to guarantee positivity of the pressure. We present the underlying theory and implementation of the new solver into the multi-physics, multi-scale adaptive mesh refinement (AMR) simulation code FLASH (http://flash.uchicago.edu). The accuracy, robustness and computational efficiency is demonstrated with a number of tests, including comparisons to available MHD implementations in FLASH.
|
|
8. |
- Derigs, Dominik, et al.
(författare)
-
Entropy Stable Finite Volume Approximations for Ideal Magnetohydrodynamics
- 2018
-
Ingår i: Jahresbericht der Deutschen Mathematiker-Vereinigung (Teubner). - : Springer Berlin/Heidelberg. - 0012-0456 .- 1869-7135. ; 120:3, s. 153-219
-
Tidskriftsartikel (refereegranskat)abstract
- This article serves as a summary outlining the mathematical entropy analysis of the ideal magnetohydrodynamic (MHD) equations. We select the ideal MHD equations as they are particularly useful for mathematically modeling a wide variety of magnetized fluids. In order to be self-contained we first motivate the physical properties of a magnetic fluid and how it should behave under the laws of thermodynamics. Next, we introduce a mathematical model built from hyperbolic partial differential equations (PDEs) that translate physical laws into mathematical equations. After an overview of the continuous analysis, we thoroughly describe the derivation of a numerical approximation of the ideal MHD system that remains consistent to the continuous thermodynamic principles. The derivation of the method and the theorems contained within serve as the bulk of the review article. We demonstrate that the derived numerical approximation retains the correct entropic properties of the continuous model and show its applicability to a variety of standard numerical test cases for MHD schemes. We close with our conclusions and a brief discussion on future work in the area of entropy consistent numerical methods and the modeling of plasmas.
|
|
9. |
- Derigs, Dominik, et al.
(författare)
-
Ideal GLM-MHD : About the entropy consistent nine-wave magnetic field divergence diminishing ideal magnetohydrodynamics equations
- 2018
-
Ingår i: Journal of Computational Physics. - : Elsevier. - 0021-9991 .- 1090-2716. ; 364, s. 420-467
-
Tidskriftsartikel (refereegranskat)abstract
- The paper presents two contributions in the context of the numerical simulation of magnetized fluid dynamics. First, we show how to extend the ideal magnetohydrodynamics (MHD) equations with an inbuilt magnetic field divergence cleaning mechanism in such a way that the resulting model is consistent with the second law of thermodynamics. As a byproduct of these derivations, we show that not all of the commonly used divergence cleaning extensions of the ideal MHD equations are thermodynamically consistent. Secondly, we present a numerical scheme obtained by constructing a specific finite volume discretization that is consistent with the discrete thermodynamic entropy. It includes a mechanism to control the discrete divergence error of the magnetic field by construction and is Galilean invariant. We implement the new high-order MHD solver in the adaptive mesh refinement code FLASH where we compare the divergence cleaning efficiency to the constrained transport solver available in FLASH (unsplit staggered mesh scheme).
|
|
10. |
- Fernandez, David C. Del Rey, et al.
(författare)
-
Entropy-stable p-nonconforming discretizations with the summation-by-parts property for the compressible Navier-Stokes equations
- 2020
-
Ingår i: Computers & Fluids. - : PERGAMON-ELSEVIER SCIENCE LTD. - 0045-7930 .- 1879-0747. ; 210
-
Tidskriftsartikel (refereegranskat)abstract
- The entropy-conservative/stable, curvilinear, nonconforming, p-refinement algorithm for hyperbolic conservation laws of Del Rey Fernandez et al. (2019) is extended from the compressible Euler equations to the compressible Navier-Stokes equations. A simple and flexible coupling procedure with planar interpolation operators between adjoining nonconforming elements is used. Curvilinear volume metric terms are numerically approximated via a minimization procedure and satisfy the discrete geometric conservation law conditions. Distinct curvilinear surface metrics are used on the adjoining interfaces to construct the interface coupling terms, thereby localizing the discrete geometric conservation law constraints to each individual element. The resulting scheme is entropy conservative/stable, element-wise conservative, and freestream preserving. Viscous interface dissipation operators that retain the entropy stability of the base scheme are developed. The accuracy and stability of the resulting numerical scheme are shown to be comparable to those of the original conforming scheme in Carpenter et al. (2014) and Parsani et al. (2016), i.e., this scheme achieves similar to p 1/2 convergence on geometrically high-order distorted element grids; this is demonstrated in the context of the viscous shock problem, the Taylor-Green vortex problem at a Reynolds number of Re = 1, 600, and a subsonic turbulent flow past a sphere at Re = 2, 000. (C) 2020 Published by Elsevier Ltd.
|
|