Update Gravitational Potential including indirect terms authored by Elena Lega's avatar Elena Lega
...@@ -43,17 +43,17 @@ The choice of the smoothing length depends mainly on the resolution: as a rule o ...@@ -43,17 +43,17 @@ The choice of the smoothing length depends mainly on the resolution: as a rule o
- **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 R_H)^2}}}`$ $`\Phi _p^{\epsilon} = -{m_pG\over {\sqrt {d^2+(\epsilon H_p)^2}}}`$
In 2 dimensional models we certainly have the problem of avoiding the singularity as well as in the 3D case; however the main constraint is to 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 certainly have the problem of avoiding the singularity as well as in the 3D case; however the main constraint is to 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 lenght is $`0.7H`$ where $`H`$ is the disc scale height. We remark that the smoothing parameter $`\epsilon`$ is considered in the code in units of the Hill radius and not in disc scale height. A value often used for the smoothing lenght is $`0.6H`$ where $`H`$ is the disc scale height. We remark that the smoothing parameter $`\epsilon`$ is considered in the code in units of the disk scale height at the planet location when selecting ``` Flat true```
To select this choice for a simulation in the config file: Precisely, select this choice for a simulation in the config file:
``` ```
Smoothing Smoothing
{ {
Flat true # select the potential for 2D (flat) case Flat true # select the potential for 2D (flat) case
Size 0.4 # smoothing parameter in units of the planet Hill radius Size 0.4 # smoothing parameter in units of the scale height at the planet position
} }
``` ```
... ...
......