Fully three-dimensional time-resolved MHD simulations in wall bounded geometry

We present a new method for computing 3D viscous-resistive MHD turbulence in wall bounded geometries of arbitrary shape [1]. The numerical scheme is based on a classical Fourier pseudo-spectral solver combined with a volume penalization method to impose the boundary condition of the velocity and magnetic field. The new code is validated using different test cases, such as three-dimensional Taylor-Couette flow and MHD in cylindrical geometry. Imposing helical magnetic boundary conditions in the latter geometry, the flow shows a self-organization to a chaotic state for elevated Hartmann number. In toroidal geometry, this self-organization leads to toroidal velocity fields. \\[4pt] [1] Morales, Leroy, Bos and Schneider (submitted).