Gravitational potential is an inseparable part of the earth’s free oscillation simulation. It is also the object of calculation in the study of gravity anomalies. Since the gravitational potential satisfies the unbounded Poisson-Laplace equation, which is zero at infinity, this is frustrating for numerical simulations. The purpose of numerical simulation is usually achieved by limiting the solution domain and approximating its boundary conditions. We attempt to use the infinite-spectral hybrid method to fit the infinity boundary directly and no longer approximate the boundary conditions. Considering the computational efficiency of the 3D Earth model, this study uses a 2.5-dimensional governing equation. Finally, numerical experiments verify the factual accuracy of this method.
REN Junsheng
,
ZHANG Huai
. Application of infinite-spectral hybrid method in 2.5-dimensional gravitational potential calculation[J]. Journal of University of Chinese Academy of Sciences, 2024
, 41(1)
: 65
-69
.
DOI: 10.7523/j.ucas.2022.041
[1] Woodhouse J H, Deuss A. Theory and observations-earth’s free oscillations [M]//Treatise on Geophysics. Amsterdam: Elsevier, 2015: 79-115 DOI:10.1016/b978-0-444-53802-4.00002-6.
[2] Cowling T G. The non-radial oscillations of polytropic stars[J]. Monthly Notices of the Royal Astronomical Society, 1941, 101(8): 367-375. DOI:10.1093/mnras/101.8.367.
[3] Yan R, Woith H, Wang R J, et al. Earth’s free oscillations excited by the 2011 Tohoku Mw 9.0 earthquake detected with a groundwater level array in mainland China[J]. Geophysical Journal International, 2016, 206(3): 1457-1466. DOI:10.1093/gji/ggw213.
[4] Dahlen F A, Tromp J. Theoretical global seismology[M]. Princeton: Princeton University Press, 2021. DOI:10.1515/9780691216157.
[5] Cai Y E, Wang C Y. Fast finite-element calculation of gravity anomaly in complex geological regions[J]. Geophysical Journal International, 2005, 162(3): 696-708. DOI:10.1111/j.1365-246X.2005.02711.x.
[6] Chaljub E, Valette B. Spectral element modelling of three-dimensional wave propagation in a self-gravitating Earth with an arbitrarily stratified outer core[J]. Geophysical Journal International, 2004, 158(1): 131-141. DOI:10.1111/j.1365-246X.2004.02267.x.
[7] Xu C Y, Huo Z J. Numerical simulation of gravity anomaly based on the unstructured element grid and finite element method[J]. Mathematical Problems in Engineering, 2020, 2020: 3604084. DOI:10.1155/2020/3604084.
[8] Hejlesen M M, Rasmussen J T, Chatelain P, et al. A high order solver for the unbounded Poisson equation[J]. Journal of Computational Physics, 2013, 252: 458-467. DOI:10.1016/j.jcp.2013.05.050.
[9] O.C Z, R.L T, Zhu J Z. The finite element method: Its basis and fundamentals[M]. Sixth. Burlington (Mass.): Elsevier, 2005.
[10] Chaljub E, Komatitsch D, Vilotte J P, et al. Spectral-element analysis in seismology[J]. Advances in Geophysics. Elsevier, 2007, 48: 365-419. DOI:10.1016/S0065-2687(06)48007-9.
[11] Ritsema J, Deuss A, van Heijst H J, et al. S40RTS: a degree-40 shear-velocity model for the mantle from new Rayleigh wave dispersion, teleseismic traveltime and normal-mode splitting function measurements[J]. Geophysical Journal International, 2011, 184(3): 1223-1236. DOI:10.1111/j.1365-246X.2010.04884.x.
[12] Nissen-Meyer T, Fournier A, Dahlen F A. A two-dimensional spectral-element method for computing spherical-earth seismograms-I. Moment-tensor source[J]. Geophysical Journal International, 2007, 168(3): 1067-1092. DOI:10.1111/j.1365-246X.2006.03121.x.
[13] Nissen-Meyer T, Fournier A, Dahlen F A. A 2-D spectral-element method for computing spherical-earth seismograms-II. Waves in solid-fluid media[J]. Geophysical Journal International, 2008, 174(3): 873-888. DOI:10.1111/j.1365-246X.2008.03813.x.
[14] Gerritsma M I, Phillips T N. Spectral element methods for axisymmetric stokes problems[J]. Journal of Computational Physics, 2000, 164(1): 81-103. DOI:10.1006/jcph.2000.6574.
[15] Shen J, Tang T. Spectral and high-order methods with applications[M]. Beijing: Science Press, 2006.
[16] 李瑞浩. 重力学引论[M]. 北京: 地震出版社, 1988.
