liu.seSearch for publications in DiVA
Change search
ReferencesLink to record
Permanent link

Direct link
Discretizing singular point sources in hyperbolic wave propagation problems
Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, Livermore, USA.
Department of Geophysics, Stanford University, Stanford, USA.
Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, Livermore, USA.
Department of Geophysics, Stanford University, Stanford, USA.
2016 (English)In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 321, 532-555 p.Article in journal (Refereed) Published
Abstract [en]

We develop high order accurate source discretizations for hyperbolic wave propagation problems in first order formulation that are discretized by finite difference schemes. By studying the Fourier series expansions of the source discretization and the finite difference operator, we derive sufficient conditions for achieving design accuracy in the numerical solution. Only half of the conditions in Fourier space can be satisfied through moment conditions on the source discretization, and we develop smoothness conditions for satisfying the remaining accuracy conditions. The resulting source discretization has compact support in physical space, and is spread over as many grid points as the number of moment and smoothness conditions. In numerical experiments we demonstrate high order of accuracy in the numerical solution of the 1-D advection equation (both in the interior and near a boundary), the 3-D elastic wave equation, and the 3-D linearized Euler equations.

Place, publisher, year, edition, pages
2016. Vol. 321, 532-555 p.
Keyword [en]
Singular sources, Hyperbolic wave propagation, Moment conditions, Smoothness conditions, Summation by parts
National Category
Computational Mathematics Mathematical Analysis
Identifiers
URN: urn:nbn:se:liu:diva-132630DOI: 10.1016/j.jcp.2016.05.060OAI: oai:DiVA.org:liu-132630DiVA: diva2:1047470
Available from: 2016-11-17 Created: 2016-11-17 Last updated: 2016-11-17
In thesis
1. Numerical methods for wave propagation in solids containing faults and fluid-filled fractures
Open this publication in new window or tab >>Numerical methods for wave propagation in solids containing faults and fluid-filled fractures
2016 (English)Doctoral thesis, comprehensive summary (Other academic)
Abstract [en]

This thesis develops numerical methods for the simulation of wave propagation in solids containing faults and fluid-filled fractures. These techniques have applications in earthquake hazard analysis, seismic imaging of reservoirs, and volcano seismology. A central component of this work is the coupling of mechanical systems. This aspect involves the coupling of both ordinary differential equations (ODE)(s) and partial differential equations (PDE)(s) along curved interfaces.  All of these problems satisfy a mechanical energy balance. This mechanical energy balance is mimicked by the numerical scheme using high-order accurate difference approximations that satisfy the principle of summation by parts, and by weakly enforcing the coupling conditions. 

The first part of the thesis considers the simulation of dynamic earthquake ruptures along non-planar fault geometries and the simulation of seismic wave radiation from earthquakes, when the earthquakes are idealized as point moment tensor sources. The dynamic earthquake rupture process is simulated by coupling the elastic wave equation at a fault interface to nonlinear ODEs that describe the fault mechanics. The fault geometry is complex and treated by combining structured and unstructured grid techniques. In other applications, when the earthquake source dimension is smaller than wavelengths of interest, the earthquake can be accurately described by a point moment tensor source localized at a single point. The numerical challenge is to discretize the point source with high-order accuracy and without producing spurious oscillations.

The second part of the thesis presents a numerical method for wave propagation in and around fluid-filled fractures. This problem requires the coupling of the elastic wave equation to a fluid inside curved and branching fractures in the solid. The fluid model is a lubrication approximation that incorporates fluid inertia, compressibility, and viscosity. The fracture geometry can have local irregularities such as constrictions and tapered tips. The numerical method discretizes the fracture geometry by using curvilinear multiblock grids and applies implicit-explicit time stepping to isolate and overcome stiffness arising in the semi-discrete equations from viscous diffusion terms, fluid compressibility, and the particular enforcement of the fluid-solid coupling conditions. This numerical method is applied to study the interaction of waves in a fracture-conduit system. A methodology to constrain fracture geometry for oil and gas (hydraulic fracturing) and volcano seismology applications is proposed.

The third part of the thesis extends the summation-by-parts methodology to staggered grids. This extension reduces numerical dispersion and enables the formulation of stable and high-order accurate multiblock discretizations for wave equations in first order form on staggered grids. Finally, the summation-by-parts methodology on staggered grids is further extended to second derivatives and used for the treatment of coordinate singularities in axisymmetric wave propagation.

Place, publisher, year, edition, pages
Linköping: Linköping University Electronic Press, 2016. 27 p.
Series
Linköping Studies in Science and Technology. Dissertations, ISSN 0345-7524 ; 1806
National Category
Computational Mathematics Mathematical Analysis
Identifiers
urn:nbn:se:liu:diva-132550 (URN)10.3384/diss.diva-132550 (DOI)9789176856352 (Print) (ISBN)
Public defence
2016-12-13, Visionen, Hus B, Campus Valla, Linköping, 13:15 (English)
Opponent
Supervisors
Available from: 2016-11-14 Created: 2016-11-14 Last updated: 2016-11-17Bibliographically approved

Open Access in DiVA

No full text

Other links

Publisher's full text
In the same journal
Journal of Computational Physics
Computational MathematicsMathematical Analysis

Search outside of DiVA

GoogleGoogle Scholar

Altmetric score

Total: 8 hits
ReferencesLink to record
Permanent link

Direct link