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

by   Zhichao Peng, et al.

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.



page 14

page 15

page 18

page 20


An adaptive stabilized finite element method based on residual minimization

We devise and analyze a new adaptive stabilized finite element method. W...

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

We introduce an unconditionally stable finite element (FE) method, the a...

A subcell-enriched Galerkin method for advection problems

In this work, we introduce a generalization of the enriched Galerkin (EG...

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...

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

This work is focused on the extension and assessment of the monotonicity...

A class of adaptive multiresolution ultra-weak discontinuous Galerkin methods for some nonlinear dispersive wave equations

In this paper, we propose a class of adaptive multiresolution (also call...

An Angular Multigrid Preconditioner for the Radiation Transport Equation with Forward-Peaked Scatter

In a previous paper (Lathouwers and Perkó, 2019) we have developed an ef...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.