drazkowska2016

dustpylib.planetesimals.formation.drazkowska2016(OmegaK, rho_dust, rho_gas, Sigma_dust, St, p2g_crit=1.0, St_crit=0.01, zeta=0.01)[source]

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

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

  • p2g_crit (float, optional, default: 1.) – Critical midplane pebbles-to-gas ratio of particles above St_crit above which planetesimal formation is triggered

  • St_crit (float, optional, default: 0.01) – Critical Stokes number above which dust particles contribute to trigger planetesimal formation

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

Returns:

S – Dust source terms due to planetesimal formation

Return type:

array-like, (Nr, Nm)