Related papers: A positivity preserving numerical scheme for the a…
In this paper, we focus on numerical methods for the genetic drift problems, which is governed by a degenerated convection-dominated parabolic equation. Due to the degeneration and convection, Dirac singularities will always be developed at…
In this paper, we derive two bound-preserving and mass-conserving schemes based on the fractional-step method and high-order compact (HOC) finite difference method for nonlinear convection-dominated diffusion equations. We split the…
Classical solvable stochastic volatility models (SVM) use a CEV process for instantaneous variance where the CEV parameter $\gamma$ takes just few values: 0 - the Ornstein-Uhlenbeck process, 1/2 - the Heston (or square root) process, 1-…
This paper presents a general positivity-preserving algorithm for implicit high-order finite volume schemes solving Euler and Navier-Stokes equations. Previous positivity-preserving algorithms are mainly based on mathematical analyses,…
In this work a simple method to enforce the positivity-preserving property for general high-order conservative schemes is proposed. The method keeps the original scheme unchanged and detects critical numerical fluxes which may lead to…
In this paper, we consider a system of one-dimensional parabolic PDEs, known as the KWC system, as a phase-field model for grain boundary motion. A key feature of this system is that the equation for the crystalline orientation angle is…
It is known in \cite{beccari} that the standard explicit Euler-type scheme (such as the exponential Euler and the linear-implicit Euler schemes) with a uniform timestep, though computationally efficient, may diverge for the stochastic…
In this paper, we introduce and analyze a class of numerical schemes that demonstrate remarkable superiority in terms of efficiency, the preservation of positivity, energy stability, and high-order precision to solve the time-dependent…
Simple modifications for higher-order Godunov-type difference schemes are presented which allow for accurate advection of multi-fluid flows in hydrodynamic simulations. The constraint that the sum of all mass fractions has to be equal to…
We give a probabilistic numerical method for solving a partial differential equation with fractional diffusion and nonlinear drift. The probabilistic interpretation of this equation uses a system of particles driven by L\'evy alpha-stable…
The continuous observation of the financial markets has identified some stylized facts which challenge the conventional assumptions, promoting the born of new approaches. On the one hand, the long-range dependence has been faced replacing…
We investigate which jump-diffusion models are convexity preserving. The study of convexity preserving models is motivated by monotonicity results for such models in the volatility and in the jump parameters. We give a necessary condition…
The {\alpha}-stable L\'evy process, commonly used to describe L\'evy flight, is characterized by discontinuous jumps and is widely used to model anomalous transport phenomena. In this study, we investigate the associated exit problem and…
The dispersion error is often the dominant error for computed solutions of wave propagation problems with high-frequency components. In this paper, we define and give explicit examples of $\alpha$-dispersion-relation-preserving schemes.…
Kinetic equations model distributions of particles in position-velocity phase space. Often, one is interested in studying the long-time behavior of particles in high-collisional regimes in which an approximate (advection)-diffusion model…
This article presents a new finite element method for convection-diffusion equations by enhancing the continuous finite element space with a flux space for flux approximations that preserve the important mass conservation locally on each…
Stable laws can be tempered by modifying the L\'evy measure to cool the probability of large jumps. Tempered stable laws retain their signature power law behavior at infinity, and infinite divisibility. This paper develops random walk…
We study the strong approximation of the solutions to singular stochastic kinetic equations (also referred to as second-order SDEs) driven by $\alpha$-stable processes, using an Euler-type scheme inspired by [11]. For these equations, the…
We propose an accurate numerical scheme for approximating the solution of the two dimensional acoustic wave problem. We use machine learning to find a stencil suitable even in the presence of high wavenumbers. The proposed scheme…
In this paper, we propose and analyze a second order accurate (in both time and space) numerical scheme for the Poisson-Nernst-Planck-Navier-Stokes system, which describes the ion electro-diffusion in fluids. In particular, the…