gsPartitionedFSI.h

Go to the documentation of this file.

 
#pragma once
 
#include <gsIO/gsOptionList.h>
 
namespace gismo
{
 
template <class T>
class gsNsTimeIntegrator;
template <class T>
class gsElTimeIntegrator;
template <class T>
class gsALE;
template <class T>
class gsMultiPatch;
 
template <class T>
class gsPartitionedFSI
{
public:
 
    gsPartitionedFSI(gsNsTimeIntegrator<T> & nsSolver,
                     gsMultiPatch<T> & velocity, gsMultiPatch<T> & pressure,
                     gsElTimeIntegrator<T> & elSolver,
                     gsMultiPatch<T> & displacement,
                     gsALE<T> & aleSolver,
                     gsMultiPatch<T> & aleDisplacement, gsMultiPatch<T> & aleVelocity);
 
    static gsOptionList defaultOptions();
 
    gsOptionList & options() { return m_options; }
 
    bool makeTimeStep(T timeStep);
 
    void formVector(const gsMultiPatch<T> & disp, gsMatrix<T> & vector);
 
    void aitken(gsMultiPatch<T> & dispA, gsMultiPatch<T> & dispB,
                gsMultiPatch<T> & dispB2, gsMultiPatch<T> & dispC);
 
    index_t numberIterations() { return numIter; }
    T timeNS() { return nsTime; }
    T timeEL() { return elTime; }
    T timeALE() { return aleTime; }
    T aitkenOmega() { return omega;}
    T residualNormAbs() { return absResNorm;}
    T residualNormRel() { return absResNorm/initResNorm; }
 
protected:
    gsNsTimeIntegrator<T> & m_nsSolver;
    gsMultiPatch<T> & m_velocity;
    gsMultiPatch<T> & m_pressure;
    gsElTimeIntegrator<T> & m_elSolver;
    gsMultiPatch<T> & m_displacement;
    gsALE<T> & m_aleSolver;
    gsMultiPatch<T> & m_ALEdisplacment;
    gsMultiPatch<T> & m_ALEvelocity;
    gsOptionList m_options;
    index_t numIter; // number of iterations at the last time step
    bool converged; // convergence flag
    T nsTime, elTime, aleTime; // component computational times
    T omega; // aitken relaxation parameter
    T absResNorm, initResNorm; // residual norms for convergence cretirion
 
};
 
} // namespace ends
 
#ifndef GISMO_BUILD_LIB
#include GISMO_HPP_HEADER(gsPartitionedFSI.hpp)
#endif