A semi-analytical pseudo-spectral method for 3D Boussinesq equations of rotating, stratified flows in unbounded cylindrical domains
Fuente:
arXiv
Guardado en:
| Autores principales: | , |
|---|---|
| Formato: | Preprint |
| Publicado: |
2026
|
| Materias: | |
| Acceso en línea: | |
| Etiquetas: |
Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
|
| _version_ | 1866910044955607040 |
|---|---|
| author | Wang, Jinge Marcus, Philip S. |
| author_facet | Wang, Jinge Marcus, Philip S. |
| contents | We present a pseudo-spectral method for solving the three-dimensional Boussinesq equations in unbounded cylindrical domains, specifically tailored for rotating, stably stratified flows subject to strong azimuthal shear. To effectively capture the global geometry without sacrificing spectral accuracy, the spatial discretization employs Fourier expansions in the azimuthal and axial directions alongside mapped associated Legendre polynomials in the radial direction. This basis spans the semi-infinite domain while analytically resolving the coordinate singularity at the origin. While this spectral framework ensures high spatial fidelity, the temporal integration of these rotating shear flows presents a formidable computational challenge due to the numerical stiffness driven by fast restorative wave forces and rapid background advection. To circumvent this, we develop an exponential time differencing (ETD) scheme that analytically integrates the fully coupled linear operator, including the radially dependent advective cross terms. By encoding the physical resonance characteristics and stability limits of the background flow directly into the integration operators, the proposed ETD formulation removes the numerical stability constraints imposed by the background shear and stratification. This permits integration time steps scaled by the slow macroscopic evolution of the physical instabilities rather than the fast background kinematics, offering significant performance gains over standard mixed implicit-explicit schemes. The method's accuracy and stability are validated through the precise conservation of energy and angular momentum, establishing a robust framework for simulating instabilities in astrophysical and geophysical vortices. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2603_07046 |
| institution | arXiv |
| publishDate | 2026 |
| record_format | arxiv |
| spellingShingle | A semi-analytical pseudo-spectral method for 3D Boussinesq equations of rotating, stratified flows in unbounded cylindrical domains Wang, Jinge Marcus, Philip S. Fluid Dynamics Computational Physics Geophysics We present a pseudo-spectral method for solving the three-dimensional Boussinesq equations in unbounded cylindrical domains, specifically tailored for rotating, stably stratified flows subject to strong azimuthal shear. To effectively capture the global geometry without sacrificing spectral accuracy, the spatial discretization employs Fourier expansions in the azimuthal and axial directions alongside mapped associated Legendre polynomials in the radial direction. This basis spans the semi-infinite domain while analytically resolving the coordinate singularity at the origin. While this spectral framework ensures high spatial fidelity, the temporal integration of these rotating shear flows presents a formidable computational challenge due to the numerical stiffness driven by fast restorative wave forces and rapid background advection. To circumvent this, we develop an exponential time differencing (ETD) scheme that analytically integrates the fully coupled linear operator, including the radially dependent advective cross terms. By encoding the physical resonance characteristics and stability limits of the background flow directly into the integration operators, the proposed ETD formulation removes the numerical stability constraints imposed by the background shear and stratification. This permits integration time steps scaled by the slow macroscopic evolution of the physical instabilities rather than the fast background kinematics, offering significant performance gains over standard mixed implicit-explicit schemes. The method's accuracy and stability are validated through the precise conservation of energy and angular momentum, establishing a robust framework for simulating instabilities in astrophysical and geophysical vortices. |
| title | A semi-analytical pseudo-spectral method for 3D Boussinesq equations of rotating, stratified flows in unbounded cylindrical domains |
| topic | Fluid Dynamics Computational Physics Geophysics |
| url | https://arxiv.org/abs/2603.07046 |