Showing 1–25 of 25 results for author: Kopriva, D A

Mimetic Metrics for the DGSEM

Authors: Daniel Bach, Andrés Rueda-Ramírez, David A. Kopriva, Gregor J. Gassner

Abstract: Free-stream preservation is an essential property for numerical solvers on curvilinear grids. Key to this property is that the metric terms of the curvilinear mapping satisfy discrete metric identities, i.e., have zero divergence. Divergence-free metric terms are furthermore essential for entropy stability on curvilinear grids. We present a new way to compute the metric terms for discontinuous Gal… ▽ More Free-stream preservation is an essential property for numerical solvers on curvilinear grids. Key to this property is that the metric terms of the curvilinear mapping satisfy discrete metric identities, i.e., have zero divergence. Divergence-free metric terms are furthermore essential for entropy stability on curvilinear grids. We present a new way to compute the metric terms for discontinuous Galerkin spectral element methods (DGSEMs) that guarantees they are divergence-free. Our proposed mimetic approach uses projections that fit within the de Rham Cohomology. △ Less

Submitted 18 October, 2024; originally announced October 2024.

Energy Bounds for Discontinuous Galerkin Spectral Element Approximations of Well-Posed Overset Grid Problems for Hyperbolic Systems

Authors: David A. Kopriva, Andrew R. Winters, Jan Nordström

Abstract: We show that even though the Discontinuous Galerkin Spectral Element Method is stable for hyperbolic boundary-value problems, and the overset domain problem is well-posed in an appropriate norm, the energy of the approximation of the latter is bounded by data only for fixed polynomial order, mesh, and time. In the absence of dissipation, coupling of the overlapping domains is destabilizing by allo… ▽ More We show that even though the Discontinuous Galerkin Spectral Element Method is stable for hyperbolic boundary-value problems, and the overset domain problem is well-posed in an appropriate norm, the energy of the approximation of the latter is bounded by data only for fixed polynomial order, mesh, and time. In the absence of dissipation, coupling of the overlapping domains is destabilizing by allowing positive eigenvalues in the system to be integrated in time. This coupling can be stabilized in one space dimension by using the upwind numerical flux. To help provide additional dissipation, we introduce a novel penalty method that applies dissipation at arbitrary points within the overlap region and depends only on the difference between the solutions. We present numerical experiments in one space dimension to illustrate the implementation of the well-posed penalty formulation, and show spectral convergence of the approximations when sufficient dissipation is applied. △ Less

Submitted 15 October, 2024; v1 submitted 7 May, 2024; originally announced May 2024.

Analysis of an Explicit, High-Order Semi-Lagrangian Nodal Method

Authors: Gustaaf B. Jacobs, Hareshram Natarajan, Pavel Popov, David A. Kopriva

Abstract: A discrete analysis of the phase and dissipation errors of an explicit, semi-Lagrangian spectral element method is performed. The semi-Lagrangian method advects the Lagrange interpolant according the Lagrangian form of the transport equations and uses a least-square fit to correct the update for interface constraints of neighbouring elements. By assuming a monomial representation instead of the La… ▽ More A discrete analysis of the phase and dissipation errors of an explicit, semi-Lagrangian spectral element method is performed. The semi-Lagrangian method advects the Lagrange interpolant according the Lagrangian form of the transport equations and uses a least-square fit to correct the update for interface constraints of neighbouring elements. By assuming a monomial representation instead of the Lagrange form, a discrete version of the algorithm on a single element is derived. The resulting algebraic system lends itself to both a Modified Equation analysis and an eigenvalue analysis. The Modified Equation analysis, which Taylor expands the stencil at a single space location and time instance, shows that the semi-Lagrangian method is consistent with the PDE form of the transport equation in the limit that the element size goes to zero. The leading order truncation term of the Modified Equation is of the order of the degree of the interpolant which is consistent with numerical tests reported in the literature. The dispersion relations show that the method is negligibly dispersive, as is common for semi-Lagrangian methods. An eigenvalue analysis shows that the semi-Lagrangian method with a nodal Chebyshev interpolant is stable for a Courant-Friedrichs-Lewy condition based on the minimum collocation node spacing within an element that is greater than unity. △ Less

Submitted 6 February, 2023; v1 submitted 21 December, 2022; originally announced December 2022.

HORSES3D: a high-order discontinuous Galerkin solver for flow simulations and multi-physics applications

