Related papers: Implicit finite difference schemes for the magneti…
The theory part of this paper is sketched as follows. Based on column stochastic average matrix $T_n$ selected as a basic substitution matrix, the method of advanced successive difference substitution is established. Then, a set of…
Finite difference schemes, using Backward Differentiation Formula (BDF), are studied for the approximation of one-dimensional diffusion equations with an obstacle term, of the form $$\min(v_t - a(t,x) v_{xx} + b(t,x) v_x + r(t,x) v, v-…
We present difference schemes for stochastic transport equations with low-regularity velocity fields. We establish $L^2$ stability and convergence of the difference approximations under conditions that are less strict than those required…
We are developing a framework for multiscale computation which enables models at a ``microscopic'' level of description, for example Lattice Boltzmann, Monte Carlo or Molecular Dynamics simulators, to perform modelling tasks at the…
We use high order finite difference methods to solve the wave equation in the second order form. The spatial discretization is performed by finite difference operators satisfying a summation-by-parts property. The focus of this work is on…
We demonstrate that we can carry over the strategy of Finite Element Exterior Calculus (FEEC) to Summation-by-Parts (SBP) Finite Difference (FD) methods to achieve divergence- and curl-free discretizations. This is not obvious at first…
In this paper, we present a numerical strategy to check the strong stability (or GKS-stability) of one-step explicit finite difference schemes for the one-dimensional advection equation with an inflow boundary condition. The strong…
In this paper, we will use the interior functions of an hierarchical basis for high order $BDM_p$ elements to enforce the divergence-free condition of a magnetic field $B$ approximated by the H(div) $BDM_p$ basis. The resulting constrained…
We consider a system of coupled free boundary problems for pricing American put options with regime-switching. To solve this system, we first employ the logarithmic transformation to map the free boundary for each regime to multi-fixed…
High-order spatial discretisations and full discretisations of parabolic partial differential equations on evolving surfaces are studied. We prove convergence of the high-order evolving surface finite element method, by showing high-order…
A Finite Volume (FV) scheme is developed for solving the extended magnetohydrodynamic (XMHD) equations, yielding accurate results in the ideal, resistive, and Hall MHD limits. This is accomplished by first re-writing the XMHD equations such…
A simple and reliable finite difference approach is presented for solution of the Dirac equation eigenproblem for states confined in rotationally symmetric systems. The method sets the boundary condition for the spinor wave function…
A scheme is presented for accurately propagating the gravitational field constraints in finite difference implementations of numerical relativity. The method is based on similar techniques used in astrophysical magnetohydrodynamics and…
To enhance the scalability and performance of the traditional finite-difference time-domain (FDTD) methods, a three-dimensional summation-by-parts simultaneous approximation term (SBP-SAT) FDTD method is developed to solve complex…
The aim of this work is to apply a semi-implicit (SI) strategy within a Rosenbrock-type and IMEX linear multistep (LM) framework to a sequence of 1D time-dependent partial differential equations (PDEs) with high order spatial derivatives.…
We introduce a fluid dynamics algorithm that performs with nearly spectral accuracy, but uses finite-differences instead of FFTs to compute gradients and thus executes 10 times faster. The finite differencing is not based on a high-order…
In this paper the numerical approximation of solutions of Liouville-Master Equations for time-dependent distribution functions of Piecewise Deterministic Processes with memory is considered. These equations are linear hyperbolic PDEs with…
Summation-by-parts (SBP) finite-difference discretizations share many attractive properties with Galerkin finite-element methods (FEMs), including time stability and superconvergent functionals; however, unlike FEMs, SBP operators are not…
An initial-boundary value problem for the 1D self-adjoint parabolic equation on the half-axis is solved. We study a broad family of two-level finite-difference schemes with two parameters related to averagings both in time and space.…
We present a hybridization technique for summation-by-parts finite difference methods with weak enforcement of interface and boundary conditions for second order, linear elliptic partial differential equations. The method is based on…