Technologies for supporting high-order geodesic mesh frameworks for computational astrophysics and space sciences
Tóm tắt
Many important problems in astrophysics, space physics, and geophysics involve flows of (possibly ionized) gases in the vicinity of a spherical object, such as a star or planet. The geometry of such a system naturally favors numerical schemes based on a spherical mesh. Despite its orthogonality property, the polar (latitude-longitude) mesh is ill suited for computation because of the singularity on the polar axis, leading to a highly non-uniform distribution of zone sizes. The consequences are (a) loss of accuracy due to large variations in zone aspect ratios, and (b) poor computational efficiency from a severe limitations on the time stepping. Geodesic meshes, based on a central projection using a Platonic solid as a template, solve the anisotropy problem, but increase the complexity of the resulting computer code. We describe a new finite volume implementation of Euler and MHD systems of equations on a triangular geodesic mesh (TGM) that is accurate up to fourth order in space and time and conserves the divergence of magnetic field to machine precision. The paper discusses in detail the generation of a TGM, the domain decomposition techniques, three-dimensional conservative reconstruction, and time stepping.
Tài liệu tham khảo
Balsara, D.S.: Multidimensional HLLE Riemann solver: application to Euler and magnetohydrodynamic flows. J. Comput. Phys. 229, 1970–1993 (2010). https://doi.org/10.1016/j.jcp.2009.11.018
Balsara, D.S.: A two-dimensional HLLC Riemann solver for conservation laws: application to Euler and magnetohydrodynamic flows. J. Comput. Phys. 231, 7476–7503 (2012). https://doi.org/10.1016/j.jcp.2011.12.025
Balsara, D.S.: Multidimensional Riemann problem with self-similar internal structure. Part I—application to hyperbolic conservation laws on structured meshes. J. Comput. Phys. 277, 163–200 (2014). https://doi.org/10.1016/j.jcp.2014.07.053
Balsara, D.S., Dumbser, M.: Divergence-free MHD on unstructured meshes using high order finite volume schemes based on multidimensional Riemann solvers. J. Comput. Phys. 299, 687–715 (2015a). https://doi.org/10.1016/j.jcp.2015.07.012
Balsara, D.S., Dumbser, M.: Multidimensional Riemann problem with self-similar internal structure. Part II—application to hyperbolic conservation laws on unstructured meshes. J. Comput. Phys. 287, 269–292 (2015b). https://doi.org/10.1016/j.jcp.2014.11.004
Balsara, D.S., Florinski, V., Garain, S., Subramanian, S., Gurski, K.F.: Efficient, divergence-free, high order MHD on 3D spherical meshes with optimal geodesic meshing. Mon. Not. R. Astron. Soc. 487, 1283–1314 (2019). https://doi.org/10.1093/mnras/stz1263
Balsara, D.S., Garain, S., Florinski, V., Boscheri, W.: An efficient class of WENO schemes with adaptive order for unstructured meshes. J. Comput. Phys. 404, 109062 (2020). https://doi.org/10.1016/j.jcp.2019.109062
Balsara, D.S., Garain, S., Shu, C.-W.: An efficient class of WENO schemes with adaptive order. J. Comput. Phys. 326, 780–804 (2016). https://doi.org/10.1016/j.jcp.2016.09.009
Balsara, D.S., Rumpf, T., Dumbser, M., Munz, C.-D.: Efficient, high accuracy ADER-WENO schemes for hydrodynamics and divergence-free magnetohydrodynamics. J. Comput. Phys. 228, 2480–2516 (2009). https://doi.org/10.1016/j.jcp.2008.12.003
Balsara, D.S., Shu, C.-W.: Monotonicity peserving weighted essentially non-oscillatory schemes with increasingly high order of accuracy. J. Comput. Phys. 160, 405–452 (2000). https://doi.org/10.1006/jcph.2000.6443
Balsara, D.S., Spicer, D.S.: A staggered mesh algorithm using high order Godunov fluxes to ensure solenoidal magnetic fields in magnetohydrodynamic simulations. J. Comput. Phys. 149, 270–292 (1999)
Batten, P., Clarke, N., Lambert, C., Causon, D.M.: On the choice of wavespeeds for the HLLC Riemann solver. SIAM J. Sci. Comput. 18, 1553–1570 (1997)
Beckmann, J., Mhaskar, H.N., Prestin, J.: Quadrature formulas for integration of multivariate trigonometric polynomials on spherical triangles. GEM Int. J. Geomath. 3, 119–138 (2012). https://doi.org/10.1007/s13137-012-0035-4
Bernard, P.-E., Remacle, J.-F., Comblen, R., Legat, V., Hillewaert, K.: High-order discontinuous Galerkin schemes on general 2D manifolds applied to the shallow water equations. J. Comput. Phys. 228, 6514–6535 (2009). https://doi.org/10.1016/j.jcp.2009.05.046
Browning, G.L., Hack, J.J., Swarztrauber, P.N.: A comparison of three numerical methods for solving differential equations on the sphere. Mon. Weather Rev. 117, 1058–1075 (1989)
Choblet, G., Čadek, O., Couturier, F., Dumoulin, C.: ŒDIPUS: a new tool to study the dynamics of planetary interiors. Geophys. J. Int. 170, 9–30 (2007). https://doi.org/10.1111/j.1365-246X.2007.03419.x
DeVore, C.R.: Flux-corrected transport trchniques for multidimensional compressible magnetohydrodynamics. J. Comput. Phys. 92, 142–160 (1991)
Du, Q., Faber, V., Gunzburger, M.: Centroidal Voronoi tessellations: applications and algorithms. SIAM Rev. 41, 637–676 (1999)
Du, Q., Gunzburger, M.D., Ju, L.: Voronoi-based finite volume methods, optimal Voronoi meshes, and PDEs on the sphere. Comput. Methods Appl. Mech. Eng. 192, 3933–3957 (2003). https://doi.org/10.1016/S0045-7825(03)00394-3
Dumbser, M., Balsara, D.S.: A new efficient formulation of the HLLEM Riemann solver for general conservative and non-conservative hyperbolic systems. J. Comput. Phys. 304, 275–319 (2016). https://doi.org/10.1016/j.jcp.2015.10.014
Dumbser, M., Balsara, D.S., Toro, E.F., Munz, C.-D.: A unified framework for the construction of one-step finite volume and discontinuous Galerkin schemes on unstructured meshes. J. Comput. Phys. 227, 8209–8253 (2008). https://doi.org/10.1016/j.jcp.2008.05.025
Dumbser, M., Käser, M.: Arbitrary high order non-oscillatory finite volume schemes on unstructured meshes for linear hyperbolic systems. J. Comput. Phys. 221, 693–723 (2007). https://doi.org/10.1016/j.jcp.2006.06.043
Dunavant, D.A.: High degree efficient symmetrical Gaussian quadrature rules for the triangle. Int. J. Numer. Methods Eng. 21, 1129–1148 (1985)
Einfeldt, B.: On Godunov-type methods for gas dynamics. SIAM J. Numer. Anal. 25, 294–318 (1988)
Evans, C.R., Hawley, J.F.: Simulation of magnetohydrodynamic flows—a constrained transport method. Astrophys. J. 332, 659–677 (1988). https://doi.org/10.1086/166684
Feng, X., Yang, L., Xiang, C., Wu, S.T., Zhou, Y., Zhong, D.: Three-dimensional solar wind modeling from the Sun to Earth by a SIP-CESE MHD model with a six-component grid. Astrophys. J. 723, 300–319 (2010). https://doi.org/10.1088/0004-637X/723/1/300
Feng, X., Zhou, Y., Wu, S.T.: A novel numerical implementation for solar wind modeling by the modified conservation element/solution wlement method. Astrophys. J. 655, 1110–1126 (2007). https://doi.org/10.1086/510121
Florinski, V., Guo, X., Balsara, D.S., Meyer, C.: Magnetohydrodynamic modeling of solar system processes on geodesic grids. Astrophys. J. Suppl. Ser. 205, 19 (2013). https://doi.org/10.1088/0067-0049/205/2/19
Friedrich, O.: Weighted essentially non-oscillatory schemes for the interpolation of mean values on unstructured grids. J. Comput. Phys. 144, 194–212 (1998)
Giraldo, F.X.: Lagrange–Galerkin methods on spherical geodesic grids. J. Comput. Phys. 136, 197–213 (1997)
Gorski, K.M., Hivon, E., Banday, A.J., Wandelt, B.D., Hansen, F.K., Reinecke, M., Bartelmann, M.: HEALPix: a framework for high-resolution discretization and fast analysis of data distributed on the sphere. Astrophys. J. 622, 759–771 (2005). https://doi.org/10.1086/427976
Guo, X., Florinski, V.: Galactic cosmic-ray intensity modulation by corotating interaction region stream interfaces at 1 AU. Astrophys. J. 826, 65 (2016). https://doi.org/10.3847/0004-637X/826/1/65
Guo, X., Florinski, V., Wang, C.: Effects of anomalous cosmic rays on the structure of the outer heliosphere. Astrophys. J. 859, 157 (2018). https://doi.org/10.3847/1538-4357/aabf42
Gurski, K.F.: An HLLC-type approximate Riemann solver for ideal magnetohydrodynamics. SIAM J. Sci. Comput. 25, 2165–2187 (2004). https://doi.org/10.1137/S1064827502407962
Harten, A., Lax, P.D., van Leer, B.: On upstream differencing and Godunov-type schemes for hyperbolic conservation laws. SIAM Rev. 25, 35–61 (1983)
Harten, A., Osher, S.: Uniformly high-order accurate nonoscillatory schemes. I. SIAM J. Numer. Anal. 2, 279–309 (1987)
Heikes, R., Randall, D.A.: Numerical integration of the shallow-water equations on a twisted icosahedral grid. Part I. Basic design and results of tests. Mon. Weather Rev. 123, 1862–1880 (1995a)
Heikes, R., Randall, D.A.: Numerical integration of the shallow-water equations on a twisted icosahedral grid. Part II. A detailed description of the grid and an analysis of numerical accuracy. Mon. Weather Rev. 123, 1881–1887 (1995b)
Ivan, L., De Sterck, H., Northrup, S.A., Groth, C.P.T.: Multi-dimensional finite-volume scheme for hyperbolic conservation laws on three-dimensional solution-adaptive cubed-sphere grids. J. Comput. Phys. 255, 205–227 (2013). https://doi.org/10.1016/j.jcp.2013.08.008
Ivan, L., De Sterck, H., Susanto, A., Groth, C.P.T.: High-order central ENO finite-volume scheme for hyperbolic conservation laws on three-dimensional cubed-sphere grids. J. Comput. Phys. 282, 157–182 (2015). https://doi.org/10.1016/j.jcp.2014.11.002
Jiang, G.-S., Shu, C.-W.: Efficient implementation of weighted ENO schemes. J. Comput. Phys. 126, 202–228 (1995)
Kageyama, A., Sato, T.: “Yin–Yang grid”: an overset grid in spherical geometry. Geochem. Geophys. Geosyst. 5, Q09005 (2004). https://doi.org/10.1029/2004GC000734
Käser, M., Iske, A.: ADER schemes on adaptive triangular meshes for scalar conservation laws. J. Comput. Phys. 205, 486–508 (2005). https://doi.org/10.1016/j.jcp.2004.11.015
Koldoba, A.V., Romanova, M.M., Ustyugova, G.V., Lovelace, R.V.E.: Three-dimensional magnetohydrodynamic simulations of accretion to an inclined rotator: the “cubed sphere” method. Astrophys. J. 576, L53–L56 (2002)
Li, S.: An HLLC Riemann solver for magneto-hydrodynamics. J. Comput. Phys. 203, 344–357 (2005). https://doi.org/10.1016/j.jcp.2004.08.020
Liu, X.-D., Osher, S., Chan, T.: Weighted essentially non-oscillatory schemes. J. Comput. Phys. 115, 200–212 (1994)
Miura, H.: An upwind-biased conservative advection scheme for spherical hexagonal–pentagonal grids. Mon. Weather Rev. 135, 4038–4044 (2007). https://doi.org/10.1175/2007MWR2101.1
Miyoshi, T., Kusano, K.: A multi-state HLL approximate Riemann solver for ideal magnetohydrodynamics. J. Comput. Phys. 208, 315–344 (2005). https://doi.org/10.1016/j.jcp.2005.02.017
Nakamizo, A., Tanaka, T., Kubo, Y., Kamei, S., Shimazu, H., Shinagawa, H.: Development of the 3-D MHD model of the solar corona-solar wind combining system. J. Geophys. Res. 114, A07109 (2009). https://doi.org/10.1029/2008JA013844
Ollivier-Gooch, C., Van Altena, M.: A high-order-accurate unstructured mesh finite-volume scheme for the advection-diffusion equation. J. Comput. Phys. 181, 729–752 (2002). https://doi.org/10.1006/jcph.2002.7159
Peixoto, P.S., Barros, S.R.M.: Analysis of grid imprinting on geodesic spherical icosahedral grids. J. Comput. Phys. 237, 61–78 (2013). https://doi.org/10.1016/j.jcp.2012.11.041
Phillips, N.A.: Numerical integration of the primitive equations on the hemisphere. Mon. Weather Rev. 87, 333–345 (1959)
Pudykiewicz, J.A.: Numerical solution of the reaction-advection-diffusion equation on the sphere. J. Comput. Phys. 213, 358–390 (2006). https://doi.org/10.1016/j.jcp.2005.08.021
Putman, W.M., Lin, S.-J.: Finite-volume transport on various cubed-sphere grids. J. Comput. Phys. 227, 55–78 (2007). https://doi.org/10.1016/j.jcp.2007.07.022
Ronchi, C., Iacono, R., Paolucci, P.S.: The “cubed sphere”: a new method for the solution of partial differential equations in spherical geometry. J. Comput. Phys. 124, 93–114 (1996)
Ryu, D., Miniati, F., Jones, T.W., Frank, A.: A divergence-free upwind code for multidimensional magnetohydrodynamic flows. Astrophys. J. 509, 244–255 (1998)
Sadourny, R., Arakawa, A., Mintz, Y.: Integration of the nondivergent barotropic vorticity equation with an icosahedral-hexagonal grid for the sphere. Mon. Weather Rev. 96, 351–356 (1968)
Shu, C.-W., Osher, S.: Efficient implementation of essentially non-oscillatory shock-capturing schemes. J. Comput. Phys. 77, 439–471 (1988). https://doi.org/10.1016/0021-9991(88)90177-5
Tomita, H., Tsugawa, M., Satoh, M., Goto, K.: Shallow water model on a modified icosahedral geodesic grid by using spring dynamics. J. Comput. Phys. 174, 579–613 (2001). https://doi.org/10.1006/jcph.2001.6897
Ullrich, P.A., Taylor, M.A.: Arbitrary-order conservative and consistent remapping and a theory of linear maps: part I. Mon. Weather Rev. 143, 2419–2440 (2015). https://doi.org/10.1175/MWR-D-14-00343.1
Usmanov, A.V., Goldstein, M.L., Matthaeus, W.H.: Three-dimensional magnetohydrodynamic modeling of the solar wind including pickup protons and turbulence transport. Astrophys. J. 754, 40 (2012). https://doi.org/10.1088/0004-637X/754/1/40
Weller, H., Thuburn, J., Cotter, C.J.: Computational modes and grid imprinting on five quasi-uniform spherical C grids. Mon. Weather Rev. 140, 2734–2755 (2012). https://doi.org/10.1175/MWR-D-11-00193.1
Zienkiewicz, O.C., Taylor, R.L., Zhu, J.Z.: The Finite Element Method: Its Basis and Fundamentals, vol. 1, 5th edn. Butterworth, Oxford (2013)