Update Gravitational Potential including indirect terms authored by Elena Lega's avatar Elena Lega
......@@ -29,9 +29,9 @@ To select this choice for a simulation in the config file:
```
<!--![SmoothingPotWiki](uploads/73ba8941570d358cbad84004cbd41548/SmoothingPotWiki.png)
<!--![SmoothingPotWiki](uploads/6b7135f781b548cc70e40c0d7e37bc2b/SmoothingPotWiki.png)
-->
<img src="uploads/73ba8941570d358cbad84004cbd41548/SmoothingPotWiki.png"
<img src="uploads/6b7135f781b548cc70e40c0d7e37bc2b/SmoothingPotWiki.png"
width="320" height="240">
......@@ -40,7 +40,7 @@ $`d -> 0`$. The choice of the smoothing length depends mainly on the resolution:
- In the 2 Dimensional case
we use a potential of the form:
$`\Phi _p^{\epsilon} = -{m_pG\over {\sqrt {d^2+\epsilon^2}}}`$
$`\Phi _p^{\epsilon} = -{m_pG\over {\sqrt {d^2+(\epsilon R_H)^2}}}`$
In 2 dimensional models we do not only have the problem of avoiding the singularity but we choose a potential that allows to mimic the average influence that the planet would have on the vertical gas column.
A value often used for the smoothing is $`\epsilon = 0.7H`$ where $`H`$ is the disc scale height.
To select this choice for a simulation in the config file:
......
......