Efficient Spectral Element Method for the Euler Equations on Unbounded Domains (2401.05624v3)
Abstract: Mitigating the impact of waves leaving a numerical domain has been a persistent challenge in numerical modeling. Reducing wave reflection at the domain boundary is crucial for accurate simulations. Absorbing layers, while common, often incur significant computational costs. This paper introduces an efficient application of a Legendre-Laguerre basis for absorbing layers for two-dimensional non-linear compressible Euler equations. The method couples a spectral-element bounded domain with a semi-infinite region, employing a tensor product of Lagrange and scaled Laguerre basis functions. Semi-infinite elements are used in the absorbing layer with Rayleigh damping. In comparison to existing methods with similar absorbing layer extensions, this approach, a pioneering application to the Euler equations of compressible and stratified flows, demonstrates substantial computational savings. The study marks the first application of semi-infinite elements to mitigate wave reflection in the solution of the Euler equations, particularly in nonhydrostatic atmospheric modeling. A comprehensive set of tests demonstrates the method's versatility for general systems of conservation laws, with a focus on its effectiveness in damping vertically propagating mountain gravity waves, a benchmark for atmospheric models. Across all tests, the model presented in this paper consistently exhibits notable performance improvements compared to a traditional Rayleigh damping approach.
- Well-posed perfectly matched layers for advective acoustics. Journal of Computational Physics, 154(2):266–283, 1999.
- T. Benacchio and L. Bonaventura. Absorbing boundary conditions: a spectral collocation approach. Int. J. Numer. Methods Fluids, 72(9):913–936, 2013.
- T. Benacchio and L. Bonaventura. An extension of dg methods for hyperbolic problems to one-dimensional semi-infinite domains. Appl. Math. Comput., 350:266–282, 2019.
- J.-P. Berenger. A perfectly matched layer for the absorption of electromagnetic waves. J. Comput. Phys., 114(2):185–200, 1994.
- J.P. Boyd. The erfc-log filter and the asymptotics of the Euler and Vandeven sequence accelerations. In Proceedings of the Third Internatioal Conference on Spectral and High Order Methods, pages 267–276. Houston Math. J, 1996.
- Adaptive hermite spectral methods in unbounded domains. Applied Numerical Mathematics, 183:201–220, 2023.
- J. R. Dea. An experimental adaptation of higdon-type non-reflecting boundary conditions to linear first-order systems. Journal of computational and applied mathematics, 235(5):1354–1366, 2011.
- High-order non-reflecting boundary conditions for the linearized 2-D Euler equations: No mean flow case. Wave Motion, 46:210–220, 2009.
- D. Durran and J. Klemp. A compressible model for the simulation of moist mountain waves. Mon. Wea. Rev., 111(12):2341–2361, 1983.
- B. Engquist and A. Majda. Absorbing boundary conditions for the numerical simulation of waves. Math Comput, 31:629–651, 1977.
- B. Engquist and A. Majda. Radiation boundary conditions for acoustic and elastic wave calculations. Comm. Pure Appl. Math., 32:313–357, 1979.
- M. B. Giles. Nonreflecting boundary conditions for Euler equation calculations. AIAA Journal, 28:2050, 1990.
- F. X. Giraldo. An Introduction to Element-based Galerkin Methods on Tensor-Product Bases - Analysis, Algorithms, and Applications, volume 24. Springer, 2020.
- A scalable spectral element Eulerian atmospheric model (SEE-AM) for NWP: Dynamical core tests. Mon. Weath. Rev., 132(1):133–153, 2004.
- Calculation of Gauss quadrature rules. Mathematics of computation, 23(106):221–230, 1969.
- J. K. Hargreaves. The solar-terrestrial environment: An introduction to geospace-the science of the terrestrial upper atmosphere, ionosphere, and magnetosphere. Cambridge university press, 1992.
- J. S. Hesthaven. On the analysis and construction of perfectly matched layers for the linearized euler equations. J. Comput. Phys., 142(1):129–147, 1998.
- R. L. Higdon. Absorbing boundary conditions for difference approximations to the multidimensional wave equation. Math. Comput., 47(176):437–459, 1986.
- R. L. Higdon. Radiation boundary conditions for dispersive waves. SIAM J. Numer. Anal., 31(1):64–100, 1994.
- M. Israeli and S. A. Orszag. Approximation of radiation boundary conditions. Journal of computational physics, 41(1):115–135, 1981.
- G. Karniadakis and S. Sherwin. Spectral/hp element methods for computational fluid dynamics. Oxford University Press, 2nd Ed., 2005.
- Continuous and discontinuous Galerkin methods for a scalable three-dimensional nonhydrostatic atmospheric model: limited-area mode. J. Comput. Phys., 231:7988–8008, 2012.
- A physics-based open atmosphere boundary condition for height-coordinate atmospheric models. J. Comput. Phys., 482:112044, 2023.
- An upper gravity-wave absorbing layer for nwp applications. Mon. Wea. Rev., 136:3987–4004, 2008.
- J. B. Klemp and D. Durran. An upper boundary condition permitting integral gravity wave radiation in numerical mesoscale models. Mon. Wea. Rev., 111:430–444, 1983.
- Numerical simulation of hydrostatic mountain waves. Journal of Atmospheric Sciences, 35(1):78–107, 1978.
- A constant pressure upper boundary formulation for models employing height-based vertical coordinates. Mon. Weath. Rev., 150(8):2175–2186, 2022.
- A pretty good sponge: Dealing with open boundaries in limited-area ocean models. Ocean Modelling, 20(3):270–292, 2008.
- A Variational Multiscale Stabilized finite element method for the solution of the Euler equations of nonhydrostatic stratified flows. Arch. Comput. Methods Eng., 23:673–722, 2015.
- S. Marras and Y. Tissaoui. Jexpresso: https://github.com/smarras79/jexpresso, 2024.
- On the parameters of absorbing layers for shallow water models. Ocean Dynamics, 60:65–79, 2010.
- A perfectly matched layer approach to the linearized shallow water equations models. Monthly Weather Review, 132(6):1369–1378, 2004.
- Application of high-order higdon non-reflecting boundary conditions to linear shallow water models. Communications in Numerical Methods in Engineering, 24(11):1459–1466, 2008.
- C. Rackauckas and Q. Nie. DifferentialEquations.jl–a performant and feature-rich ecosystem for solving differential equations in Julia. JORS, 5(1), 2017.
- Development and testing of a next generation spectral element model for the US navy. 17th Workshop on High Performance Computing (ECMWF), 24-28 October, Reading, UK, 2016.
- J. Shen. Stable and efficient spectral methods in unbounded domains using laguerre functions. SIAM Journal on Numerical Analysis, 38(4):1113–1133, 2000.
- J. Shen and L.-L. Wang. Some recent advances on spectral methods for unbounded domains. Communications in computational physics, 5(2-4):195–241, 2009.
- R. B. Smith. Linear theory of stratified hydrostatic flow past an isolated mountain. Tellus, 32(4):348–364, 1980.
- Large-eddy simulations with ClimateMachine v0.2.0: a new open-source code for atmospheric simulations on GPUs and CPUs. Geoscient. Model Develop., gmd-2021-335, 2022.
- M. A. Taylor and A. Fournier. A compatible and conservative spectral element method on unstructured grids. J. Comput. Phys., 229:5879–5895, 2010.
- J. Thuburn. Some conservation issues for the dynamical cores of NWP and climate models. J. Comput. Phys., 227(7):3715–3730, 2008.
- A non-column based, fully unstructured implementation of Kessler microphysics with warm rain using continuous and discontinuous spectral elements. J. Adv. Modeling Earth Sys., 2022MS003283, 2023.
- F. Vismara and T. Benacchio. A seamless, extended DG approach for advection–diffusion problems on unbounded domains. J. Scient. Computing, 90(1):1–27, 2023.
- F. Vismara and T. Benacchio. Efficient hyperbolic–parabolic models on multi-dimensional unbounded domains using an extended DG approach. Int. J. Numer. Methods Fluids, 96(2):125–230, 2023.
- Spectrally adapted physics-informed neural networks for solving unbounded domain problems. Machine Learning: Science and Technology, 4, 2023.
- Efficient scaling and moving techniques for spectral methods in unbounded domains. SIAM Journal for Scientific Computing, 43(5):A3244–A3268, 2021.
- A coupled Legendre-Laguerre spectral-element method for the Navier-Stokes equations in unbounded domains. Journal of Scientific Computing, 42(1):1, 2010.