Related papers: Accurate simulations of nonlinear dynamic shear ru…
We present an energy/entropy stable and high order accurate finite difference (FD) method for solving the nonlinear (rotating) shallow water equations (SWEs) in vector invariant form using the newly developed dual-pairing and…
High-order accurate summation-by-parts (SBP) finite difference (FD) methods constitute efficient numerical methods for simulating large-scale hyperbolic wave propagation problems. Traditional SBP FD operators that approximate first-order…
A recent study by Gassner et al. [J. Sci. Comput. 90:79 (2022)] demonstrates that local energy stability--that is, ensuring the asymptotic numerical growth rate does not exceed the continuous growth rate--is crucial for achieving accurate…
Robust and convergent high-order numerical methods for solving partial differential equations are highly attractive due to their efficiency on modern and next-generation hardware architectures. However, designing such methods for nonlinear…
The Finite Difference (FD) and the Spectral Boundary Integral (SBI) methods have been used extensively to model spontaneously propagating shear cracks in a variety of engineering and geophysical applications. In this paper, we propose a new…
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…
We study numerical methods for dissipative particle dynamics (DPD), which is a system of stochastic differential equations and a popular stochastic momentum-conserving thermostat for simulating complex hydrodynamic behavior at mesoscales.…
This paper focuses on the denoising and enhancing of 3-D reflection seismic data. We propose a pre-processing step based on a non linear diffusion filtering leading to a better detection of seismic faults. The non linear diffusion…
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…
Earthquake fault zones are more complex, both geometrically and rheologically, than an idealised infinitely thin plane embedded in linear elastic material. To incorporate nonlinear material behaviour, natural complexities, and multi-physics…
Simulating long-term, fully dynamic sequences of earthquakes and aseismic slip (SEAS) on geometrically complex fault networks remains computationally demanding due to the cost of resolving elastodynamic interactions. Although…
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…
Flexible elastic structures, such as beams, rods, ribbons, plates, and shells, exhibit complex nonlinear dynamical behaviors that are central to a wide range of engineering and scientific applications, including soft robotics, deployable…
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…
Shear dilation based hydraulic stimulations enable exploitation of geothermal energy from reservoirs with inadequate initial permeability. While contributing to enhancing the reservoir's permeability, hydraulic stimulation processes may…
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.…
Simulating dynamic rupture propagation is challenging due to the uncertainties involved in the underlying physics of fault slip, stress conditions, and frictional properties of the fault. A trial and error approach is often used to…
The fracture simulation of random particle reinforced composite structures remains a challenge. Current techniques either assumed a homogeneous model, ignoring the microstructure characteristics of composite structures, or considered a…
High order entropy stable discontinuous Galerkin (DG) methods for nonlinear conservation laws satisfy an inherent discrete entropy inequality. The construction of such schemes has relied on the use of carefully chosen nodal points or volume…