-
Unified $hp$-HDG Frameworks for Friedrichs' PDE systems
Authors:
Jau-Uei Chen,
Shinhoo Kang,
Tan Bui-Thanh,
John N. Shadid
Abstract:
This work proposes a unified $hp$-adaptivity framework for hybridized discontinuous Galerkin (HDG) method for a large class of partial differential equations (PDEs) of Friedrichs' type. In particular, we present unified $hp$-HDG formulations for abstract one-field and two-field structures and prove their well-posedness. In order to handle non-conforming interfaces we simply take advantage of HDG b…
▽ More
This work proposes a unified $hp$-adaptivity framework for hybridized discontinuous Galerkin (HDG) method for a large class of partial differential equations (PDEs) of Friedrichs' type. In particular, we present unified $hp$-HDG formulations for abstract one-field and two-field structures and prove their well-posedness. In order to handle non-conforming interfaces we simply take advantage of HDG built-in mortar structures. With split-type mortars and the approximation space of trace, a numerical flux can be derived via Godunov approach and be naturally employed without any additional treatment. As a consequence, the proposed formulations are parameter-free. We perform several numerical experiments for time-independent and linear PDEs including elliptic, hyperbolic, and mixed-type to verify the proposed unified $hp$-formulations and demonstrate the effectiveness of $hp$-adaptation. Two adaptivity criteria are considered: one is based on a simple and fast error indicator, while the other is rigorous but more expensive using an adjoint-based error estimate. The numerical results show that these two approaches are comparable in terms of convergence rate even for problems with strong gradients, discontinuities, or singularities.
△ Less
Submitted 7 April, 2023;
originally announced April 2023.
-
Structure-preserving finite-element schemes for the Euler-Poisson equations
Authors:
Matthias Maier,
John N. Shadid,
Ignacio Tomas
Abstract:
We discuss structure-preserving numerical discretizations for repulsive and attractive Euler-Poisson equations that find applications in fluid-plasma and self-gravitation modeling. The scheme is fully discrete and structure preserving in the sense that it maintains a discrete energy law, as well as hyperbolic invariant domain properties, such as positivity of the density and a minimum principle of…
▽ More
We discuss structure-preserving numerical discretizations for repulsive and attractive Euler-Poisson equations that find applications in fluid-plasma and self-gravitation modeling. The scheme is fully discrete and structure preserving in the sense that it maintains a discrete energy law, as well as hyperbolic invariant domain properties, such as positivity of the density and a minimum principle of the specific entropy. A detailed discussion of algorithmic details is given, as well as proofs of the claimed properties. We present computational experiments corroborating our analytical findings and demonstrating the computational capabilities of the scheme.
△ Less
Submitted 27 April, 2023; v1 submitted 16 July, 2022;
originally announced July 2022.
-
An adaptive scalable fully implicit algorithm based on stabilized finite element for reduced visco-resistive MHD
Authors:
Qi Tang,
Luis Chacon,
Tzanio V. Kolev,
John N. Shadid,
Xian-Zhu Tang
Abstract:
The magnetohydrodynamics (MHD) equations are continuum models used in the study of a wide range of plasma physics systems, including the evolution of complex plasma dynamics in tokamak disruptions. However, efficient numerical solution methods for MHD are extremely challenging due to disparate time and length scales, strong hyperbolic phenomena, and nonlinearity. Therefore the development of scala…
▽ More
The magnetohydrodynamics (MHD) equations are continuum models used in the study of a wide range of plasma physics systems, including the evolution of complex plasma dynamics in tokamak disruptions. However, efficient numerical solution methods for MHD are extremely challenging due to disparate time and length scales, strong hyperbolic phenomena, and nonlinearity. Therefore the development of scalable, implicit MHD algorithms and high-resolution adaptive mesh refinement strategies is of considerable importance. In this work, we develop a high-order stabilized finite-element algorithm for the reduced visco-resistive MHD equations based on the MFEM finite element library (mfem.org). The scheme is fully implicit, solved with the Jacobian-free Newton-Krylov (JFNK) method with a physics-based preconditioning strategy. Our preconditioning strategy is a generalization of the physics-based preconditioning methods in [Chacon, et al, JCP 2002] to adaptive, stabilized finite elements. Algebraic multigrid methods are used to invert sub-block operators to achieve scalability. A parallel adaptive mesh refinement scheme with dynamic load-balancing is implemented to efficiently resolve the multi-scale spatial features of the system. Our implementation uses the MFEM framework, which provides arbitrary-order polynomials and flexible adaptive conforming and non-conforming meshes capabilities. Results demonstrate the accuracy, efficiency, and scalability of the implicit scheme in the presence of large scale disparity. The potential of the AMR approach is demonstrated on an island coalescence problem in the high Lundquist-number regime ($\ge 10^7$) with the successful resolution of plasmoid instabilities and thin current sheets.
△ Less
Submitted 12 January, 2022; v1 submitted 1 June, 2021;
originally announced June 2021.
-
A Monolithic Algebraic Multigrid Framework for Multiphysics Applications with Examples from Resistive MHD
Authors:
Peter Ohm,
Tobias Wiesner,
Eric C. Cyr,
Jonathan J. Hu,
John N. Shadid,
Raymond S. Tuminaro
Abstract:
A multigrid framework is described for multiphysics applications. The framework allows one to construct, adapt, and tailor a monolithic multigrid methodology to different linear systems coming from discretized partial differential equations. The main idea centers on developing multigrid components in a blocked fashion where each block corresponds to separate sets of physical unknowns and equations…
▽ More
A multigrid framework is described for multiphysics applications. The framework allows one to construct, adapt, and tailor a monolithic multigrid methodology to different linear systems coming from discretized partial differential equations. The main idea centers on developing multigrid components in a blocked fashion where each block corresponds to separate sets of physical unknowns and equations within the larger discretization matrix. Once defined, these components are ultimately assembled into a monolithic multigrid solver for the entire system. We demonstrate the potential of the framework by applying it to representative linear solution sub-problems arising from resistive MHD.
△ Less
Submitted 22 March, 2021; v1 submitted 12 March, 2021;
originally announced March 2021.
-
A Multilevel Block Preconditioner for the HDG Trace System Applied to Incompressible Resistive MHD
Authors:
Sriramkrishnan Muralikrishnan,
Stephen Shannon,
Tan Bui-Thanh,
John N. Shadid
Abstract:
We present a scalable block preconditioning strategy for the trace system coming from the high-order hybridized discontinuous Galerkin (HDG) discretization of incompressible resistive magnetohydrodynamics (MHD). We construct the block preconditioner with a least squares commutator (BFBT) approximation for the inverse of the Schur complement that segregates out the pressure unknowns of the trace sy…
▽ More
We present a scalable block preconditioning strategy for the trace system coming from the high-order hybridized discontinuous Galerkin (HDG) discretization of incompressible resistive magnetohydrodynamics (MHD). We construct the block preconditioner with a least squares commutator (BFBT) approximation for the inverse of the Schur complement that segregates out the pressure unknowns of the trace system. The remaining velocity, magnetic field, and Lagrange multiplier unknowns form a coupled nodal unknown block (the upper block), for which a system algebraic multigrid (AMG) is used for the approximate inverse. The complexity of the MHD equations together with the algebraic nature of the statically condensed HDG trace system makes the choice of smoother in the system AMG part critical for the convergence and performance of the block preconditioner. Our numerical experiments show GMRES preconditioned by ILU(0) of overlap zero as a smoother inside system AMG performs best in terms of robustness, time per nonlinear iteration and memory requirements. With several transient test cases in 2D and 3D including the island coalescence problem at high Lundquist number we demonstrate the robustness and parallel scalability of the block preconditioner. Additionally for the upper block a preliminary study of an alternate nodal block system solver based on a multilevel approximate nested dissection is presented. On a 2D island coalescence problem the multilevel approximate nested dissection preconditioner shows better scalability with respect to mesh refinement than the system AMG, but is relatively less robust with respect to Lundquist number scaling.
△ Less
Submitted 14 December, 2020;
originally announced December 2020.
-
High Rayleigh number variational multiscale large eddy simulations of Rayleigh-Bénard Convection
Authors:
David Sondak,
Thomas M. Smith,
Roger P. Pawlowski,
Sidafa Conde,
John N. Shadid
Abstract:
The variational multiscale (VMS) formulation is used to develop residual-based VMS large eddy simulation (LES) models for Rayleigh-Bénard convection. The resulting model is a mixed model that incorporates the VMS model and an eddy viscosity model. The Wall-Adapting Local Eddy-viscosity (WALE) model is used as the eddy viscosity model in this work. The new LES models were implemented in the finite…
▽ More
The variational multiscale (VMS) formulation is used to develop residual-based VMS large eddy simulation (LES) models for Rayleigh-Bénard convection. The resulting model is a mixed model that incorporates the VMS model and an eddy viscosity model. The Wall-Adapting Local Eddy-viscosity (WALE) model is used as the eddy viscosity model in this work. The new LES models were implemented in the finite element code Drekar. Simulations are performed using continuous, piecewise linear finite elements. The simulations ranged from $Ra = 10^6$ to $Ra = 10^{14}$ and were conducted at $Pr = 1$ and $Pr = 7$. Two domains were considered: a two-dimensional domain of aspect ratio 2 with a fluid confined between two parallel plates and a three-dimensional cylinder of aspect ratio $1/4$. The Nusselt number from the VMS results is compared against three dimensional direct numerical simulations and experiments. In all cases, the VMS results are in good agreement with existing literature.
△ Less
Submitted 20 May, 2020;
originally announced May 2020.
-
An a posteriori error analysis for the equations of stationary incompressible magnetohydrodynamics
Authors:
J. H. Chaudhry,
A. E. Rappaport,
J. N. Shadid
Abstract:
Magnetohydrodynamics (MHD) is a continuum level model for conducting fluids subject to external magnetic fields, e.g. plasmas and liquid metals. The efficient and robust solution of the MHD system poses many challenges due to it's nonlinear, non self-adjoint, and highly coupled nature. In this paper, we develop a robust and accurate a posteriori error estimate for the numerical solution of the MHD…
▽ More
Magnetohydrodynamics (MHD) is a continuum level model for conducting fluids subject to external magnetic fields, e.g. plasmas and liquid metals. The efficient and robust solution of the MHD system poses many challenges due to it's nonlinear, non self-adjoint, and highly coupled nature. In this paper, we develop a robust and accurate a posteriori error estimate for the numerical solution of the MHD equations based on the exact penalty method. The error estimate also isolates particular contributions of error in a quantity of interest (QoI) to inform discretization choices to arrive at accurate solutions. The tools required for these estimates involve duality arguments and computable residuals.
△ Less
Submitted 17 December, 2020; v1 submitted 21 April, 2020;
originally announced April 2020.
-
On differentiable local bounds preserving stabilization for Euler equations
Authors:
Santiago Badia,
Jesús Bonilla,
Sibusiso Mabuza,
John N. Shadid
Abstract:
This work presents the design of nonlinear stabilization techniques for the finite element discretization of Euler equations in both steady and transient form. Implicit time integration is used in the case of the transient form. A differentiable local bounds preserving method has been developed, which combines a Rusanov artificial diffusion operator and a differentiable shock detector. Nonlinear s…
▽ More
This work presents the design of nonlinear stabilization techniques for the finite element discretization of Euler equations in both steady and transient form. Implicit time integration is used in the case of the transient form. A differentiable local bounds preserving method has been developed, which combines a Rusanov artificial diffusion operator and a differentiable shock detector. Nonlinear stabilization schemes are usually stiff and highly nonlinear. This issue is mitigated by the differentiability properties of the proposed method. Moreover, in order to further improve the nonlinear convergence, we also propose a continuation method for a subset of the stabilization parameters. The resulting method has been successfully applied to steady and transient problems with complex shock patterns. Numerical experiments show that it is able to provide sharp and well resolved shocks. The importance of the differentiability is assessed by comparing the new scheme with its non-differentiable counterpart. Numerical experiments suggest that, for up to moderate nonlinear tolerances, the method exhibits improved robustness and nonlinear convergence behavior for steady problems. In the case of transient problem, we also observe a reduction in the computational cost.
△ Less
Submitted 7 July, 2020; v1 submitted 15 November, 2019;
originally announced November 2019.
-
A Multilevel Approach for Trace System in HDG Discretizations
Authors:
Sriramkrishnan Muralikrishnan,
Tan Bui-Thanh,
John N. Shadid
Abstract:
We propose a multilevel approach for trace systems resulting from hybridized discontinuous Galerkin (HDG) methods. The key is to blend ideas from nested dissection, domain decomposition, and high-order characteristic of HDG discretizations. Specifically, we first create a coarse solver by eliminating and/or limiting the front growth in nested dissection. This is accomplished by projecting the trac…
▽ More
We propose a multilevel approach for trace systems resulting from hybridized discontinuous Galerkin (HDG) methods. The key is to blend ideas from nested dissection, domain decomposition, and high-order characteristic of HDG discretizations. Specifically, we first create a coarse solver by eliminating and/or limiting the front growth in nested dissection. This is accomplished by projecting the trace data into a sequence of same or high-order polynomials on a set of increasingly $h-$coarser edges/faces. We then combine the coarse solver with a block-Jacobi fine scale solver to form a two-level solver/preconditioner. Numerical experiments indicate that the performance of the resulting two-level solver/preconditioner depends only on the smoothness of the solution and is independent of the nature of the PDE under consideration. While the proposed algorithms are developed within the HDG framework, they are applicable to other hybrid(ized) high-order finite element methods. Moreover, we show that our multilevel algorithms can be interpreted as a multigrid method with specific intergrid transfer and smoothing operators. With several numerical examples from Poisson, pure transport, and convection-diffusion equations we demonstrate the robustness and scalability of the algorithms.
△ Less
Submitted 26 March, 2019;
originally announced March 2019.
-
Embedded error estimation and adaptive step-size control for optimal explicit strong stability preserving Runge--Kutta methods
Authors:
Sidafa Conde,
Imre Fekete,
John N. Shadid
Abstract:
We construct a family of embedded pairs for optimal strong stability preserving explicit Runge-Kutta methods of order $2 \leq p \leq 4$ to be used to obtain numerical solution of spatially discretized hyperbolic PDEs. In this construction, the goals include non-defective methods, large region of absolute stability, and optimal error measurement as defined in [5,19]. The new family of embedded pair…
▽ More
We construct a family of embedded pairs for optimal strong stability preserving explicit Runge-Kutta methods of order $2 \leq p \leq 4$ to be used to obtain numerical solution of spatially discretized hyperbolic PDEs. In this construction, the goals include non-defective methods, large region of absolute stability, and optimal error measurement as defined in [5,19]. The new family of embedded pairs offer the ability for strong stability preserving (SSP) methods to adapt by varying the step-size based on the local error estimation while maintaining their inherent nonlinear stability properties. Through several numerical experiments, we assess the overall effectiveness in terms of precision versus work while also taking into consideration accuracy and stability.
△ Less
Submitted 22 June, 2018;
originally announced June 2018.
-
Analysis of an HDG Method for Linearized Incompressible Resistive MHD Equations
Authors:
Jeonghun J. Lee,
Stephen Shannon,
Tan Bui-Thanh,
John N. Shadid
Abstract:
We present a hybridized discontinuous Galerkin (HDG) method for stationary linearized incompressible magnetohydrodynamics (MHD) equations. At the heart of the paper is the introduction of an HDG flux of the dual saddle-point form of the MHD equations that facilitates the hybridization of discontinuous Galerkin (DG) method. We carry out the $\textit{a priori}$ error estimates for the proposed HDG m…
▽ More
We present a hybridized discontinuous Galerkin (HDG) method for stationary linearized incompressible magnetohydrodynamics (MHD) equations. At the heart of the paper is the introduction of an HDG flux of the dual saddle-point form of the MHD equations that facilitates the hybridization of discontinuous Galerkin (DG) method. We carry out the $\textit{a priori}$ error estimates for the proposed HDG method on simplicial meshes in both two- and three-dimensions. The analysis provides optimal convergence for the fluid velocity and the magnetic variables, and quasi-optimal convergence for the remaining quantities. Numerical examples are presented to verify the theoretical findings.
△ Less
Submitted 12 January, 2019; v1 submitted 16 February, 2017;
originally announced February 2017.
-
Implicit and Implicit-Explicit Strong Stability Preserving Runge-Kutta Methods with High Linear Order
Authors:
Sidafa Conde,
Sigal Gottlieb,
Zachary J. Grant,
John N. Shadid
Abstract:
When evolving in time the solution of a hyperbolic partial differential equation, it is often desirable to use high order strong stability preserving (SSP) time discretizations. These time discretizations preserve the monotonicity properties satisfied by the spatial discretization when coupled with the first order forward Euler, under a certain time-step restriction. While the allowable time-step…
▽ More
When evolving in time the solution of a hyperbolic partial differential equation, it is often desirable to use high order strong stability preserving (SSP) time discretizations. These time discretizations preserve the monotonicity properties satisfied by the spatial discretization when coupled with the first order forward Euler, under a certain time-step restriction. While the allowable time-step depends on both the spatial and temporal discretizations, the contribution of the temporal discretization can be isolated by taking the ratio of the allowable time-step of the high order method to the forward Euler time-step. This ratio is called the strong stability coefficient. The search for high order strong stability time-stepping methods with high order and large allowable time-step had been an active area of research. It is known that implicit SSP Runge-Kutta methods exist only up to sixth order. However, if we restrict ourselves to solving only linear autonomous problems, the order conditions simplify and we can find implicit SSP Runge-Kutta methods of any linear order. In the current work we aim to find very high linear order implicit SSP Runge-Kutta methods that are optimal in terms of allowable time-step. Next, we formulate an optimization problem for implicit-explicit (IMEX) SSP Runge-Kutta methods and find implicit methods with large linear stability regions that pair with known explicit SSP Runge-Kutta methods of orders plin=3,4,6 as well as optimized IMEX SSP Runge-Kutta pairs that have high linear order and nonlinear orders p=2,3,4. These methods are then tested on sample problems to verify order of convergence and to demonstrate the sharpness of the SSP coefficient and the typical behavior of these methods on test problems.
△ Less
Submitted 15 February, 2017;
originally announced February 2017.
-
Error Estimation for Multi-Stage Runge-Kutta IMEX Schemes
Authors:
Jehanzeb H. Chaudhry,
J. B. Collins,
John N. Shadid
Abstract:
Implicit-Explicit (IMEX) schemes are widely used for time integration methods for approximating solutions to a large class of problems. In this work, we develop accurate a posteriori error estimates of a quantity of interest for approximations obtained from multi-stage IMEX schemes. This is done by first defining a finite element method that is nodally equivalent to an IMEX scheme, then using typi…
▽ More
Implicit-Explicit (IMEX) schemes are widely used for time integration methods for approximating solutions to a large class of problems. In this work, we develop accurate a posteriori error estimates of a quantity of interest for approximations obtained from multi-stage IMEX schemes. This is done by first defining a finite element method that is nodally equivalent to an IMEX scheme, then using typical methods for adjoint-based error estimation. The use of a nodally equivalent finite element method allows a decomposition of the error into multiple components, each describing the effect of a different portion of the method on the total error in a quantity of interest.
△ Less
Submitted 18 October, 2016; v1 submitted 28 September, 2015;
originally announced September 2015.
-
A new class of finite element variational multiscale turbulence models for incompressible magnetohydrodynamics
Authors:
David Sondak,
John N. Shadid,
Assad A. Oberai,
Roger P. Pawlowski,
Eric C. Cyr,
Tom M. Smith
Abstract:
New large eddy simulation (LES) turbulence models for incompressible magnetohydrodynamics (MHD) derived from the variational multiscale (VMS) formulation for finite element simulations are introduced. The new models include the variational multiscale formulation, a residual-based eddy viscosity model, and a mixed model that combines both of these component models. Each model contains terms that ar…
▽ More
New large eddy simulation (LES) turbulence models for incompressible magnetohydrodynamics (MHD) derived from the variational multiscale (VMS) formulation for finite element simulations are introduced. The new models include the variational multiscale formulation, a residual-based eddy viscosity model, and a mixed model that combines both of these component models. Each model contains terms that are proportional to the residual of the incompressible MHD equations and is therefore numerically consistent. Moreover, each model is also dynamic, in that its effect vanishes when this residual is small. The new models are tested on the decaying MHD Taylor Green vortex at low and high Reynolds numbers. The evaluation of the models is based on comparisons with available data from direct numerical simulations (DNS) of the time evolution of energies as well as energy spectra at various discrete times. A numerical study, on a sequence of meshes, is presented that demonstrates that the large eddy simulation approaches the DNS solution for these quantities with spatial mesh refinement.
△ Less
Submitted 2 December, 2014;
originally announced December 2014.