Authors: E. Ferrer, G. Rubio, G. Ntoukas, W. Laskowski, O. A. Mariño, S. Colombo, A. Mateo-Gabín, F. Manrique de Lara, D. Huergo, J. Manzanero, A. M. Rueda-Ramírez, D. A. Kopriva, E. Valero

Abstract: We present the latest developments of our High-Order Spectral Element Solver (HORSES3D), an open source high-order discontinuous Galerkin framework, capable of solving a variety of flow applications, including compressible flows (with or without shocks), incompressible flows, various RANS and LES turbulence models, particle dynamics, multiphase flows, and aeroacoustics. We provide an overview of t… ▽ More We present the latest developments of our High-Order Spectral Element Solver (HORSES3D), an open source high-order discontinuous Galerkin framework, capable of solving a variety of flow applications, including compressible flows (with or without shocks), incompressible flows, various RANS and LES turbulence models, particle dynamics, multiphase flows, and aeroacoustics. We provide an overview of the high-order spatial discretisation (including energy/entropy stable schemes) and anisotropic p-adaptation capabilities. The solver is parallelised using MPI and OpenMP showing good scalability for up to 1000 processors. Temporal discretisations include explicit, implicit, multigrid, and dual time-stepping schemes with efficient preconditioners. Additionally, we facilitate meshing and simulating complex geometries through a mesh-free immersed boundary technique. We detail the available documentation and the test cases included in the GitHub repository. △ Less

Submitted 20 June, 2022; originally announced June 2022.

On the Theoretical Foundation of Overset Grid Methods for Hyperbolic Problems II: Entropy Bounded Formulations for Nonlinear Conservation Laws

Authors: David A. Kopriva, Gregor J. Gassner, Jan Nordstrom

Abstract: We derive entropy conserving and entropy dissipative overlapping domain formulations for systems of nonlinear hyperbolic equations in conservation form, such as would be approximated by overset mesh methods. The entropy conserving formulation imposes two-way coupling at the artificial interface boundaries through nonlinear penalty functions that vanish when the solutions coincide. The penalty func… ▽ More We derive entropy conserving and entropy dissipative overlapping domain formulations for systems of nonlinear hyperbolic equations in conservation form, such as would be approximated by overset mesh methods. The entropy conserving formulation imposes two-way coupling at the artificial interface boundaries through nonlinear penalty functions that vanish when the solutions coincide. The penalty functions are expressed in terms of entropy conserving fluxes originally introduced for finite volume schemes. Entropy dissipation and additional coupling in the overlap region are added through the use of linear penalties. △ Less

Submitted 12 September, 2022; v1 submitted 21 March, 2022; originally announced March 2022.

MSC Class: 65N85

On the Theoretical Foundation of Overset Grid Methods for Hyperbolic Problems: Well-Posedness and Conservation

Authors: David A. Kopriva, Jan Nordström, Gregor J. Gassner

Abstract: We use the energy method to study the well-posedness of initial-boundary value problems approximated by overset mesh methods in one and two space dimensions for linear constant-coefficient hyperbolic systems. We show that in one space dimension, for both scalar equations and systems of equations, the problem where one domain partially oversets another is well-posed when characteristic coupling con… ▽ More We use the energy method to study the well-posedness of initial-boundary value problems approximated by overset mesh methods in one and two space dimensions for linear constant-coefficient hyperbolic systems. We show that in one space dimension, for both scalar equations and systems of equations, the problem where one domain partially oversets another is well-posed when characteristic coupling conditions are used. If a system cannot be diagonalized, as is ususally the case in multiple space dimensions, then the energy method does not give proper bounds in terms of initial and boundary data. For those problems, we propose a novel penalty approach. We show, by using a global energy that accounts for the energy in the overlap region of the domains, that under well-defined conditions on the coupling matrices the penalized overset domain problems are energy bounded, conservative, well-posed and have solutions equivalent to the original single domain problem. △ Less

Submitted 10 May, 2021; originally announced May 2021.

A Split-Form, Stable CG/DG-SEM for Wave Propagation Modeled by Linear Hyperbolic Systems

Authors: David A. Kopriva, Gregor J. Gassner

