Simulations of seismic wave propagation using a spectral-element method in a Lagrangian framework with logarithmic strain

Uno Vaaland, Hom Nath Gharti, Jeroen Tromp

Research output: Contribution to journalArticle

Abstract

We show how the linearized equations that govern the motion of a body that undergoes deformation can be generalized to capture geometrical non-linearities in a spectral-element formulation. Generalizing the equations adds little complexity, the main addition being that we have to track the deformation gradient. Geometrical changes due to deformation are captured using the logarithmic strain. We test the geometrically non-linear formulation by considering numerical experiments in seismic wave propagation and cantilever beam bending and compare the results with the linearized formulation. In cases where finite deformation occurs, the effect of solving the geometrically non-linear equations can be significant while in cases where deformation is smaller the result is similar to solving the linearized equations. We find that the time it takes to run the geometrically non-linear simulations is on the same order of magnitude as running the linearized simulations. The limited amount of added cost and complexity suggests that we might as well solve the geometrically non-linear equations since it does not assume anything about the size of deformations.

Original languageEnglish (US)
Pages (from-to)2148-2157
Number of pages10
JournalGeophysical Journal International
Volume216
Issue number3
DOIs
StatePublished - Mar 1 2019

Fingerprint

Seismic waves
seismic waves
seismic wave
Wave propagation
wave propagation
simulation
formulations
Nonlinear equations
nonlinear equations
cantilever beams
Cantilever beams
nonlinearity
method
costs
gradients
cost
Costs
experiment
Experiments

All Science Journal Classification (ASJC) codes

  • Geochemistry and Petrology
  • Geophysics

Cite this

@article{5c2620fc1f414711853ba4e5a4ef4d7e,
title = "Simulations of seismic wave propagation using a spectral-element method in a Lagrangian framework with logarithmic strain",
abstract = "We show how the linearized equations that govern the motion of a body that undergoes deformation can be generalized to capture geometrical non-linearities in a spectral-element formulation. Generalizing the equations adds little complexity, the main addition being that we have to track the deformation gradient. Geometrical changes due to deformation are captured using the logarithmic strain. We test the geometrically non-linear formulation by considering numerical experiments in seismic wave propagation and cantilever beam bending and compare the results with the linearized formulation. In cases where finite deformation occurs, the effect of solving the geometrically non-linear equations can be significant while in cases where deformation is smaller the result is similar to solving the linearized equations. We find that the time it takes to run the geometrically non-linear simulations is on the same order of magnitude as running the linearized simulations. The limited amount of added cost and complexity suggests that we might as well solve the geometrically non-linear equations since it does not assume anything about the size of deformations.",
author = "Uno Vaaland and Gharti, {Hom Nath} and Jeroen Tromp",
year = "2019",
month = "3",
day = "1",
doi = "https://doi.org/10.1093/gji/ggy546",
language = "English (US)",
volume = "216",
pages = "2148--2157",
journal = "Geophysical Journal International",
issn = "0956-540X",
publisher = "Wiley-Blackwell",
number = "3",

}

Simulations of seismic wave propagation using a spectral-element method in a Lagrangian framework with logarithmic strain. / Vaaland, Uno; Gharti, Hom Nath; Tromp, Jeroen.

In: Geophysical Journal International, Vol. 216, No. 3, 01.03.2019, p. 2148-2157.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Simulations of seismic wave propagation using a spectral-element method in a Lagrangian framework with logarithmic strain

AU - Vaaland, Uno

AU - Gharti, Hom Nath

AU - Tromp, Jeroen

PY - 2019/3/1

Y1 - 2019/3/1

N2 - We show how the linearized equations that govern the motion of a body that undergoes deformation can be generalized to capture geometrical non-linearities in a spectral-element formulation. Generalizing the equations adds little complexity, the main addition being that we have to track the deformation gradient. Geometrical changes due to deformation are captured using the logarithmic strain. We test the geometrically non-linear formulation by considering numerical experiments in seismic wave propagation and cantilever beam bending and compare the results with the linearized formulation. In cases where finite deformation occurs, the effect of solving the geometrically non-linear equations can be significant while in cases where deformation is smaller the result is similar to solving the linearized equations. We find that the time it takes to run the geometrically non-linear simulations is on the same order of magnitude as running the linearized simulations. The limited amount of added cost and complexity suggests that we might as well solve the geometrically non-linear equations since it does not assume anything about the size of deformations.

AB - We show how the linearized equations that govern the motion of a body that undergoes deformation can be generalized to capture geometrical non-linearities in a spectral-element formulation. Generalizing the equations adds little complexity, the main addition being that we have to track the deformation gradient. Geometrical changes due to deformation are captured using the logarithmic strain. We test the geometrically non-linear formulation by considering numerical experiments in seismic wave propagation and cantilever beam bending and compare the results with the linearized formulation. In cases where finite deformation occurs, the effect of solving the geometrically non-linear equations can be significant while in cases where deformation is smaller the result is similar to solving the linearized equations. We find that the time it takes to run the geometrically non-linear simulations is on the same order of magnitude as running the linearized simulations. The limited amount of added cost and complexity suggests that we might as well solve the geometrically non-linear equations since it does not assume anything about the size of deformations.

UR - http://www.scopus.com/inward/record.url?scp=85062236464&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85062236464&partnerID=8YFLogxK

U2 - https://doi.org/10.1093/gji/ggy546

DO - https://doi.org/10.1093/gji/ggy546

M3 - Article

VL - 216

SP - 2148

EP - 2157

JO - Geophysical Journal International

JF - Geophysical Journal International

SN - 0956-540X

IS - 3

ER -