... | @@ -35,14 +35,14 @@ To select this choice for a simulation in the config file: |
... | @@ -35,14 +35,14 @@ To select this choice for a simulation in the config file: |
|
width="320" height="240">
|
|
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.
|
|
$`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
|
|
- 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 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 $`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:
|
|
To select this choice for a simulation in the config file:
|
|
```
|
|
```
|
|
Smoothing
|
|
Smoothing
|
... | @@ -52,7 +52,7 @@ To select this choice for a simulation in the config file: |
... | @@ -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.
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|