TurbDiffusionHyperreduction.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
    TurbDiffusionHyperreduction
Description
    Class to compute and hyperreduce the non-polynomial diffusion for a NON-Stationary turbulent NS problem
SourceFiles
    TurbDiffusionHyperreduction.C
 
\*---------------------------------------------------------------------------*/

 
#ifndef TurbDiffusionHyperreduction_H
#define TurbDiffusionHyperreduction_H
 
#include "UnsteadyNSTurb.H"
#include "PODTemplate.H"
#include "PODTemplateH1.H"
#include "hyperReduction.templates.H"
 
 
template<typename volF, typename T, typename S>
class TurbDiffusionHyperreduction
{
  protected :
    StoredParameters* m_parameters;

    autoPtr<UnsteadyNSTurb> ithacaUnsteadyNSTurb;

    autoPtr<HyperReduction<PtrList<S>>> ithacaHyperreduction;

    label l_nSnapshot;

    label l_startTime;

    label l_nSnapshotSimulation;

    word snapshotsFolder;
    label POD_nSnapshot = -1;
    label POD_startTime = -1;
    label POD_endTime = -1;

    T template_HRInterpField;

    S template_HRSnapshotsField;

    PtrList<volF> f_spatialModesU; // TO DO : need to be duplicated (local => global) if //

    volF* f_meanU; // TO DO : need to be duplicated (local => global) if //
 
    // Spatial POD-HR modes
    PtrList<S> f_spatialModesHR;
 
    // Temporal modes of HR field
    Eigen::MatrixXd m_temporalModesHR;
 
    // Test temporal modes of HR field
    Eigen::MatrixXd m_temporalModesHRSimulation;
 
    // Covariance matrix the resolved modes of HR field
    Eigen::MatrixXd covMatrixHR;

    word folder_HR;

    Foam::Time runTime2;

    bool mgPoints_0_Neighborhoods_1;
 
    // Name of the HR method: DEIM or ECP
    word HRMethod;
 
    // Name of the interpolated field: fullStressFunction or nut or their reduced version
    word HRInterpField;
 
    // Name of the snapshots' field if ECP
    word HRSnapshotsField;
 
    // Name of the algorithm choice if ECP: Global or EachMode
    word ECPAlgo;
 
    // Centered field to perform the hyperreduction
    bool interpFieldCenteredOrNot;
 
 
 
  public :
 
        //--------------------------------------------------------------------------
        TurbDiffusionHyperreduction(Parameters* myParameters, bool mgPoints_0_Neighborhoods_1,
                                    T& template_HRInterpField, S& template_HRSnapshotsField,
                                    volF*& meanU, PtrList<volF>& spatialModesU);
 
        //--------------------------------------------------------------------------
        ~TurbDiffusionHyperreduction();
 
 
        //--------------------------------------------------------------------------
        void precomputeTurbDiffusionFunctions(word& fieldToCompute);
 
        //--------------------------------------------------------------------------
        template<typename volGradF>
        void computeDefTensor(volGradF tensorMeanU);
 
        //--------------------------------------------------------------------------
        void common_MeanHRInterpField();
       
        //--------------------------------------------------------------------------
        void computeTurbDiffusionHyperreduction();
 
};
 
 
#endif