Abstract: We present a hybrid continuous and discontinuous Galerkin spectral element approximation that leverages the advantages of each approach. The continuous Galerkin approximation is used on interior element faces where the equation properties are continuous. A discontinuous Galerkin approximation is used at physical boundaries and if there is a jump in properties at a face. The approximation uses a sp… ▽ More We present a hybrid continuous and discontinuous Galerkin spectral element approximation that leverages the advantages of each approach. The continuous Galerkin approximation is used on interior element faces where the equation properties are continuous. A discontinuous Galerkin approximation is used at physical boundaries and if there is a jump in properties at a face. The approximation uses a split form of the equations and two-point fluxes to ensure stability for unstructured quadrilateral/hexahedral meshes with curved elements. The approximation is also conservative and constant state preserving on such meshes. Spectral accuracy is obtained for all examples, which include wave scattering at a discontinuous medium boundary. △ Less

Submitted 11 December, 2020; originally announced December 2020.

Stability of Discontinuous Galerkin Spectral Element Schemes for Wave Propagation when the Coefficient Matrices have Jumps

Authors: David A. Kopriva, Gregor J. Gassner, Jan Nordström

Abstract: We use the behavior of the $L_{2}$ norm of the solutions of linear hyperbolic equations with discontinuous coefficient matrices as a surrogate to infer stability of discontinuous Galerkin spectral element methods (DGSEM). Although the $L_{2}$ norm is not bounded by the initial data for homogeneous and dissipative boundary conditions for such systems, the $L_{2}$ norm is easier to work with than a… ▽ More We use the behavior of the $L_{2}$ norm of the solutions of linear hyperbolic equations with discontinuous coefficient matrices as a surrogate to infer stability of discontinuous Galerkin spectral element methods (DGSEM). Although the $L_{2}$ norm is not bounded by the initial data for homogeneous and dissipative boundary conditions for such systems, the $L_{2}$ norm is easier to work with than a norm that discounts growth due to the discontinuities. We show that the DGSEM with an upwind numerical flux that satisfies the Rankine-Hugoniot (or conservation) condition has the same energy bound as the partial differential equation does in the $L_{2}$ norm, plus an added dissipation that depends on how much the approximate solution fails to satisfy the Rankine-Hugoniot jump. △ Less

Submitted 22 December, 2020; v1 submitted 23 November, 2020; originally announced November 2020.

Construction of Modern Robust Nodal Discontinuous Galerkin Spectral Element Methods for the Compressible Navier-Stokes Equations

Authors: Andrew R. Winters, David A. Kopriva, Gregor J. Gassner, Florian Hindenlang

Abstract: Discontinuous Galerkin (DG) methods have a long history in computational physics and engineering to approximate solutions of partial differential equations due to their high-order accuracy and geometric flexibility. However, DG is not perfect and there remain some issues. Concerning robustness, DG has undergone an extensive transformation over the past seven years into its modern form that provide… ▽ More Discontinuous Galerkin (DG) methods have a long history in computational physics and engineering to approximate solutions of partial differential equations due to their high-order accuracy and geometric flexibility. However, DG is not perfect and there remain some issues. Concerning robustness, DG has undergone an extensive transformation over the past seven years into its modern form that provides statements on solution boundedness for linear and nonlinear problems. This chapter takes a constructive approach to introduce a modern incarnation of the DG spectral element method for the compressible Navier-Stokes equations in a three-dimensional curvilinear context. The groundwork of the numerical scheme comes from classic principles of spectral methods including polynomial approximations and Gauss-type quadratures. We identify aliasing as one underlying cause of the robustness issues for classical DG spectral methods. Removing said aliasing errors requires a particular differentiation matrix and careful discretization of the advective flux terms in the governing equations. △ Less

Submitted 5 May, 2020; originally announced May 2020.

Comments: 85 pages, 2 figures, book chapter

MSC Class: 65M70 (primary) 35M13; 65M20 (secondary) ACM Class: G.1.4; G.1.8

A Statically Condensed Discontinuous Galerkin Spectral Element Method on Gauss-Lobatto Nodes for the Compressible Navier-Stokes Equations

Authors: Andrés M. Rueda-Ramírez, Esteban Ferrer, David A. Kopriva, Gonzalo Rubio, Eusebio Valero

