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: ...@@ -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"> width="320" height="240">
...@@ -40,7 +40,7 @@ $`d -> 0`$. The choice of the smoothing length depends mainly on the resolution: ...@@ -40,7 +40,7 @@ $`d -> 0`$. The choice of the smoothing length depends mainly on the resolution:
- In the 2 Dimensional case - In the 2 Dimensional case
we use a potential of the form: 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. 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. 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: To select this choice for a simulation in the config file:
... ...
......