-
An adaptive Newton-based free-boundary Grad-Shafranov solver
Authors:
Daniel A. Serino,
Qi Tang,
Xian-Zhu Tang,
Tzanio V. Kolev,
Konstantin Lipnikov
Abstract:
Equilibriums in magnetic confinement devices result from force balancing between the Lorentz force and the plasma pressure gradient. In an axisymmetric configuration like a tokamak, such an equilibrium is described by an elliptic equation for the poloidal magnetic flux, commonly known as the Grad--Shafranov equation. It is challenging to develop a scalable and accurate free-boundary Grad--Shafrano…
▽ More
Equilibriums in magnetic confinement devices result from force balancing between the Lorentz force and the plasma pressure gradient. In an axisymmetric configuration like a tokamak, such an equilibrium is described by an elliptic equation for the poloidal magnetic flux, commonly known as the Grad--Shafranov equation. It is challenging to develop a scalable and accurate free-boundary Grad--Shafranov solver, since it is a fully nonlinear optimization problem that simultaneouly solves for the magnetic field coil current outside the plasma to control the plasma shape. In this work, we develop a Newton-based free-boundary Grad--Shafranov solver using adaptive finite elements and preconditioning strategies. The free-boundary interaction leads to the evaluation of a domain-dependent nonlinear form of which its contribution to the Jacobian matrix is achieved through shape calculus. The optimization problem aims to minimize the distance between the plasma boundary and specified control points while satisfying two non-trivial constraints, which correspond to the nonlinear finite element discretization of the Grad--Shafranov equation and a constraint on the total plasma current involving a nonlocal coupling term. The linear system is solved by a block factorization, and AMG is called for subblock elliptic operators. The unique contributions of this work include the treatment of a global constraint, preconditioning strategies, nonlocal reformulation, and the implementation of adaptive finite elements. It is found that the resulting Newton solver is robust, successfully reducing the nonlinear residual to 1e-6 and lower in a small handful of iterations while addressing the challenging case to find a Taylor state equilibrium where conventional Picard-based solvers fail to converge.
△ Less
Submitted 3 July, 2024;
originally announced July 2024.
-
A mimetic finite difference based quasi-static magnetohydrodynamic solver for force-free plasmas in tokamak disruptions
Authors:
Zakariae Jorti,
Qi Tang,
Konstantin Lipnikov,
Xian-Zhu Tang
Abstract:
Force-free plasmas are a good approximation where the plasma pressure is tiny compared with the magnetic pressure, which is the case during the cold vertical displacement event (VDE) of a major disruption in a tokamak. On time scales long compared with the transit time of Alfven waves, the evolution of a force-free plasma is most efficiently described by the quasi-static magnetohydrodynamic (MHD)…
▽ More
Force-free plasmas are a good approximation where the plasma pressure is tiny compared with the magnetic pressure, which is the case during the cold vertical displacement event (VDE) of a major disruption in a tokamak. On time scales long compared with the transit time of Alfven waves, the evolution of a force-free plasma is most efficiently described by the quasi-static magnetohydrodynamic (MHD) model, which ignores the plasma inertia. Here we consider a regularized quasi-static MHD model for force-free plasmas in tokamak disruptions and propose a mimetic finite difference (MFD) algorithm. The full geometry of an ITER-like tokamak reactor is treated, with a blanket module region, a vacuum vessel region, and the plasma region. Specifically, we develop a parallel, fully implicit, and scalable MFD solver based on PETSc and its DMStag data structure for the discretization of the five-field quasi-static perpendicular plasma dynamics model on a 3D structured mesh. The MFD spatial discretization is coupled with a fully implicit DIRK scheme. The algorithm exactly preserves the divergence-free condition of the magnetic field under the resistive Ohm's law. The preconditioner employed is a four-level fieldsplit preconditioner, which is created by combining separate preconditioners for individual fields, that calls multigrid or direct solvers for sub-blocks or exact factorization on the separate fields. The numerical results confirm the divergence-free constraint is strongly satisfied and demonstrate the performance of the fieldsplit preconditioner and overall algorithm. The simulation of ITER VDE cases over the actual plasma current diffusion time is also presented.
△ Less
Submitted 23 March, 2023; v1 submitted 14 March, 2023;
originally announced March 2023.
-
Flow and Transport in Three-Dimensional Discrete Fracture Matrix Models using Mimetic Finite Difference on a Conforming Multi-Dimensional Mesh
Authors:
Jeffrey D. Hyman,
Matthew R. Sweeney,
Carl W. Gable,
Daniil Svyatsky,
Konstantin Lipnikov,
J. David Moulton
Abstract:
We present a comprehensive workflow to simulate single-phase flow and transport in fractured porous media using the discrete fracture matrix approach. The workflow has three primary parts: (1) a method for conforming mesh generation of and around a three-dimensional fracture network, (2) the discretization of the governing equations using a second-order mimetic finite difference method, and (3) im…
▽ More
We present a comprehensive workflow to simulate single-phase flow and transport in fractured porous media using the discrete fracture matrix approach. The workflow has three primary parts: (1) a method for conforming mesh generation of and around a three-dimensional fracture network, (2) the discretization of the governing equations using a second-order mimetic finite difference method, and (3) implementation of numerical methods for high-performance computing environments. A method to create a conforming Delaunay tetrahedralization of the volume surrounding the fracture network, where the triangular cells of the fracture mesh are faces in the volume mesh, that addresses pathological cases which commonly arise and degrade mesh quality is also provided. Our open-source subsurface simulator uses a hierarchy of process kernels (one kernel per physical process) that allows for both strong and weak coupling of the fracture and matrix domains. We provide verification tests based on analytic solutions for flow and transport, as well as numerical convergence. We also provide multiple expositions of the method in complex fracture networks. In the first example, we demonstrate that the method is robust by considering two scenarios where the fracture network acts as a barrier to flow, as the primary pathway, or offers the same resistance as the surrounding matrix. In the second test, flow and transport through a three-dimensional stochastically generated network containing 257 fractures is presented.
△ Less
Submitted 20 December, 2021;
originally announced December 2021.
-
Verification benchmarks for single-phase flow in three-dimensional fractured porous media
Authors:
Inga Berre,
Wietse M. Boon,
Bernd Flemisch,
Alessio Fumagalli,
Dennis Gläser,
Eirik Keilegavlen,
Anna Scotti,
Ivar Stefansson,
Alexandru Tatomir,
Konstantin Brenner,
Samuel Burbulla,
Philippe Devloo,
Omar Duran,
Marco Favino,
Julian Hennicker,
I-Hsien Lee,
Konstantin Lipnikov,
Roland Masson,
Klaus Mosthaf,
Maria Giuseppina Chiara Nestola,
Chuen-Fa Ni,
Kirill Nikitin,
Philipp Schädle,
Daniil Svyatskiy,
Ruslan Yanbarisov
, et al. (1 additional authors not shown)
Abstract:
Flow in fractured porous media occurs in the earth's subsurface, in biological tissues, and in man-made materials. Fractures have a dominating influence on flow processes, and the last decade has seen an extensive development of models and numerical methods that explicitly account for their presence. To support these developments, we present a portfolio of four benchmark cases for single-phase flo…
▽ More
Flow in fractured porous media occurs in the earth's subsurface, in biological tissues, and in man-made materials. Fractures have a dominating influence on flow processes, and the last decade has seen an extensive development of models and numerical methods that explicitly account for their presence. To support these developments, we present a portfolio of four benchmark cases for single-phase flow in three-dimensional fractured porous media. The cases are specifically designed to test the methods' capabilities in handling various complexities common to the geometrical structures of fracture networks. Based on an open call for participation, results obtained with 17 numerical methods were collected. This paper presents the underlying mathematical model, an overview of the features of the participating numerical methods, and their performance in solving the benchmark cases.
△ Less
Submitted 17 February, 2020;
originally announced February 2020.
-
Multidimensional staggered grid residual distribution scheme for Lagrangian hydrodynamics
Authors:
R. Abgrall,
K. Lipnikov,
N. Morgan,
S. Tokareva
Abstract:
We present the second-order multidimensional Staggered Grid Hydrodynamics Residual Distribution (SGH RD) scheme for Lagrangian hydrodynamics. The SGH RD scheme is based on the staggered finite element discretizations as in [Dobrev et al., SISC, 2012]. However, the advantage of the residual formulation over classical FEM approaches consists in the natural mass matrix diagonalization which allows on…
▽ More
We present the second-order multidimensional Staggered Grid Hydrodynamics Residual Distribution (SGH RD) scheme for Lagrangian hydrodynamics. The SGH RD scheme is based on the staggered finite element discretizations as in [Dobrev et al., SISC, 2012]. However, the advantage of the residual formulation over classical FEM approaches consists in the natural mass matrix diagonalization which allows one to avoid the solution of the linear system with the global sparse mass matrix while retaining the desired order of accuracy. This is achieved by using Bernstein polynomials as finite element shape functions and coupling the space discretization with the deferred correction type timestepping method. Moreover, it can be shown that for the Lagrangian formulation written in non-conservative form, our residual distribution scheme ensures the exact conservation of the mass, momentum and total energy. In this paper we also discuss construction of numerical viscosity approximations for the SGH RD scheme allowing to reduce the dissipation of the numerical solution. Thanks to the generic formulation of the staggered grid residual distribution scheme, it can be directly applied to both single- and multimaterial and multiphase models. Finally, we demonstrate computational results obtained with the proposed residual distribution scheme for several challenging test problems.
△ Less
Submitted 31 October, 2018;
originally announced November 2018.
-
Convergence analysis of the mimetic finite difference method for elliptic problems with staggered discretization of the diffusion coefficients
Authors:
G. Manzini,
K. Lipnikov,
J. D. Moulton,
M. Shashkov
Abstract:
We study the convergence of the new family of mimetic finite difference schemes for linear diffusion problems recently proposed in [38]. In contrast to the conventional approach, the diffusion coefficient enters both the primary mimetic operator, i.e., the discrete divergence, and the inner product in the space of gradients. The diffusion coefficient is therefore evaluated on different mesh locati…
▽ More
We study the convergence of the new family of mimetic finite difference schemes for linear diffusion problems recently proposed in [38]. In contrast to the conventional approach, the diffusion coefficient enters both the primary mimetic operator, i.e., the discrete divergence, and the inner product in the space of gradients. The diffusion coefficient is therefore evaluated on different mesh locations, i.e., inside mesh cells and on mesh faces. Such a staggered discretization may provide the exibility necessary for future development of efficient numerical schemes for nonlinear problems, especially for problems with degenerate coefficients. These new mimetic schemes preserve symmetry and positive-definiteness of the continuum problem, which allow us to use efficient algebraic solvers such as the preconditioned Conjugate Gradient method. We show that these schemes are inf-sup stable and establish a priori error estimates for the approximation of the scalar and vector solution fields. Numerical examples confirm the convergence analysis and the effectiveness of the method in providing accurate approximations.
△ Less
Submitted 6 December, 2016;
originally announced December 2016.
-
Discretization of mixed formulations of elliptic problems on polyhedral meshes
Authors:
Konstantin Lipnikov,
Gianmarco Manzini
Abstract:
We review basic design principles underpinning the construction of mimetic finite difference and a few finite volume and finite element schemes for mixed formulations of elliptic problems. For a class of low-order mixed-hybrid schemes, we show connections between these principles and prove that the consistency and stability conditions must lead to a member of the mimetic family of schemes regardle…
▽ More
We review basic design principles underpinning the construction of mimetic finite difference and a few finite volume and finite element schemes for mixed formulations of elliptic problems. For a class of low-order mixed-hybrid schemes, we show connections between these principles and prove that the consistency and stability conditions must lead to a member of the mimetic family of schemes regardless of the selected discretization framework. Finally, we give two examples of using flexibility of the mimetic framework: derivation of higher-order schemes and convergent schemes for nonlinear problems with small diffusion coefficients.
△ Less
Submitted 18 October, 2016;
originally announced October 2016.
-
The mimetic finite difference method for the Landau-Lifshitz equation
Authors:
Eugenia Kim,
Konstantin Lipnikov
Abstract:
The Landau-Lifshitz equation describes the dynamics of the magnetization inside ferromagnetic materials. This equation is highly nonlinear and has a non-convex constraint (the magnitude of the magnetization is constant) which pose interesting challenges in developing numerical methods. We develop and analyze explicit and implicit mimetic finite difference schemes for this equation. These schemes w…
▽ More
The Landau-Lifshitz equation describes the dynamics of the magnetization inside ferromagnetic materials. This equation is highly nonlinear and has a non-convex constraint (the magnitude of the magnetization is constant) which pose interesting challenges in developing numerical methods. We develop and analyze explicit and implicit mimetic finite difference schemes for this equation. These schemes work on general polytopal meshes which provide enormous flexibility to model magnetic devices with various shapes. A projection on the unit sphere is used to preserve the magnitude of the magnetization. We also provide a proof that shows the exchange energy is decreasing in certain conditions. The developed schemes are tested on general meshes that include distorted and randomized meshes. The numerical experiments include a test proposed by the National Institute of Standard and Technology and a test showing formation of domain wall structures in a thin film.
△ Less
Submitted 25 August, 2016;
originally announced August 2016.
-
The nonconforming virtual element method
Authors:
B. Ayuso de Dios,
K. Lipnikov,
G. Manzini
Abstract:
We introduce the nonconforming Virtual Element Method (VEM) for the approximation of second order elliptic problems.
We present the construction of the new element in two and three dimensions, highlighting the main differences with the conforming VEM and the classical nonconforming finite element methods.
We provide the error analysis and establish the equivalence with a family of mimetic fini…
▽ More
We introduce the nonconforming Virtual Element Method (VEM) for the approximation of second order elliptic problems.
We present the construction of the new element in two and three dimensions, highlighting the main differences with the conforming VEM and the classical nonconforming finite element methods.
We provide the error analysis and establish the equivalence with a family of mimetic finite difference methods.
△ Less
Submitted 16 May, 2014; v1 submitted 15 May, 2014;
originally announced May 2014.