TY - JOUR
T1 - Efficient spectral element method for the Euler equations on unbounded domains
AU - Tissaoui, Yassine
AU - Kelly, James F.
AU - Marras, Simone
N1 - Publisher Copyright:
© 2024 Elsevier Inc.
PY - 2025/2/15
Y1 - 2025/2/15
N2 - 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.
AB - 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.
UR - http://www.scopus.com/inward/record.url?scp=85204789208&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85204789208&partnerID=8YFLogxK
U2 - 10.1016/j.amc.2024.129080
DO - 10.1016/j.amc.2024.129080
M3 - Article
AN - SCOPUS:85204789208
SN - 0096-3003
VL - 487
JO - Applied Mathematics and Computation
JF - Applied Mathematics and Computation
M1 - 129080
ER -