Related papers: Study Notes on Numerical Solutions of the Wave Equ…
We present an efficient second-order finite difference scheme for solving the 2D sine-Gordon equation, which can inherit the discrete energy conservation for the undamped model theoretically. Due to the semi-implicit treatment for the…
We consider the numerical solution of partial differential equations with coefficients that are strongly heterogeneous in space. We provide an overview of higher-order localized orthogonal decomposition (LOD) methods for the elliptic…
To avoid instabilities in the continuum semi-classical limit of loop quantum cosmology models, refinement of the underlying lattice is necessary. The lattice refinement leads to new dynamical difference equations which, in general, do not…
In this work, we present a semi-numerical solution of a fractal telegraphic dual-porosity fluid flow model. It combines Laplace transform and finite difference schemes. The Laplace transform handles the time variable whereas the finite…
This paper presents a stable finite element approximation for the acoustic wave equation on second-order form, with perfectly matched layers (PML) at the boundaries. Energy estimates are derived for varying PML damping for both the discrete…
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…
This paper develops a robust three-level time split high-order Leapfrog/Crank-Nicolson technique for solving the two-dimensional unsteady sobolev and regularized long wave equations arising in fluid mechanics. A deep analysis of the…
Due to its highly oscillating solution, the Helmholtz equation is numerically challenging to solve. To obtain a reasonable solution, a mesh size that is much smaller than the reciprocal of the wavenumber is typically required (known as the…
The 3-level leapfrog time integration algorithm is an attractive choice for numerical relativity simulations since it is time-symmetric and avoids non-physical damping. In Newtonian problems without velocity dependent forces, this method…
The present study concerns the numerical homogenization of second order hyperbolic equations in non-divergence form, where the model problem includes a rapidly oscillating coefficient function. These small scales influence the large scale…
The analysis of wave propagation problems in linear damped media must take into account both propagation features and attenuation process. To perform accurate numerical investigations by the finite differences or finite element method, one…
A new modified Galerkin / Finite Element Method is proposed for the numerical solution of the fully nonlinear shallow water wave equations. The new numerical method allows the use of low-order Lagrange finite element spaces, despite the…
Stiff and chaotic differential equations are challenging for time-stepping numerical methods. For explicit methods, the required time step resolution significantly exceeds the resolution associated with the smoothness of the exact solution…
We extend and analyze the energy-based discontinuous Galerkin method for second order wave equations on staggered and structured meshes. By combining spatial staggering with local time-stepping near boundaries, the method overcomes the…
A stochastic leap-frog algorithm for the numerical integration of Brownian motion stochastic differential equations with multiplicative noise is proposed and tested. The algorithm has a second-order convergence of moments in a finite time…
We propose two stable and one conditionally stable finite difference schemes of second-order in both time and space for the time-fractional diffusion-wave equation. In the first scheme, we apply the fractional trapezoidal rule in time and…
This paper introduces filtered finite difference methods for numerically solving a dispersive evolution equation with solutions that are highly oscillatory in both space and time. We consider a semiclassically scaled nonlinear Schr\"odinger…
Langevin simulation provides an effective way to study collisional effects in beams by reducing the six-dimensional Fokker-Planck equation to a group of stochastic ordinary differential equations. These resulting equations usually have…
We begin with a treatment of the Caputo time-fractional diffusion equation, by using the Laplace transform, to obtain a Volterra intego-differential equation where we may examine the weakly singular nature of this convolution…
In this chapter, we aim at presenting the basic techniques necessary to go beyond the widely accepted paradigm of second-order numerics. We specifically focus on finite-volume schemes for hyperbolic conservation laws occuring in fluid…