Muhammad Mohebujjaman

Muhammad Mohebujjaman

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ässer 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.

Birgul Koc, Muhammad Mohebujjaman, Changhong Mou et al.

For reduced order models (ROMs) of fluid flows, we investigate theoretically and computationally whether differentiation and ROM spatial filtering commute, i.e., whether the commutation error (CE) is nonzero. We study the CE for the Laplacian and two ROM filters: the ROM projection and the ROM differential filter. Furthermore, when the CE is nonzero, we investigate whether it has any significant effect on ROMs that are constructed by using spatial filtering. As numerical tests, we use the Burgers equation with viscosities $ν=10^{-1}$ and $ν=10^{-3}$ and a 2D flow past a circular cylinder at Reynolds numbers $Re=1$ and $Re=100$. Our investigation shows that: (i) the CE exists, and (ii) the CE has a significant effect on ROM development for low Reynolds numbers, but not so much for higher Reynolds numbers.