{"title":"Compressible and anelastic governing-equation solution methods for thermospheric gravity waves with realistic background parameters","authors":"Harold Knight, Dave Broutman, Stephen Eckermann","doi":"10.1007/s00162-024-00709-x","DOIUrl":null,"url":null,"abstract":"<div><h3>abstract</h3><p>A previously developed numerical-multilayer modeling approach for systems of governing equations is extended so that unwanted terms, resulting from vertical variations in certain background parameters, can be removed from the dispersion-relation polynomial associated with the system. The new approach is applied to linearized anelastic and compressible systems of governing equations for gravity waves including molecular viscosity and thermal diffusion. The ability to remove unwanted terms from the dispersion-relation polynomial is crucial for solving the governing equations when realistic background parameters, such as horizontal velocity and temperature, with strong vertical gradients, are included. With the unwanted terms removed, previously studied dispersion-relation polynomials, for which methods for defining upgoing and downgoing vertical wavenumber roots already exist, are obtained. The new methods are applied to a comprehensive set of medium-scale time-wavepacket examples, with realistic background parameters, lower boundary conditions at 30 km altitude, and modeled wavefields extending up to 500 km altitude. Results from the compressible and anelastic model versions are compared, with compressible governing-equation solutions understood as the more physically accurate of the two. The new methods provide significantly less computationally expensive alternatives to nonlinear time-step methods, which makes them useful for comprehensive studies of the behavior of viscous/diffusive gravity waves and also for large studies of cases based on observational data. Additionally, they generalize previously existing Fourier methods that have been applied to inviscid problems while providing a theoretical framework for the study of viscous/diffusive gravity waves.</p><h3>Graphic abstract</h3><div><figure><div><div><picture><source><img></source></picture></div></div></figure></div></div>","PeriodicalId":795,"journal":{"name":"Theoretical and Computational Fluid Dynamics","volume":"38 4","pages":"479 - 509"},"PeriodicalIF":2.2000,"publicationDate":"2024-07-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s00162-024-00709-x.pdf","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Theoretical and Computational Fluid Dynamics","FirstCategoryId":"5","ListUrlMain":"https://link.springer.com/article/10.1007/s00162-024-00709-x","RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"MECHANICS","Score":null,"Total":0}
引用次数: 0
Abstract
abstract
A previously developed numerical-multilayer modeling approach for systems of governing equations is extended so that unwanted terms, resulting from vertical variations in certain background parameters, can be removed from the dispersion-relation polynomial associated with the system. The new approach is applied to linearized anelastic and compressible systems of governing equations for gravity waves including molecular viscosity and thermal diffusion. The ability to remove unwanted terms from the dispersion-relation polynomial is crucial for solving the governing equations when realistic background parameters, such as horizontal velocity and temperature, with strong vertical gradients, are included. With the unwanted terms removed, previously studied dispersion-relation polynomials, for which methods for defining upgoing and downgoing vertical wavenumber roots already exist, are obtained. The new methods are applied to a comprehensive set of medium-scale time-wavepacket examples, with realistic background parameters, lower boundary conditions at 30 km altitude, and modeled wavefields extending up to 500 km altitude. Results from the compressible and anelastic model versions are compared, with compressible governing-equation solutions understood as the more physically accurate of the two. The new methods provide significantly less computationally expensive alternatives to nonlinear time-step methods, which makes them useful for comprehensive studies of the behavior of viscous/diffusive gravity waves and also for large studies of cases based on observational data. Additionally, they generalize previously existing Fourier methods that have been applied to inviscid problems while providing a theoretical framework for the study of viscous/diffusive gravity waves.
期刊介绍:
Theoretical and Computational Fluid Dynamics provides a forum for the cross fertilization of ideas, tools and techniques across all disciplines in which fluid flow plays a role. The focus is on aspects of fluid dynamics where theory and computation are used to provide insights and data upon which solid physical understanding is revealed. We seek research papers, invited review articles, brief communications, letters and comments addressing flow phenomena of relevance to aeronautical, geophysical, environmental, material, mechanical and life sciences. Papers of a purely algorithmic, experimental or engineering application nature, and papers without significant new physical insights, are outside the scope of this journal. For computational work, authors are responsible for ensuring that any artifacts of discretization and/or implementation are sufficiently controlled such that the numerical results unambiguously support the conclusions drawn. Where appropriate, and to the extent possible, such papers should either include or reference supporting documentation in the form of verification and validation studies.