Related papers: Upwind summation by parts finite difference method…
We construct new, efficient, and accurate high-order finite differencing operators which satisfy summation by parts. Since these operators are not uniquely defined, we consider several optimization criteria: minimizing the bandwidth, the…
We consider the numerical simulation of the acoustic wave equations arising from seismic applications, for which staggered grid finite difference methods are popular choices due to their simplicity and efficiency. We relax the uniform grid…
In this paper, we design high order accurate and stable finite difference schemes for the initial-boundary value problem, associated with the magnetic induction equation with resistivity. We use Summation-By-Parts (SBP) finite difference…
We present a new class of efficient and robust discontinuous spectral-element methods of arbitrary order for nonlinear hyperbolic systems of conservation laws on curved triangular and tetrahedral unstructured grids. Such discretizations…
We employ the summation-by-parts (SBP) framework to extend the recent domain-of-dependence (DoD) stabilization for cut cells to linear kinetic models in diffusion scaling. Numerical methods for these models are challenged by increased…
Summation-by-parts (SBP) operators are finite-difference operators that mimic integration by parts. This property can be useful in constructing energy-stable discretizations of partial differential vequations. SBP operators are defined by a…
High-order entropy stable summation-by-parts (SBP) schemes are a class of robust and accurate numerical methods for hyperbolic conservation laws that are numerically stable at arbitrary order without the need for artificial stabilization.…
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…
Curvilinear, multiblock summation-by-parts finite difference operators with the simultaneous approximation term method provide a stable and accurate framework for solving the wave equation in second order form. That said, the standard…
Context: Calculating stellar pulsations requires a sufficient accuracy to match the quality of the observations. Many current pulsation codes apply a second order finite-difference scheme, combined with Richardson extrapolation to reach…
There has been an increasing interest in developing efficient immersed boundary method (IBM) based on Cartesian grids, recently in the context of high-order methods. IBM based on volume penalization is a robust and easy to implement method…
Non-conforming numerical approximations offer increased flexibility for applications that require high resolution in a localized area of the computational domain or near complex geometries. Two key properties for non-conforming methods to…
Finite difference schemes for the simulation of elastic waves in materi- als with jump discontinuities are presented. The key feature is the highly accurate treatment of interfaces where media discontinuities arise. The schemes are…
We investigate numerical methods for wave equations in $n+2$ spacetime dimensions, written in spherical coordinates, decomposed in spherical harmonics on $S^n$, and finite-differenced in the remaining coordinates $r$ and $t$. Such an…
Many applications rely on solving time-dependent partial differential equations (PDEs) that include second derivatives. Summation-by-parts (SBP) operators are crucial for developing stable, high-order accurate numerical methodologies for…
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…
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 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, we consider finite difference approximations of the second order wave equation. We use finite difference operators satisfying the summation-by-parts property to discretize the equation in space. Boundary conditions and grid…
We discuss the order, efficiency, stability and positivity of several meshless schemes for linear scalar hyperbolic equations. Meshless schemes are Generalised Finite Difference Methods (GFDMs) for arbitrary irregular grids in which there…