A 3D Spectral Anelastic Hydrodynamic Code for Shearing, Stratified Flows

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.

Submitted to JCP, 36 pages, 9 figures (1 color), For higher quality figures, see http://www.cfa.harvard.edu/~jbarranco

References in corpus (1)

• Three-Dimensional Vortices in Stratified Protoplanetary Disks

Cited by in corpus (17)

• Vortex generation in protoplanetary disks with an embedded giant planet
• Stable a posteriori LES of 2D turbulence using convolutional neural networks: Backscattering analysis and generalization to higher Re via transfer learning
• Self-Replicating Three-Dimensional Vortices in Neutrally-Stable Stratified Rotating Shear Flows
• The Universal Aspect Ratio of Vortices in Rotating Stratified Flows: Theory and Simulation
• Three-Dimensional Simulations of Kelvin-Helmholtz Instability in Settled Dust Layers in Protoplanetary Disks
• Forming Planetesimals by Gravitational Instability: I. The Role of the Richardson Number in Triggering the Kelvin-Helmholtz Instability
• Nonlinear Evolution of Hydrodynamical Shear Flows in Two Dimensions
• Zombie Vortex Instability. III. Persistence with Nonuniform Stratification and Radiative Damping
• Zombie Vortex Instability. II. Thresholds to Trigger Instability and the Properties of Zombie Turbulence in the Dead Zones of Protoplanetary Disks
• Formation, Survival, and Destruction of Vortices in Accretion Disks
• Forming Planetesimals by Gravitational Instability: II. How Dust Settles to its Marginally Stable State
• Stability of 3D Gaussian vortices in an unbounded, rotating, vertically-stratified, Boussinesq flow: Linear analysis
• First Results From New 3D Spectral Simulations Of Anelastic Turbulent Convection
• Reduced-order modeling of fully turbulent buoyancy-driven flows using the Green's function method
• Solitary zonal structures in subcritical drift waves: a minimum model
• Mini-conference on angular momentum transport in laboratory and nature
• Spectral Difference method with a posteriori limiting: II- Application to low Mach number flows