ReducedFsi.H

/*---------------------------------------------------------------------------*\
 ██╗████████╗██╗  ██╗ █████╗  ██████╗ █████╗       ███████╗██╗   ██╗
 ██║╚══██╔══╝██║  ██║██╔══██╗██╔════╝██╔══██╗      ██╔════╝██║   ██║
 ██║   ██║   ███████║███████║██║     ███████║█████╗█████╗  ██║   ██║
 ██║   ██║   ██╔══██║██╔══██║██║     ██╔══██║╚════╝██╔══╝  ╚██╗ ██╔╝
 ██║   ██║   ██║  ██║██║  ██║╚██████╗██║  ██║      ██║      ╚████╔╝
 ╚═╝   ╚═╝   ╚═╝  ╚═╝╚═╝  ╚═╝ ╚═════╝╚═╝  ╚═╝      ╚═╝       ╚═══╝
 
* In real Time Highly Advanced Computational Applications for Finite Volumes
* Copyright (C) 2017 by the ITHACA-FV authors
-------------------------------------------------------------------------------
License
This file is part of ITHACA-FV
ITHACA-FV is free software: you can redistribute it and/or modify
it under the terms of the GNU Lesser General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.
ITHACA-FV is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU Lesser General Public License for more details.
You should have received a copy of the GNU Lesser General Public License
along with ITHACA-FV. If not, see <http://www.gnu.org/licenses/>.
Class
ReducedFsi
Description
Reduced classs for fluid structure interaction  problem with dynamic mesh
SourceFiles
ReducedFsi.C
\*---------------------------------------------------------------------------*/
#ifndef ReducedFsi_H
#define ReducedFsi_H
 
#include "fsiBasic.H"
#include "ITHACAPOD.H"
#include "ReducedProblem.H"
#include "pointConstraints.H"
#include "mathematicalConstants.H"
#include "ReducedSimpleSteadyNS.H"
 
class ReducedFsi: public reducedSimpleSteadyNS
{
    public:
        // Constructors
        ReducedFsi();
        explicit ReducedFsi(fsiBasic& problem);
        ~ReducedFsi() {};

        fsiBasic* problem;
        volScalarModes Pmodes;
        volVectorModes Umodes;
        pointVectorModes Dmodes;
        PtrList<dynamicFvMesh> OnlineMeshes;

        scalar startTime = 0.0;
        scalar finalTime = 0.0;
        scalar timeStep = 0.0;
        scalar writeEvery = timeStep;
        scalar nextWrite = 0.0;
        List<scalar> centerofmassx;
        List<scalar> centerofmassy;
        List<scalar> centerofmassz;
        List<scalar> romforcey;
        List<scalar> romforcex;
        // List of scaler for saving interpolation  coefficients using RBF
        List<scalar> pdcoeffrbf;
        List<Eigen::MatrixXd> CoeffU;
        List<Eigen::MatrixXd> CoeffP;
 
        PtrList<volScalarField> PredFields;
        PtrList<volVectorField> UredFields;
        PtrList<pointVectorField> Dfield;
        PtrList<volScalarField> Pfield;
        PtrList<volVectorField> Ufield;
        PtrList<pointField> ListOfpoints;
        label counter = problem->counter;
        //std::ofstream Errors;
 
        // Functions
 

        void solveOnline_Pimple(int NmodesUproj,
                                int NmodesPproj,
                                int NmodesDproj,
                                fileName folder = "./ITHACAoutput/Online/");
        void restart();
        bool checkWrite(Time& timeObject);
        void PODI(Eigen::MatrixXd coeffL2, Eigen::MatrixXd muu, label NPdModes);
        void prepareRomData(const word& outputPath);
};
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
#endif