[HTML][HTML] Boundary-optimized summation-by-parts operators for finite difference approximations of second derivatives with variable coefficients

V Stiernström, M Almquist, K Mattsson - Journal of Computational Physics, 2023 - Elsevier
Boundary-optimized summation-by-parts (SBP) finite difference operators for second
derivatives with variable coefficients are presented. The operators achieve increased …

[HTML][HTML] A residual-based artificial viscosity finite difference method for scalar conservation laws

V Stiernström, L Lundgren, M Nazarov… - Journal of Computational …, 2021 - Elsevier
In this paper, we present an accurate, stable and robust shock-capturing finite difference
method for solving scalar non-linear conservation laws. The spatial discretization uses high …

Infrasound radiation from impulsive volcanic eruptions: Nonlinear aeroacoustic 2D simulations

LM Watson, EM Dunham, D Mohaddes… - Journal of …, 2021 - Wiley Online Library
Infrasound observations are increasingly used to constrain properties of volcanic eruptions.
In order to better interpret infrasound observations, however, there is a need to better …

Modeling and inversion in acoustic-elastic coupled media using energy-stable summation-by-parts operators

M Bader, M Almquist, EM Dunham - Geophysics, 2023 - library.seg.org
Seismic data acquired at the seafloor are valuable in characterizing the subsurface and
monitoring producing hydrocarbon fields. To fully use such data, a stable, accurate, and …

[HTML][HTML] Adjoint-based inversion for stress and frictional parameters in earthquake modeling

V Stiernström, M Almquist, EM Dunham - Journal of Computational Physics, 2024 - Elsevier
We present an adjoint-based optimization method to invert for stress and frictional
parameters used in earthquake modeling. The forward problem is linear elastodynamics …

Elastic wave propagation in anisotropic solids using energy-stable finite differences with weakly enforced boundary and interface conditions

M Almquist, EM Dunham - Journal of Computational Physics, 2021 - Elsevier
Abstract Summation-by-parts (SBP) finite difference methods have several desirable
properties for second-order wave equations. They combine the computational efficiency of …

A high‐order accurate summation‐by‐parts finite difference method for fully‐dynamic earthquake sequence simulations within sedimentary basins

TW Harvey, BA Erickson… - Journal of Geophysical …, 2023 - Wiley Online Library
We present an efficient numerical method for earthquake sequences in 2D antiplane shear
that incorporates wave propagation. A vertical strike‐slip fault governed by rate‐and‐state …

Simulation of flexural-gravity wave propagation for elastic plates in shallow water using an energy-stable finite difference method with weakly enforced boundary and …

N Tazhimbetov, M Almquist, J Werpers… - Journal of Computational …, 2023 - Elsevier
We introduce an energy stable, high-order-accurate finite difference approximation of the
dynamic, pure bending Kirchhoff plate equations for complex geometries and spatially …

Programming approaches for scalability, performance, and portability of combustion physics codes

KA Schau, D Purushotham, JC Oefelein - Proceedings of the Combustion …, 2023 - Elsevier
This paper presents the process, strategy, and results associated with porting a typical
combustion physics flow solver to current state-of-the-art and future massively-parallel …

[HTML][HTML] A hybrid numerical method for elastic wave propagation in discontinuous media with complex geometry

A Granath, S Wang - Journal of Computational Physics, 2025 - Elsevier
We develop a high order accurate numerical method for solving the elastic wave equation in
second-order form. We hybridize the computationally efficient Cartesian grid formulation of …