High order efficient algorithm for computation of MHD flow ensembles

DOI10.4208/AAMM.OA-2021-0154MaRDI QIDQ6358483arXivFDO

Authors Muhammad Mohebujjaman

Publication date 17 January 2021

Published in arXiv

Full work available at URL https://doi.org/10.4208/AAMM.OA-2021-0154

Mathematics Subject Classification ID

PDEs in connection with fluid mechanics (35Q35) Finite difference methods for initial value and initial-boundary value problems involving PDEs (65M06) Stability and convergence of numerical methods for initial value and initial-boundary value problems involving PDEs (65M12) Numerical solution of discretized equations for initial value and initial-boundary value problems involving PDEs (65M22) Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs (65M60) Error bounds for boundary value problems involving PDEs (65N15) Magnetohydrodynamics and electrohydrodynamics (76W05)

arXiv classification

cs.NAmath.NA

Abstract: In this paper, we propose, analyze, and test a new fully discrete, efficient, decoupled, stable, and practically second-order time-stepping algorithm for computing MHD ensemble flow averages under uncertainties in the initial conditions and forcing. For each viscosity and magnetic diffusivity pair, the algorithm picks the largest possible parameter $hetain[0,1]$ to avoid the instability that arises due to the presence of some explicit viscous terms. At each time step, the algorithm shares the same system matrix with all $J$ realizations but with different right-hand-side vectors. That saves assembling time and computer memory, allows the reuse of the same preconditioner, and can take the advantage of block linear solvers. For the proposed algorithm, we prove stability and convergence rigorously. To illustrate the predicted convergence rates of our analysis, numerical experiments with manufactured solutions are given on a unit square domain. Finally, we test the scheme on a benchmark channel flow over a step problem and it performs well.

This page was built for publication: High order efficient algorithm for computation of MHD flow ensembles

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6358483)