Abstract: We present a static-condensation method for time-implicit discretizations of the Discontinuous Galerkin Spectral Element Method on Gauss-Lobatto points (GL-DGSEM). We show that, when solving the compressible Navier-Stokes equations, it is possible to reorganize the linear system that results from the implicit time-integration of the GL-DGSEM as a Schur complement problem, which can be efficiently… ▽ More We present a static-condensation method for time-implicit discretizations of the Discontinuous Galerkin Spectral Element Method on Gauss-Lobatto points (GL-DGSEM). We show that, when solving the compressible Navier-Stokes equations, it is possible to reorganize the linear system that results from the implicit time-integration of the GL-DGSEM as a Schur complement problem, which can be efficiently solved using static condensation. The use of static condensation reduces the linear system size and improves the condition number of the system matrix, which translates into shorter computational times when using direct and iterative solvers. The statically condensed GL-DGSEM presented here can be applied to linear and nonlinear advection-diffusion partial differential equations in conservation form. To test it we solve the compressible Navier-Stokes equations with direct and Krylov subspace solvers, and we show for a selected problem that using the statically condensed GL-DGSEM leads to speed-ups of up to $200$ when compared to the time-explicit GL-DGSEM, and speed-ups of up to three when compared with the time-implicit GL-DGSEM that solves the global system. The GL-DGSEM has gained increasing popularity in recent years because it satisfies the summation-by-parts property, which enables the construction of provably entropy stable schemes, and because it is computationally very efficient. In this paper, we show that the GL-DGSEM has an additional advantage: It can be statically condensed. △ Less

Submitted 12 December, 2019; v1 submitted 3 November, 2019; originally announced November 2019.

Assessing Standard and Kinetic Energy Conserving Discontinuous Galerkin Formulations for Marginally Resolved Navier-Stokes Flows

Authors: B. F. Klose, G. B. Jacobs, D. A Kopriva

Abstract: The robustness and accuracy of marginally resolved discontinuous Galerkin spectral element computations are evaluated for the standard formulation and a kinetic energy conserving split form on complex flow problems of physical and engineering interest, including the flow over a square cylinder, an airfoil and a plane jet. It is shown that the kinetic energy conserving formulation is significantly… ▽ More The robustness and accuracy of marginally resolved discontinuous Galerkin spectral element computations are evaluated for the standard formulation and a kinetic energy conserving split form on complex flow problems of physical and engineering interest, including the flow over a square cylinder, an airfoil and a plane jet. It is shown that the kinetic energy conserving formulation is significantly more robust than the standard scheme for under-resolved simulations. A disadvantage of the split form is the restriction to Gauss-Lobatto nodes with the inherent underintegration and lower accuracy as compared to Gauss quadrature used with the standard scheme. While the results support the higher accuracy of the standard Gauss form, lower numerical robustness and spurious oscillations are evident in some cases, giving the advantage to the kinetic energy conserving scheme for marginally resolved numerical simulations. △ Less

Submitted 30 October, 2019; originally announced October 2019.

Entropy-stable discontinuous Galerkin approximation with summation-by-parts property for the incompressible Navier-Stokes/Cahn-Hilliard system

Authors: Juan Manzanero, Gonzalo Rubio, David A. Kopriva, Esteban Ferrer, Eusebio Valero

Abstract: We develop an entropy stable two-phase incompressible Navier--Stokes/Cahn--Hilliard discontinuous Galerkin (DG) flow solver method. The model poses the Cahn-Hilliard equation as the phase field method, a skew-symmetric form of the momentum equation, and an artificial compressibility method to compute the pressure. We design the model so that it satisfies an entropy law, including free- and no-slip… ▽ More We develop an entropy stable two-phase incompressible Navier--Stokes/Cahn--Hilliard discontinuous Galerkin (DG) flow solver method. The model poses the Cahn-Hilliard equation as the phase field method, a skew-symmetric form of the momentum equation, and an artificial compressibility method to compute the pressure. We design the model so that it satisfies an entropy law, including free- and no-slip wall boundary conditions with non-zero wall contact angle. We then construct a high-order DG approximation of the model that satisfies the SBP-SAT property. With the help of a discrete stability analysis, the scheme has two modes: an entropy conserving approximation with central advective fluxes and the Bassi-Rebay 1 (BR1) method for diffusion, and an entropy stable approximation with an exact Riemann solver for advection and interface stabilization added to the BR1 method. The scheme is applicable to, and the stability proofs hold for, three-dimensional unstructured meshes with curvilinear hexahedral elements. We test the convergence of the schemes on a manufactured solution, and their robustness by solving a flow initialized from random numbers. In the latter, we find that a similar scheme that does not satisfy an entropy inequality had 30% probability to fail, while the entropy stable scheme never does. We also solve the static and rising bubble test problems, and to challenge the solver capabilities we compute a three-dimensional pipe flow in the annular regime. △ Less

Submitted 25 October, 2019; v1 submitted 24 October, 2019; originally announced October 2019.

Naturally curved quadrilateral mesh generation using an adaptive spectral element solver

