Skip main navigationJournal menuClose Drawer MenuOpen Drawer Menu
Skip main navigationJournal menuClose Drawer MenuOpen Drawer Menu
Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences
Restricted accessArticles

Space–time adaptive numerical methods for geophysical applications

Published:28 November 2009https://doi.org/10.1098/rsta.2009.0158

    Abstract

    In this paper we present high-order formulations of the finite volume and discontinuous Galerkin finite-element methods for wave propagation problems with a space–time adaptation technique using unstructured meshes in order to reduce computational cost without reducing accuracy. Both methods can be derived in a similar mathematical framework and are identical in their first-order version. In their extension to higher order accuracy in space and time, both methods use spatial polynomials of higher degree inside each element, a high-order solution of the generalized Riemann problem and a high-order time integration method based on the Taylor series expansion. The static adaptation strategy uses locally refined high-resolution meshes in areas with low wave speeds to improve the approximation quality. Furthermore, the time step length is chosen locally adaptive such that the solution is evolved explicitly in time by an optimal time step determined by a local stability criterion. After validating the numerical approach, both schemes are applied to geophysical wave propagation problems such as tsunami waves and seismic waves comparing the new approach with the classical global time-stepping technique. The problem of mesh partitioning for large-scale applications on multi-processor architectures is discussed and a new mesh partition approach is proposed and tested to further reduce computational cost.

    References

    • Amante C.& Eakins B. W.. 2008Etopo1 1 arc-minute global relief model: procedures, data sources and analysis. National Geophysical Data Center, NESDIS, NOAA, US Department of Commerce, Boulder, CO. Google Scholar
    • Babuska I., Henshaw W., Oliger J., Flaherty J., Hopcroft J.& Tezduyar T.. 1995Modeling, mesh generation, and adaptive numerical methods for partial differential equationsIMA Volumes in Mathematics and its Applications 75, 417-430Berlin, GermanySpringer. Crossref, Google Scholar
    • Castro C. E.. 2007High order ADER FV/DG numerical methods for hyperbolic equationsMonographs of the School of Doctoral Studies in Environmental EngineeringTrento, ItalyUniversity of Trento. Google Scholar
    • Castro C. E.& Toro E. F.. 2008Solvers for the high-order Riemann problem for hyperbolic balance laws. J. Comput. Phys. 227, 2481-2513(doi:10.1016/j.jcp.2007.11.013). Crossref, ISI, Google Scholar
    • Castro C. E., Toro E. F.& Käser M.SubmittedADER schemes on unstructured meshes for shallow water: simulation of tsunami waves. Google Scholar
    • Courant R., Friedrichs K. O.& Lewy H.. 1928Über die partiellen Differenzialgleichungen der mathematischen Physik. Mathemat. Annal. 100, 32-74(doi:10.1007/BF01448839). Crossref, Google Scholar
    • Dumbser M.& Käser M.. 2006An arbitrary high order discontinuous Galerkin method for elastic waves on unstructured meshes: II. The three-dimensional isotropic case. Geophys. J. Int. 167, 319-336(doi:10.1111/j.1365-246X.2006.03120.x). Crossref, ISI, Google Scholar
    • Dumbser M.& Käser M.. 2007Arbitrary high order non-oscillatory finite volume schemes on unstructured meshes for linear hyperbolic systems. J. Comput. Phys. 221, 693-723(doi:10.1016/j.jcp.2006.06.043). Crossref, ISI, Google Scholar
    • Dumbser M., Käser M.& Toro E. F.. 2007An arbitrary high order discontinuous Galerkin method for elastic waves on unstructured meshes: V. Local time stepping and p-adaptivity. Geophys. J. Int. 171, 695-717(doi:10.1111/j.1365-246X.2007.03427.x). Crossref, ISI, Google Scholar
    • Dumbser M., Enaux C.& Toro E. F.. 2008Finite volume schemes of very high order of accuracy for stiff hyperbolic balance laws. J. Comput. Phys. 227, 3971-4001(doi:10.1016/j.jcp.2007.12.005). Crossref, ISI, Google Scholar
    • Flaherty J., Loy R., Shephard M., Szymanski B., Teresco J.& Ziantz L.. 1997Adaptive local refinement with octree load-balancing for the parallel solution of three-dimensional conservation laws. J. Parallel Distrib. Comput. 47, 139-152(doi:10.1006/jpdc.1997.1412). Crossref, ISI, Google Scholar
    • Fumeaux C., Baumann D., Leuchtmann P.& Vahldieck R.. 2004A generalized local time-step scheme for efficient FVTD simulations in strongly inhomogeneous meshes. IEEE Trans. Microwave Theory Tech. 52, 1067-1076(doi:10.1109/TMTT.2004.823595). Crossref, ISI, Google Scholar
    • Gassner G., Lörcher F.& Munz C. D.. 2008A discontinuous Galerkin scheme based on a space-time expansion II. Viscous flow equations in multi dimensions. J. Sci. Comput. 34, 260-286(doi:10.1007/s10915-007-9169-1). Crossref, ISI, Google Scholar
    • Karypis G.& Kumar V.. 1998Multilevel k-way partitioning scheme for irregular graphs. J. Parallel Distrib. Comput. 48, 96-129(doi:10.1006/jpdc.1997.1404). Crossref, ISI, Google Scholar
    • Käser M.& Dumbser M.. 2006An arbitrary high order discontinuous Galerkin method for elastic waves on unstructured meshes: I. The two-dimensional isotropic case with external source terms. Geophys. J. Int. 166, 855-877(doi:10.1111/j.1365-246X.2006.03051.x). Crossref, ISI, Google Scholar
    • Käser M.& Iske A.. 2005ADER schemes on adaptive triangular meshes for scalar conservation laws. J. Comput. Phys. 205, 486-508(doi:10.1016/j.jcp.2004.11.015). Crossref, ISI, Google Scholar
    • Lörcher F., Gassner G.& Munz C. D.. 2007A discontinuous Galerkin scheme based on a space-time expansion I. Inviscid compressible flow in one space dimension. J. Sci. Comput. 32, 175-199(doi:10.1007/s10915-007-9128-x). Crossref, ISI, Google Scholar
    • Oeser J., Bunge H.-P.& Mohr M.. 2006Cluster design in the Earth sciences: TETHYS. High performance computing and communications, Gerndt M.& Kranzlmüller D.Springer Lecture Notes in Computer Science 4208, 31–40Berlin, GermanySpringer. Google Scholar
    • Plewa T., Linde T.& Weirs V. G.. 2000Adaptive mesh refinement: theory and applicationsSpringer Lecture notes in Computational Science and Engineering 41Berlin, GermanySpringer. Google Scholar
    • Shi Z.-C., Chen Z., Tang T.& Yu E.. 2006Recent advances in adaptive computation. Contemporary Mathematics 383Oxford, UKOxford University Press. Google Scholar
    • Tessmer E.. 2000Seismic finite-difference modeling with spatially varying time steps. Geophysics 65, 1290-1293(doi:10.1190/1.1444820). Crossref, ISI, Google Scholar
    • Titarev V. A.& Toro E. F.. 2002ADER: arbitrary high order Godunov approach. J. Sci. Comput. 17, 609-618(doi:10.1023/A:1015126814947). Crossref, ISI, Google Scholar
    • Toro E. F.. 2001Shock-capturing methods for free-surface shallow flowsChichester, UKJohn Wiley & Sons Ltd. Google Scholar
    • Toro E. F., Millington R. C.& Nejad L. A. M.. 2001Towards very high order Godunov schemes. Godunov methods: theory and applicationsed. & Toro E. F.907–940New York, NYKluwer Academic/Plenum Publishers. Crossref, Google Scholar
    Previous Article
    Next Article