Date of Award
Doctor of Philosophy (PhD)
Dr. Leo Rebholz, Committee Chair
Dr. Timo Heister, Co-chair
Dr. Hyesuk Lee
Dr. Qingshan Chen
This dissertation studies eﬃcient numerical methods for approximating solu-tions to viscous, incompressible, time-dependent magnetohydrodynamic (MHD) ﬂows and computing MHD ﬂows ensembles. Chapter 3 presents and analyzes a fully discrete, decoupled eﬃcient algorithm for MHD ﬂow that is based on the Els¨asser variable formulation, proves its uncondi-tional stability with respect to the timestep size, and proves its unconditional con-vergence. Numerical experiments are given which verify all predicted convergence rates of our analysis, show the results of the scheme on a set of channel ﬂow problems match well the results found when the computation is done with MHD in primitive variables, and ﬁnally illustrate that the scheme performs well for channel ﬂow over a step. In chapter 4, we propose, analyze, and test a new MHD discretization which decouples the system into two Oseen problems at each timestep, yet maintains un-conditional stability with respect to timestep size. The scheme is optimally accu-rate in space, and behaves like second order in time in practice. The proposed method chooses θ ∈ [0, 1], dependent on the viscosity ν and magnetic diﬀusiv-ity νm, so that unconditionally stability is achieved, and gives temporal accuracy O(∆t2 + (1 − θ)|ν − νm|∆t). In practice, ν and νm are small, and so the method be-haves like second order. We show the θ-method provides excellent accuracy in cases
where usual BDF2 is unstable. Chapter 5 proposes an eﬃcient algorithm and studies for computing ﬂow en-sembles of incompressible MHD ﬂows under uncertainties in initial or boundary data. The ensemble average of J realizations is approximated through an eﬃcient algo-rithm that, at each time step, uses the same coeﬃcient matrix for each of the J system 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¨asser 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. This work concludes with chapter 6, which proposes, analyzes and tests high order algebraic splitting methods for MHD ﬂows. The key idea is to applying Yosida-type algebraic splitting to the incremental part of the unknowns at each time step. This reduces the block Schur complement by decoupling it into two Navier-Stokes-type Schur complements, each of which is symmetric positive deﬁnite and the same at each time step. We prove the splitting is third order in ∆t, and if used together with (block-)pressure correction, is fourth order. A full analysis of the solver is given, both as a linear algebraic approximation, and as a ﬁnite element discretization of an approximation to the un-split discrete system. Numerical tests are given to illustrate the theory and show the eﬀectiveness of the method. Finally, conclusions and future works are discussed in the ﬁnal chapter.
Mohebujjaman, Muhammad, "Efficient Numerical Methods for Magnetohydrodynamic Flow" (2017). All Dissertations. 2027.