Authors: Julian Marcon, David A. Kopriva, Spencer J. Sherwin, Joaquim Peiró

Abstract: We describe an adaptive version of a method for generating valid naturally curved quadrilateral meshes. The method uses a guiding field, derived from the concept of a cross field, to create block decompositions of multiply connected two dimensional domains. The a priori curved quadrilateral blocks can be further split into a finer high-order mesh as needed. The guiding field is computed by a Lapla… ▽ More We describe an adaptive version of a method for generating valid naturally curved quadrilateral meshes. The method uses a guiding field, derived from the concept of a cross field, to create block decompositions of multiply connected two dimensional domains. The a priori curved quadrilateral blocks can be further split into a finer high-order mesh as needed. The guiding field is computed by a Laplace equation solver using a continuous Galerkin or discontinuous Galerkin spectral element formulation. This operation is aided by using $p$-adaptation to achieve faster convergence of the solution with respect to the computational cost. From the guiding field, irregular nodes and separatrices can be accurately located. A first version of the code is implemented in the open source spectral element framework Nektar++ and its dedicated high order mesh generation platform NekMesh. △ Less

Submitted 18 September, 2019; v1 submitted 12 August, 2019; originally announced August 2019.

Comments: 13 pages, 14 figures, accepted for publication in the Proceedings of 28th International Meshing Roundtable

MSC Class: 65N50

Journal ref: Proceedings of the 28th International Meshing Roundtable (2020) 254-266

Entropy-stable discontinuous Galerkin approximation with summation-by-parts property for the incompressible Navier-Stokes equations with variable density and artificial compressibility

Authors: Juan Manzanero, Gonzalo Rubio, David A Kopriva, Esteban Ferrer, Eusebio Valero

Abstract: We present a provably stable discontinuous Galerkin spectral element method for the incompressible Navier-Stokes equations with artificial compressibility and variable density. Stability proofs, which include boundary conditions, that follow a continuous entropy analysis are provided. We define a mathematical entropy function that combines the traditional kinetic energy and an additional energy te… ▽ More We present a provably stable discontinuous Galerkin spectral element method for the incompressible Navier-Stokes equations with artificial compressibility and variable density. Stability proofs, which include boundary conditions, that follow a continuous entropy analysis are provided. We define a mathematical entropy function that combines the traditional kinetic energy and an additional energy term for the artificial compressiblity, and derive its associated entropy conservation law. The latter allows us to construct a provably stable split-form nodal Discontinuous Galerkin (DG) approximation that satisfies the summation-by-parts simultaneous-approximation-term (SBP-SAT) property. The scheme and the stability proof are presented for general curvilinear three-dimensional hexahedral meshes. We use the exact Riemann solver and the Bassi-Rebay 1 (BR1) scheme at the inter-element boundaries for inviscid and viscous fluxes respectively, and an explicit low storage Runge-Kutta RK3 scheme to integrate in time. We assess the accuracy and robustness of the method by solving the Kovasznay flow, the inviscid Taylor-Green vortex, and the Rayleigh-Taylor instability. △ Less

Submitted 16 July, 2019; v1 submitted 12 July, 2019; originally announced July 2019.

A free-energy stable nodal discontinuous Galerkin approximation with summation-by-parts property for the Cahn-Hilliard equation

Authors: Juan Manzanero, Gonzalo Rubio, David A. Kopriva, Esteban Ferrer, Eusebio Valero

Abstract: We present a nodal Discontinuous Galerkin (DG) scheme for the Cahn-Hilliard equation that satisfies the summation-by-parts simultaneous-approximation-term (SBP-SAT) property. The latter permits us to show that the discrete free-energy is bounded, and as a result, the scheme is provably stable. The scheme and the stability proof are presented for general curvilinear three-dimensional hexahedral mes… ▽ More We present a nodal Discontinuous Galerkin (DG) scheme for the Cahn-Hilliard equation that satisfies the summation-by-parts simultaneous-approximation-term (SBP-SAT) property. The latter permits us to show that the discrete free-energy is bounded, and as a result, the scheme is provably stable. The scheme and the stability proof are presented for general curvilinear three-dimensional hexahedral meshes. We use the Bassi-Rebay 1 (BR1) scheme to compute interface fluxes, and an IMplicit-EXplicit (IMEX) scheme to integrate in time. Lastly, we test the theoretical findings numerically and present examples for two and three-dimensional problems. △ Less

Submitted 21 February, 2019; originally announced February 2019.

