Related papers: A Non-compact Positivity-Preserving Numerical Sche…
We propose a novel non-compact, positivity-preserving scheme for linear non-divergence form parabolic equations. Based on the Feynman-Kac formula, the solution is expressed as a conditional expectation of an associated diffusion process.…
The Feynman-Kac formula provides a way to understand solutions to elliptic partial differential equations in terms of expectations of continuous time Markov processes. This connection allows for the creation of numerical schemes for…
This work presents a probabilistic scheme for solving semilinear nonlocal diffusion equations with volume constraints and integrable kernels. The nonlocal model of interest is defined by a time-dependent semilinear partial…
We present a positive and asymptotic preserving numerical scheme for solving linear kinetic, transport equations that relax to a diffusive equation in the limit of infinite scattering. The proposed scheme is developed using a standard…
We prove Feynman-Kac formulas for solutions to elliptic and parabolic boundary value and obstacle problems associated with a general Markov diffusion process. Our diffusion model covers several popular stochastic volatility models, such as…
In this work, we develop novel structure-preserving numerical schemes for a class of nonlinear Fokker--Planck equations with nonlocal interactions. Such equations can cover many cases of importance, such as porous medium equations with…
The concern of the present work is the introduction of a very efficient Asymptotic Preserving scheme for the resolution of highly anisotropic diffusion equations. The characteristic features of this scheme are the uniform convergence with…
We present a computational alternative to probabilistic simulations for non-smooth stochastic dynamical systems that are prevalent in engineering mechanics. As examples, we target (1) stochastic elasto-plastic problems, which involve…
The fractional Feynman-Kac equations describe the distribution of functionals of non-Brownian motion, or anomalous diffusion, including two types called the forward and backward fractional Feynman-Kac equations, where the fractional…
The Feynman-Kac equations are a type of partial differential equations describing the distribution of functionals of diffusive motion. The probability density function (PDF) of Brownian functionals satisfies the Feynman-Kac formula, being a…
In this paper, we study numerical methods for the homogenization of linear second-order elliptic equations in nondivergence-form with periodic diffusion coefficients and large drift terms. Upon noting that the effective diffusion matrix can…
In this work we consider an extension of a recently proposed structure preserving numerical scheme for nonlinear Fokker-Planck-type equations to the case of nonconstant full diffusion matrices. While in existing works the schemes are…
In this paper we study numerical positivity and contractivity in the infinite norm of Crank-Nicolson method when it is applied to the diffusion equation with homogeneous Dirichlet boundary conditions. For this purpose, the amplification…
We present a novel structure-preserving numerical scheme for discontinuous finite element approximations of nonlinear hyperbolic systems. The method can be understood as a generalization of the Lax-Friedrichs flux to a high-order staggered…
In this article we introduce an asymptotic preserving scheme designed to compute the solution of a two dimensional elliptic equation presenting large anisotropies. We focus on an anisotropy aligned with one direction, the dominant part of…
In this paper, we use the theory of symmetric Dirichlet forms to derive Feynman-Kac formulae for the forward problem of electrical impedance tomography with possibly anisotropic, merely measurable conductivities corresponding to different…
The present paper introduces an efficient and accurate numerical scheme for the solution of a highly anisotropic elliptic equation, the anisotropy direction being given by a variable vector field. This scheme is based on an asymptotic…
In this paper, a high-order exponential scheme is developed to solve the 1D unsteady convection-diffusion equation with Neumann boundary conditions. The present method applies fourth-order compact exponential difference scheme in spatial…
We propose a locally conservative enriched Galerkin scheme that preserves the physical bounds for an elliptic problem. To this end, we use a substantial over-penalization of the discrete solution's jumps to obtain optimal convergence. To…
We study the uniqueness of reaction-diffusion steady states in general domains with Dirichlet boundary data. Here we consider "positive" (monostable) reactions. We describe geometric conditions on the domain that ensure uniqueness and we…