fsiBasic.H

Go to the documentation of this file.

/*---------------------------------------------------------------------------*\
     ██╗████████╗██╗  ██╗ █████╗  ██████╗ █████╗       ███████╗██╗   ██╗
     ██║╚══██╔══╝██║  ██║██╔══██╗██╔════╝██╔══██╗      ██╔════╝██║   ██║
     ██║   ██║   ███████║███████║██║     ███████║█████╗█████╗  ██║   ██║
     ██║   ██║   ██╔══██║██╔══██║██║     ██╔══██║╚════╝██╔══╝  ╚██╗ ██╔╝
     ██║   ██║   ██║  ██║██║  ██║╚██████╗██║  ██║      ██║      ╚████╔╝
     ╚═╝   ╚═╝   ╚═╝  ╚═╝╚═╝  ╚═╝ ╚═════╝╚═╝  ╚═╝      ╚═╝       ╚═══╝
 
 * 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
    fsiBasic
Description
    fluid structure interaction  problem with dynamic mesh
SourceFiles
    fsiBasic.C
\*---------------------------------------------------------------------------*/

 
#ifndef fsiBasic_H
#define fsiBasic_H
#include "fvCFD.H"
#include "dynamicFvMesh.H"
#include "singlePhaseTransportModel.H"
#include "turbulentTransportModel.H"
#include "pimpleControl.H"
#include "fvOptions.H"
#include "CorrectPhi.H"
#include "sixDoFRigidBodyMotionSolver.H"
#include "sixDoFRigidBodyMotion.H"
#include "unsteadyNS.H"
#include "ITHACAstream.H"
#include "ITHACAutilities.H"
#include "ITHACAPOD.H"
#include "forces.H"
//#include "MapConsistentVolFields.H"
#include "meshToMesh0.H"
//#include "IOmanip.H"
#include "steadyNS.H"
#include "UnsteadyProblem.H"
//#include "IOstreams.H"
//#include "objectRegistry.H"
//#include "localEulerDdtScheme.H"
//#include "reductionProblem.H"
#include <iostream>
 
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
/*---------------------------------------------------------------------------*\
                        Class FsiBasic Declaration
\*---------------------------------------------------------------------------*/


class fsiBasic:  public steadyNS, public UnsteadyProblem
{
    public:
        // Constructors
        fsiBasic();

        fsiBasic(int argc, char* argv[]);
        ~fsiBasic() {};
 
        autoPtr<volScalarField> _p;
        // /// Dynamic mesh field
        autoPtr<dynamicFvMesh> meshPtr;
        autoPtr<fvMesh> Mesh0;
        // PtrList of dynamicMesh
        PtrList<dynamicFvMesh> meshes;
        autoPtr<pimpleControl> _pimple;
        autoPtr<sixDoFRigidBodyMotionSolver> sDRBMS;
        pointField point0;
        faceList faces0;
        cellList celllist0;
 
        autoPtr<surfaceVectorField> _Uf;
        PtrList<pointVectorField> DFields;
        List<scalar> centerofmassx;
        List<scalar> centerofmassy;
        List<scalar> centerofmassz;
        List<scalar> omegaz;
        PtrList<pointVectorField> Dfield;
        PtrList<volVectorField> Ufield;
        PtrList<surfaceVectorField> NormalFields;
        pointVectorModes Dmodes;
        autoPtr<pointVectorField> _pointDisplacement;
        List<scalar> fomforcey;
        List<scalar> fomforcex;
        Eigen::MatrixXd coeffL2;
 
        Eigen::MatrixXd CylDispl;
        Eigen::MatrixXd CylRot;
        Eigen::VectorXd pdCoeff;

        std::vector<SPLINTER::DataTable*> samples;

        std::vector<SPLINTER::RBFSpline*> rbfSplines;
 
        bool  correctPhi;
        bool  checkMeshCourantNo;
        bool  moveMeshOuterCorrectors;
        //IOdictionary* dyndict;
        autoPtr<IOdictionary> dyndict;
        autoPtr<IOobject> oMesh;
        //const DimensionedField<scalar, volMesh>* V0Ptr_;
 
        // Functions
        //--------------------------------------------------------------------------
        void truthSolve(label paramIndex,  fileName folder = "./ITHACAoutput/Offline/");
        void restart();

        void change_viscosity(double mu);
        void change_stiffness(scalar& mu);
        void updateStiffnessAndRebuildSolver(scalar& newMu,
                                             word param_name = "stiffness");
        void exportFoamFieldToNpy(const word& outputDir,
                                  const word& fileName,
                                  const List<scalar>& foamField);

        void prepareFoamData(const word& outputPath);
        void loadCentreOfMassY(const fileName& baseDir = "ITHACAoutput");
 
        /*
         mapFieldToMesh0( const volScalarField& sourceField,
                  const Foam::fvMesh& mesh0,
                  volScalarField& mappedField);
 
         mapFieldToMesh0( const volVectorField& sourceField,
                  const Foam::fvMesh& mesh0,
                  volVectorField& mappedField);*/
};
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
#endif