Related papers: Numerical treatment of interfaces for second-order…
Current spectral simulations of Einstein's equations require writing the equations in first-order form, potentially introducing instabilities and inefficiencies. We present a new penalty method for pseudo-spectral evolutions of second order…
Synchronizations of processing elements (PEs) in massively parallel simulations, which arise due to communication or load imbalances between PEs, significantly affect the scalability of scientific applications. We have recently proposed a…
Explicit Runge-Kutta schemes with large stable step sizes are developed for integration of high order spectral difference spatial discretization on quadrilateral grids. The new schemes permit an effective time step that is substantially…
Obtaining exact solutions to the Schr\"odinger equation in complex quantum systems poses significant challenges. In this context, numerical methods emerge as valuable tools for analyzing such systems. This article proposes a numerical…
In this paper, we develop a high order finite difference boundary treatment method for the implicit-explicit (IMEX) Runge-Kutta (RK) schemes solving hyperbolic systems with possibly stiff source terms on a Cartesian mesh. The main challenge…
We study a class of general purpose linear multisymplectic integrators for Hamiltonian wave equations based on a diamond-shaped mesh. On each diamond, the PDE is discretized by a symplectic Runge--Kutta method. The scheme advances in time…
A numerical method is described for studying how elastic waves interact with imperfect contacts such as fractures or glue layers existing between elastic solids. These contacts have been classicaly modeled by interfaces, using a simple…
In this paper, a class of finite difference numerical techniques is presented to solve the second-order linear inhomogeneous damped wave equation. The consistency, stability, and convergences of these numerical schemes are discussed. The…
In this paper we present two semi-implicit-type second order Compact Approximate Taylor (CAT2) numerical schemes and blend them with a local a posteriori Multi-dimensional Optimal Order Detection (MOOD) paradigm to solve hyperbolic systems…
We present a new method for two-material Lagrangian hydrodynamics, which combines the Shifted Interface Method (SIM) with a high-order Finite Element Method. Our approach relies on an exact (or sharp) material interface representation, that…
A numerical scheme is described for accurately accommodating oblique, non-aligned, boundaries, on a three-dimensional cartesian grid. The scheme gives second-order accuracy in the solution for potential of Poisson's equation using compact…
In this paper, a symmetrized two-scale finite element method is proposed for a class of partial differential equations with symmetric solutions. With this method, the finite element approximation on a fine tensor product grid is reduced to…
Implicit Runge--Kutta (IRK) methods are highly effective for solving stiff ordinary differential equations (ODEs) but can be computationally expensive for large-scale problems due to the need of solving coupled algebraic equations at each…
In this work we present an extension of the Virtual Element Method with curved edges for the numerical approximation of the second order wave equation in a bidimensional setting. Curved elements are used to describe the domain boundary, as…
In this paper, we develop a higher order symmetric partitioned Runge-Kutta method for a coupled system of differential equations on Lie groups. We start with a discussion on partitioned Runge-Kutta methods on Lie groups of arbitrary order.…
We consider the finite difference discretization of isotropic elastic wave equations on nonuniform grids. The intended applications are seismic studies, where heterogeneity of the earth media can lead to severe oversampling for simulations…
In this paper, we develop a computational multiscale to solve the parabolic wave approximation with heterogeneous and variable media. Parabolic wave approximation is a technique to approximate the full wave equation. One benefit of the…
We present two approaches for enhancing the accuracy of second order finite difference approximations of two-dimensional semilinear parabolic systems. These are the fourth order compact difference scheme and the fourth order scheme based on…
The goal of the present paper is to understand the impact of numerical schemes for the reconstruction of data at cell faces in finite-volume methods, and to assess their interaction with the quadrature rule used to compute the average over…
Summation-by-parts (SBP) finite difference methods have several desirable properties for second-order wave equations. They combine the computational efficiency of narrow-stencil finite difference operators with provable stability on…