Stability of Wall Boundary Condition Procedures for Discontinuous Galerkin Spectral Element Approximations of the Compressible Euler Equations

Authors: Florian J. Hindenlang, Gregor J. Gassner, David A. Kopriva

Abstract: We perform a linear and entropy stability analysis for wall boundary condition procedures for discontinuous Galerkin spectral element approximations of the compressible Euler equations. Two types of boundary procedures are examined. The first defines a special wall boundary flux that incorporates the boundary condition. The other is the commonly used reflection condition where an external state is… ▽ More We perform a linear and entropy stability analysis for wall boundary condition procedures for discontinuous Galerkin spectral element approximations of the compressible Euler equations. Two types of boundary procedures are examined. The first defines a special wall boundary flux that incorporates the boundary condition. The other is the commonly used reflection condition where an external state is specified that has an equal and opposite normal velocity. The internal and external states are then combined through an approximate Riemann solver to weakly impose the boundary condition. We show that with the exact upwind and Lax-Friedrichs solvers the approximations are energy dissipative, with the amount of dissipation proportional to the square of the normal Mach number. Standard approximate Riemann solvers, namely Lax-Friedrichs, HLL, HLLC are entropy stable. The Roe flux is entropy stable under certain conditions. An entropy conserving flux with an entropy stable dissipation term (EC-ES) is also presented. The analysis gives insight into why these boundary conditions are robust in that they introduce large amounts of energy or entropy dissipation when the boundary condition is not accurately satisfied, e.g. due to an impulsive start or under resolution. △ Less

Submitted 23 April, 2019; v1 submitted 15 January, 2019; originally announced January 2019.

Comments: ICOSAHOM 2018 conference proceedings, 14 pages, 2 Figures

A High Resolution PDE Approach to Quadrilateral Mesh Generation

Authors: Julian Marcon, David A. Kopriva, Spencer J. Sherwin, Joaquim Peiró

Abstract: We describe a high order technique to generate quadrilateral decompositions and meshes for complex two dimensional domains using spectral elements in a field guided procedure. Inspired by cross field methods, we never actually compute crosses. Instead, we compute a high order accurate guiding field using a continuous Galerkin (CG) or discontinuous Galerkin (DG) spectral element method to solve a L… ▽ More We describe a high order technique to generate quadrilateral decompositions and meshes for complex two dimensional domains using spectral elements in a field guided procedure. Inspired by cross field methods, we never actually compute crosses. Instead, we compute a high order accurate guiding field using a continuous Galerkin (CG) or discontinuous Galerkin (DG) spectral element method to solve a Laplace equation for each of the field variables using the open source code Nektar++. The spectral method provides spectral convergence and sub-element resolution of the fields. The DG approximation allows meshing of corners that are not multiples of $π/2$ in a discretization consistent manner, when needed. The high order field can then be exploited to accurately find irregular nodes, and can be accurately integrated using a high order separatrix integration method to avoid features like limit cycles. The result is a mesh with naturally curved quadrilateral elements that do not need to be curved a posteriori to eliminate invalid elements. The mesh generation procedure is implemented in the open source mesh generation program NekMesh. △ Less

Submitted 30 August, 2019; v1 submitted 8 January, 2019; originally announced January 2019.

Comments: 31 pages, 21 figures, accepted for publication in Journal of Computational Physics

MSC Class: 65N50

Journal ref: Journal of Computational Physics 399C (2019) 108918

Free-Stream Preservation for Curved Geometrically Non-Conforming Discontinuous Galerkin Spectral Elements

Authors: David A. Kopriva, Florian Hindenlang, Thomas Boleman, Gregor J. Gassner

Abstract: The under integration of the volume terms in the discontinuous Galerkin spectral element approximation introduces errors at non-conforming element faces that do not cancel and lead to free-stream preservation errors. We derive volume and face conditions on the geometry under which a constant state is preserved. From those, we catalog six special cases on the geometry that preserve a constant state… ▽ More The under integration of the volume terms in the discontinuous Galerkin spectral element approximation introduces errors at non-conforming element faces that do not cancel and lead to free-stream preservation errors. We derive volume and face conditions on the geometry under which a constant state is preserved. From those, we catalog six special cases on the geometry that preserve a constant state, the most general being to approximate the geometry sub-parametrically to one half the order of the solution. Numerical examples are presented to illustrate the results. △ Less

Submitted 5 July, 2019; v1 submitted 13 September, 2018; originally announced September 2018.

Journal ref: J Sci Comput (2019) 79: 1389

