33 namespace tetherPotentials
59 tetherPotentialProperties.
subDict(typeName +
"Coeffs")
61 mu_(
readScalar(pitchForkRingCoeffs_.lookup(
"mu"))),
62 alpha_(
readScalar(pitchForkRingCoeffs_.lookup(
"alpha"))),
63 rOrbit_(
readScalar(pitchForkRingCoeffs_.lookup(
"rOrbit")))
71 scalar
p =
sqrt(r.
x()*r.
x() + r.
y()*r.
y());
73 scalar pMinusRSqr = (p - rOrbit_)*(p - rOrbit_);
75 return -0.5 * mu_ * pMinusRSqr
76 + 0.25 * pMinusRSqr * pMinusRSqr
77 + 0.5 * alpha_ * r.
z() * r.
z();
83 scalar
p =
sqrt(r.
x()*r.
x() + r.
y()*r.
y());
85 scalar pMinusR = (p - rOrbit_);
89 (mu_ - pMinusR * pMinusR) * pMinusR * r.
x()/
p,
90 (mu_ - pMinusR * pMinusR) * pMinusR * r.
y()/
p,
100 pitchForkRingCoeffs_ =
101 tetherPotentialProperties.
subDict(typeName +
"Coeffs");
103 pitchForkRingCoeffs_.
lookup(
"mu") >> mu_;
104 pitchForkRingCoeffs_.
lookup(
"alpha") >> alpha_;
105 pitchForkRingCoeffs_.
lookup(
"rOrbit") >> rOrbit_;