Second order ensemble simulation for MHD flow in Els\"asser variable with noisy input data

TL;DR

This paper introduces a second-order, decoupled ensemble simulation algorithm for MHD flows using Els"asser variables, which reduces computational costs and is validated through stability analysis and numerical experiments.

Contribution

The paper presents a novel second-order ensemble algorithm for MHD flows that is efficient, decoupled, and based on Els"asser variables, with proven stability and convergence.

Findings

Algorithm achieves optimal convergence rates.

Reduces storage and computational costs.

Performs well on benchmark channel flow.

Abstract

We propose, analyze and test a fully discrete, efficient second-order algorithm for computing flow ensembles average of viscous, incompressible, and time-dependent magnetohydrodynamic (MHD) flows under uncertainties in initial conditions. The scheme is decoupled and based on Els\"asser variable formulation. The algorithm uses the breakthrough idea of Jiang and Layton, 2014 to approximate the ensemble average of $J$ realizations. That is, at each time step, each of the $J$ realization shares the same coefficient matrix for different right-hand side matrices. Thus, storage requirements and computational time are reduced by building preconditioners once per time step and reuse them. We prove stability and optimal convergence with respect to the time step restriction. On some manufactured solutions, numerical experiments are given to verify the predicted convergence rates of our analysis. Finally, we test the scheme on a benchmark channel flow over a step and it performs well.