Stability of Overintegration Methods for Nodal Discontinuous Galerkin Spectral Element Methods

Authors: David A. Kopriva

Abstract: We perform stability analyses for discontinuous Galerkin spectral element approximations of linear variable coefficient hyperbolic systems in three dimensional domains with curved elements. Although high order, the precision of the quadratures used are typically too low with respect to polynomial order associated with their arguments, which introduces aliasing errors that can destabilize an approx… ▽ More We perform stability analyses for discontinuous Galerkin spectral element approximations of linear variable coefficient hyperbolic systems in three dimensional domains with curved elements. Although high order, the precision of the quadratures used are typically too low with respect to polynomial order associated with their arguments, which introduces aliasing errors that can destabilize an approximation, especially when the solution is underresolved. We show that using a larger number of points in the volume quadrature, often called "overintegration", can eliminate the aliasing term associated with the volume, but introduces new aliasing errors at the surfaces that can destabilize the solution. Increased quadrature precision on both the volume and surface terms, on the other hand, leads to a stable approximation. The results support the findings of Mengaldo et al. [Dealiasing techniques for high-order spectral element methods on regular and irregular grids. Journal of Computational Physics, 299:56 -- 81, 2015] who found that fully consistent integration was more robust for the solution of compressible flows than the volume only version. △ Less

Submitted 16 July, 2017; originally announced July 2017.

MSC Class: 65M70; 65M12

Journal ref: J Sci Comput (2019) 79: 1389

Insights on aliasing driven instabilities for advection equations with application to Gauss-Lobatto discontinuous Galerkin methods

Authors: Juan Manzanero, Gonzalo Rubio, Esteban Ferrer, Eusebio Valero, David A. Kopriva

Abstract: We analyse instabilities due to aliasing errors when solving one dimensional non-constant advection speed equations and discuss means to alleviate these types of errors when using high order discontinuous Galerkin (DG) schemes. First, we compare analytical bounds for the continuous and discrete version of the PDEs. Whilst traditional $L^2$ norm energy bounds applied to the discrete PDE do not alwa… ▽ More We analyse instabilities due to aliasing errors when solving one dimensional non-constant advection speed equations and discuss means to alleviate these types of errors when using high order discontinuous Galerkin (DG) schemes. First, we compare analytical bounds for the continuous and discrete version of the PDEs. Whilst traditional $L^2$ norm energy bounds applied to the discrete PDE do not always predict the physical behaviour of the continuous version of the equation, more strict elliptic norm bounds correctly bound the behaviour of the continuous PDE. Having derived consistent bounds, we analyse the effectiveness of two stabilising techniques: over-integration and split form variations (conservative, non-conservative and skew-symmetric). Whilst the former is shown to not alleviate aliasing in general, the latter ensures an aliasing-free solution if the splitting form of the discrete PDE is consistent with the continuous equation. The success of split form de-aliasing is restricted to DG schemes with the summation-by-parts simultaneous-approximation-term (SBP-SAT) properties (e.g. DG with Gauss-Lobatto points). Numerical experiments are included to illustrate the theoretical findings. △ Less

Submitted 3 May, 2017; originally announced May 2017.

The BR1 Scheme is Stable for the Compressible Navier-Stokes Equations

Authors: Gregor J. Gassner, Andrew R. Winters, Florian J. Hindenlang, David A. Kopriva

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 th… ▽ More 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. △ Less

Submitted 9 April, 2018; v1 submitted 12 April, 2017; originally announced April 2017.

A Polynomial Spectral Calculus for Analysis of DG Spectral Element Methods

Authors: David A. Kopriva

Abstract: We introduce a polynomial spectral calculus that follows from the summation by parts property of the Legendre-Gauss-Lobatto quadrature. We use the calculus to simplify the analysis of two multidimensional discontinuous Galerkin spectral element approximations. We introduce a polynomial spectral calculus that follows from the summation by parts property of the Legendre-Gauss-Lobatto quadrature. We use the calculus to simplify the analysis of two multidimensional discontinuous Galerkin spectral element approximations. △ Less

Submitted 3 April, 2017; originally announced April 2017.

MSC Class: 65M70

Split Form Nodal Discontinuous Galerkin Schemes with Summation-By-Parts Property for the Compressible Euler Equations

Authors: Gregor J. Gassner, Andrew R. Winters, David A. Kopriva

