A Priori Error Analysis of a Discontinuous Galerkin Scheme for the Magnetic Induction Equation

2020 ◽  
Vol 20 (1) ◽  
pp. 121-140 ◽  
Author(s):  
Tanmay Sarkar

AbstractWe perform the error analysis of a stabilized discontinuous Galerkin scheme for the initial boundary value problem associated with the magnetic induction equations using standard discontinuous Lagrange basis functions. In order to obtain the quasi-optimal convergence incorporating second-order Runge–Kutta schemes for time discretization, we need a strengthened {4/3}-CFL condition ({\Delta t\sim h^{4/3}}). To overcome this unusual restriction on the CFL condition, we consider the explicit third-order Runge–Kutta scheme for time discretization. We demonstrate the error estimates in {L^{2}}-sense and obtain quasi-optimal convergence for smooth solution in space and time for piecewise polynomials with any degree {l\geq 1} under the standard CFL condition.

Author(s):  
Manas Vijay Upadhyay ◽  
Jérémy Bleyer

Abstract A time-explicit Runge-Kutta discontinuous Galerkin (RKDG) finite element scheme is proposed to solve the dislocation transport initial boundary value problem in 3D. The dislocation density transport equation, which lies at the core of this problem, is a first-order unsteady-state advection-reaction-type hyperbolic partial differential equation; the DG approach is well suited to solve such equations that lack any diffusion terms. The development of the RKDG scheme follows the method of lines approach. First, a space semi-discretization is performed using the DG approach with upwinding to obtain a system of ordinary differential equations in time. Then, time discretization is performed using explicit RK schemes to solve this system. The 3D numerical implementation of the RKDG scheme is performed for the first-order (forward Euler), second-order and third-order RK methods using the strong stability preserving approach. These implementations provide (quasi-)optimal convergence rates for smooth solutions. A slope limiter is used to prevent spurious Gibbs oscillations arising from high-order space approximations (polynomial degree ≥ 1) of rough solutions. A parametric study is performed to understand the influence of key parameters of the RKDG scheme on the stability of the solution predicted during a screw dislocation transport simulation. Then, annihilation of two oppositely signed screw dislocations and the expansion of a polygonal dislocation loop are simulated. The RKDG scheme is able to resolve the shock generated during dislocation annihilation without any spurious oscillations and predict the prismatic loop expansion with very low numerical diffusion. These results demonstrate the robustness of the scheme.


Sign in / Sign up

Export Citation Format

Share Document