A 3D spectral anelastic hydrodynamic code for shearing, stratified flows

DOI10.1016/J.JCP.2006.03.015MaRDI QIDQ856365zbMATHOpenAlexFDO

Authors J. Barranco, Philip S. Marcus

Publication date 7 December 2006

Published in Journal of Computational Physics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/astro-ph/0509063

zbMATH Keywords

spectral methodshydrodynamicsvortex dynamicscoordinate mappinginfinite domainrotating flowsshear flowsstratified flowsanelastic approximation

Mathematics Subject Classification ID

Vortex flows for incompressible inviscid fluids (76B47) Spectral methods applied to problems in fluid mechanics (76M22) General theory of rotating fluids (76U05) Internal waves for incompressible inviscid fluids (76B55) Galactic and stellar structure (85A15)

Abstract: We have developed a three-dimensional (3D) spectral hydrodynamic code to study vortex dynamics in rotating, shearing, stratified systems (e.g. the atmosphere of gas giant planets, protoplanetary disks around newly forming protostars). The time-independent background state is stably stratified in the vertical direction and has a unidirectional linear shear flow aligned with one horizontal axis. Superposed on this background state is an unsteady, subsonic flow that is evolved with the Euler equations subject to the anelastic approximation to filter acoustic phenomena. A Fourier-Fourier basis in a set of quasi-Lagrangian coordinates that advect with the background shear is used for spectral expansions in the two horizontal directions. For the vertical direction, two different sets of basis functions have been implemented: (1) Chebyshev polynomials on a truncated, finite domain, and (2) rational Chebyshev functions on an infinite domain. Use of this latter set is equivalent to transforming the infinite domain to a finite one with a cotangent mapping, and using cosine and sine expansions in the mapped coordinate. The nonlinear advection terms are time integrated explicitly, whereas the Coriolis force, buoyancy terms, and pressure/enthalpy gradient are integrated semi-implicitly. We show that internal gravity waves can be damped by adding new terms to the Euler equations. The code exhibits excellent parallel performance with the Message Passing Interface (MPI). As a demonstration of the code, we simulate the merger of two 3D vortices in the midplane of a protoplanetary disk.

Recommendations

• Computational vortex dynamics in two and three dimensions
• A compressible high-order unstructured spectral difference code for stratified convection in rotating spherical shells
• Spectral simulation for the development of vortex arrays in stratified fluids
• scientific article; zbMATH DE number 2232697
• A numerical algorithm for geophysical and astrophysical inhomogeneous fluid flows

Cited in

(7)

• Stable \textit{a posteriori} LES of 2D turbulence using convolutional neural networks: backscattering analysis and generalization to higher \(Re\) via transfer learning
• A spectral anelastic Navier-Stokes solver for a stratified two-miscible-layer system in infinite horizontal channel
• Fast, accurate, and small-scale direct trajectory optimization using a Gegenbauer transcription method
• A compressible high-order unstructured spectral difference code for stratified convection in rotating spherical shells
• The universal aspect ratio of vortices in rotating stratified flows: theory and simulation
• Asymptotic analysis of Rayleigh-Taylor flow for Newtonian miscible fluids
• Stability of three-dimensional Gaussian vortices in an unbounded, rotating, vertically stratified, Boussinesq flow: linear analysis

Describes a project that uses

Uses Software

• ZEUS

This page was built for publication: A 3D spectral anelastic hydrodynamic code for shearing, stratified flows

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q856365)