Stochastics and Partial Differential Equations-Analysis and Computations

A solution theory for a general class of SPDEs
Süß A and Waurick M
In this article we present a way of treating stochastic partial differential equations with multiplicative noise by rewriting them as stochastically perturbed evolutionary equations in the sense of Picard and McGhee (Partial differential equations: a unified Hilbert space approach, DeGruyter, Berlin, 2011), where a general solution theory for deterministic evolutionary equations has been developed. This allows us to present a unified solution theory for a general class of stochastic partial differential equations (SPDEs) which we believe has great potential for further generalizations. We will show that many standard stochastic PDEs fit into this class as well as many other SPDEs such as the stochastic Maxwell equation and time-fractional stochastic PDEs with multiplicative noise on sub-domains of . The approach is in spirit similar to the approach in DaPrato and Zabczyk (Stochastic equations in infinite dimensions, Cambridge University Press, Cambridge, 2008), but complementing it in the sense that it does not involve semi-group theory and allows for an effective treatment of coupled systems of SPDEs. In particular, the existence of a (regular) fundamental solution or Green's function is not required.
Eikonal equations and pathwise solutions to fully non-linear SPDEs
Friz PK, Gassiat P, Lions PL and Souganidis PE
We study the existence and uniqueness of the stochastic viscosity solutions of fully nonlinear, possibly degenerate, second order stochastic pde with quadratic Hamiltonians associated to a Riemannian geometry. The results are new and extend the class of equations studied so far by the last two authors.
A pedestrian approach to the invariant Gibbs measures for the 2- defocusing nonlinear Schrödinger equations
Oh T and Thomann L
We consider the defocusing nonlinear Schrödinger equations on the two-dimensional compact Riemannian manifold without boundary or a bounded domain in . Our aim is to give a pedagogic and self-contained presentation on the Wick renormalization in terms of the Hermite polynomials and the Laguerre polynomials and construct the Gibbs measures corresponding to the Wick ordered Hamiltonian. Then, we construct global-in-time solutions with initial data distributed according to the Gibbs measure and show that the law of the random solutions, at any time, is again given by the Gibbs measure.
Non-stationary phase of the MALA algorithm
Kuntz J, Ottobre M and Stuart AM
The Metropolis-Adjusted Langevin Algorithm (MALA) is a Markov Chain Monte Carlo method which creates a Markov chain reversible with respect to a given target distribution, , with Lebesgue density on ; it can hence be used to approximately sample the target distribution. When the dimension is large a key question is to determine the computational cost of the algorithm as a function of . The measure of efficiency that we consider in this paper is the (ESJD), introduced in Roberts et al. (Ann Appl Probab 7(1):110-120, 1997). To determine how the cost of the algorithm (in terms of ESJD) increases with dimension , we adopt the widely used approach of deriving a diffusion limit for the Markov chain produced by the MALA algorithm. We study this problem for a class of target measures which is in product form and we address the situation of practical relevance in which the algorithm is started out of stationarity. We thereby significantly extend previous works which consider either measures of product form, when the Markov chain is started out of stationarity, or non-product measures (defined via a density with respect to a Gaussian), when the Markov chain is started in stationarity. In order to work in this non-stationary and non-product setting, significant new analysis is required. In particular, our diffusion limit comprises a stochastic PDE coupled to a scalar ordinary differential equation which gives a measure of how far from stationarity the process is. The family of non-product target measures that we consider in this paper are found from discretization of a measure on an infinite dimensional Hilbert space; the discretised measure is defined by its density with respect to a Gaussian random field. The results of this paper demonstrate that, in the non-stationary regime, the cost of the algorithm is of in contrast to the stationary regime, where it is of .
An order approach to SPDEs with antimonotone terms
Scarpa L and Stefanelli U
We consider a class of parabolic stochastic partial differential equations featuring an antimonotone nonlinearity. The existence of unique maximal and minimal variational solutions is proved via a fixed-point argument for nondecreasing mappings in ordered spaces. This relies on the validity of a comparison principle.
Some results on the penalised nematic liquid crystals driven by multiplicative noise: weak solution and maximum principle
Brzeźniak Z, Hausenblas E and Razafimandimby PA
In this paper, we prove several mathematical results related to a system of highly nonlinear stochastic partial differential equations (PDEs). These stochastic equations describe the dynamics of penalised nematic liquid crystals under the influence of stochastic external forces. Firstly, we prove the existence of a global weak solution (in the sense of both stochastic analysis and PDEs). Secondly, we show the pathwise uniqueness of the solution in a 2D domain. In contrast to several works in the deterministic setting we replace the Ginzburg-Landau function by an appropriate polynomial and we give sufficient conditions on the polynomial for these two results to hold. Our third result is a maximum principle type theorem. More precisely, if we consider and if the initial condition satisfies , then the solution also remains in the unit ball.
The hyperbolic Anderson model: moment estimates of the Malliavin derivatives and applications
Balan RM, Nualart D, Quer-Sardanyons L and Zheng G
In this article, we study the hyperbolic Anderson model driven by a space-time Gaussian homogeneous noise with spatial dimension . Under mild assumptions, we provide -estimates of the iterated Malliavin derivative of the solution in terms of the fundamental solution of the wave solution. To achieve this goal, we rely heavily on the of the solution. Our first application are for spatial averages of the solution to the hyperbolic Anderson model, where the rates of convergence are described by the total variation distance. These quantitative results have been elusive so far due to the temporal correlation of the noise blocking us from using the Itô calculus. A ingredient to overcome this difficulty is the coupled with the application of the aforementioned -estimates of the first two Malliavin derivatives. Besides, we provide the corresponding functional central limit theorems. As a second application, we establish the absolute continuity of the law for the hyperbolic Anderson model. The -estimates of Malliavin derivatives are crucial ingredients to verify a local version of Bouleau-Hirsch criterion for absolute continuity. Our approach substantially simplifies the arguments for the one-dimensional case, which has been studied in the recent work by [2].
Existence of dynamical low rank approximations for random semi-linear evolutionary equations on the maximal interval
Kazashi Y and Nobile F
An existence result is presented for the dynamical low rank (DLR) approximation for random semi-linear evolutionary equations. The DLR solution approximates the true solution at each time instant by a linear combination of products of deterministic and stochastic basis functions, both of which evolve over time. A key to our proof is to find a suitable equivalent formulation of the original problem. The so-called Dual Dynamically Orthogonal formulation turns out to be convenient. Based on this formulation, the DLR approximation is recast to an abstract Cauchy problem in a suitable linear space, for which existence and uniqueness of the solution in the maximal interval are established.
Multilevel quadrature for elliptic problems on random domains by the coupling of FEM and BEM
Harbrecht H and Schmidlin M
Elliptic boundary value problems which are posed on a random domain can be mapped to a fixed, nominal domain. The randomness is thus transferred to the diffusion matrix and the loading. While this domain mapping method is quite efficient for theory and practice, since only a single domain discretisation is needed, it also requires the knowledge of the domain mapping. However, in certain applications, the random domain is only described by its random boundary, while the quantity of interest is defined on a fixed, deterministic subdomain. In this setting, it thus becomes necessary to compute a random domain mapping on the whole domain, such that the domain mapping is the identity on the fixed subdomain and maps the boundary of the chosen fixed, nominal domain on to the random boundary. To overcome the necessity of computing such a mapping, we therefore couple the finite element method on the fixed subdomain with the boundary element method on the random boundary. We verify on one hand the regularity of the solution with respect to the random domain mapping required for many multilevel quadrature methods, such as the multilevel quasi-Monte Carlo quadrature using Halton points, the multilevel sparse anisotropic Gauss-Legendre and Clenshaw-Curtis quadratures and multilevel interlaced polynomial lattice rules. On the other hand, we derive the coupling formulation and show by numerical results that the approach is feasible.
Well-posedness for a stochastic 2D Euler equation with transport noise
Lang O and Crisan D
We prove the existence of a unique global strong solution for a stochastic two-dimensional Euler vorticity equation for incompressible flows with noise of transport type. In particular, we show that the initial smoothness of the solution is preserved. The arguments are based on approximating the solution of the Euler equation with a family of viscous solutions which is proved to be relatively compact using a tightness criterion by Kurtz.
Some approximation results for mild solutions of stochastic fractional order evolution equations driven by Gaussian noise
Fahim K, Hausenblas E and Kovács M
We investigate the quality of space approximation of a class of stochastic integral equations of convolution type with Gaussian noise. Such equations arise, for example, when considering mild solutions of stochastic fractional order partial differential equations but also when considering mild solutions of classical stochastic partial differential equations. The key requirement for the equations is a smoothing property of the deterministic evolution operator which is typical in parabolic type problems. We show that if one has access to nonsmooth data estimates for the deterministic error operator together with its derivative of a space discretization procedure, then one obtains error estimates in pathwise Hölder norms with rates that can be read off the deterministic error rates. We illustrate the main result by considering a class of stochastic fractional order partial differential equations and space approximations performed by spectral Galerkin methods and finite elements. We also improve an existing result on the stochastic heat equation.
Delayed blow-up and enhanced diffusion by transport noise for systems of reaction-diffusion equations
Agresti A
This paper is concerned with the problem of regularization by noise of systems of reaction-diffusion equations with mass control. It is known that solutions to such systems of PDEs may blow-up in finite time. Moreover, for many systems of practical interest, establishing whether the blow-up occurs or not is an open question. Here we prove that a suitable multiplicative noise of transport type has a regularizing effect. More precisely, for both a sufficiently noise intensity and a high spectrum, the blow-up of strong solutions is delayed up to an arbitrary large time. Global existence is shown for the case of exponentially decreasing mass. The proofs combine and extend recent developments in regularization by noise and in the -approach to stochastic PDEs, highlighting new connections between the two areas.
Importance sampling for stochastic reaction-diffusion equations in the moderate deviation regime
Gasteratos I, Salins M and Spiliopoulos K
We develop a provably efficient importance sampling scheme that estimates exit probabilities of solutions to small-noise stochastic reaction-diffusion equations from scaled neighborhoods of a stable equilibrium. The moderate deviation scaling allows for a local approximation of the nonlinear dynamics by their linearized version. In addition, we identify a finite-dimensional subspace where exits take place with high probability. Using stochastic control and variational methods we show that our scheme performs well both in the zero noise limit and pre-asymptotically. Simulation studies for stochastically perturbed bistable dynamics illustrate the theoretical results.