Update Gravitational Potential including indirect terms authored by Elena Lega's avatar Elena Lega
......@@ -35,14 +35,14 @@ To select this choice for a simulation in the config file:
width="320" height="240">
The plot shows the effect of smoothing the potential for various values of the smoothing length. In the 3 dimensional case we smooth the potential in order to avoid the singularity in the case
The plot shows the effect of smoothing the potential for various values of the smoothing parameter. In the 3 dimensional case we smooth the potential in order to avoid the singularity in the case
$`d -> 0`$. The choice of the smoothing length depends mainly on the resolution: as a rule of the thumb the minimum smoothing length should correspond to the size of about 4 grid cells.
- In the 2 Dimensional case
we use a potential of the form:
$`\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.
A value often used for the smoothing is $`0.7H`$ where $`H`$ is the disc scale height. We remark that the parameter $`\epsilon`$ is considered in the code in units of the Hill radius and not in disc scale height.
To select this choice for a simulation in the config file:
```
Smoothing
......@@ -52,7 +52,7 @@ To select this choice for a simulation in the config file:
}
```
We remark that the parameter $`\epsilon`$ is considered in the code in units of the Hill radius and not in disc scale height.