... | ... | @@ -36,4 +36,10 @@ $`\Sigma_0`$=Density.Start |
|
|
|
|
|
We recall that in this case the initial azimuthal velocity profile
|
|
|
is given by centrifugal balance for the above density profile, i.e.:
|
|
|
$` v_{\theta}(r) = \Omega_K r\sqrt{1-h_0^2(2-\frac{r}{\sigma_c}(1-\tanh(\frac {r-r_c} {\sigma_c})))}`$ |
|
|
\ No newline at end of file |
|
|
$` v_{\theta}(r) = \Omega_K r\sqrt{1-h_0^2(2-\frac{r}{\sigma_c}(1-\tanh(\frac {r-r_c} {\sigma_c})))}`$
|
|
|
|
|
|
In order to trigger the instability in the region $r\sim r_c$ we have introduced small radial velocity perturbations as:
|
|
|
$` v_r (r,\theta) = v_{\theta}(r) \psi(r,\theta)\exp-\frac {(r-r_c)^2}{10\sigma_c^2}`$
|
|
|
|
|
|
where $`\psi(r,\theta) \in [-10^{-2},10^{-2}]`$ is a uniformly distributed random variable drawn for each grid cell.
|
|
|
|