Explicit discontinuous Galerkin methods for unsteady problems

En esta tesis se presentan varias técnicas de solución para métodos de Galerkin discontinuo (DG, p.s.s.e.i.1) de alto orden. Si bien el objetivo principal es el desarrollo de técnicas precisas y computacionalmente eficientes para resolver las ecuaciones de Navier-Stokes de flujo compresible en 3D por medio del método Discontinuous Galerkin Spectral Element Method (DGSEM), todos los métodos desarrollados en esta tesis se pueden aplicar a otras ecuaciones de conservación no lineales y a otros esquemas DG. Las técnicas de solución acá presentadas pueden ser clasificadas en dos grupos: (i) las que están orientadas a acelarar la integración temporal, y (ii) las que mejoran la discretización espacial con adaptación-p local. El primer grupo de técnicas disminuye el tiempo de cálculo dada una discretización espacial con un número fijo de grados de libertad (NDOF, p.s.s.e.i.) mediante la mejora de los algoritmos de solución. Por ese motivo, se presta especial atención al desarrollo de algoritmos eficientes para integración temporal implícita y de métodos multimalla. También se detalla una metodología para obtener el Jacobiano analítico DG para ecuaciónes de advección-difusión. Además, se demuestra que la discretización implícita del DGSEM en nodos de Gauss-Lobatto puede ser resuelta eficientemente mediante el método static condensation. Lo anterior permite resolver el problema con sistemas lineales más pequeños y mejor condicionados, lo cual disminuye los tiempos de cálculo para órdenes moderados. El segundo grupo de técnicas consiste en modificar localmente la discretización espacial para reducir NDOF sin deteriorar la precisión del método. Para lograr este obetivo, en esta tesis se desarrolla un nuevo estimador de error para el método DGSEM p-anisitrópico basado en el método _-estimation, el cual es más preciso y requiere menos recursos computacionales que enfoques previos. Posteriormente, este estimador se usa con éxito para realizar adaptación-p anisotrópica en casos estacionaros y no estacionarios, y se demuestra que éste se puede acoplar con métodos multimalla para acelerar más las simulaciones. Adicionalmente, se desarrollan dos técnicas de adaptación-p para casos no estacionarios: un algoritmo dinámico y un algoritmo estático, que es computacionalmente más eficiente en algunos casos estadísticamente estacionarios. ----------ABSTRACT---------- In this thesis, we develop space and time solution techniques for high-order Discontinuous Galerkin (DG) methods. Even though the main focus is the accurate and computationally efficient solution of the 3D compressible Navier-Stokes equations using the Discontinuous Galerkin Spectral Element Method (DGSEM), the methods developed here can be applied to other systems of nonlinear conservation laws and other DG methods. The solution techniques can be classified into two groups: (i) improved solution techniques for temporal discretizations, and (ii) local p-adaptation techniques that enhance the spatial discretization. The first group of techniques decreases the computational time for a given spatial discretization with a fixed number of degrees of freedom (NDOF) by improving the linear and nonlinear solver algorithms. In that regard, special attention is paid to the development of efficient implicit time-integration solvers and to multigrid methods (both explicit and implicit in time). We derive a methodology to obtain the analytical Jacobian of DG discretizations of nonlinear advectiondiffusion equations. Furthermore, we show that time-implicit DGSEM discretizations on Gauss-Lobatto points can be formulated as Schur complement problems and solved using the static-condensation method. This enables one to solve smaller, better-conditioned problems, which yields shorter computation times for moderate polynomial orders. The second group of techniques modifies the spatial discretization by reducing the NDOF while maintaining similar accuracy. To do that, we develop a novel truncation error estimator that uses the _-estimation method for the p-anisotropic DGSEM. The new method is more accurate and computationally cheaper than previous approaches. The error estimator is then successfully applied to perform local anisotropic p-adaptation for steady and unsteady flows. We show that the truncation error estimator can be readily coupled with multigrid techniques to obtain further speed-ups. Moreover, we develop two p-adaptation strategies for unsteady flows. One of them is a dynamic _-based anisotropic p-adaptation method. The other one is a static _-based anisotropic p-adaptation method, which may be computationally cheaper in statistically steady flows.

High-order methods are well-suited for the numerical simulation of complex compressible turbulent flows, but require additional stabilization techniques to capture instabilities arising from the underlying non-linear hyperbolic equations. This paper provides a detailed comparison of the effectiveness of entropy stable discontinuous Galerkin methods for the stabilization of compressible (wall-bounded) turbulent flows. For this investigation, an entropy stable discontinuous Galerkin spectral element method is applied on Gauss–Legendre and Gauss–Lobatto nodes. In the compressible regime, an additional stabilization technique for shock capturing based on a convex blending of a low-order finite volume with the high-order discontinuous Galerkin operator is utilized. The present investigation provides a systematic study from convergence tests, to the Taylor–Green vortex and finally to a more intricate turbulent wall-bounded 3D diffuser flow, encompassing both weakly compressible and compressible regimes. The comparison demonstrates that the DGSEM on Gauss–Lobatto nodes is generally less accurate for an equal amount of degrees of freedom. Conversely, it is faster than the DGSEM on Gauss—Legendre nodes due to a less severe time step restriction and simpler numerical operator. A performance comparison reveals that the DGSEM on Gauss–Lobatto nodes generally outperforms the DGSEM on Gauss nodes for underresolved turbulence in the subsonic regime on a periodic domain. Conversely, the opposite effect can be observed for wall-bounded flows as well as the supersonic regime, the latter depending of course on the chosen shock-capturing scheme. To the author’s knowledge, this is the first time for which a comparison of entropy stable DGSEM on Gauss–Lobatto and Gauss–Legendre has been performed for compressible, wall-bounded turbulent flows with separation.

