Symplectic 2nd-order explicit time-integrator for 6DoF solid-body motion. More...
Public Member Functions | |
TypeName ("symplectic") | |
Runtime type information. More... | |
symplectic (const dictionary &dict, sixDoFRigidBodyMotion &body) | |
Construct from a dictionary and the body. More... | |
virtual | ~symplectic () |
Destructor. More... | |
virtual void | solve (bool firstIter, const vector &fGlobal, const vector &tauGlobal, scalar deltaT, scalar deltaT0) |
Drag coefficient. More... | |
Public Member Functions inherited from sixDoFSolver | |
TypeName ("sixDoFSolver") | |
Runtime type information. More... | |
declareRunTimeSelectionTable (autoPtr, sixDoFSolver, dictionary,(const dictionary &dict, sixDoFRigidBodyMotion &body),(dict, body)) | |
sixDoFSolver (sixDoFRigidBodyMotion &body) | |
virtual | ~sixDoFSolver () |
Destructor. More... | |
Additional Inherited Members | |
Static Public Member Functions inherited from sixDoFSolver | |
static autoPtr< sixDoFSolver > | New (const dictionary &dict, sixDoFRigidBodyMotion &body) |
Protected Member Functions inherited from sixDoFSolver | |
point & | centreOfRotation () |
Return the current centre of rotation. More... | |
tensor & | Q () |
Return the orientation. More... | |
vector & | v () |
Return non-const access to vector. More... | |
vector & | a () |
Return non-const access to acceleration. More... | |
vector & | pi () |
Return non-const access to angular momentum. More... | |
vector & | tau () |
Return non-const access to torque. More... | |
const point & | centreOfRotation0 () const |
Return the centre of rotation at previous time-step. More... | |
const tensor & | Q0 () const |
Return the orientation at previous time-step. More... | |
const vector & | v0 () const |
Return the velocity at previous time-step. More... | |
const vector & | a0 () const |
Return the acceleration at previous time-step. More... | |
const vector & | pi0 () const |
Return the angular momentum at previous time-step. More... | |
const vector & | tau0 () const |
Return the torque at previous time-step. More... | |
scalar | aDamp () const |
Acceleration damping coefficient (for steady-state simulations) More... | |
tensor | tConstraints () const |
Translational constraint tensor. More... | |
tensor | rConstraints () const |
Rotational constraint tensor. More... | |
Tuple2< tensor, vector > | rotate (const tensor &Q0, const vector &pi, const scalar deltaT) const |
Apply rotation tensors to Q0 for the given torque (pi) and deltaT. More... | |
void | updateAcceleration (const vector &fGlobal, const vector &tauGlobal) |
Update and relax accelerations from the force and torque. More... | |
Protected Attributes inherited from sixDoFSolver | |
sixDoFRigidBodyMotion & | body_ |
The rigid body. More... | |
Symplectic 2nd-order explicit time-integrator for 6DoF solid-body motion.
Reference:
Dullweber, A., Leimkuhler, B., & McLachlan, R. (1997). Symplectic splitting methods for rigid body molecular dynamics. The Journal of chemical physics, 107(15), 5840-5851.
Can only be used for explicit integration of the motion of the body, i.e. may only be called once per time-step, no outer-correctors may be applied. For implicit integration with outer-correctors choose either CrankNicolson or Newmark schemes.
Example specification in dynamicMeshDict:
solver { type symplectic; }
Definition at line 74 of file symplectic.H.
symplectic | ( | const dictionary & | dict, |
sixDoFRigidBodyMotion & | body | ||
) |
Construct from a dictionary and the body.
Definition at line 43 of file symplectic.C.
|
virtual |
Destructor.
Definition at line 55 of file symplectic.C.
TypeName | ( | "symplectic" | ) |
Runtime type information.
|
virtual |
Drag coefficient.
Implements sixDoFSolver.
Definition at line 61 of file symplectic.C.
References Foam::constant::atomic::a0, Tuple2< Type1, Type2 >::first(), Foam::constant::mathematical::pi(), and Tuple2< Type1, Type2 >::second().