symplectic Class Reference

Symplectic 2nd-order explicit time-integrator for 6DoF solid-body motion. More...

Inheritance diagram for symplectic:
Collaboration diagram for symplectic:

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< sixDoFSolverNew (const dictionary &dict, sixDoFRigidBodyMotion &body)
 
- Protected Member Functions inherited from sixDoFSolver
pointcentreOfRotation ()
 Return the current centre of rotation. More...
 
tensorQ ()
 Return the orientation. More...
 
vectorv ()
 Return non-const access to vector. More...
 
vectora ()
 Return non-const access to acceleration. More...
 
vectorpi ()
 Return non-const access to angular momentum. More...
 
vectortau ()
 Return non-const access to torque. More...
 
const pointcentreOfRotation0 () const
 Return the centre of rotation at previous time-step. More...
 
const tensorQ0 () const
 Return the orientation at previous time-step. More...
 
const vectorv0 () const
 Return the velocity at previous time-step. More...
 
const vectora0 () const
 Return the acceleration at previous time-step. More...
 
const vectorpi0 () const
 Return the angular momentum at previous time-step. More...
 
const vectortau0 () 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, vectorrotate (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
sixDoFRigidBodyMotionbody_
 The rigid body. More...
 

Detailed Description

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;
}
See also
Foam::sixDoFSolvers::CrankNicolson Foam::sixDoFSolvers::Newmark
Source files

Definition at line 74 of file symplectic.H.

Constructor & Destructor Documentation

◆ symplectic()

symplectic ( const dictionary dict,
sixDoFRigidBodyMotion body 
)

Construct from a dictionary and the body.

Definition at line 43 of file symplectic.C.

◆ ~symplectic()

~symplectic ( )
virtual

Destructor.

Definition at line 55 of file symplectic.C.

Member Function Documentation

◆ TypeName()

TypeName ( "symplectic"  )

Runtime type information.

◆ solve()

void solve ( bool  firstIter,
const vector fGlobal,
const vector tauGlobal,
scalar  deltaT,
scalar  deltaT0 
)
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().

Here is the call graph for this function:

The documentation for this class was generated from the following files: