An efficient algorithm for simulating ensembles of parameterized MHD flow problems

by   Muhammad Mohebujjaman, et al.

In this paper, we propose, analyze, and test an efficient algorithm for computing ensemble average of incompressible magnetohydrodynamics (MHD) flows, where instances/members correspond to varying kinematic viscosity, magnetic diffusivity, body forces, and initial conditions. The algorithm is decoupled in Elsässer variables and permits a shared coefficient matrix for all members at each time-step. Thus, the algorithm is much more computationally efficient than separately computing simulations for each member using usual MHD algorithms. We prove the proposed algorithm is unconditionally stable and convergent. Several numerical tests are given to support the predicted convergence rates. Finally, we test the proposed scheme and observe how the physical behavior changes as the coupling number increases in a lid-driven cavity problem with mean Reynolds number Re≈ 15000, and as the deviation of uncertainties in the initial and boundary conditions increases in a channel flow past a step problem.



There are no comments yet.


page 18

page 19

page 20


High order efficient algorithm for computation of MHD flow ensembles

In this paper, we propose, analyze, and test a new fully discrete, effic...

A gPAV-Based Unconditionally Energy-Stable Scheme for Incompressible Flows with Outflow/Open Boundaries

We present an unconditionally energy-stable scheme for approximating the...

Numerical Approximations and Error Analysis of the Cahn-Hilliard Equation with Dynamic Boundary Conditions

We consider the numerical approximations of the Cahn-Hilliard equation w...

Numerical approximations and error analysis of the Cahn-Hilliard equation with reaction rate dependent dynamic boundary conditions

We consider numerical approximations and error analysis for the Cahn-Hil...

A Two-Level Fourth-Order Approach For Time-Fractional Convection-Diffusion-Reaction Equation With Variable Coefficients

This paper develops a two-level fourth-order scheme for solving time-fra...

Industrial scale large eddy simulations (LES) with adaptive octree meshes using immersogeometric analysis

We present a variant of the immersed boundary method integrated with oct...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.