Geodynamic Simulations using the Fast Multipole Boundary Element Method

Simple item page
dc.contributor.advisorHier-Majumder, Saswataen_US
dc.contributor.authorDrombosky, Tyler Williamen_US
dc.contributor.departmentApplied Mathematics and Scientific Computationen_US
dc.contributor.publisherDigital Repository at the University of Marylanden_US
dc.contributor.publisherUniversity of Maryland (College Park, Md.)en_US
dc.date.accessioned2014-06-24T06:05:41Z
dc.date.available2014-06-24T06:05:41Z
dc.date.issued2014en_US
dc.description.abstractInteraction between viscous fluids models two important phenomena in geophysics: (i) the evolution of partially molten rocks, and (ii) the dynamics of Ultralow-Velocity Zones. Previous attempts to numerically model these behaviors have been plagued either by poor resolution at the fluid interfaces or high computational costs. We employ the Fast Multipole Boundary Element Method, which tracks the evolution of the fluid interfaces explicitly and is scalable to large problems, to model these systems. The microstructure of partially molten rocks strongly influences the macroscopic physical properties. The fractional area of intergranular contact, contiguity, is a key parameter that controls the elastic strength of the grain network in the partially molten aggregate. We study the influence of matrix deformation on the contiguity of an aggregate by carrying out pure shear and simple shear deformations of an aggregate. We observe that the differential shortening, the normalized difference between the major and minor axes of grains is inversely related to the ratio between the principal components of the contiguity tensor. From the numerical results, we calculate the seismic anisotropy resulting from melt redistribution during pure and simple shear deformation. During deformation, the melt is expelled from tubules along three grain corners to films along grain edges. The initially isotropic fractional area of intergranular contact, contiguity, becomes anisotropic due to deformation. Consequently, the component of contiguity evaluated on the plane parallel to the axis of maximum compressive stress decreases. We demonstrate that the observed global shear wave anisotropy and shear wave speed reduction of the Lithosphere-Asthenosphere Boundary are best explained by 0.1 vol\% partial melt distributed in horizontal films created by deformation. We use our microsimulation in conjunction with a large scale mantle deep Earth simulation to gain insight into the formation of anisotropy within dense and transient regions known as Ultralow-Velocity Zones, where there is an observed slowdown of shear waves. The results demonstrate a geometric steady state of the dynamic reservoirs by mechanical processes. Within the steady state Ultralow-Velocity Zones, we find significant anisotropy that can explain the speed reduction in shear waves passing through the region.en_US
dc.identifier.urihttp://hdl.handle.net/1903/15302
dc.language.isoenen_US
dc.subject.pqcontrolledApplied mathematicsen_US
dc.subject.pqcontrolledGeophysicsen_US
dc.subject.pquncontrolledBoundary Element Methoden_US
dc.subject.pquncontrolledFast Multipole Methoden_US
dc.subject.pquncontrolledLithosphere-Asthenosphere Boundaryen_US
dc.subject.pquncontrolledMicrostructureen_US
dc.subject.pquncontrolledStokes Flowen_US
dc.subject.pquncontrolledUltralow-Velocity Zoneen_US
dc.titleGeodynamic Simulations using the Fast Multipole Boundary Element Methoden_US
dc.typeDissertationen_US

Files

Original bundle
Now showing 1 - 1 of 1
Thumbnail Image
Name:
Drombosky_umd_0117E_15109.pdf
Size:
14.88 MB
Format:
Adobe Portable Document Format
Download

Collections

UMD Theses and Dissertations
Computer Science Theses and Dissertations
Mathematics Theses and Dissertations