Related papers: Complete monotonicity-preserving numerical methods…
We develop a rapid and accurate contour method for the solution of time-fractional PDEs. The method inverts the Laplace transform via an optimised stable quadrature rule, suitable for infinite-dimensional operators, whose error decreases…
A conforming finite element scheme with mixed explicit-implicit time discretization for quasi-incompressible Navier-Stokes-Maxwell-Stefan systems in a bounded domain with periodic boundary conditions is presented. The system consists of the…
In this paper we study the stability of explicit finite difference discretizations of linear advection-diffusion equations (ADE) with arbitrary order of accuracy in the context of method of lines. The analysis first focuses on the stability…
Ordinary differential equations (ODEs) and ordinary difference systems (O$\Delta$Ss) invariant under the actions of the Lie groups $\mathrm{SL}_x(2)$, $\mathrm{SL}_y(2)$ and $\mathrm{SL}_x(2)\times\mathrm{SL}_y(2)$ of projective…
We consider linear scalar wave equations with a hereditary integral term of the kind used to model viscoelastic solids. The kernel in this Volterra integral is a sum of decaying exponentials (The so-called Maxwell, or Zener model) and this…
Many Hamiltonian systems can be recast in multi-symplectic form. We develop a reduced-order model (ROM) for multi-symplectic Hamiltonian partial differential equations (PDEs) that preserves the global energy. The full-order solutions are…
The method and the formula of variation of constants for ordinary differential equations (ODEs) is a fundamental tool to analyze the dynamics of an ODE near an equilibrium. It is natural to expect that such a formula works for delay…
Computationally efficient, structure-preserving reduced-order methods are developed for the Korteweg de Vries (KdV) equations in Hamiltonian form. The KdV equation is discretized in space by finite differences. The resulting skew-gradient…
We investigate a local modification of a variable-order fractional wave equation, which describes the propagation of diffusive wave in viscoelastic media with evolving physical property. We incorporate an equivalent formulation to prove the…
We propose and analyse a boundary-preserving numerical scheme for the weak approximation for some stochastic partial differential equations (SPDEs) with bounded state-space. We impose regularity assumptions on the drift and diffusion…
We propose some finite element schemes to solve a class of fourth-order nonlinear PDEs, which include the vector-valued Landau--Lifshitz--Baryakhtar equation, the Swift--Hohenberg equation, and various Cahn--Hilliard-type equations with…
A new explicit stochastic scheme of order 1 is proposed for solving commutative stochastic differential equations (SDEs) with non-globally Lipschitz continuous coefficients. The proposed method is a semi-tamed version of Milstein scheme to…
This papers deals with a construction and convergence analysis of a finite difference scheme for solving time-fractional porous medium equation. The governing equation exhibits both nonlocal and nonlinear behaviour making the numerical…
In this paper, we study 1D autonomous fractional ODEs $D_c^{\gamma}u=f(u), 0< \gamma <1$, where $u: [0,\infty)\mapsto\mathbb{R}$ is the unknown function and $D_c^{\gamma}$ is the generalized Caputo derivative introduced by Li and Liu (…
We compare the performance of several discretizations of the simple pendulum equation in a series of numerical experiments. The stress is put on the long-time behaviour. We choose for the comparison numerical schemes which preserve the…
We develop a general framework for numerically solving differential equations while preserving invariants. As in standard projection methods, we project an arbitrary base integrator onto an invariant-preserving manifold, however, our method…
Compatible Discrete Operator schemes preserve basic properties of the continuous model at the discrete level. They combine discrete differential operators that discretize exactly topological laws and discrete Hodge operators that…
In this paper, we develop a second-order accurate time-stepping scheme for the tempered time-fractional advection-dispersion equation based on a sum-of-exponentials (SOE) approximation to the convolution kernel involved in the fractional…
We develop a unified and easy to use framework to study robust fully discrete numerical methods for nonlinear degenerate diffusion equations $$ \partial_t u-\mathfrak{L}^{\sigma,\mu}[\varphi(u)]=f \quad\quad\text{in}\quad\quad…
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…