|Authors||D. G. E. Grigoriadis and S. C. Kassinos|
|Title||Extension of the Immersed Boundary Method for Liquid-Metal Magnetohydrodynamics|
|Afilliation||, Scientific Computing|
|Project(s)||Center for Biomedical Computing (SFF)|
|Publication Type||Proceedings, refereed|
|Year of Publication||2009|
|Conference Name||Immersed Boundary Methods: Current status and Future Research Directions|
Liquid-metal magnetohydrodynamics (MHD hereafter) are met in a variety of challenging industrial applications (metalheating, pumping, stirring, blankets of fusion reactors) as well as in geophysics whenever a magnetic field interacts with the motion of a conductive fluid. The present study forms the first attempt to simulate these flows using a non-boundary conforming methodology such as the immersed boundary method (IB hereafter). It is motivated by the belief that the IB method could bring to MHD simulations the same advantages it offers for hydrodynamics. We have recently presented  an extension of the IB methodology for efficient MHD simulations in geometrically complicated domains using a Cartesian flow solver and an identical, direct Poisson solver for both the hydrodynamic pressure and the electric potential fields. A projection scheme and a suitable forcing scheme for electric density currents in the vicinity of non-conducting immersed surfaces were suggested. The methodology has been extensively tested for several MHD flows at Hartmann numbers Ha=5-2000, reaching an excellent agreement with numerical or analytical solutions. The projection and forcing schemes for current densities were found capable of satisfying the charge conservation law in the presence of immersed boundaries. Finally, we demonstrate how the proposed extension can be used to extend the applicability of existing flow solvers using the IB concept, and extend the range of computable wall-bounded MHD flows.