arXiv:2610.33385v1 Announce Type: new 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.