Abstract: Fisher and Carpenter (\textit{High-order entropy stable finite difference schemes for non-linear conservation laws: Finite domains, Journal of Computational Physics, 252:518--557, 2013}) found a remarkable equivalence of general diagonal norm high-order summation-by-parts operators to a subcell based high-order finite volume formulation. This equivalence enables the construction of provably entrop… ▽ More Fisher and Carpenter (\textit{High-order entropy stable finite difference schemes for non-linear conservation laws: Finite domains, Journal of Computational Physics, 252:518--557, 2013}) found a remarkable equivalence of general diagonal norm high-order summation-by-parts operators to a subcell based high-order finite volume formulation. This equivalence enables the construction of provably entropy stable schemes by a specific choice of the subcell finite volume flux. We show that besides the construction of entropy stable high order schemes, a careful choice of subcell finite volume fluxes generates split formulations of quadratic or cubic terms. Thus, by changing the subcell finite volume flux to a specific choice, we are able to generate, in a systematic way, all common split forms of the compressible Euler advection terms, such as the Ducros splitting and the Kennedy and Gruber splitting. Although these split forms are not entropy stable, we present a systematic way to prove which of those split forms are at least kinetic energy preserving. With this, we show we construct a unified high-order split form DG framework. We investigate with three dimensional numerical simulations of the inviscid Taylor-Green vortex and show that the new split forms enhance the robustness of high order simulations in comparison to the standard scheme when solving turbulent vortex dominated flows. In fact, we show that for certain test cases, the novel split form discontinuous Galerkin schemes are more robust than the discontinuous Galerkin scheme with over-integration. △ Less

Submitted 23 September, 2016; v1 submitted 22 April, 2016; originally announced April 2016.

A Provably Stable Discontinuous Galerkin Spectral Element Approximation for Moving Hexahedral Meshes

Authors: David A. Kopriva, Andrew R. Winters, Marvin Bohm, Gregor J. Gassner

Abstract: We design a novel provably stable discontinuous Galerkin spectral element (DGSEM) approximation to solve systems of conservation laws on moving domains. To incorporate the motion of the domain, we use an arbitrary Lagrangian-Eulerian formulation to map the governing equations to a fixed reference domain. The approximation is made stable by a discretization of a skew-symmetric formulation of the pr… ▽ More We design a novel provably stable discontinuous Galerkin spectral element (DGSEM) approximation to solve systems of conservation laws on moving domains. To incorporate the motion of the domain, we use an arbitrary Lagrangian-Eulerian formulation to map the governing equations to a fixed reference domain. The approximation is made stable by a discretization of a skew-symmetric formulation of the problem. We prove that the discrete approximation is stable, conservative and, for constant coefficient problems, maintains the free-stream preservation property. We also provide details on how to add the new skew-symmetric ALE approximation to an existing discontinuous Galerkin spectral element code. Lastly, we provide numerical support of the theoretical results. △ Less

Submitted 30 October, 2015; originally announced October 2015.

An Entropy Stable Nodal Discontinuous Galerkin Method for the Two Dimensional Shallow Water Equations on Unstructured Curvilinear Meshes with Discontinuous Bathymetry

Authors: Niklas Wintermeyer, Andrew R. Winters, Gregor J. Gassner, David A. Kopriva

Abstract: We design an arbitrary high-order accurate nodal discontinuous Galerkin spectral element approximation for the nonlinear two dimensional shallow water equations with non-constant, possibly discontinuous, bathymetry on unstructured, possibly curved, quadrilateral meshes. The scheme is derived from an equivalent flux differencing formulation of the split form of the equations. We prove that this dis… ▽ More We design an arbitrary high-order accurate nodal discontinuous Galerkin spectral element approximation for the nonlinear two dimensional shallow water equations with non-constant, possibly discontinuous, bathymetry on unstructured, possibly curved, quadrilateral meshes. The scheme is derived from an equivalent flux differencing formulation of the split form of the equations. We prove that this discretisation exactly preserves the local mass and momentum. Furthermore, combined with a special numerical interface flux function, the method exactly preserves the mathematical entropy, which is the total energy for the shallow water equations. By adding a specific form of interface dissipation to the baseline entropy conserving scheme we create a provably entropy stable scheme. That is, the numerical scheme discretely satisfies the second law of thermodynamics. Finally, with a particular discretisation of the bathymetry source term we prove that the numerical approximation is well-balanced. We provide numerical examples that verify the theoretical findings and furthermore provide an application of the scheme for a partial break of a curved dam test problem. △ Less

Submitted 22 June, 2016; v1 submitted 23 September, 2015; originally announced September 2015.