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