Improved lattice Boltzmann method for conjugate magnetohydrodynamic simulations

arXiv:2609.18097v1 Announce Type: new Abstract: In the simulation of conjugate magnetohydrodynamic (MHD) flows, where a moving conducting fluid dynamically interacts with bounding conducting walls under the influence of an externally imposed magnetic field, resolving the electrical conductivity transition across the fluid-solid interface introduces significant numerical challenges. In problems with high magnetic Reynolds numbers, it is required to solve the full magnetic induction formulation that has a curl-of-curl term for the magnetic diffusion, which is simplified to a divergence term in the original vector-valued lattice Boltzmann method (LBM). The resulting LBM scheme is valid for simulating fluid domains with a constant conductivity where its bounding solid walls are modelled using appropriate boundary conditions. The current study shows that this scheme is also valid for two-dimensional conjugate MHD simulations where the magnetic field component with non-zero gradients is perpendicular to the conductivity gradient. For general conjugate MHD simulations, we improve the LBM scheme by considering the difference between the curl-of-curl term and the simplified divergence form, which is treated as a source term in the LBM evolution algorithm and computed using an efficient LBM discretisation scheme. The improved LBM scheme automatically satisfies the required conjugate constraints at the fluid-solid interface as a volume integral inside two adjacent boundary layers. The accuracy of the improved LBM scheme is verified by piecewise analytical solutions in benchmark problems with an abrupt conductivity jump across the fluid-solid interface.