Related papers: Monotone finite difference schemes for anisotropic…
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.…
We present a robust and accurate numerical method for the anisotropic diffusion equation in curvilinear coordinates. This study extends the recent work [Muir et al., Computer Physics Communications, 2025] for solving the anisotropic…
We design a monotone meshfree finite difference method for linear elliptic equations in the non-divergence form on point clouds via a nonlocal relaxation method. The key idea is a novel combination of a nonlocal integral relaxation of the…
On the example of the Poynting-Thomson-Zener rheological model for solids, which exhibits both dissipation and wave propagation - with nonlinear dispersion relation -, we introduce and investigate a finite difference numerical scheme. Our…
We study a finite volume scheme for the approximation of the solution to convection diffusion equations with nonlinear convection and Robin boundary conditions. The scheme builds on the interpretation of such a continuous equation as the…
Our goal is to develop a flux limiter of the Flux-Corrected Transport method for a nonconservative convection-diffusion equation. For this, we consider a hybrid difference scheme that is a linear combination of a monotone scheme and a…
Sticky diffusion processes on bounded domains spend finite time (and finite mean time) on the lower-dimensional space given by the boundary. Once the process hits the boundary, then it starts again after a random amount of time. While on…
We propose a finite volume scheme for a class of nonlinear parabolic equations endowed with non-homogeneous Dirichlet boundary conditions and which admit relative en-tropy functionals. For this kind of models including porous media…
Anisotropic mesh adaptation is studied for linear finite element solution of 3D anisotropic diffusion problems. The M-uniform mesh approach is used, where an anisotropic adaptive mesh is generated as a uniform one in the metric specified by…
We investigate the connections between several recent methods for the discretization of anisotropic heterogeneous diffusion operators on general grids. We prove that the Mimetic Finite Difference scheme, the Hybrid Finite Volume scheme and…
A common numerical task is to represent functions which are highly spatially anisotropic, and to solve differential equations related to these functions. One way such anisotropy arises is that information transfer along one spatial…
We study an inverse design problem for the linear multiple fragmentation equation arising in particle dynamics. Our objective is to reconstruct an unknown initial size distribution that evolves, under a prescribed fragmentation law, into a…
Explicit non-oscillatory central difference schemes become excessively diffusive when applied to highly nonlinear advection problems where small time steps are necessary to maintain stability. Here, we present a correction to reduce such…
This paper provides a full characterization of the value function and solution(s) of an optimal stopping problem for a one-dimensional diffusion with an integral criterion. The results hold under very weak assumptions, namely, the diffusion…
The multiple-Dirichlet-to-Neumann (multiple-DtN) non-reflecting boundary condition is adapted to acoustic scattering from obstacles embedded in the half-plane. The multiple-DtN map is coupled with the method of images as an alternative…
We propose a class of nonlocal diffusion systems on time-varying domains, and fully characterize their asymptotic dynamics in the asymptotically fixed, time-periodic and unbounded cases. The kernel is not necessarily symmetric or compactly…
The standard method to generate dynamical models with a finite extent is to apply a truncation in binding energy to the distribution function. This approach has the disadvantages that one cannot choose the density to start with, that the…
Preserving scalar boundedness is important for numerical schemes used in turbulent compressible multi-component flow simulations to prevent unphysical results and unstable simulations. However, ensuring scalar boundedness for high-order,…
We report numerical simulations of a strongly biased diffusion process on a one-dimensional substrate with directed shortcuts between randomly chosen sites, i.e. with a small-world-like structure. We find that, unlike many other dynamical…
A paradigmatic nonhyperbolic dynamical system exhibiting deterministic diffusion is the smooth nonlinear climbing sine map. We find that this map generates fractal hierarchies of normal and anomalous diffusive regions as functions of the…