This work is focused on the entropy analysis of a semi-discrete nodal discontinuous Galerkin spectral element method (DGSEM) on moving meshes for hyperbolic conservation laws. The DGSEM is constructed with a local tensor-product Lagrange-polynomial basis computed from Legendre–Gauss–Lobatto points. Furthermore, the collocation of interpolation and quadrature nodes is used in the spatial discretization. This approach leads to discrete derivative approximations in space that are summation-by-parts (SBP) operators. On a static mesh, the SBP property and suitable two-point flux functions, which satisfy the entropy condition from Tadmor, allow to mimic results from the continuous entropy analysis, if it is ensured that properties such as positivity preservation (of the water height, density or pressure) are satisfied on the discrete level. In this paper, Tadmor’s condition is extended to the moving mesh framework. We show that the volume terms in the semi-discrete moving mesh DGSEM do not contribute to the discrete entropy evolution when a two-point flux function that satisfies the moving mesh entropy condition is applied in the split form DG framework. The discrete entropy behavior then depends solely on the interface contributions and on the domain boundary contribution. The interface contributions are directly controlled by proper choice of the numerical element interface fluxes. If an entropy conserving two-point flux is chosen, the interface contributions vanish. To increase the robustness of the discretization we use so-called entropy stable two-point fluxes at the interfaces that are guaranteed entropy dissipative and thus give a bound on the interface contributions in the discrete entropy balance. The remaining boundary condition contributions depend on the type of the considered boundary condition. E.g. for periodic boundary conditions that are of entropy conserving type, our methodology with the entropy conserving interface fluxes is fully entropy conservative and with the entropy stable interface fluxes is guaranteed entropy stable. The presented proof does not require any exactness of quadrature in the spatial integrals of the variational forms. As it is the case for static meshes, these results rely on the assumption that additional properties like positivity preservation are satisfied on the discrete level. Besides the entropy stability, the time discretization of the moving mesh DGSEM will be investigated and it will be proven that the moving mesh DGSEM satisfies the free stream preservation property for an arbitrary s-stage Runge–Kutta method, when periodic boundary conditions are used. The theoretical properties of the moving mesh DGSEM will be validated by numerical experiments for the compressible Euler equations with periodic boundary conditions.

Entropy-stable Gauss collocation methods for ideal magneto-hydrodynamics

Spatial discontinuous Galerkin spectral element method for a family of chromatography models in CADET

Simulation of underresolved turbulent flows by adaptive filtering using the high order discontinuous Galerkin spectral element method

This paper summarizes our progress in the application of a high-order discontinuous Galerkin (DG) method for scale resolving fluid dynamics simulations on the Cray XC40 Hazel Hen cluster at HLRS. We present the large eddy simulation (LES) of flow around a wall mounted cylinder, a LES of flow around an airfoil at realistic Reynolds number using a recently introduced kinetic energy preserving flux formulation and a simulation of transitional flow in a low pressure turbine. Furthermore, it provides an overview over the parallel efficiency reached by our code when using up to 49, 152 CPUs and the latest developments of our DG framework.

The aim of this work is to provide a solver for viscoelastic multi-phase flows within the Bounded Support Spectral Solver (BoSSS) currently under development at the Chair of Fluid Dynamics at the Technical University of Darmstadt. The discretisation in BoSSS consists of a high-order Discontinuous Galerkin (DG) method for single-phase flow and a high-order eXtended Discontinuous Galerkin (XDG) method for the multi-phase purpose. The solver shall be used to investigate numerically the behaviour of viscoelastic droplets. The macroscopic Oldroyd B model which is used in a wide range of applications is chosen as the constitutive model. A detailed derivation of the system of equations including the modeling principles for the Oldroyd B model is presented. A DG discretisation of the system of equations including the Local Discontinuous Galerkin (LDG) method is presented after introducing the field of the DG method. The derivation of appropriate flux functions for the constitutive equations and the extra stress tensor are one of the key derivations of this scientific work. Difficulties arising in the numerical solution of viscoelastic flow problems for higher Weissenberg numbers for different discretisation methods are due to the convection dominated, mixed hyperbolic-elliptic-parabolic nature of the system of equations. Several strategies are presented which overcome these problems and are known from the literature. A key achievement of this scientific work is the application of the LDG method, originally developed for a hyperbolic system of equations for a Newtonian fluid, on the viscoelastic system of equations which renders methods for preserving ellipticity unnecessary. Furthermore, various strategiesn to enhance and to support convergence of the solution of the DG discretised system are presented. These are the Newton method with different approaches determining the Jacobian of the system, a homotopy continuation method based on the Weissenberg number for a better initial guess for the Newton method, and a troubled cell indicator combined with an artificial diffusion approach or an adaptive mesh refinement strategy, respectively. For the completeness of this work the XDG method is presented using a sharp interface approach with a signed distance level-set function as it is implemented in BoSSS. The single-phase solver is combined with these methods and appropriate flux functions for the interface are implemented to enable multi-phase applications for viscoelastic fluid. Several numerical experiments are conducted to verify and to validate the viscoelastic singlephase solver and to show the capability of the viscoelastic multi-phase solver to simulate viscoelastic droplets. Advantages and disadvantages of the implementation and an outlook for future research can be found in the conclusion.

