Numerical analysis for time-dependent advection-diffusion problems with random discontinuous coefficients

As an extension to the well-established stationary elliptic partial differential equation (PDE) with random continuous coefficients we study a time-dependent advection-diffusion problem, where the coefficients may have random spatial discontinuities. In a subsurface flow model, the randomness in a parabolic equation may account for insufficient measurements or uncertain material procurement, while the discontinuities could represent transitions in heterogeneous media. Specifically, a scenario with coupled advection and diffusion coefficients that are modeled as sums of continuous random fields and discontinuous jump components are considered. The respective coefficient functions allow a very flexible modeling, however, they also complicate the analysis and numerical approximation of the corresponding random parabolic PDE. We show that the model problem is indeed well-posed under mild assumptions and show measurability of the pathwise solution. For the numerical approximation we employ a sample-adapted, pathwise discretization scheme based on a finite element approach. This semi-discrete method accounts for the discontinuities in each sample, but leads to stochastic, finite-dimensional approximation spaces. We ensure measurability of the semi-discrete solution, which in turn enables us to derive moments bounds on the mean-squared approximation error. By coupling this semi-discrete approach with suitable coefficient approximation and a stable time stepping, we obtain a fully discrete algorithm to solve the random parabolic PDE. We provide an overall error bound for this scheme and illustrate our results with several numerical experiments.