DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: An efficient time-domain perfectly matched layers formulation for elastodynamics on spherical domains

Abstract

Many practical applications require the analysis of elastic wave propagation in a homogeneous isotropic media in an unbounded domain. One widely used approach for truncating the infinite domain is the so called method of perfectly matched layers (PMLs). Most existing PML formulations are developed for finite difference methods based on the first-order velocity-stress form of the elasticity equations, and they are not straight-forward to implement using standard finite element methods (FEMs) on unstructured meshes. Some of the problems with these formulations include the application of boundary conditions in half-space problems and in the treatment of edges and/or corners for time-domain problems. Several PML formulations, which do work with FEMs have been proposed, although most of them still have some of these problems and/or they require a large number of auxiliary nodal history/memory variables. Here in this work, we develop a new PML formulation for time-domain elastodynamics on a spherical domain, which reduces to a two-dimensional formulation under the assumption of axisymmetry. Our formulation is well-suited for implementation using FEMs, where it requires lower memory than existing formulations, and it allows for natural application of boundary conditions. We solve example problems on two-dimensional and three-dimensional domains using a high-order discontinuous Galerkinmore » (DG) discretization on unstructured meshes and explicit time-stepping. We also study an approach for stabilization of the discrete equations, and we show several practical applications for quality factor predictions of micromechanical resonators along with verifying the accuracy and versatility of our formulation.« less

Authors:
 [1];  [1];  [1]
  1. Univ. of California, Berkeley, CA (United States). Dept. of Civil Engineering
Publication Date:
Research Org.:
Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States)
Sponsoring Org.:
USDOE Office of Science (SC); National Science Foundation (NSF)
OSTI Identifier:
1523916
Grant/Contract Number:  
AC02-05CH11231; FA9550-10-1- 0229
Resource Type:
Accepted Manuscript
Journal Name:
International Journal for Numerical Methods in Engineering
Additional Journal Information:
Journal Volume: 100; Journal Issue: 6; Journal ID: ISSN 0029-5981
Publisher:
Wiley
Country of Publication:
United States
Language:
English
Subject:
42 ENGINEERING; PML; perfectly matched layer; elastodynamics; Discontinuous Galerkin; stabilization

Citation Formats

