An Efficient Algorithm for Computation of MHD Flow Ensembles

Abstract An efficient algorithm is proposed and studied for computing flow ensembles of incompressible magnetohydrodynamic (MHD) flows under uncertainties in initial or boundary data. The ensemble average of J realizations is approximated through a clever algorithm (adapted from a breakthrough idea of Jiang and Layton [23]) that, at each time step, uses the same matrix for each of the J systems solves. Hence, preconditioners need to be built only once per time step, and the algorithm can take advantage of block linear solvers. Additionally, an Elsässer variable formulation is used, which allows for a stable decoupling of each MHD system at each time step. We prove stability and convergence of the algorithm, and test it with two numerical experiments.

[1]  Nan Jiang,et al.  AN ALGORITHM FOR FAST CALCULATION OF FLOW ENSEMBLES , 2014 .

[2]  Shangyou Zhang,et al.  A new family of stable mixed finite elements for the 3D Stokes equations , 2004, Math. Comput..

[3]  Leo G. Rebholz,et al.  Decoupled, Unconditionally Stable, Higher Order Discretizations for MHD Flow Simulation , 2017, J. Sci. Comput..

[4]  L. R. Scott,et al.  The Mathematical Theory of Finite Element Methods , 1994 .

[5]  Nan Jiang,et al.  Numerical analysis of two ensemble eddy viscosity numerical regularizations of fluid motion , 2015 .

[6]  Friedrich Kupka,et al.  Interdisciplinary aspects of turbulence , 2009 .

[7]  Christopher J. Freitas,et al.  The issue of numerical uncertainty , 2002 .

[8]  Frédéric Hecht,et al.  New development in freefem++ , 2012, J. Num. Math..

[9]  J. D. Barrow,et al.  Cosmology with inhomogeneous magnetic fields , 2007 .

[10]  Nan Jiang,et al.  A Higher Order Ensemble Simulation Algorithm for Fluid Flows , 2015, J. Sci. Comput..

[11]  José A. Font,et al.  General Relativistic Hydrodynamics and Magnetohydrodynamics: Hyperbolic Systems in Relativistic Astrophysics , 2008 .

[12]  Tim N. Palmer,et al.  Ensemble forecasting , 2008, J. Comput. Phys..

[13]  R. Rannacher,et al.  Finite-element approximations of the nonstationary Navier-Stokes problem. Part IV: error estimates for second-order time discretization , 1990 .

[14]  Douglas N. Arnold,et al.  Quadratic velocity/linear pressure Stokes elements , 1992 .

[15]  William Layton,et al.  Introduction to the Numerical Analysis of Incompressible Viscous Flows , 2008 .

[16]  Peter Bodenheimer.  Numerical methods in astrophysics , 2007 .

[17]  John M. Lewis,et al.  Roots of Ensemble Forecasting , 2005 .

[18]  Volker John,et al.  On the Divergence Constraint in Mixed Finite Element Methods for Incompressible Flows , 2015, SIAM Rev..

[19]  Vivette Girault,et al.  Finite Element Methods for Navier-Stokes Equations - Theory and Algorithms , 1986, Springer Series in Computational Mathematics.

[20]  R. Ghanem,et al.  Stochastic Finite Elements: A Spectral Approach , 1990 .

[21]  Max D. Gunzburger,et al.  An Ensemble-Proper Orthogonal Decomposition Method for the Nonstationary Navier-Stokes Equations , 2016, SIAM J. Numer. Anal..

[22]  Jiajia Waters,et al.  Time relaxation algorithm for flow ensembles , 2016 .

[23]  Hidetoshi Hashizume,et al.  Numerical and experimental research to solve MHD problem in liquid blanket system , 2006 .

[24]  B. Punsly,et al.  Black hole gravitohydromagnetics , 2001 .

[25]  Omar M. Knio,et al.  Spectral Methods for Uncertainty Quantification , 2010 .

[26]  Sergey Smolentsev,et al.  MHD thermofluid issues of liquid-metal blankets: Phenomena and advances , 2010 .

[27]  M. Akbas,et al.  Numerical analysis and testing of a fully discrete, decoupled penalty-projection algorithm for mhd in elsässer variable , 2016 .

[28]  Jim Dowling,et al.  Predicting probability distributions for surf height using an ensemble of mixture density networks , 2005, ICML.

[29]  Catalin Trenchea,et al.  Unconditional stability of a partitioned IMEX method for magnetohydrodynamic flows , 2014, Appl. Math. Lett..

[30]  Ming Xue,et al.  Sensitivity Analysis of Convection of the 24 May 2002 IHOP Case Using Very Large Ensembles , 2006 .

[31]  L. Barleon,et al.  MHD flow in liquid-metal-cooled blankets , 1991 .

[32]  J. G. Osorio,et al.  Building hazard maps of extreme daily rainy events from PDF ensemble, via REA method, on Senegal River Basin , 2011 .

[33]  L. Driel-Gesztelyi.  An Introduction to Magnetohydrodynamics , 2004 .