K. B. Winters, Mariona Claret, M.-Pascale Lelong, Yann Ourmières
{"title":"A Pressure Projection Scheme With Near-Spectral Accuracy for Nonhydrostatic Flow in Domains With Open Boundaries","authors":"K. B. Winters, Mariona Claret, M.-Pascale Lelong, Yann Ourmières","doi":"10.1029/2023MS004040","DOIUrl":null,"url":null,"abstract":"<p>We describe a pressure projection scheme for the simulation of incompressible flow in cubic domains with open boundaries based on fast Fourier transforms. The scheme is implemented in <b>flow_solve</b>, a numerical code designed for process studies of rotating, density-stratified flow. The main algorithmic features of the open-boundary code are the near-spectral accuracy of the discrete differentiation and a dynamic two-dimensional domain decomposition that scales efficiently to large numbers of processors. The simulated flows are not required to be periodic or to satisfy symmetry conditions at the open boundaries owing to the use of mixed series expansions combining cosine and singular Bernoulli polynomial basis functions. These expansions facilitate the imposition of inhomogeneous boundary conditions and allow the code to be used for offline, one-way nesting within an arbitrarily embedded subdomain of a larger scale simulation. The projection scheme is designed to exploit a simple and powerful numerical engine: inversion of Poisson's equation with homogeneous Neumann boundary conditions using fast cosine transforms. Here, we describe the mathematical transformations used to accommodate the imposition of space- and time-varying boundary conditions. The utility of the approach for process studies and for nesting within submesoscale-resolving ocean models is demonstrated with simulations of wind-driven near-inertial waves in the upper ocean.</p>","PeriodicalId":14881,"journal":{"name":"Journal of Advances in Modeling Earth Systems","volume":null,"pages":null},"PeriodicalIF":4.4000,"publicationDate":"2024-06-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1029/2023MS004040","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Advances in Modeling Earth Systems","FirstCategoryId":"89","ListUrlMain":"https://onlinelibrary.wiley.com/doi/10.1029/2023MS004040","RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"METEOROLOGY & ATMOSPHERIC SCIENCES","Score":null,"Total":0}
引用次数: 0
Abstract
We describe a pressure projection scheme for the simulation of incompressible flow in cubic domains with open boundaries based on fast Fourier transforms. The scheme is implemented in flow_solve, a numerical code designed for process studies of rotating, density-stratified flow. The main algorithmic features of the open-boundary code are the near-spectral accuracy of the discrete differentiation and a dynamic two-dimensional domain decomposition that scales efficiently to large numbers of processors. The simulated flows are not required to be periodic or to satisfy symmetry conditions at the open boundaries owing to the use of mixed series expansions combining cosine and singular Bernoulli polynomial basis functions. These expansions facilitate the imposition of inhomogeneous boundary conditions and allow the code to be used for offline, one-way nesting within an arbitrarily embedded subdomain of a larger scale simulation. The projection scheme is designed to exploit a simple and powerful numerical engine: inversion of Poisson's equation with homogeneous Neumann boundary conditions using fast cosine transforms. Here, we describe the mathematical transformations used to accommodate the imposition of space- and time-varying boundary conditions. The utility of the approach for process studies and for nesting within submesoscale-resolving ocean models is demonstrated with simulations of wind-driven near-inertial waves in the upper ocean.
期刊介绍:
The Journal of Advances in Modeling Earth Systems (JAMES) is committed to advancing the science of Earth systems modeling by offering high-quality scientific research through online availability and open access licensing. JAMES invites authors and readers from the international Earth systems modeling community.
Open access. Articles are available free of charge for everyone with Internet access to view and download.
Formal peer review.
Supplemental material, such as code samples, images, and visualizations, is published at no additional charge.
No additional charge for color figures.
Modest page charges to cover production costs.
Articles published in high-quality full text PDF, HTML, and XML.
Internal and external reference linking, DOI registration, and forward linking via CrossRef.