Show simple item record

dc.contributor.advisor de Hoop, Maarten
dc.creatorYe, Jingchen
dc.date.accessioned 2019-01-02T22:11:47Z
dc.date.created 2017-12
dc.date.issued 2018-02-21
dc.date.submitted December 2017
dc.identifier.urihttps://hdl.handle.net/1911/104942
dc.description.abstract We develop a robust algorithm to compute seismic normal modes in a spherically symmetric, non-rotating Earth. A well-known problem is the cross-contamination of modes near "intersections" of dispersion curves for separate waveguides. Our novel computational approach completely avoids artificial degeneracies by guaranteeing orthonormality among the eigenfunctions. We extend Buland’s work, and reformulate the Sturm-Liouville problem as a generalized eigenvalue problem with the Rayleigh-Ritz Galerkin method. A special projection operator incorporating the gravity terms proposed by de Hoop and a displacement/pressure formulation are utilized in the fluid outer core to project out the essential spectrum. Moreover, the weak variational form enables us to achieve high accuracy across the solid-fluid boundary, especially for Stoneley modes, which have exponentially decaying behavior. We also employ the mixed finite element technique to avoid spurious pressure modes arising from discretization schemes and a numerical inf-sup test is performed following Bathe’s work. In addition, the self-gravitation terms are reformulated to avoid computations outside the Earth, thanks to the domain decomposition technique. Our package enables us to study the physical properties of intersection points of waveguides. According to Okal's classification theory, the group velocities should be continuous within a branch of the same mode family. However, we have found that there will be a small “bump” near intersection points, which is consistent with Miropol'sky’s observation. In fact, we can loosely regard Earth’s surface and the CMB as independent waveguides. For those modes that are far from the intersection points, their eigenfunctions are localized in the corresponding waveguides. However, those that are close to intersections will have physical features of both waveguides, which means they cannot be classified in either family. Our results improve on Okal’s classification, demonstrating that dispersion curves from independent waveguides should be considered to break at intersection points. Moreover, intersection points have close relations with Stoneley-related modes observed from earthquakes data, which casts light on studying deep Earth's structures.
dc.format.mimetype application/pdf
dc.language.iso en
dc.subjectNormal modes computation
Cross-contamination
Intersection points
dc.title REVISITING THE COMPUTATION OF NORMAL MODES IN SNREI MODELS OF PLANETS - close eigenfrequencies
dc.type Thesis
dc.date.updated 2019-01-02T22:11:48Z
dc.type.material Text
thesis.degree.department Applied Physics
thesis.degree.discipline Appl Physics/Comp & Appl Math
thesis.degree.grantor Rice University
thesis.degree.level Masters
thesis.degree.name Master of Science
dc.embargo.terms 2019-12-01
dc.embargo.lift 2019-12-01
dc.identifier.citation Ye, Jingchen. "REVISITING THE COMPUTATION OF NORMAL MODES IN SNREI MODELS OF PLANETS - close eigenfrequencies." (2018) Master’s Thesis, Rice University. https://hdl.handle.net/1911/104942.


Files in this item

Thumbnail

This item appears in the following Collection(s)

Show simple item record