Related papers: Higher order A-stable schemes for the wave equatio…
This paper is concerned with developing and analyzing two novel implicit temporal discretization methods for the stochastic semilinear wave equations with multiplicative noise. The proposed methods are natural extensions of well-known…
Direct methods to obtain global stability modes are restricted by the daunting sizes and complexity of Jacobians encountered in general three-dimensional flows. Jacobian-free iterative approaches such as Arnoldi methods have greatly…
A system of high-order adaptive multiresolution wavelet collocation upwind schemes are developed for the solution of hyperbolic conservation laws. A couple of asymmetrical wavelet bases with interpolation property are built to realize the…
In this paper, we propose an adaptive high-order method for hyperbolic systems of conservation laws. The proposed method is based on a dual formulation approach: Two numerical solutions, corresponding to conservative and nonconservative…
The Fokker-Planck (FP) model is one of the commonly used methods for studies of the dynamical evolution of dense spherical stellar systems such as globular clusters and galactic nuclei. The FP model is numerically stable in most cases, but…
Certain solutions of autonomous PDEs without any boundary conditions describing the spatiotemporal evolution of a dependent variable in an unbounded spatial domain can be characterised as a travelling wave moving with constant speed. In the…
A new method is formulated and analyzed for the approximate solution of a two-dimensional time-fractional diffusion-wave equation. In this method, orthogonal spline collocation is used for the spatial discretization and, for the…
We present fast, spatially dispersionless and unconditionally stable high-order solvers for Partial Differential Equations (PDEs) with variable coefficients in general smooth domains. Our solvers, which are based on (i) A certain "Fourier…
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…
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…
This paper introduces an adaptive time splitting technique for the solution of stiff evolutionary PDEs that guarantees an effective error control of the simulation, independent of the fastest physical time scale for highly unsteady…
We propose a solution for linear inverse problems based on higher-order Langevin diffusion. More precisely, we propose pre-conditioned second-order and third-order Langevin dynamics that provably sample from the posterior distribution of…
Efficient and accurate numerical simulation of 3D acoustic wave propagation in heterogeneous media plays an important role in the success of seismic full waveform inversion (FWI) problem. In this work, we employed the combined scheme and…
In this paper, a compact alternating direction implicit (ADI) finite difference scheme for the two-dimensional time fractional diffusion-wave equation is developed, with temporal and spatial accuracy order equal to two and four…
The manuscript presents a new technique for computing the exponential of skew-Hermitian operators. Principal advantages of the proposed method include: stability even for large time-steps, the possibility to parallelize in time over many…
In this paper, we consider the variable-order time fractional mobile/immobile diffusion (TF-MID) equation in two-dimensional spatial domain, where the fractional order $\alpha(t)$ satisfies $0<\alpha_{*}\leq \alpha(t)\leq \alpha^{*}<1$. We…
Steady state simulations} of magnetized electron fluid equations with strong anisotropic diffusion based on the first-order hyperbolic approach is carried out using cell-centered higher order upwind schemes, linear and weighted essentially…
The propagation of electromagnetic waves in general media is modeled by the time-dependent Maxwell's partial differential equations (PDEs), coupled with constitutive laws that describe the response of the media. In this work, we focus on…
We introduce a class of high order accurate, semi-implicit Runge-Kutta schemes in the general setting of evolution equations that arise as gradient flow for a cost function, possibly with respect to an inner product that depends on the…
We consider a space-time variational formulation of the second-order wave equation, where integration by parts is also applied with respect to the time variable. Conforming tensor-product finite element discretisations with piecewise…