... | ... | @@ -47,13 +47,13 @@ $`\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.
|
|
|
|
|
|
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```
|
|
|
Precisely, select this choice for a simulation in the config file:
|
|
|
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```.
|
|
|
Precisely, for a 2D simulation the config file will have:
|
|
|
```
|
|
|
Smoothing
|
|
|
{
|
|
|
Flat true # select the potential for 2D (flat) case
|
|
|
Size 0.4 # smoothing parameter in units of the scale height at the planet position
|
|
|
Size 0.6 # smoothing parameter in units of the scale height at the planet position
|
|
|
}
|
|
|
```
|
|
|
|
... | ... | |