Fractional magneto-hydrodynamics: Algorithms and applications

• Published in: Journal of computational physics Vol. 378; pp. 44 - 62
• Main Authors: Song, Fangying, Karniadakis, George Em
• Format: Journal Article
• Language: English
• Published: Cambridge Elsevier Inc 01.02.2019; Elsevier Science Ltd
• Subjects: Algorithms; Boundary conditions; Cavity flow; Computational fluid dynamics; Computational physics; Convergence; Discretization; Divergence; Domains; Fluid flow; Fluid mechanics; FMHD; Fractional Laplacian; Fractional Navier–Stokes equation; Lid-driven cavity flow; Magnetic fields; Magnetohydrodynamics; Navier-Stokes equations; Time-splitting; Lid-driven cavity flow; Time-splitting; FMHD; Fractional Laplacian; Fractional Navier–Stokes equation
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2018.10.047

•We present the first ever three discretization schemes for solving the fractional magneto-hydrodynamic (FMHD) equations in bounded domains.•We use the lifting method for extending our methods to solve the FMHD equations with inhomogeneous boundary conditions.•The divergence-free constraint of the magnetic field is preserved for all the three schemes.•We applied these schemes for solving the lid-driven cavity flow with and without magnetic field, and observing that new flow patterns emerge as the fractional orders vary. We present two discretization methods for solving the fractional magneto-hydrodynamic (FMHD) equations with the fractional Laplacian defined in bounded domains. In the first method, we add a pseudo-pressure in the magnetic field equation to enforce that the magnetic field is divergence-free. In the second method, the magnetic field satisfies the divergence-free condition automatically with the no-slip boundary condition for the velocity and the perfectly conducting boundary condition for the magnetic field equation. In addition, we propose a discretization method for solving the modified fractional magneto-hydrodynamic (M-FMHD) equation. The divergence-free condition is preserved by employing the stream function in the M-FMHD system. Numerical experiments with fabricated solutions show that all three methods exhibit second-order accuracy in time for the velocity and magnetic fields. The pseudo-pressure and pressure also exhibit second-order convergence in time for small values of time step, however the pressure exhibits 1.5-order for convergence for large values of time step size when the fractional order α=1. We use the spectral decomposition method for spatial discretization, and we demonstrate that exponential convergence is achieved for all fields and for any fractional order. We also applied these methods to solve the lid-driven cavity flow with and without magnetic field, and we observe new flow patterns emerging as the fractional order varies.