Source author record

Andrea Barth

Andrea Barth appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

8works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2026arXiv

Deep learning methods for stochastic Galerkin approximations of elliptic random PDEs

This work considers stochastic Galerkin approximations of linear elliptic partial differential equations (PDEs) with stochastic forcing terms and stochastic diffusion coefficients, that cannot be bounded uniformly away from zero and infinity. A traditional numerical method for solving the resulting high-dimensional coupled system of PDEs is replaced by deep learning techniques. In order to achieve this, physics-informed neural networks (PINNs), which typically operate on the strong residual of the PDE and can therefore be applied in a wide range of settings, are considered. As a second approach, the Deep Ritz method, which is a neural network that minimizes the Ritz energy functional to find the weak solution, is employed. While the second approach only works in special cases, it overcomes the necessity of testing in variational problems while maintaining mathematical rigor and ensuring the existence of a unique solution. Furthermore, the residual is of a lower differentiation order, reducing the training cost considerably. The efficiency of the method is demonstrated on several model problems.

preprint2022arXiv

On properties and applications of Gaussian subordinated Lévy fields

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 $\mathcal{D}\subset \mathbb{R}^d$, $d\geq 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.

preprint2021arXiv

Multilevel Monte Carlo estimators for elliptic PDEs with Lévy-type diffusion coefficient

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.

preprint2021arXiv

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

Subsurface flows are commonly modeled by advection-diffusion equations. Insufficient measurements or uncertain material procurement may be accounted for by random coefficients. To represent, for example, transitions in heterogeneous media, the parameters of the equation are spatially discontinuous. Specifically, a scenario with coupled advection- and diffusion coefficients that are modeled as sums of continuous random fields and discontinuous jump components are considered. For the numerical approximation of the solution, an adaptive, pathwise discretization scheme based on a Finite Element approach is introduced. To stabilize the numerical approximation and accelerate convergence, the discrete space-time grid is chosen with respect to the varying discontinuities in each sample of the coefficients, leading to a stochastic formulation of the Galerkin projection and the Finite Element basis.

preprint2020arXiv

A Multilevel Monte Carlo Algorithm for Parabolic Advection-Diffusion Problems with Discontinuous Coefficients

The Richards' equation is a model for flow of water in unsaturated soils. The coefficients of this (nonlinear) partial differential equation describe the permeability of the medium. Insufficient or uncertain measurements are commonly modeled by random coefficients. For flows in heterogeneous\textbackslash fractured\textbackslash porous media, the coefficients are modeled as discontinuous random fields, where the interfaces along the stochastic discontinuities represent transitions in the media. More precisely, the random coefficient is given by the sum of a (continuous) Gaussian random field and a (discontinuous) jump part. In this work moments of the solution to the random partial differential equation are calculated using a path-wise numerical approximation combined with multilevel Monte Carlo sampling. The discontinuities dictate the spatial discretization, which leads to a stochastic grid. Hence, the refinement parameter and problem-dependent constants in the error analysis are random variables and we derive (optimal) a-priori convergence rates in a mean-square sense.

preprint2020arXiv

Subordinated Gaussian Random Fields

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 validate our theoretical findings in various numerical examples.

preprint2020arXiv

Subordinated Gaussian Random Fields in Elliptic Partial Differential Equations

To model subsurface flow in uncertain heterogeneous\ 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. In this paper a new subordination approach is employed 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.

preprint2016arXiv

Weak convergence of Galerkin approximations of stochastic partial differential equations driven by additive Lévy noise

This work considers weak approximations of stochastic partial differential equations (SPDEs) driven by Lévy noise. The SPDEs at hand are parabolic with additive noise processes. A weak-convergence rate for the corresponding Galerkin approximation is derived. The convergence result is derived by use of the Malliavin derivative rather then the common approach via the Kolmogorov backward equation.