Abstract:We develop new multilevel Monte Carlo (MLMC) methods to estimate the expectation of the smallest eigenvalue of a stochastic convection-diffusion operator with random coefficients. The MLMC method is based on a sequence of finite element (FE) discretizations of the eigenvalue problem on a hierarchy of increasingly finer meshes. For the discretized, algebraic eigenproblems we use both the Rayleigh quotient (RQ) iteration and implicitly restarted Arnoldi (IRA), providing an analysis of the cost in each case. By studying the variance on each level and adapting classical FE error bounds to the stochastic setting, we are able to bound the total error of our MLMC estimator and provide a complexity analysis. As expected, the complexity bound for our MLMC estimator is superior to plain Monte Carlo. To improve the efficiency of the MLMC further, we exploit the hierarchy of meshes and use coarser approximations as starting values for the eigensolvers on finer ones. To improve the stability of the MLMC method for convection-dominated problems, we employ two additional strategies. First, we consider the streamline upwind Petrov--Galerkin formulation of the discrete eigenvalue problem, which allows us to start the MLMC method on coarser meshes than is possible with standard FEs. Second, we apply a homotopy method to add stability to the eigensolver for each sample. Finally, we present a multilevel quasi-Monte Carlo method that replaces Monte Carlo with a quasi-Monte Carlo (QMC) rule on each level. Due to the faster convergence of QMC, this improves the overall complexity. We provide detailed numerical results comparing our different strategies to demonstrate the practical feasibility of the MLMC method in different use cases. The results support our complexity analysis and further demonstrate the superiority over plain Monte Carlo in all cases.

A variational MAX ensemble numerical algorism for a transient heat model with random inputs

Improved error estimates of ensemble Monte Carlo methods for random transient heat equations with uncertain inputs

Heat Equation Stein Variational Ensemble: Rethinking and Advancing Uncertainty-Aware Soft Sensor Modeling

Monte-Carlo Simulation of a Stochastic Differential Equation

Explicit numerical approximations for McKean-Vlasov stochastic differential equations in finite and infinite time

Multilevel and Multi-index Monte Carlo methods for the McKean-Vlasov equation

Ensemble Timestepping Algorithms for the Heat Equation with Uncertain Conductivity

Moving Mesh Method with Variational Multiscale Finite Element Method for Convection–diffusion–reaction Equations

Stochastic simulation of magnetically controlled convection in porous media with random porosity via Karhunen-Loève expansion and intrusive polynomial chaos expansion

An ensemble algorithm for numerical solutions to deterministic and random parabolic PDEs

Multilevel Monte Carlo methods for stochastic convection-diffusion eigenvalue problems

Lattice Models for Magnetic Fluids: Correlations Between Order Parameters

New Ensemble Domain Decomposition Method for the Steady-state Random Stokes-Darcy Coupled Problems with Uncertain Parameters

Variational multiscale virtual element method for the convection-dominated diffusion problem

Multi-Group Maximum Entropy Method: Modeling Translational Non-Equilibrium

A semi-implicit stochastic multiscale method for radiative heat transfer problem

Extended Ensemble Monte Carlo

Constrained Ensemble Langevin Monte Carlo

Variational Microcanonical Estimator

Unconditionally energy stable and first-order accurate numerical schemes for the heat equation with uncertain temperature-dependent conductivity

Numerical Analysis of Penalty-based Ensemble Methods