inverseLaplacianProblem.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
    inverseLaplacianProblem
Description
    A general full order implementation of an inverse boundary condition estimation
    problem
SourceFiles
    inverseLaplacianProblem.C
\*---------------------------------------------------------------------------*/
 
 
 
#ifndef inverseLaplacianProblem_H
#define inverseLaplacianProblem_H
#include <iostream>
#include "fvCFD.H"
#include "interpolation.H"
#include "fvOptions.H"
#include "simpleControl.H"
#include "IOmanip.H"
#include "Time.H"
#include "laplacianProblem.H"
#include "ITHACAPOD.H"
#include "ITHACAutilities.H"
#include "ITHACAregularization.H"
#include <Eigen/Dense>
#include <cmath>
#include "Foam2Eigen.H"
#include "cnpy.H"
#include "mixedFvPatchFields.H"
#include "cellDistFuncs.H"
#include "reductionProblem.H"
#include "Modes.H"
#define _USE_MATH_DEFINES
 
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
/*---------------------------------------------------------------------------*\
                 Class inverseLaplacianProblem Declaration
\*---------------------------------------------------------------------------*/
 
 
class inverseLaplacianProblem: public laplacianProblem
{
 
    public:
        // Constructors
        inverseLaplacianProblem();
        inverseLaplacianProblem(int argc, char* argv[]);
        virtual ~inverseLaplacianProblem() {};
 
        ITHACAparameters* para;
 
        autoPtr<volScalarField> _T;
 
        mutable autoPtr<fvMesh> _mesh;
 
        autoPtr<simpleControl> _simple;
 
        autoPtr<fv::options> _fvOptions;
 
        autoPtr<Time> _runTime;
 
        dimensionedScalar DT;
 
        double k;
 
        double H;
 
        bool thermocouplesRead = 0;
 
        int thermocouplesNum;
 
        double J;
        double L2norm;
        double LinfNorm;
 
        scalar homogeneousBC = 0.0;
 
        List<scalar> homogeneousBCcoldSide;
 
        List<scalar> Tf;
 
        List<scalar> refGrad;
 
        List<scalar> valueFraction;
 
        label hotSide_ind;
 
        label coldSide_ind;
 
        List<scalar> g;
 
        List<List<scalar>> gList;
 
        List<scalar> gTrue;
 
        List<scalar> faceCellArea;
 
        List<vector> thermocouplesPos;
 
        List<int> thermocouplesCellID;
 
        List<int> thermocouplesCellProc;
 
        Eigen::VectorXd Tmeas;
 
        Eigen::VectorXd Tdirect;
 
        Eigen::VectorXd Tdiff;
 
        label nProcs;
 
 
        // Functions
 
        //--------------------------------------------------------------------------
        void set_g();
 
        //--------------------------------------------------------------------------
        volScalarField list2Field(List<scalar> list, scalar innerField = 0.0);
 
        //--------------------------------------------------------------------------
        void set_valueFraction();
 
        //--------------------------------------------------------------------------
        virtual void solveTrue() {};
 
        //--------------------------------------------------------------------------
        virtual void assignDirectBC();
 
        //--------------------------------------------------------------------------
        void solveDirect();
 
        //--------------------------------------------------------------------------
        void solve(const char* problem);
 
        //--------------------------------------------------------------------------
        virtual void readThermocouples();
 
        //--------------------------------------------------------------------------
        Eigen::VectorXd fieldValueAtThermocouples(volScalarField& field);
 
        //--------------------------------------------------------------------------
        void differenceBetweenDirectAndMeasure();
 
        //--------------------------------------------------------------------------
        Foam::vector cellDim(const faceList& ff, const pointField& pp,
                             const cell& cc, labelList pLabels, pointField pLocal);
 
        //--------------------------------------------------------------------------
        void restart();
};
 
#endif