To model subsurface flow in uncertain heterogeneous or fractured media an elliptic equation with a discontinuous stochastic diffusion coefficient—also called random field—may be used. In case of a one-dimensional parameter space, Lévy processes allow for jumps and display great flexibility in the distributions used. However, in various situations (e.g. microstructure modeling), a one-dimensional parameter space is not sufficient. Classical extensions of Lévy processes on two parameter dimensions suffer from the fact that they do not allow for spatial discontinuities [see for example Barth and Stein (Stoch Part Differ Equ Anal Comput 6(2):286–334, 2018)]. In this paper a new subordination approach is employed [see also Barth and Merkle (Subordinated gaussian random fields. ArXiv e-prints, arXiv:2012.06353 [math.PR], 2020)] to generate Lévy-type discontinuous random fields on a two-dimensional spatial parameter domain. Existence and uniqueness of a (pathwise) solution to a general elliptic partial differential equation is proved and an approximation theory for the diffusion coefficient and the corresponding solution provided. Further, numerical examples using a Monte Carlo approach on a Finite Element discretization validate our theoretical results.
General elliptic equations with spatially discontinuous diffusion coefficients may be used as a simplified model for subsurface flow in heterogeneous or fractured porous media. In such a model, data sparsity and measurement errors are often taken into account by a randomization of the diffusion coefficient of the elliptic equation which reveals the necessity of the construction of flexible, spatially discontinuous random fields. Subordinated Gaussian random fields are random functions on higher dimensional parameter domains with discontinuous sample paths and great distributional flexibility. In the present work, we consider a random elliptic partial differential equation (PDE) where the discontinuous subordinated Gaussian random fields occur in the diffusion coefficient. Problem specific multilevel Monte Carlo (MLMC) Finite Element methods are constructed to approximate the mean of the solution to the random elliptic PDE. We prove a-priori convergence of a standard MLMC estimator and a modified MLMC—control variate estimator and validate our results in various numerical examples.
Motivated by the subordinated Brownian motion, we define a new class of (in general discontinuous) random fields on higher-dimensional parameter domains: the subordinated Gaussian random field. We investigate the pointwise marginal distribution of the constructed random fields, derive a Lévy-Khinchin-type formula and semi-explicit formulas for the covariance function. Further, we study the pointwise stochastic regularity and present various numerical examples.
We consider Gaussian subordinated Lévy fields (GSLFs) that arise by subordinating Lévy processes with positive transformations of Gaussian random fields on some spatial domain D ⊂ R d , d ≥ 1. The resulting random fields are distributionally flexible and have in general discontinuous sample paths. Theoretical investigations of the random fields include pointwise distributions, possible approximations and their covariance function. As an application, a random elliptic PDE is considered, where the constructed random fields occur in the diffusion coefficient. Further, we present various numerical examples to illustrate our theoretical findings.
We consider Gaussian subordinated Lévy fields (GSLFs) that arise by subordinating Lévy processes with positive transformations of Gaussian random fields on some spatial domain. The resulting random fields are distributionally flexible and have in general discontinuous sample paths. Theoretical investigations of the random fields include pointwise distributions, possible approximations and their covariance function. As an application, a random elliptic PDE is considered, where the constructed random fields occur in the diffusion coefficient. Further, we present various numerical examples to illustrate our theoretical findings.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
customersupport@researchsolutions.com
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.