“…This approach turned out to be successful in various coupled multiphysics PDE simulations -e.g. in poroelasticity [1], geophysical electromagnetics [8], Cahn-Hilliard Navier-Stokes systems [9], multiphase poromechanics of heterogeneous media [11], linear elasticity in mixed form [13], incompressible (reduced) resistive magnetohydrodynamic [15], incompressible flow simulations [21], elliptic optimal control problems [22], fluid-structure interaction problems in hemodynamics [32], models of coupled magma/mantle dynamics [39] and incompressible Navier-Stokes problems [44]. For fully implicit discretisations of multidimensional radiation diffusion equations, An et al [2], Feng et al [18], and Mousseau et al [33] proposed various operator-based preconditioners in the Jacobian-free Newton-Krylov framework.…”