On a fully-implicit VMS-stabilized FE formulation for low Mach number compressible resistive MHD with application to MCF

被引：1

作者：

Bonillaa, J. ^{[1
]}

Shadidb, J. N. ^{[2
,3
]}

Tanga, X. -Z. ^{[1
]}

Crockattb, M. M. ^{[2
]}

Ohmd, P. ^{[4
]}

Phillipsb, E. G. ^{[2
]}

Pawlowskib, R. P. ^{[2
]}

Condeb, S. ^{[2
]}

Beznosova, O. ^{[1
]}

机构：

[1] Los Alamos Natl Lab, Los Alamos, NM 87545 USA

[2] Sandia Natl Labs, Albuquerque, NM USA

[3] Univ New Mexico, Albuquerque, NM USA

[4] RIKEN, Ctr Computat Sci, Wako, Japan

来源：

COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING | 2023年 / 417卷

关键词：

Compressible resistive MHD; VMS FE; Tokamak VDE; Fully-implicit Newton-Krylov; Block preconditioning; AMG; FINITE-ELEMENT APPROXIMATION; HEAT-TRANSPORT; NAVIER-STOKES; MAGNETOHYDRODYNAMICS; PRECONDITIONERS; ALGORITHM; SYSTEMS; FLOW;

D O I：

10.1016/j.cma.2023.116359

中图分类号：

T [工业技术];

学科分类号：

08 ;

摘要：

This study presents the development and evaluation of a fully-implicit variational multiscale (VMS) stabilized unstructured finite element (FE) formulation for compressible magnetohydrodynamics (MHD) model, at low Mach number regime. The model describes the dynamics of a compressible conducting fluid in the low Mach number limit in the presence of electromagnetic fields and can be used to study aspects of astrophysical phenomena, important science and technology applications, and basic plasma physics phenomena. The specific applications that motivate this study are macroscopic simulations of the longer time-scale stability and disruptions of magnetic confinement fusion (MCF) devices, specifically the ITER tokamak. The discussion considers the development of the VMS FE representation, the structure of the stabilizing terms that deal with significant convective flows, the stabilization of the nearly incompressible response of the fluid flow, and the stabilization of the constraint that enforces the solenoidal involution on the magnetic field. The nonlinear discretized system is solved with scalable preconditioned Newton-Krylov iterative methods, which employs a multiphysics block preconditioning method based on approximate block factorizations and Schur complements. The study presents an evaluation of the VMS method on a 2D cartesian tearing mode instability, and illustrates the scalability of the solvers on MCF relevant problems. A set of results are also presented for longer time-scale stability and disruptions for the ITER tokamak. These include a vertical displacement event (VDE), and a (1,1) internal kink mode. The formulation is demonstrated to be scalable and also reasonably robust with respect to the Lundquist number scaling. (c) 2023 Elsevier B.V. All rights reserved.

引用

页数：23

共 61 条

[1] Unconditionally stable operator splitting algorithms for the incompressible magnetohydrodynamics system discretized by a stabilized finite element formulation based on projections [J].