HyperReducedCompressibleUnSteadyNS.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
    reducedSteadyNS
Description
    A reduced problem for the stationary NS equations
SourceFiles
    reducedSteadyNS.C
\*---------------------------------------------------------------------------*/

 
#ifndef HyperReducedCompressibleUnSteadyNS_H
#define HyperReducedCompressibleUnSteadyNS_H
 
#include "fvCFD.H"
#include "IOmanip.H"
#include "ReducedProblem.H"
#include "UnsteadyProblem.H"
#include "ReducedCompressibleSteadyNS.H"
#include "ITHACAutilities.H"
#include <Eigen/Eigen>
#include "Modes.H"
#include "CompressibleUnSteadyRhoPimple.H"
 
/*---------------------------------------------------------------------------*\
              Class HyperReducedCompressibleUnSteadyNS Declaration
\*---------------------------------------------------------------------------*/


class HyperReducedCompressibleUnSteadyNS: public ReducedCompressibleSteadyNS
{
    private:
 
    public:
        // Constructors
        HyperReducedCompressibleUnSteadyNS();

        explicit HyperReducedCompressibleUnSteadyNS(CompressibleUnSteadyRhoPimple&
                problem);
 
        ~HyperReducedCompressibleUnSteadyNS() {};
 
        // Functions

        void SolveHyperReducedSys(int NmodesUproj,
                                  int NmodesPproj,
                                  int NmodesEproj,
                                  fileName folder = "./ITHACAoutput/Online/");

        void setOnlineVelocity(Eigen::MatrixXd vel);

        void projectReducedOperators(int NmodesUproj, int NmodesPproj, int NmodesEproj);

        std::tuple<Eigen::MatrixXd, Eigen::VectorXd> HyperReducedSys(
            Eigen::SparseMatrix<double>& S,
            Eigen::VectorXd& se,
            List<label>& MagicPoints,
            Eigen::MatrixXd& Modes);

        bool checkWrite(Time& timeObject);
        // Variables

        volVectorModes Umodes;
        volScalarModes Emodes, Pmodes;
        scalar startTime = 0.0;
        scalar finalTime = 0.0;
        scalar timeStep = 0.0;
        scalar writeEvery = timeStep;
        scalar nextWrite = 0.0;
        PtrList<volScalarField> Pfield, Efield;
        PtrList<volVectorField> Ufield;
        // /// List of snapshots for the solution for eddy viscosity
        // PtrList<volScalarField> nutFields;
 
        // /// List of POD modes for eddy viscosity
        volScalarModes nutModes;

        CompressibleUnSteadyRhoPimple* problem;

        Eigen::MatrixXd vel_now;

        int counter = 0;
 
        //Reduced operators

        Eigen::MatrixXd projGradModP;
};
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
 
 
#endif