Coupling local discontinuous and continuous galerkin methods for flow problems

Split form ALE discontinuous Galerkin methods with applications to under-resolved turbulent low-Mach number flows

An entropy stable high-order discontinuous Galerkin spectral element method for the Baer-Nunziato two-phase flow model

The discontinuous Galerkin Spectral Element Method (DG-SEM) is highly attractive for both DNS and LES of turbulent flows due to its low dispersion and dissipation errors, but also because of its good parallel scaling property. We show that especially for underresolved simulations the method has highly beneficial properties for LES, as well as for the direct treatment of the acoustic propagation. We also discuss different approaches for non-reflecting boundary conditions for DG-SEM and show the behavior of the methods for airfoil flows. Our main intend is to directly simulate trailing edge noise for airfoil flows at medium Reynolds numbers.

We further investigate the high order positivity-preserving discontinuous Galerkin (DG) methods for linear hyperbolic and radiative transfer equations developed in Yuan et al. (SIAM J Sci Comput 38:A2987–A3019, 2016). The DG methods in Yuan et al. (2016) can maintain positivity and high order accuracy, but they rely both on the scaling limiter in Zhang and Shu (J Comput Phys 229:8918–8934, 2010) and a rotational limiter, the latter may alter cell averages of the unmodulated DG scheme, thereby affecting conservation. Even though a Lax–Wendroff type theorem is proved in Yuan et al. (2016), guaranteeing convergence to weak solutions with correct shock speed when such rotational limiter is applied, it would still be desirable if a conservative DG method without changing the cell averages can be obtained which has both high order accuracy and positivity-preserving capability. In this paper, we develop and analyze such a DG method for both linear hyperbolic equations and radiative transfer equations. In the one-dimensional case, the method uses traditional DG space $$P^k$$ of piecewise polynomials of degree at most k. A key result is proved that the unmodulated DG solver in this case can maintain positivity of the cell average if the inflow boundary value and the source term are both positive, therefore the positivity-preserving framework in Zhang and Shu (2010) can be used to obtain a high order conservative positivity-preserving DG scheme. Unfortunately, in two-dimensions this is no longer the case. We show that the unmodulated DG solver based either on $$P^k$$ or $$Q^k$$ spaces (piecewise kth degree polynomials or piecewise tensor-product kth degree polynomials) could generate negative cell averages. We augment the DG space with additional functions so that the positivity of cell averages from the unmodulated DG solver can be restored, thereby leading to high order conservative positivity-preserving DG scheme based on these augmented DG spaces following the framework in Zhang and Shu (2010). Computational results are provided to demonstrate the good performance of our DG schemes.

The discontinuous Galerkin (dG) method was introduced by Reed and Hill [73] in 1973 for steady-state neutron transport as an hyperbolic problem. This was followed by other studies; by Bassi and Rebay [12] for the compressible Navier-Stokes equations, Cockburn and Shu [31] developed the local discontinuous Galerkin (ldG) method for advection-diffusion equations, and Peraire and Persson [69] introduced the compact discontinuous Galerkin (cdG) method. Independent of the dG methods, interior penalty (IP) methods have been developed for elliptic and parabolic problems by Douglas and Dupont [40] and Wheeler [95]. Then, in the 1980’s, Arnold et al. [6] proposed a unified classification and analysis of various kinds of dG methods. Later on, the dG methods were developed for elliptic problems [8, 24, 77] and for problems with advection [7, 14, 51, 56].

We describe an efficient parallelization strategy for the discontinuous Galerkin spectral element method, illustrated by a structured grid framework. Target applications are large scale DNS and LES calculations on massively parallel systems. Due to the simple and efficient formulation of the method, a parallelization aiming at one-element-per-processor calculations is feasible; a highly desired feature for emerging multi- and many-core architectures. We show scale-up tests on up to 131,000 processors.KeywordsDiscontinuous GalerkinSpectral ElementMPIParallelHPC