DAMotionTranslationCoupled.H

Go to the documentation of this file.

/*---------------------------------------------------------------------------*\
 
    DAFoam  : Discrete Adjoint with OpenFOAM
    Version : v3
 
    Description:
        Child class for coupled translation motion (FSI)
        with the Mass-spring-damper model
        
        M * a + C * V + K * y = F
 
        where M is the mass, a is the acceleration, V is the velocity, 
        C is the damping coefficient, K is the spring constant, y is 
        the displacement, and F is the external force
 
        To solve this equation, we can break it into two ODEs
 
        y1 = y;   y2 = dy/dt
 
        dy1/dt = V = y2
        dy2/dt = a = (F - C * y2 - K * y1) / M
 
        We can then use the RK2 method to solve the above two ODEs
 
\*---------------------------------------------------------------------------*/
 
#ifndef DAMotionTranslationCoupled_H
#define DAMotionTranslationCoupled_H
 
#include "DAMotion.H"
#include "addToRunTimeSelectionTable.H"
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
namespace Foam
{
 
/*---------------------------------------------------------------------------*\
      Class DAMotionTranslationCoupled Declaration
\*---------------------------------------------------------------------------*/
 
class DAMotionTranslationCoupled
    : public DAMotion
{
 
protected:
    scalar M_;
 
    scalar C_;
 
    scalar K_;
    
    scalar forceScale_;
 
    vector direction_ = vector::zero;
 
    scalar y_;
 
    scalar V_;
 
    scalar y0_;
 
    scalar V0_;
 
    scalar t_;
 
public:
    TypeName("translationCoupled");
    // Constructors
 
    //- Construct from components
    DAMotionTranslationCoupled(
        const dynamicFvMesh& mesh,
        const DAOption& daOption);
 
    //- Destructor
    virtual ~DAMotionTranslationCoupled()
    {
    }
 
    virtual void correct();
};
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
} // End namespace Foam
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
#endif
 
// ************************************************************************* //