schoonenberg2018

dustpylib.planetesimals.formation.schoonenberg2018(OmegaK, rho_dust, rho_gas, Sigma_dust, St, d2g_crit=1.0, zeta=0.1)[source]

Function calculates the dust source term due to planetesimal formation of Schoonenberg et al. (2018).

Parameters:

  • OmegaK (array-like, (Nr,)) – Keplerian frequency

  • rho_dust (array-like, (Nr, Nm)) – Midplane dust volume density

  • rho_gas (array-like, (Nr,)) – Midplane gas volume density

  • Sigma_dust (array-like, (Nr, Nm)) – Dust surface density

  • St (array-like, (Nr, Nm)) – Stokes numbers

  • d2g_crit (float, optional, default: 1.) – Critical midplane dust-to-gas ratio above which planetesimal formation is triggered

  • zeta (float, optional, default: 0.1) – Planetesimal formation efficiency

Returns:

S – Dust source terms due to planetesimal formation

Return type:

array-like, (Nr, Nm)