Length and Velocity Scales in Protoplanetary Disk Turbulence

Debanjan Sengupta,Jeffrey N. Cuzzi,Orkan M. Umurhan,Wladimir Lyra

2024-02-24

Abstract:In the theory of protoplanetary disk turbulence, a widely adopted \emph{ansatz}, or assumption, is that the turnover frequency of the largest turbulent eddy, $\Omega_L$, is the local Keplerian frequency $\Omega_K$. In terms of the standard dimensionless Shakura-Sunyaev $\alpha$ parameter that quantifies turbulent viscosity or diffusivity, this assumption leads to characteristic length and velocity scales given respectively by $\sqrt{\alpha}H$ and $\sqrt{\alpha}c$, in which $H$ and $c$ are the local gas scale height and sound speed. However, this assumption is not applicable in cases when turbulence is forced numerically or driven by some natural processes such as Vertical Shear Instability. Here we explore the more general case where $\Omega_L\ge\Omega_K$ and show that under these conditions, the characteristic length and velocity scales are respectively $\sqrt{\alpha/R'}H$ and $\sqrt{\alpha R'}c$, where $R'\equiv \Omega_L/\Omega_K$ is twice the Rossby number. It follows that $\alpha=\alphat/R'$, where $\sqrt{\alphat} c$ is the root-mean-square average of the turbulent velocities. Properly allowing for this effect naturally explains the reduced particle scale heights produced in shearing box simulations of particles in forced turbulence, and may help with interpreting recent edge-on disk observations; more general implications for observations are also presented. For $R'>1$ the effective particle Stokes numbers are increased, which has implications for particle collision dynamics and growth, as well as for planetesimal formation.

Earth and Planetary Astrophysics

What problem does this paper attempt to address?

The paper primarily discusses the issues of length and velocity scales in the theory of protoplanetary disk turbulence. In the theory of protoplanetary disk turbulence, a widely accepted assumption is that the turnover frequency of the largest turbulent eddies, Ω_L, is equal to the local Keplerian frequency, Ω_K. This leads to the characteristic length and velocity scales of turbulence viscosity or diffusivity being quantified by the standard dimensionless Shakura-Sunyaev α parameter, represented as √αH and √αc, respectively, where H is the local gas scale height, and c is the sound speed.

However, when turbulence is driven by numerical forcing or natural processes such as vertical shear instability, this assumption does not hold. This paper explores the more general case where Ω_L≥Ω_K, and demonstrates that under this condition, the characteristic length and velocity scales become √α/R' * H and √αR' * c, respectively, where R'≡Ω_L/Ω_K is twice the Rossby number. This implies that α= ˜α/R', where √˜αc is the root mean square of the turbulent velocity.

Properly considering this effect can naturally explain the phenomenon of reduced particle scale height produced in forced turbulence shear box simulations, which may help to interpret recent observations of edge-on disks; the paper also proposes broader implications for observational results. For cases where R'>1, the effective particle Stokes number increases, which has significant implications for particle collision dynamics, growth, and planet formation processes.

In summary, the paper attempts to address the issue of how to more accurately describe the length and velocity scales of turbulence in protoplanetary disks when the turnover frequency of the largest eddies is greater than or equal to the Keplerian frequency, as well as the impact of this phenomenon on particle behavior and the planet formation process.