Show simple item record

dc.contributor.authorGharti, Hom Nath
dc.contributor.authorTromp, Jeroen
dc.contributor.authorZampini, Stefano
dc.date.accessioned2018-11-21T13:12:54Z
dc.date.available2018-11-21T13:12:54Z
dc.date.issued2018-08-07
dc.identifier.citationGharti HN, Tromp J, Zampini S (2018) Spectral-infinite-element simulations of gravity anomalies. Geophysical Journal International 215: 1098–1117. Available: http://dx.doi.org/10.1093/GJI/GGY324.
dc.identifier.issn0956-540X
dc.identifier.issn1365-246X
dc.identifier.doi10.1093/GJI/GGY324
dc.identifier.urihttp://hdl.handle.net/10754/629951
dc.description.abstractGravity anomalies induced by density heterogeneities are governed by Poisson's equation. Most existing methods for modelling such anomalies rely on its integral solution. In this approach, for each observation point, an integral over the entire density distribution needs to be carried out, and the computational cost is proportional to the number of observation points. Frequently, such methods are sensitive to high density contrasts due to inaccurate resolution of the volume integral. We introduce a new approach which directly solves a discretized form of the Poisson/Laplace equation. The main challenge in our approach involves the unbounded nature of the problem, because the potential exists in all of space. To circumvent this challenge, we combine a mapped infinite-element approach with a spectral-element method. Spectral elements represent the domain of interest, and a single layer of infinite elements captures outer space. To solve the weak form of the Poisson/Laplace equation, we use Gauss-Legendre- Lobatto (GLL) quadrature in spectral elements inside the domain of interest. Outside the domain, we use Gauss-Radau quadrature in the infinite direction, and GLL quadrature in the other directions. We illustrate the efficiency and accuracy of our method by comparing calculated gravity anomalies for various density heterogeneities with corresponding analytical solutions. Finally, we consider a complex 3-D model of an ore mine, which consists of both positive and negative density anomalies.
dc.description.sponsorshipWe thank Volker Oye, Katja Sahala, and ISS for access to the mine model. We thank Frederik J. Simons for helpful discussions. Parallel programs were run on computers provided by the Princeton Institute for Computational Science and Engineering (PICSciE). 3D data were visualized using the open-source parallel visualization software ParaView/VTK (www.paraview.org). This research was partially supported by NSF grants 1644826 and 1550901. Our software is open source and freely available via the Computational Infrastructure for Geodynamics (CIG; geodynamics.org). We thank the editor Prof. Bert Vermeersen, Bernhard Steinberger, David Al-Attar, and an anonymous reviewer for their insightful comments which helped to improve the manuscript..
dc.publisherOxford University Press (OUP)
dc.relation.urlhttps://academic.oup.com/gji/article/215/2/1098/5067876
dc.rightsThis is a pre-copyedited, author-produced PDF of an article accepted for publication in Geophysical Journal International following peer review. The version of record is available online at: https://academic.oup.com/gji/article/215/2/1098/5067876.
dc.subjectGravity anomaly
dc.subjectPoisson's equation
dc.subjectSpectral-infiniteelement method
dc.subjectUnbounded domain
dc.titleSpectral-infinite-element simulations of gravity anomalies
dc.typeArticle
dc.contributor.departmentExtreme Computing Research Center
dc.contributor.departmentComputer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division
dc.identifier.journalGeophysical Journal International
dc.eprint.versionPublisher's Version/PDF
dc.contributor.institutionDepartment of Geosciences, Princeton University, Princeton, New Jersey, USA
dc.contributor.institutionProgram in Applied & Computational Mathematics, Princeton University, Princeton, New Jersey, USA
kaust.personZampini, Stefano
refterms.dateFOA2018-11-22T07:29:06Z
dc.date.published-online2018-08-07
dc.date.published-print2018-11-01


Files in this item

Thumbnail
Name:
ggy324.pdf
Size:
7.492Mb
Format:
PDF
Description:
Published version

This item appears in the following Collection(s)

Show simple item record