Robust augmented Block Preconditioning for two-phase Incompressible Inductionless MHD
Abstract
In this paper, we develop a robust augmented block preconditioner for three-dimensional, two-phase, incompressible inductionless MHD equations. We employ a high-order monolithic and charge-conservative arbitrary Lagrangian--Eulerian (ALE) finite element method as base scheme, for which the solution variables are current density, velocity, electric potential and pressure. Grad-div stabilization is added in both the current density equation and discrete momentum equation. Backward Euler method is used for the time discretization. At each time step, fully coupled solution strategy is adopted. This leads to large-scale, ill-conditioned, nonsymmetric, and indefinite dual saddle-point system. To solve this system efficiently, we construct an augmented block preconditioner within the constraint preconditioning framework. The ultimate preconditioner combines the augmented Lagrangian technique with approximate block factorization for current density sub-system and fluid sub-system. The fully coupled four-field system is then solved using preconditioned FGMRES iteration. Numerical experiments confirm the high-order accuracy and discrete charge conservation of the ALE-FEM. Drop of heavy fluid in light fluid is taken as benchmark to demonstrate the robustness of the preconditioner with respect to the discontinuity of physical parameters and the mesh size.
Related
- High-order Energy-stable and Charge-conservative Lagrangian FEM for 3D Incompressible Inductionless MHD equations with Variable Density
- A unified variational framework for augmented Lagrangian preconditioning of saddle point systems: applications to mixed-dimensional and fictitious domain methods
- A parameterized block LU preconditioning for double saddle point problems
- Block Preconditioning for All-at-Once Variable-Coefficient Fractional Evolution Equations via the GLT Analysis
- Dobb: Shard-Local Preconditioning in Block-Diagonal Eigenbases
BibTeX
@misc{song2026robust,
title = {{Robust augmented Block Preconditioning for two-phase Incompressible Inductionless MHD}},
author = {Bijie Song and Lingxiao Li and Jie Zhang},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.33385},
eprint = {2610.33385},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.33385},
}
RIS
TY - UNPB
TI - Robust augmented Block Preconditioning for two-phase Incompressible Inductionless MHD
AU - Bijie Song
AU - Lingxiao Li
AU - Jie Zhang
PY - 2026
DA - 2026/10/24
PB - arXiv
DO - 10.48550/arxiv.2610.33385
UR - https://arxiv.org/abs/2610.33385
AB - In this paper, we develop a robust augmented block preconditioner for three-dimensional, two-phase, incompressible inductionless MHD equations. We employ a high-order monolithic and charge-conservative arbitrary Lagrangian--Eulerian (ALE) finite element method as base scheme, for which the solution variables are current density, velocity, electric potential and pressure. Grad-div stabilization is added in both the current density equation and discrete momentum equation. Backward Euler method is used for the time discretization. At each time step, fully coupled solution strategy is adopted. This leads to large-scale, ill-conditioned, nonsymmetric, and indefinite dual saddle-point system. To solve this system efficiently, we construct an augmented block preconditioner within the constraint preconditioning framework. The ultimate preconditioner combines the augmented Lagrangian technique with approximate block factorization for current density sub-system and fluid sub-system. The fully coupled four-field system is then solved using preconditioned FGMRES iteration. Numerical experiments confirm the high-order accuracy and discrete charge conservation of the ALE-FEM. Drop of heavy fluid in light fluid is taken as benchmark to demonstrate the robustness of the preconditioner with respect to the discontinuity of physical parameters and the mesh size.
ER -