Newton type methods for solving a Hasegawa–Mima plasma model

• Published in: Applied mathematics and computation 2023-12, Vol.459, p.128141, Article 128141
• Main Authors: Moufawad, Sophie M., Nassif, Nabil R.
• Format: Article
• Language: English
• Subjects: Finite-element method; Hasegawa–Mima; Implicit finite-differences; Newton-type methods; Periodic Sobolev spaces; Petrov–Galerkin approximations
• ISSN: 0096-3003; 1873-5649
• DOI: 10.1016/j.amc.2023.128141

•Hasegawa–Mima equation as a system of two linear PDE’s, elliptic and hyperbolic.•Weak solutions to the Hasegawa–Mima equation.•New Newton-type methods for solving finite element in space, finite difference in time discretization.•Theoretical and numerical validations of the proposed methods. In Karakazian and Nassif (2019), the non-linear space-time Hasegawa–Mima plasma equation is formulated as a coupled system of two linear PDEs, a solution of which is a pair (u,w), with w=(I−Δ)u. The first equation is of hyperbolic type and the second of elliptic type. Variational frames for obtaining weak solutions to the initial value Hasegawa–Mima problem with periodic boundary conditions were also derived. In a more recent work Karakazian et al. (2022), a numerical approach consisting of a finite element space-domain combined with an Euler-implicit time scheme was used to discretize the coupled variational Hasegawa–Mima model. A semi-linear version of this implicit nonlinear scheme was tested for several types of initial conditions. This semi-linear scheme proved to lack efficiency for long time, which necessitates imposing a cap on the magnitude of the solution. To circumvent this difficulty, in this paper, we use Newton-type methods (Newton, Chord and an introduced Modified Newton method) to solve numerically the fully-implicit non-linear scheme. Testing these methods in FreeFEM++ indicates significant improvements as no cap needs to be imposed for long time. In the sequel, we demonstrate the validity of these methods by proving several results, in particular the convergence of the implemented methods.