skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Tutorial on wave propagation: Spectral Element Method for modeling at Earth scale

Authors:
 [1]
  1. Los Alamos National Laboratory
Publication Date:
Research Org.:
Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1261801
Report Number(s):
LA-UR-16-24741
DOE Contract Number:
AC52-06NA25396
Resource Type:
Conference
Resource Relation:
Conference: SIMNM Summer School ; 2016-07-05 - 2016-07-08 ; Reno, Nevada, United States
Country of Publication:
United States
Language:
English
Subject:
Earth Sciences; wave propagation, seismic waves, Earth scale, SPECFEM3D, spectral element method

Citation Formats

Larmat, Carene. Tutorial on wave propagation: Spectral Element Method for modeling at Earth scale. United States: N. p., 2016. Web.
Larmat, Carene. Tutorial on wave propagation: Spectral Element Method for modeling at Earth scale. United States.
Larmat, Carene. 2016. "Tutorial on wave propagation: Spectral Element Method for modeling at Earth scale". United States. doi:. https://www.osti.gov/servlets/purl/1261801.
@article{osti_1261801,
title = {Tutorial on wave propagation: Spectral Element Method for modeling at Earth scale},
author = {Larmat, Carene},
abstractNote = {},
doi = {},
journal = {},
number = ,
volume = ,
place = {United States},
year = 2016,
month = 7
}

Conference:
Other availability
Please see Document Availability for additional information on obtaining the full-text document. Library patrons may search WorldCat to identify libraries that hold this conference proceeding.

Save / Share:
  • A wide variety of plasma waves play an important role in the energization and loss of particles in the inner magnetosphere. Our ability to understand and model wave-particle interactions in this region requires improved knowledge of the spatial distribution and properties of these waves as well as improved understanding of how the waves depend on changes in solar wind forcing and/or geomagnetic activity. To this end, we have developed a two-dimensional, finite element code that solves the full wave equations in global magnetospheric geometry. The code describes three-dimensional wave structure including mode conversion when ULF, EMIC, and whistler waves aremore » launched in a two-dimensional axisymmetric background plasma with general magnetic field topology. We illustrate the capabilities of the code by examining the role of plasmaspheric plumes on magnetosonic wave propagation; mode conversion at the ion-ion and Alfven resonances resulting from external, solar wind compressions; and wave structure and mode conversion of electromagnetic ion cyclotron waves launched in the equatorial magnetosphere, which propagate along the magnetic field lines toward the ionosphere. We also discuss advantages of the finite element method for resolving resonant structures, and how the model may be adapted to include nonlocal kinetic effects.« less
  • An important challenge for seismic monitoring of nuclear explosions at low magnitude to verify a nuclear-test-ban treaty is the development of techniques that use regional phases for detection, location, and identification. In order to use such phases, region-specific earth models and tools are needed that accurately predict features such as travel times, amplitudes, and spectral characteristics. In this paper, we present our efforts to use two-dimensional finite-difference modeling to help develop and validate regional earth models for the Middle East and North Africa and to develop predictive algorithms for identifying anomalous regional phases. To help develop and validate a modelmore » for the Middle East and North Africa, we compare data and finite-difference simulations for selected regions. We show that the proposed three-dimensional regional model is a significant improvement over standard one-dimensional models by comparing features of broadband data and simulations and differences between observed and predicted features such as narrow-band group velocities. We show how a potential trade-off between source and structure can be avoided by constraining source parameters such as depth, mechanism, and moment/source-time function with independent data. We also present numerous observations of anomalous timing and amplitude of regional phases and show how incorporation of two-dimensional structure can explain many of these observations. Based on these observations, and the predictive capability of our simulations, we develop a simple model that can accurately predict the timing of such phases.« less
  • In this research, the authors present the numerical results of three-dimensional finite element modeling of elastic wave propagation in a half space of a transversely isotropic medium. The input physical parameters are taken from the physical model used by Chang and Gardner (1992). The synthetic seismograms of the three displacement components along the axis at an azimuth of 45{degree} to the fracture orientation on the surface are shown. Snapshots taken on the free surface are illustrated which are the displacement fields with polarization in the direction perpendicular to the fracture orientation. These snapshots show clearly the elliptical patterns demonstrating themore » effects of anisotropy on wave propagation. Four waves, namely, the fast and slow S-waves and the fast and slow P-waves, can be clearly observed.« less
  • A numerical method is described for approximately solving stationary wave propagation models for which the propagation is, for the most part, occurring in one principal direction. This occurs, for example, in connection with the exterior scattering problem at distances sufficiently far away from the obstacle. This also frequently occurs in wave propagation problems in a duct or waveguide (including underwater acoustic wave propagation). The numerical method is based on first factoring out the main propagation direction is the solution and then applying the finite element method. Since the solution being approximated is much smoother, it is seen that there ismore » a great reduction in the number of grid points compared to standard finite element and finite difference methods. This is validated by means of numerical experiments. 132 references.« less