Sagiyama, K., Govindjee, S., and Persson, P. -O. An efficient time-domain perfectly matched layers formulation for elastodynamics on spherical domains. United States: N. p., 2014. Web. doi:10.1002/nme.4740.
Sagiyama, K., Govindjee, S., & Persson, P. -O. An efficient time-domain perfectly matched layers formulation for elastodynamics on spherical domains. United States. https://doi.org/10.1002/nme.4740
Sagiyama, K., Govindjee, S., and Persson, P. -O. Tue . "An efficient time-domain perfectly matched layers formulation for elastodynamics on spherical domains". United States. https://doi.org/10.1002/nme.4740. https://www.osti.gov/servlets/purl/1523916.
@article{osti_1523916,
title = {An efficient time-domain perfectly matched layers formulation for elastodynamics on spherical domains},
author = {Sagiyama, K. and Govindjee, S. and Persson, P. -O.},
abstractNote = {Many practical applications require the analysis of elastic wave propagation in a homogeneous isotropic media in an unbounded domain. One widely used approach for truncating the infinite domain is the so called method of perfectly matched layers (PMLs). Most existing PML formulations are developed for finite difference methods based on the first-order velocity-stress form of the elasticity equations, and they are not straight-forward to implement using standard finite element methods (FEMs) on unstructured meshes. Some of the problems with these formulations include the application of boundary conditions in half-space problems and in the treatment of edges and/or corners for time-domain problems. Several PML formulations, which do work with FEMs have been proposed, although most of them still have some of these problems and/or they require a large number of auxiliary nodal history/memory variables. Here in this work, we develop a new PML formulation for time-domain elastodynamics on a spherical domain, which reduces to a two-dimensional formulation under the assumption of axisymmetry. Our formulation is well-suited for implementation using FEMs, where it requires lower memory than existing formulations, and it allows for natural application of boundary conditions. We solve example problems on two-dimensional and three-dimensional domains using a high-order discontinuous Galerkin (DG) discretization on unstructured meshes and explicit time-stepping. We also study an approach for stabilization of the discrete equations, and we show several practical applications for quality factor predictions of micromechanical resonators along with verifying the accuracy and versatility of our formulation.},
doi = {10.1002/nme.4740},
journal = {International Journal for Numerical Methods in Engineering},
number = 6,
volume = 100,
place = {United States},
year = {Tue Jul 15 00:00:00 EDT 2014},
month = {Tue Jul 15 00:00:00 EDT 2014}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 5 works
Citation information provided by
Web of Science

Figures / Tables:

Figure 1 Figure 1: (a) Problem setup for a stability study. (b) Plots of computed $\bar{u}$$R$ versus time on a semi-log scale.

Save / Share:

Works referenced in this record:

Mixed Spectral Finite Elements for the Linear Elasticity System in Unbounded Domains
journal, January 2005


A Simple Mesh Generator in MATLAB
journal, January 2004


A perfectly matched layer for the absorption of electromagnetic waves
journal, October 1994


Perfectly matched layers for elastic waves in cylindrical and spherical coordinates
journal, April 1999

  • Liu, Q. H.
  • The Journal of the Acoustical Society of America, Vol. 105, Issue 4
  • DOI: 10.1121/1.426812

The Newmark scheme as velocity-stress time-staggering: an efficient PML implementation for spectral element simulations of elastodynamics
journal, June 2005


MEMS technology for timing and frequency control
journal, January 2007

  • Nguyen, Clark
  • IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, Vol. 54, Issue 2, p. 251-270
  • DOI: 10.1109/TUFFC.2007.240

A time-domain Discontinuous Galerkin method for mechanical resonator quality factor computations
journal, August 2012


Mixed perfectly-matched-layers for direct transient analysis in 2D elastic heterogeneous media
journal, January 2011

  • Kucukcoban, S.; Kallivokas, L. F.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 200, Issue 1-4
  • DOI: 10.1016/j.cma.2010.07.013

Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations
journal, September 1982

  • Brooks, Alexander N.; Hughes, Thomas J. R.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 32, Issue 1-3
  • DOI: 10.1016/0045-7825(82)90071-8

Explicit finite element perfectly matched layer for transient three-dimensional elastic waves
journal, January 2009

  • Basu, Ushnish
  • International Journal for Numerical Methods in Engineering, Vol. 77, Issue 2
  • DOI: 10.1002/nme.2397

On the implementation of perfectly matched layers in a three-dimensional fourth-order velocity-stress finite difference scheme: IMPLEMENTING PERFECTLY MATCHED LAYERS
journal, May 2003

  • Marcinkovich, Carey; Olsen, Kim
  • Journal of Geophysical Research: Solid Earth, Vol. 108, Issue B5
  • DOI: 10.1029/2002JB002235

The implementation of an improved NPML absorbing boundary condition in elastic wave modeling
journal, June 2009


A 3D perfectly matched medium from modified maxwell's equations with stretched coordinates
journal, September 1994

  • Chew, Weng Cho; Weedon, William H.
  • Microwave and Optical Technology Letters, Vol. 7, Issue 13
  • DOI: 10.1002/mop.4650071304

A perfectly matched layer absorbing boundary condition for the second-order seismic wave equation
journal, July 2003


An unsplit convolutional perfectly matched layer improved at grazing incidence for the seismic wave equation
journal, September 2007

  • Komatitsch, Dimitri; Martin, Roland
  • GEOPHYSICS, Vol. 72, Issue 5
  • DOI: 10.1190/1.2757586

Perfectly matched layers for time-harmonic elastodynamics of unbounded domains: theory and finite-element implementation
journal, March 2003


Auxiliary differential equation formulation: an efficient implementation of the perfectly matched layer
journal, February 2003


A symmetric hybrid formulation for transient wave simulations in PML-truncated heterogeneous media
journal, January 2013


Complex frequency shifted convolution PML for FDTD modelling of elastic waves
journal, August 2007


Elastic PMLs for resonator anchor loss simulation
journal, January 2005

  • Bindel, David S.; Govindjee, Sanjay
  • International Journal for Numerical Methods in Engineering, Vol. 64, Issue 6
  • DOI: 10.1002/nme.1394

Convolutional perfectly matched layer for elastic second-order wave equation
journal, March 2010

  • Li, YiFeng; Bou Matar, Olivier
  • The Journal of the Acoustical Society of America, Vol. 127, Issue 3
  • DOI: 10.1121/1.3290999

Runge-Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems
journal, September 2001

  • Cockburn, Bernardo; Shu, Chi-Wang
  • Journal of Scientific Computing, Vol. 16, Issue 3, p. 173-261
  • DOI: 10.1023/A:1012873910884

Perfectly matched layers for transient elastodynamics of unbounded domains
journal, January 2004

  • Basu, Ushnish; Chopra, Anil K.
  • International Journal for Numerical Methods in Engineering, Vol. 59, Issue 8
  • DOI: 10.1002/nme.896

An unsplit convolutional perfectly matched layer technique improved at grazing incidence for the viscoelastic wave equation
journal, October 2009


The Compact Discontinuous Galerkin (CDG) Method for Elliptic Problems
journal, January 2008

  • Peraire, J.; Persson, P. -O.
  • SIAM Journal on Scientific Computing, Vol. 30, Issue 4
  • DOI: 10.1137/070685518

Frequency dependence of the constitutive parameters of causal perfectly matched anisotropic absorbers
journal, January 1996

  • Kuzuoglu, M.; Mittra, R.
  • IEEE Microwave and Guided Wave Letters, Vol. 6, Issue 12
  • DOI: 10.1109/75.544545

Convolution PML (CPML): An efficient FDTD implementation of the CFS-PML for arbitrary media
journal, January 2000


Perfectly Matched Layers for Elastodynamics: a new Absorbing Boundary Condition
journal, December 1996


An MSI Micromechanical Differential Disk-Array Filter
conference, June 2007

  • Li, Sheng-Shian; Lin, Yu-Wei; Ren, Zeying
  • TRANSDUCERS 2007 - 2007 International Solid-State Sensors, Actuators and Microsystems Conference
  • DOI: 10.1109/SENSOR.2007.4300130

Harmonic inversion of time signals and its applications
journal, November 1997

  • Mandelshtam, Vladimir A.; Taylor, Howard S.
  • The Journal of Chemical Physics, Vol. 107, Issue 17
  • DOI: 10.1063/1.475324

An Auxiliary Differential Equation Formulation for the Complex-Frequency Shifted PML
journal, March 2010

  • Gedney, Stephen D.; Zhao, Bo
  • IEEE Transactions on Antennas and Propagation, Vol. 58, Issue 3
  • DOI: 10.1109/TAP.2009.2037765

FDTD for Nth-order dispersive media
journal, January 1992

  • Luebbers, R. J.; Hunsberger, F.
  • IEEE Transactions on Antennas and Propagation, Vol. 40, Issue 11
  • DOI: 10.1109/8.202707