An adaptive discontinuous Petrov-Galerkin method for the Grad-Shafranov equation

01/13/2020
by   Zhichao Peng, et al.
0

In this work, we propose and develop an arbitrary-order adaptive discontinuous Petrov-Galerkin (DPG) method for the nonlinear Grad-Shafranov equation. An ultraweak formulation of the DPG scheme for the equation is given based on a minimal residual method. The DPG scheme has the advantage of providing more accurate gradients compared to conventional finite element methods, which is desired for numerical solutions to the Grad-Shafranov equation. The numerical scheme is augmented with an adaptive mesh refinement approach, and a criterion based on the residual norm in the minimal residual method is developed to achieve dynamic refinement. Nonlinear solvers for the resulting system are explored and a Picard iteration with Anderson acceleration is found to be efficient to solve the system. Finally, the proposed algorithm is implemented in parallel on MFEM using a domain-decomposition approach, and our implementation is general, supporting arbitrary order of accuracy and general meshes. Numerical results are presented to demonstrate the efficiency and accuracy of the proposed algorithm.

READ FULL TEXT

page 14

page 15

page 18

page 20

research
07/29/2019

An adaptive stabilized finite element method based on residual minimization

We devise and analyze a new adaptive stabilized finite element method. W...
research
06/11/2020

An Unconditionally Stable Space-Time FE Method for the Korteweg-de Vries Equation

We introduce an unconditionally stable finite element (FE) method, the a...
research
03/13/2023

Adaptive mesh refinement for the Landau-Lifshitz-Gilbert equation

We propose a new adaptive algorithm for the approximation of the Landau-...
research
12/01/2020

Minimal residual Hermitian and skew-Hermitian splitting iteration method for the continuous Sylvester equation

By applying the minimal residual technique to the Hermitian and skew-Her...
research
12/24/2019

Monotonicity-preserving finite element schemes with adaptive mesh refinement for hyperbolic problems

This work is focused on the extension and assessment of the monotonicity...
research
06/04/2020

A Memory-efficient Implementation of Perfectly Matched Layer with Smoothly-varying Coefficients in Discontinuous Galerkin Time-Domain Method

Wrapping a computation domain with a perfectly matched layer (PML) is on...
research
11/25/2021

An efficient IMEX-DG solver for the compressible Navier-Stokes equations with a general equation of state

We propose an efficient, accurate and robust IMEX solver for the compres...

Please sign up or login with your details

Forgot password? Click here to reset