36 namespace sixDoFSolvers
53 aoc_(
dict.getOrDefault<scalar>(
"aoc", 0.5)),
54 voc_(
dict.getOrDefault<scalar>(
"voc", 0.5))
76 updateAcceleration(fGlobal, tauGlobal);
83 & (v0() + aDamp()*deltaT*(aoc_*a() + (1 - aoc_)*
a0()));
87 & (pi0() + aDamp()*deltaT*(aoc_*tau() + (1 - aoc_)*tau0()));
91 centreOfRotation0() + deltaT*(voc_*v() + (1 - voc_)*v0());
95 rotate(Q0(), (voc_*
pi() + (1 - voc_)*pi0()), deltaT);
Six degree of freedom motion for a rigid body.
CrankNicolson(const dictionary &dict, sixDoFRigidBodyMotion &body)
Construct from a dictionary and the body.
A list of keyword definitions, which are a keyword followed by a number of values (eg...
A 2-tuple for storing two objects of dissimilar types. The container is similar in purpose to std::pa...
virtual void solve(bool firstIter, const vector &fGlobal, const vector &tauGlobal, scalar deltaT, scalar deltaT0)
Drag coefficient.
defineTypeNameAndDebug(CrankNicolson, 0)
Macros for easy insertion into run-time selection tables.
addToRunTimeSelectionTable(sixDoFSolver, CrankNicolson, dictionary)
constexpr scalar pi(M_PI)
Crank-Nicolson 2nd-order time-integrator for 6DoF solid-body motion.
A Vector of values with scalar precision, where scalar is float/double depending on the compilation f...
virtual ~CrankNicolson()
Destructor.
const dimensionedScalar a0
Bohr radius: default SI units: [m].