Update Gravitational Potential including indirect terms authored by Elena Lega's avatar Elena Lega
......@@ -8,14 +8,14 @@ The gravitational potential entering in the source terms of the Navier-Stokes eq
we use a cubic-potential of the form:
$`\Phi _p = \left\lbrace \begin{array}{ll}
-{m_pG\over d} & d > r_{\rm sm} \\
-{m_pG\over d} & d > r_{\rm sm}} \\
-{m_pG\over d}f({d\over r_{\rm sm}}) & d\leq r_{\rm sm}
\end{array} \right.`$
with $`f({d\over r_{\rm sm}}) = \left [ \left( {d\over r_{\rm sm}}\right)^4-2\left( {d\over r_{\rm sm}}\right)^3+2{d\over r_{\rm sm}} \right]`$;
$`d`$ is the distance from the disc element to the planet, and $`r_{\rm sm}`$ the
smoothing length:
$`r_{\rm sm} = \alpha _{sm} R_H`$, where $`\alpha _{sm}`$ is a smoothing parameter and
$`r_{\rm sm} = \epsilon R_H`$, where $`\epsilon`$ is a smoothing parameter and
$`R_H`$ is the Hill radius of the planet.
To select this choice for a simulation in the config file:
......
......