skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


  • NSF Public Access
  • Search Results
  • Higher-order Hamilton–Jacobi perturbation theory for anisotropic heterogeneous media: transformation between Cartesian and ray-centred coordinates
Title: Higher-order Hamilton–Jacobi perturbation theory for anisotropic heterogeneous media: transformation between Cartesian and ray-centred coordinates
SUMMARY Within the field of seismic modelling in anisotropic media, dynamic ray tracing is a powerful technique for computation of amplitude and phase properties of the high-frequency Green’s function. Dynamic ray tracing is based on solving a system of Hamilton–Jacobi perturbation equations, which may be expressed in different 3-D coordinate systems. We consider two particular coordinate systems; a Cartesian coordinate system with a fixed origin and a curvilinear ray-centred coordinate system associated with a reference ray. For each system we form the corresponding 6-D phase spaces, which encapsulate six degrees of freedom in the variation of position and momentum. The formulation of (conventional) dynamic ray tracing in ray-centred coordinates is based on specific knowledge of the first-order transformation between Cartesian and ray-centred phase-space perturbations. Such transformation can also be used for defining initial conditions for dynamic ray tracing in Cartesian coordinates and for obtaining the coefficients involved in two-point traveltime extrapolation. As a step towards extending dynamic ray tracing in ray-centred coordinates to higher orders we establish detailed information about the higher-order properties of the transformation between the Cartesian and ray-centred phase-space perturbations. By numerical examples, we (1) visualize the validity limits of the ray-centred coordinate system, (2) demonstrate the transformation of higher-order derivatives of traveltime from Cartesian to ray-centred coordinates and (3) address the stability of function value and derivatives of volumetric parameters in a higher-order representation of the subsurface model.  more » « less
Award ID(s):
1815143
PAR ID:
10298817
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Geophysical Journal International
Volume:
226
Issue:
2
ISSN:
0956-540X
Page Range / eLocation ID:
893 to 927
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; (, Geophysical Journal International)
    null (Ed.)
    SUMMARY Dynamic ray tracing is a robust and efficient method for computation of amplitude and phase attributes of the high-frequency Green’s function. A formulation of dynamic ray tracing in Cartesian coordinates was recently extended to higher orders. Extrapolation of traveltime and geometrical spreading was demonstrated to yield significantly higher accuracy—for isotropic as well as anisotropic heterogeneous 3-D models of an elastic medium. This is of value in mapping, modelling and imaging, where kernel operations are based on extrapolation or interpolation of Green’s function attributes to densely sampled 3-D grids. We introduce higher-order dynamic ray tracing in ray-centred coordinates, which has certain advantages: (1) such coordinates fit naturally with wave propagation; (2) they lead to a reduction of the number of ordinary differential equations; (3) the initial conditions are simple and intuitive and (4) numerical errors due to redundancies are less likely to influence the computation of the Green’s function attributes. In a 3-D numerical example, we demonstrate that paraxial extrapolation based on higher-order dynamic ray tracing in ray-centred coordinates yields results highly consistent with those obtained using Cartesian coordinates. Furthermore, in a 2-D example we show that interpolation of dynamic ray tracing quantities along a wavefront can be done with much better consistency in ray-centred coordinates than in Cartesian coordinates. In both examples we measure consistency by means of constraints on the dynamic ray tracing quantities in the 3-D position space and in the 6-D phase space. 
    more » « less
  2. ; (, Analysis & PDE)
    Consider a one-dimensional simple small-amplitude solution (ϱ(bkg), v1(bkg)) to the isentropic compressible Euler equations which has smooth initial data, coincides with a constant state outside a compact set, and forms a shock in finite time. Viewing (ϱ(bkg), v1(bkg)) as a plane-symmetric solution to the full compressible Euler equations in three dimensions, we prove that the shock-formation mechanism for the solution (ϱ(bkg), v1(bkg)) is stable against all sufficiently small and compactly supported perturbations. In particular, these perturbations are allowed to break the symmetry and have nontrivial vorticity and variable entropy. Our approach reveals the full structure of the set of blowup-points at the first singular time: within the constant-time hypersurface of first blowup, the solution’s first-order Cartesian coordinate partial derivatives blow up precisely on the zero level set of a function that measures the inverse foliation density of a family of characteristic hypersurfaces. Moreover, relative to a set of geometric coordinates constructed out of an acoustic eikonal function, the fluid solution and the inverse foliation density function remain smooth up to the shock; the blowup of the solution’s Cartesian coordinate partial derivatives is caused by a degeneracy between the geometric and Cartesian coordinates, signified by the vanishing of the inverse foliation density (i.e., the intersection of the characteristics). 
    more » « less
  3. ; ; (, European Journal of Physics)
    Abstract Upper-division undergraduate physics coursework necessitates a firm grasp on and fluid use of mathematical knowledge, including an understanding of non-cartesian (specifically polar, spherical and cylindrical) coordinates and how to use them. A limited body of research into physics students’ thinking about coordinate systems suggests that even for upper-division students, understanding of coordinate system concepts is emergent. To more fully grasp upper-division physics students’ incoming understanding of non-cartesian coordinates, the prevalence of non-cartesian content in seven popular Calculus textbooks was studied. Using content analysis techniques, a coding scheme was developed to gain insight into the presentation of coordinate system content both quantitatively and qualitatively. An initial finding was that non-cartesian basis unit vectors were absent in all but one book. A deeper analysis of three of the calculus textbooks showed that cartesian coordinates comprise an overwhelming proportion of the textbooks’ content and that qualitatively the cartesian coordinate system is presented as the default coordinate system. Quantitative and qualitative results are presented with implications for how these results might impact physics teaching and research at the middle and upper-division. 
    more » « less
  4. ; ; ; ; (, American Mineralogist)
    Abstract The influence of Al substitution on the elastic properties of stishovite and its transition to post-stishovite is of great importance for interpreting the seismic wave velocities of subducted mid-ocean ridge basalt (MORB) within the mantle transition zone and the lower mantle. However, atomistic mechanisms of Al substitution effects on the transition and its associated elasticity remain debated. Here synchrotron single-crystal X-ray diffraction measurements have been performed at room temperature on Al1.3-SiO2 (1.3 mol% Al in the chemical formula of Si0.965(3)Al0.041(1)O2H0.017(4)) and Al2.1-SiO2 (2.1 mol% Al in Si0.948(2)Al0.064(1)O2H 0.018(3)) crystals in diamond anvil cells with Boehler-Almax designed anvils up to 38.0 GPa and 28.5 GPa, respectively. Refinements of the diffraction patterns show that a transformation from stishovite (space group P42/mnm; No. 136) to CaCl2-typed post-stishovite (space group Pnnm; No. 58) is accompanied by splitting of O coordinates. The Al substitution in stishovite results in a faster decrease in the O coordinate, softer apical (Si,Al)-O bonds, and a softer and less distorted (Si,Al)O6 octahedron under compression. This leads to reduced adiabatic bulk modulus (KS), shear modulus (G), shear wave velocity (VS), and compressional wave velocity (VP) in the stishovite phase, explaining seismic wave perturbations in the mantle transition zone. Together with Raman data, Landau theory modeling shows that Al substitution increases the order parameter and excess free energy, stabilizing the post-stishovite phase at lower pressures. Correlation between elasticity and octahedral distortion index (D) reveals that at certain D, the Al substitution reduces KS, G, VS, and VP of the stishovite phase while increasing G, VS, and VP of the post-stishovite phase. Importantly, the maximum shear reduction is slightly enhanced at D = 0.00620(9) at the transition point. Our results help explain the seismically observed small-scale VS anomalies beneath subduction regions in the shallow lower mantle where Al,H-bearing stishovite undergoes the post-stishovite transition. 
    more » « less
  5. ; ; (, Atoms)
    We generalize Schrödinger’s factorization method for Hydrogen from the conventional separation into angular and radial coordinates to a Cartesian-based factorization. Unique to this approach is the fact that the Hamiltonian is represented as a sum over factorizations in terms of coupled operators that depend on the coordinates and momenta in each Cartesian direction. We determine the eigenstates and energies, the wavefunctions in both coordinate and momentum space, and we also illustrate how this technique can be employed to develop the conventional confluent hypergeometric equation approach. The methodology developed here could potentially be employed for other Hamiltonians that can be represented as the sum over coupled Schrödinger factorizations. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email