ITHACA-FV/msrProblem_8H_source.html

#ifndef msrProblem_H

#define msrProblem_H

#include "fvCFD.H"

#include "singlePhaseTransportModel.H"

#include "turbulentTransportModel.H"

#include "simpleControl.H"

#include "pisoControl.H"

#include "fvOptions.H"

#include "reductionProblem.H"

#include "ITHACAstream.H"

#include "ITHACAforces.H"

#include "volFields.H"

#include <iostream>

#include "IOmanip.H"

#include "IFstream.H"

#include "primitiveFields.H"

#include "FieldFields.H"

#include "scalarMatrices.H"

#include "SortableList.H"

#include "volFieldsFwd.H"

#include "forces.H"

#include "forceCoeffs.H"

#include "volFields.H"

#include <fstream>

#include <sstream>

#include <vector>

#include <string>

#include <stdio.h>

#include "ITHACAPOD.H"

#include <math.h>


class msrProblem: public reductionProblem

{


    public:

        //Constructors

        msrProblem();

        msrProblem(int argc, char* argv[]);

        ~msrProblem() {};


        // Dummy variables to transform msrFoam into a class, createTime.H


        autoPtr<Time> _runTime;


        // Dummy variables to transform msrFoam into a class, createMesh.H


        mutable autoPtr<fvMesh> _mesh;


        // Dummy variables to transform msrFoam into a class, createFields.H


        autoPtr<volVectorField> _U;

        autoPtr<surfaceScalarField> _phi;

        autoPtr<singlePhaseTransportModel> _laminarTransport;

        autoPtr<incompressible::turbulenceModel> turbulence;

        autoPtr<IOMRFZoneList> _MRF;

        autoPtr<volScalarField> _p;

        label pRefCell;

        scalar pRefValue;

        autoPtr<dimensionedScalar> _nu;

        autoPtr<dimensionedScalar> _Pr;

        autoPtr<dimensionedScalar> _Prt;

        autoPtr<dimensionedScalar> _rhoRef;

        autoPtr<volScalarField> _v;

        autoPtr<dimensionedScalar> _CpRef;

        autoPtr<dimensionedScalar> _Tref;

        autoPtr<dimensionedScalar> _TrefXS;

        autoPtr<volScalarField> _TXS;

        autoPtr<volScalarField> _p0;

        autoPtr<volVectorField> _U0;

        autoPtr<surfaceScalarField> _phi0;

        autoPtr<volScalarField> _flux0;

        autoPtr<volScalarField> _prec10;

        autoPtr<volScalarField> _prec20;

        autoPtr<volScalarField> _prec30;

        autoPtr<volScalarField> _prec40;

        autoPtr<volScalarField> _prec50;

        autoPtr<volScalarField> _prec60;

        autoPtr<volScalarField> _prec70;

        autoPtr<volScalarField> _prec80;

        autoPtr<volScalarField> _T0;

        autoPtr<volScalarField> _dec10;

        autoPtr<volScalarField> _dec20;

        autoPtr<volScalarField> _dec30;

        autoPtr<dimensionedScalar> _K0;

        autoPtr<volScalarField> _v0;

        autoPtr<volScalarField> _NSF0;

        autoPtr<volScalarField> _A0;

        autoPtr<volScalarField> _D0;

        autoPtr<volScalarField> _SP0;

        autoPtr<volScalarField> _TXS0;


        // Dummy variables to transform msrFoam into a class, createFields_neutronics.H


        autoPtr<dimensionedScalar> _Keff;

        autoPtr<dimensionedScalar> _IV1;

        autoPtr<dimensionedScalar> _D1_0;

        autoPtr<volScalarField> _D;

        autoPtr<dimensionedScalar> _alfa_D1;

        autoPtr<dimensionedScalar> _A1_0;

        autoPtr<volScalarField> _A;

        autoPtr<dimensionedScalar> _alfa_A1;

        autoPtr<dimensionedScalar> _NSF1_0;

        autoPtr<volScalarField> _NSF;

        autoPtr<dimensionedScalar> _alfa_NSF1;

        autoPtr<dimensionedScalar> _SP1_0;

        autoPtr<volScalarField> _SP;

        autoPtr<dimensionedScalar> _alfa_SP1;

        autoPtr<dimensionedScalar> _lam1;

        autoPtr<dimensionedScalar> _lam2;

        autoPtr<dimensionedScalar> _lam3;

        autoPtr<dimensionedScalar> _lam4;

        autoPtr<dimensionedScalar> _lam5;

        autoPtr<dimensionedScalar> _lam6;

        autoPtr<dimensionedScalar> _lam7;

        autoPtr<dimensionedScalar> _lam8;

        autoPtr<dimensionedScalar> _beta1;

        autoPtr<dimensionedScalar> _beta2;

        autoPtr<dimensionedScalar> _beta3;

        autoPtr<dimensionedScalar> _beta4;

        autoPtr<dimensionedScalar> _beta5;

        autoPtr<dimensionedScalar> _beta6;

        autoPtr<dimensionedScalar> _beta7;

        autoPtr<dimensionedScalar> _beta8;

        autoPtr<dimensionedScalar> _betaTot;

        autoPtr<volScalarField> _flux;

        autoPtr<dimensionedScalar> _decLam1;

        autoPtr<dimensionedScalar> _decLam2;

        autoPtr<dimensionedScalar> _decLam3;

        autoPtr<dimensionedScalar> _decBeta1;

        autoPtr<dimensionedScalar> _decBeta2;

        autoPtr<dimensionedScalar> _decBeta3;

        autoPtr<dimensionedScalar> _decbetaTot;

        autoPtr<dimensionedScalar> _betaTE;

        autoPtr<volScalarField> _prec1;

        autoPtr<volScalarField> _prec2;

        autoPtr<volScalarField> _prec3;

        autoPtr<volScalarField> _prec4;

        autoPtr<volScalarField> _prec5;

        autoPtr<volScalarField> _prec6;

        autoPtr<volScalarField> _prec7;

        autoPtr<volScalarField> _prec8;

        autoPtr<dimensionedScalar> _Sc;

        autoPtr<dimensionedScalar> _Sct;

        autoPtr<volScalarField> _alphat;

        autoPtr<volScalarField> _difft;


        //recirculation time

        double tau;


        // Dummy variables to transform msrFoam into a class, createFields_T.H


        autoPtr<volScalarField> _T;

        autoPtr<volScalarField> _logT;

        autoPtr<volScalarField> _dec1;

        autoPtr<volScalarField> _dec2;

        autoPtr<volScalarField> _dec3;


        // Dummy variables to transform msrFoam into a class, createFvOptions.H


        autoPtr<fv::options> _fvOptions;


        // Dummy variables to transform msrFoam into a class, remaining in .C


        autoPtr<simpleControl> _simple;

        scalar tolerance;

        scalar maxIter;

        scalar cumulativeContErr = 0;


        bool precInBool = false;

        label precinIndex;

        label precoutIndex;

        Eigen::MatrixXd bc_prec;


        // Fields containing the offline solution i.e. the snapshot matrix


        PtrList<volScalarField> Pfield;


        PtrList<volVectorField> Ufield;


        PtrList<volScalarField> Fluxfield;


        PtrList<volScalarField> Prec1field;


        PtrList<volScalarField> Prec2field;


        PtrList<volScalarField> Prec3field;


        PtrList<volScalarField> Prec4field;


        PtrList<volScalarField> Prec5field;


        PtrList<volScalarField> Prec6field;


        PtrList<volScalarField> Prec7field;


        PtrList<volScalarField> Prec8field;


        PtrList<volScalarField> Tfield;


        PtrList<volScalarField> Dec1field;


        PtrList<volScalarField> Dec2field;


        PtrList<volScalarField> Dec3field;


        PtrList<volScalarField> PowerDensfield;


        PtrList<volScalarField> vFields;


        PtrList<volScalarField> DFields;


        PtrList<volScalarField> NSFFields;


        PtrList<volScalarField> AFields;


        PtrList<volScalarField> SPFields;


        PtrList<volScalarField> TXSFields;


        PtrList<volScalarField> Pmodes;


        PtrList<volVectorField> Umodes;


        PtrList<volScalarField> Fluxmodes;


        PtrList<volScalarField> Prec1modes;


        PtrList<volScalarField> Prec2modes;


        PtrList<volScalarField> Prec3modes;


        PtrList<volScalarField> Prec4modes;


        PtrList<volScalarField> Prec5modes;


        PtrList<volScalarField> Prec6modes;


        PtrList<volScalarField> Prec7modes;


        PtrList<volScalarField> Prec8modes;


        PtrList<volScalarField> Tmodes;


        PtrList<volScalarField> Dec1modes;


        PtrList<volScalarField> Dec2modes;


        PtrList<volScalarField> Dec3modes;


        PtrList<volScalarField> vmodes;


        PtrList<volScalarField> Dmodes;


        PtrList<volScalarField> NSFmodes;


        PtrList<volScalarField> Amodes;


        PtrList<volScalarField> SPmodes;


        PtrList<volScalarField> TXSmodes;


        //----------------------------------------------------------------------//

        // Coefficient matrices:


        //Fluid-dynamics matrices:


        Eigen::MatrixXd B_matrix;


        Eigen::MatrixXd M_matrix;


        Eigen::MatrixXd K_matrix;


        List <Eigen::MatrixXd> C_matrix;


        Eigen::MatrixXd P_matrix;


        Eigen::MatrixXd D_matrix;


        List <Eigen::MatrixXd> G_matrix;


        Eigen::MatrixXd BC1_matrix;


        List <Eigen::MatrixXd> BC2_matrix;


        Eigen::MatrixXd BC3_matrix;


        //Neutronics matrices:


        List<Eigen::MatrixXd> LF_matrix;

        Eigen::MatrixXd MF_matrix;

        List<Eigen::MatrixXd> PF_matrix;

        List<Eigen::MatrixXd> AF_matrix;


        ITHACAparameters* para;


        Eigen::MatrixXd PS1_matrix;

        Eigen::MatrixXd PS2_matrix;

        Eigen::MatrixXd PS3_matrix;

        Eigen::MatrixXd PS4_matrix;

        Eigen::MatrixXd PS5_matrix;

        Eigen::MatrixXd PS6_matrix;

        Eigen::MatrixXd PS7_matrix;

        Eigen::MatrixXd PS8_matrix;


        List<Eigen::MatrixXd> ST1_matrix;

        List<Eigen::MatrixXd> ST2_matrix;

        List<Eigen::MatrixXd> ST3_matrix;

        List<Eigen::MatrixXd> ST4_matrix;

        List<Eigen::MatrixXd> ST5_matrix;

        List<Eigen::MatrixXd> ST6_matrix;

        List<Eigen::MatrixXd> ST7_matrix;

        List<Eigen::MatrixXd> ST8_matrix;

        Eigen::MatrixXd MP1_matrix;

        Eigen::MatrixXd MP2_matrix;

        Eigen::MatrixXd MP3_matrix;

        Eigen::MatrixXd MP4_matrix;

        Eigen::MatrixXd MP5_matrix;

        Eigen::MatrixXd MP6_matrix;

        Eigen::MatrixXd MP7_matrix;

        Eigen::MatrixXd MP8_matrix;

        Eigen::MatrixXd LP1_matrix;

        Eigen::MatrixXd LP2_matrix;

        Eigen::MatrixXd LP3_matrix;

        Eigen::MatrixXd LP4_matrix;

        Eigen::MatrixXd LP5_matrix;

        Eigen::MatrixXd LP6_matrix;

        Eigen::MatrixXd LP7_matrix;

        Eigen::MatrixXd LP8_matrix;

        List<Eigen::MatrixXd> FS1_matrix;

        List<Eigen::MatrixXd> FS2_matrix;

        List<Eigen::MatrixXd> FS3_matrix;

        List<Eigen::MatrixXd> FS4_matrix;

        List<Eigen::MatrixXd> FS5_matrix;

        List<Eigen::MatrixXd> FS6_matrix;

        List<Eigen::MatrixXd> FS7_matrix;

        List<Eigen::MatrixXd> FS8_matrix;


        //Energy matrices:


        List<Eigen::MatrixXd> SD1_matrix;

        List<Eigen::MatrixXd> SD2_matrix;

        List<Eigen::MatrixXd> SD3_matrix;

        Eigen::MatrixXd MD1_matrix;

        Eigen::MatrixXd MD2_matrix;

        Eigen::MatrixXd MD3_matrix;

        Eigen::MatrixXd LD1_matrix;

        Eigen::MatrixXd LD2_matrix;

        Eigen::MatrixXd LD3_matrix;

        List<Eigen::MatrixXd> DFS1_matrix;

        List<Eigen::MatrixXd> DFS2_matrix;

        List<Eigen::MatrixXd> DFS3_matrix;


        Eigen::MatrixXd TM_matrix;

        List<Eigen::MatrixXd> TS_matrix;

        Eigen::MatrixXd LT_matrix;

        List<Eigen::MatrixXd> TXS_matrix;

        List<Eigen::MatrixXd> THS1_matrix;

        List<Eigen::MatrixXd> THS2_matrix;

        List<Eigen::MatrixXd> THS3_matrix;

        //--------------------------------------------------------------------------------//


        label NUmodes;

        label NPmodes;

        label NFluxmodes;

        Eigen::VectorXi NPrecmodes;

        Eigen::VectorXi NDecmodes;

        label NTmodes;

        label NCmodes;


        //--------------------------------------------------------------------------------//

        //Lift-functions and homogenization related terms:


        PtrList<volVectorField> liftfield;

        PtrList<volScalarField> liftfieldT;

        PtrList<volVectorField> Uomfield;

        PtrList<volScalarField> Tomfield;

        bool homboolU = false;

        bool homboolT = false;


        //-----------------------------------------------------------------------------//

        //RBF interpolation procedure quantities, i.e. constants changing with

        //temperature


        // Create a SAMPLES for interpolation

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

        // Create a SAMPLES for interpolation

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

        // Create a SAMPLES for interpolation

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

        // Create a SAMPLES for interpolation

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

        // Create a SAMPLES for interpolation

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

        // Create a SAMPLES for interpolation

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

        // Create a RBF splines for interpolation

        std::vector<SPLINTER::RBFSpline*> rbfsplines_v;

        // Create a RBF splines for interpolation

        std::vector<SPLINTER::RBFSpline*> rbfsplines_D;

        // Create a RBF splines for interpolation

        std::vector<SPLINTER::RBFSpline*> rbfsplines_NSF;

        // Create a RBF splines for interpolation

        std::vector<SPLINTER::RBFSpline*> rbfsplines_A;

        // Create a RBF splines for interpolation

        std::vector<SPLINTER::RBFSpline*> rbfsplines_SP;

        // Create a RBF splines for interpolation

        std::vector<SPLINTER::RBFSpline*> rbfsplines_TXS;


        //-----------------------------------------------------------------------------------------//

        // Methods:


        void truthSolve(List<scalar> mu_now);

        //----------------------------------------------------------------------------------


        void msrcoeff(label& NC);

        //----------------------------------------------------------------------------------


        void liftSolve();

        //----------------------------------------------------------------------------------


        void liftSolveT();

        //----------------------------------------------------------------------------------


        void readMSRfields();

        void readMSRfields(std::string& dir);

        //-----------------------------------------------------------------------------------


        void homogenizeU();

        //-----------------------------------------------------------------------------------


        void homogenizeT();

        //--------------------------------------------------------------------------


        void msrgetModesSVD();

        //--------------------------------------------------------------------------


        void msrgetModesEVD();

        //--------------------------------------------------------------------------


        void projectPPE(fileName folder, label NUmodes, label NPmodes, label NFluxmodes,

                        Eigen::VectorXi Nprecmodes, label NTmodes, Eigen::VectorXi Ndecmodes,

                        label NCmodes);

        //--------------------------------------------------------------------------


        Eigen::MatrixXd diffusive_term(label NUmodes, label NPmodes);

        Eigen::MatrixXd pressure_gradient_term(label NUmodes, label NPmodes);

        List < Eigen::MatrixXd > convective_term(label NUmodes, label NPmodes);

        Eigen::MatrixXd mass_term(label NUmodes, label NPmodes);


        Eigen::MatrixXd divergence_term(label NUmodes, label NPmodes);

        List < Eigen::MatrixXd > div_momentum(label NUmodes, label NPmodes);

        Eigen::MatrixXd laplacian_pressure(label NPmodes);

        Eigen::MatrixXd pressure_BC1(label NUmodes, label NPmodes);

        List < Eigen::MatrixXd > pressure_BC2(label NUmodes, label NPmodes);

        Eigen::MatrixXd pressure_BC3(label NUmodes, label NPmodes);


        List<Eigen::MatrixXd> laplacian_flux(label NFluxmodes, label NCmodes);

        Eigen::MatrixXd mass_flux(label NFluxmodes);

        List<Eigen::MatrixXd> prod_flux(label NFluxmodes, label NCmodes);

        List<Eigen::MatrixXd> abs_flux(label NFluxmodes, label NCmodes);

        Eigen::MatrixXd prec_source(label NFluxmodes, label NPrecmodes, label family);


        List<Eigen::MatrixXd> stream_term(label NUmodes, label NPrecmodes,

                                          label family);

        Eigen::MatrixXd prec_mass(label NPrecmodes, label family);

        Eigen::MatrixXd laplacian_prec(label NPrecmodes, label family);

        List<Eigen::MatrixXd> flux_source(label NFluxmodes, label NPrecmodes,

                                          label NCmodes, label family);


        List<Eigen::MatrixXd> stream_dec(label NUmodes, label NDecmodes,

                                         label decgroup);

        Eigen::MatrixXd dec_mass(label NDecmodes, label decgroup);

        Eigen::MatrixXd laplacian_dec(label NDecmodes, label decgroup);

        List<Eigen::MatrixXd> dec_fluxsource(label NFluxmodes, label NDecmodes,

                                             label NCmodes, label decgroup);


        Eigen::MatrixXd mass_temp(label NTmodes);

        List<Eigen::MatrixXd> temp_stream(label NUmodes, label NTmodes);

        Eigen::MatrixXd laplacian_temp(label NTmodes);

        List<Eigen::MatrixXd> temp_heatsource(label NTmodes, label NDecmodes,

                                              label NCmodes, label decgroup);

        List<Eigen::MatrixXd> temp_XSfluxsource(label NTmodes, label NFluxmodes,

                                                label NCmodes);


        //--------------------------------------------------------------------------

        void change_viscosity(double mu);


        //--------------------------------------------------------------------------

        void restart();


    protected:


        //--------------------------------------------------------------------------


        PtrList<volScalarField> choose_group(string field, label ith);


        //--------------------------------------------------------------------------


        template<typename M>

        void savegroupMatrix(string nome, label n, word folder, M matrice);


};


#endif


ITHACAPOD.H
Header file of the ITHACAPOD class.

ITHACAstream.H
Header file of the ITHACAstream class, it contains the implementation of several methods for input ou...

ITHACAparameters
Class for the definition of some general parameters, the parameters must be defined from the file ITH...
Definition ITHACAparameters.H:41

msrProblem::stream_term
List< Eigen::MatrixXd > stream_term(label NUmodes, label NPrecmodes, label family)
precursor eq. methods:
Definition msrProblem.C:1078

msrProblem::MP7_matrix
Eigen::MatrixXd MP7_matrix
precursor mass term-7
Definition msrProblem.H:418

msrProblem::MD1_matrix
Eigen::MatrixXd MD1_matrix
decay heat mass term-1
Definition msrProblem.H:464

msrProblem::homogenizeT
void homogenizeT()
Method to compute the homogenized temperature field, it also sets homboolT=true.
Definition msrProblem.C:1659

msrProblem::ST5_matrix
List< Eigen::MatrixXd > ST5_matrix
precursor stream term-5
Definition msrProblem.H:398

msrProblem::Prec6field
PtrList< volScalarField > Prec6field
List of pointers used to form the prec6 snapshots matrix.
Definition msrProblem.H:219

msrProblem::THS2_matrix
List< Eigen::MatrixXd > THS2_matrix
temperature decay heat source term-2
Definition msrProblem.H:493

msrProblem::LP3_matrix
Eigen::MatrixXd LP3_matrix
precursor laplacian term-3
Definition msrProblem.H:426

msrProblem::SD1_matrix
List< Eigen::MatrixXd > SD1_matrix
decay heat stream term-1
Definition msrProblem.H:458

msrProblem::BC1_matrix
Eigen::MatrixXd BC1_matrix
PPE BC1.
Definition msrProblem.H:350

msrProblem::NSFmodes
PtrList< volScalarField > NSFmodes
List of pointers used to form the NSF snapshosts matrix.
Definition msrProblem.H:312

msrProblem::Prec5modes
PtrList< volScalarField > Prec5modes
List of pointers used to form the prec1 modes.
Definition msrProblem.H:282

msrProblem::PF_matrix
List< Eigen::MatrixXd > PF_matrix
production flux
Definition msrProblem.H:366

msrProblem::choose_group
PtrList< volScalarField > choose_group(string field, label ith)
method to choose one field among precs & decs field can be "prec" or "dec" only if field==prec then i...
Definition msrProblem.C:1585

msrProblem::savegroupMatrix
void savegroupMatrix(string nome, label n, word folder, M matrice)
method to save matrices for precs and decs M can be an Eigen::MatrixXd or List<Eigen::MatrixXd> nome ...
Definition msrProblem.C:1646

msrProblem::Prec4modes
PtrList< volScalarField > Prec4modes
List of pointers used to form the prec1 modes.
Definition msrProblem.H:279

msrProblem::FS3_matrix
List< Eigen::MatrixXd > FS3_matrix
precursor flux source term-3
Definition msrProblem.H:442

msrProblem::BC2_matrix
List< Eigen::MatrixXd > BC2_matrix
PPE BC2.
Definition msrProblem.H:353

msrProblem::ST7_matrix
List< Eigen::MatrixXd > ST7_matrix
precursor stream term-7
Definition msrProblem.H:402

msrProblem::Prec8field
PtrList< volScalarField > Prec8field
List of pointers used to form the prec8 snapshots matrix.
Definition msrProblem.H:225

msrProblem::FS1_matrix
List< Eigen::MatrixXd > FS1_matrix
precursor flux source term-1
Definition msrProblem.H:438

msrProblem::Prec1modes
PtrList< volScalarField > Prec1modes
List of pointers used to form the prec1 modes.
Definition msrProblem.H:270

msrProblem::TXS_matrix
List< Eigen::MatrixXd > TXS_matrix
temperature flux source term TXS
Definition msrProblem.H:489

msrProblem::LF_matrix
List< Eigen::MatrixXd > LF_matrix
laplacian_flux
Definition msrProblem.H:362

msrProblem::Pmodes
PtrList< volScalarField > Pmodes
List of pointers used to form the pressure modes.
Definition msrProblem.H:261

msrProblem::NUmodes
label NUmodes
Number of modes adopted during Galerkin projection.
Definition msrProblem.H:500

msrProblem::MP1_matrix
Eigen::MatrixXd MP1_matrix
precursor mass term-1
Definition msrProblem.H:406

msrProblem::ST8_matrix
List< Eigen::MatrixXd > ST8_matrix
precursor stream term-8
Definition msrProblem.H:404

msrProblem::Uomfield
PtrList< volVectorField > Uomfield
List of pointers used to form the homogeneous velocity snapshots.
Definition msrProblem.H:517

msrProblem::ST4_matrix
List< Eigen::MatrixXd > ST4_matrix
precursor stream term-4
Definition msrProblem.H:396

msrProblem::Prec7modes
PtrList< volScalarField > Prec7modes
List of pointers used to form the prec1 modes.
Definition msrProblem.H:288

msrProblem::vFields
PtrList< volScalarField > vFields
List of pointers used to form the v snapshosts matrix.
Definition msrProblem.H:243

msrProblem::PS3_matrix
Eigen::MatrixXd PS3_matrix
prec_source 3
Definition msrProblem.H:377

msrProblem::Tmodes
PtrList< volScalarField > Tmodes
List of pointers used to form the temperature modes.
Definition msrProblem.H:294

msrProblem::precinIndex
label precinIndex
indexes of inlet and outlet to adopt for precursors boundary conditions
Definition msrProblem.H:184

msrProblem::K_matrix
Eigen::MatrixXd K_matrix
Gradient of pressure matrix.
Definition msrProblem.H:335

msrProblem::G_matrix
List< Eigen::MatrixXd > G_matrix
Divergence of momentum PPE.
Definition msrProblem.H:347

msrProblem::AF_matrix
List< Eigen::MatrixXd > AF_matrix
absorption flux
Definition msrProblem.H:368

msrProblem::MP2_matrix
Eigen::MatrixXd MP2_matrix
precursor mass term-2
Definition msrProblem.H:408

msrProblem::FS2_matrix
List< Eigen::MatrixXd > FS2_matrix
precursor flux source term-2
Definition msrProblem.H:440

msrProblem::LP8_matrix
Eigen::MatrixXd LP8_matrix
precursor laplacian term-8
Definition msrProblem.H:436

msrProblem::PS6_matrix
Eigen::MatrixXd PS6_matrix
prec_source 6
Definition msrProblem.H:383

msrProblem::Prec8modes
PtrList< volScalarField > Prec8modes
List of pointers used to form the prec1 modes.
Definition msrProblem.H:291

msrProblem::FS4_matrix
List< Eigen::MatrixXd > FS4_matrix
precursor flux source term-4
Definition msrProblem.H:444

msrProblem::Fluxfield
PtrList< volScalarField > Fluxfield
List of pointers used to form the flux snapshots matrix.
Definition msrProblem.H:201

msrProblem::BC3_matrix
Eigen::MatrixXd BC3_matrix
PPE BC3.
Definition msrProblem.H:356

msrProblem::LP4_matrix
Eigen::MatrixXd LP4_matrix
precursor laplacian term-4
Definition msrProblem.H:428

msrProblem::LP7_matrix
Eigen::MatrixXd LP7_matrix
precursor laplacian term-7
Definition msrProblem.H:434

msrProblem::TM_matrix
Eigen::MatrixXd TM_matrix
temperature mass term
Definition msrProblem.H:483

msrProblem::SD3_matrix
List< Eigen::MatrixXd > SD3_matrix
decay heat stream term-3
Definition msrProblem.H:462

msrProblem::Prec2modes
PtrList< volScalarField > Prec2modes
List of pointers used to form the prec1 modes.
Definition msrProblem.H:273

msrProblem::LP6_matrix
Eigen::MatrixXd LP6_matrix
precursor laplacian term-6
Definition msrProblem.H:432

msrProblem::TS_matrix
List< Eigen::MatrixXd > TS_matrix
temperature stream term
Definition msrProblem.H:485

msrProblem::MF_matrix
Eigen::MatrixXd MF_matrix
mass flux
Definition msrProblem.H:364

msrProblem::Prec4field
PtrList< volScalarField > Prec4field
List of pointers used to form the prec4 snapshots matrix.
Definition msrProblem.H:213

msrProblem::MP3_matrix
Eigen::MatrixXd MP3_matrix
precursor mass term-3
Definition msrProblem.H:410

msrProblem::DFields
PtrList< volScalarField > DFields
List of pointers used to form the D snapshosts matrix.
Definition msrProblem.H:246

msrProblem::D_matrix
Eigen::MatrixXd D_matrix
Laplacian term PPE.
Definition msrProblem.H:344

msrProblem::DFS2_matrix
List< Eigen::MatrixXd > DFS2_matrix
decay heat flux source term-2
Definition msrProblem.H:478

msrProblem::change_viscosity
void change_viscosity(double mu)
method to change the viscosity in UEqn
Definition msrProblem.C:1783

msrProblem::laplacian_flux
List< Eigen::MatrixXd > laplacian_flux(label NFluxmodes, label NCmodes)
diffusion eq. methods:
Definition msrProblem.C:938

msrProblem::AFields
PtrList< volScalarField > AFields
List of pointers used to form the A snapshosts matrix.
Definition msrProblem.H:252

msrProblem::PS7_matrix
Eigen::MatrixXd PS7_matrix
prec_source 7
Definition msrProblem.H:385

msrProblem::FS7_matrix
List< Eigen::MatrixXd > FS7_matrix
precursor flux source term-7
Definition msrProblem.H:450

msrProblem::ST6_matrix
List< Eigen::MatrixXd > ST6_matrix
precursor stream term-6
Definition msrProblem.H:400

msrProblem::readMSRfields
void readMSRfields()
Method to read all the fieds of the MSR problem in the offline folder, it also reads mu_samples matri...
Definition msrProblem.C:1665

msrProblem::LP1_matrix
Eigen::MatrixXd LP1_matrix
precursor laplacian term-1
Definition msrProblem.H:422

msrProblem::Dec1modes
PtrList< volScalarField > Dec1modes
List of pointers used to form the dec1 modes.
Definition msrProblem.H:297

msrProblem::THS1_matrix
List< Eigen::MatrixXd > THS1_matrix
temperature decay heat source term-1
Definition msrProblem.H:491

msrProblem::Dec3field
PtrList< volScalarField > Dec3field
List of pointers used to form the dec3 snapshots matrix.
Definition msrProblem.H:237

msrProblem::ST1_matrix
List< Eigen::MatrixXd > ST1_matrix
precursor stream term-1
Definition msrProblem.H:390

msrProblem::liftSolve
void liftSolve()
Perform a lift solve for the velocity field.
Definition msrProblem.C:291

msrProblem::LT_matrix
Eigen::MatrixXd LT_matrix
temperature laplacian term
Definition msrProblem.H:487

msrProblem::Prec6modes
PtrList< volScalarField > Prec6modes
List of pointers used to form the prec1 modes.
Definition msrProblem.H:285

msrProblem::liftSolveT
void liftSolveT()
Perform a lift solve for the temperature.
Definition msrProblem.C:1930

msrProblem::FS6_matrix
List< Eigen::MatrixXd > FS6_matrix
precursor flux source term-6
Definition msrProblem.H:448

msrProblem::mass_temp
Eigen::MatrixXd mass_temp(label NTmodes)
temperature eq. methods
Definition msrProblem.C:1340

msrProblem::vmodes
PtrList< volScalarField > vmodes
List of pointers used to form the v modes.
Definition msrProblem.H:306

msrProblem::P_matrix
Eigen::MatrixXd P_matrix
Div of velocity.
Definition msrProblem.H:341

msrProblem::cumulativeContErr
scalar cumulativeContErr
continuity error
Definition msrProblem.H:179

msrProblem::_fvOptions
autoPtr< fv::options > _fvOptions
fvOptions
Definition msrProblem.H:168

msrProblem::tolerance
scalar tolerance
Tolerance for the residual of the stationary problems, there is the same tolerance for velocity and p...
Definition msrProblem.H:175

msrProblem::Dec1field
PtrList< volScalarField > Dec1field
List of pointers used to form the dec1 snapshots matrix.
Definition msrProblem.H:231

msrProblem::maxIter
scalar maxIter
Number of maximum iterations to be done for the computation of the truth solution.
Definition msrProblem.H:177

msrProblem::TXSmodes
PtrList< volScalarField > TXSmodes
List of pointers used to form the SP snapshosts matrix.
Definition msrProblem.H:321

msrProblem::divergence_term
Eigen::MatrixXd divergence_term(label NUmodes, label NPmodes)
continuity eq. methods:
Definition msrProblem.C:672

msrProblem::LD1_matrix
Eigen::MatrixXd LD1_matrix
decay heat laplacian term-1
Definition msrProblem.H:470

msrProblem::SD2_matrix
List< Eigen::MatrixXd > SD2_matrix
decay heat stream term-2
Definition msrProblem.H:460

msrProblem::_simple
autoPtr< simpleControl > _simple
simpleControl
Definition msrProblem.H:173

msrProblem::MP5_matrix
Eigen::MatrixXd MP5_matrix
precursor mass term-5
Definition msrProblem.H:414

msrProblem::PowerDensfield
PtrList< volScalarField > PowerDensfield
List of pointers used to form the powerDens snapshots matrix.
Definition msrProblem.H:240

msrProblem::Prec3modes
PtrList< volScalarField > Prec3modes
List of pointers used to form the prec1 modes.
Definition msrProblem.H:276

msrProblem::LP2_matrix
Eigen::MatrixXd LP2_matrix
precursor laplacian term-2
Definition msrProblem.H:424

msrProblem::_p0
autoPtr< volScalarField > _p0
Initial fields (for restart purposes).
Definition msrProblem.H:73

msrProblem::diffusive_term
Eigen::MatrixXd diffusive_term(label NUmodes, label NPmodes)
sub-functions needed by projectPPE
Definition msrProblem.C:507

msrProblem::homboolU
bool homboolU
boolean variables to check if the homogenization of U and T is performed (true) or not (false)
Definition msrProblem.H:522

msrProblem::LD2_matrix
Eigen::MatrixXd LD2_matrix
decay heat laplacian term-2
Definition msrProblem.H:472

msrProblem::Prec7field
PtrList< volScalarField > Prec7field
List of pointers used to form the prec7 snapshots matrix.
Definition msrProblem.H:222

msrProblem::PS8_matrix
Eigen::MatrixXd PS8_matrix
prec_source 8
Definition msrProblem.H:387

msrProblem::MP8_matrix
Eigen::MatrixXd MP8_matrix
precursor mass term-8
Definition msrProblem.H:420

msrProblem::MD2_matrix
Eigen::MatrixXd MD2_matrix
decay heat mass term-2
Definition msrProblem.H:466

msrProblem::liftfieldT
PtrList< volScalarField > liftfieldT
List of pointers used to form the list of lifting functions.
Definition msrProblem.H:515

msrProblem::Dec2modes
PtrList< volScalarField > Dec2modes
List of pointers used to form the dec2 modes.
Definition msrProblem.H:300

msrProblem::restart
void restart()
method to set all fields back to values in 0 folder
Definition msrProblem.C:2016

msrProblem::liftfield
PtrList< volVectorField > liftfield
List of pointers used to form the list of lifting functions.
Definition msrProblem.H:513

msrProblem::PS4_matrix
Eigen::MatrixXd PS4_matrix
prec_source 4
Definition msrProblem.H:379

msrProblem::Pfield
PtrList< volScalarField > Pfield
List of pointers used to form the pressure snapshots matrix.
Definition msrProblem.H:195

msrProblem::MP6_matrix
Eigen::MatrixXd MP6_matrix
precursor mass term-6
Definition msrProblem.H:416

msrProblem::NSFFields
PtrList< volScalarField > NSFFields
List of pointers used to form the NSF snapshosts matrix.
Definition msrProblem.H:249

msrProblem::Prec3field
PtrList< volScalarField > Prec3field
List of pointers used to form the prec3 snapshots matrix.
Definition msrProblem.H:210

msrProblem::Umodes
PtrList< volVectorField > Umodes
List of pointers used to form the velocity modes.
Definition msrProblem.H:264

msrProblem::Dec3modes
PtrList< volScalarField > Dec3modes
List of pointers used to form the dec3 modes.
Definition msrProblem.H:303

msrProblem::Tfield
PtrList< volScalarField > Tfield
List of pointers used to form the temperature snapshots matrix.
Definition msrProblem.H:228

msrProblem::B_matrix
Eigen::MatrixXd B_matrix
Diffusion term.
Definition msrProblem.H:329

msrProblem::Tomfield
PtrList< volScalarField > Tomfield
List of pointers used to form the homogeneous velocity snapshots.
Definition msrProblem.H:519

msrProblem::Ufield
PtrList< volVectorField > Ufield
List of pointers used to form the velocity snapshots matrix.
Definition msrProblem.H:198

msrProblem::Prec1field
PtrList< volScalarField > Prec1field
List of pointers used to form the prec1 snapshots matrix.
Definition msrProblem.H:204

msrProblem::M_matrix
Eigen::MatrixXd M_matrix
Mass Matrix.
Definition msrProblem.H:332

msrProblem::MP4_matrix
Eigen::MatrixXd MP4_matrix
precursor mass term-4
Definition msrProblem.H:412

msrProblem::SPmodes
PtrList< volScalarField > SPmodes
List of pointers used to form the SP snapshosts matrix.
Definition msrProblem.H:318

msrProblem::ST2_matrix
List< Eigen::MatrixXd > ST2_matrix
precursor stream term-2
Definition msrProblem.H:392

msrProblem::FS8_matrix
List< Eigen::MatrixXd > FS8_matrix
precursor flux source term-8
Definition msrProblem.H:452

msrProblem::SPFields
PtrList< volScalarField > SPFields
List of pointers used to form the SP snapshosts matrix.
Definition msrProblem.H:255

msrProblem::msrcoeff
void msrcoeff(label &NC)
method to apply RBF interpolation procedure NC is the number of modes to adopt for construncting the ...
Definition msrProblem.C:1795

msrProblem::PS1_matrix
Eigen::MatrixXd PS1_matrix
prec_source 1
Definition msrProblem.H:373

msrProblem::DFS1_matrix
List< Eigen::MatrixXd > DFS1_matrix
decay heat flux source term-1
Definition msrProblem.H:476

msrProblem::DFS3_matrix
List< Eigen::MatrixXd > DFS3_matrix
decay heat flux source term-3
Definition msrProblem.H:480

msrProblem::homogenizeU
void homogenizeU()
Method to compute the homogenized velocity field, it also sets homboolU=true.
Definition msrProblem.C:1653

msrProblem::stream_dec
List< Eigen::MatrixXd > stream_dec(label NUmodes, label NDecmodes, label decgroup)
decay heat eq. methods:
Definition msrProblem.C:1209

msrProblem::FS5_matrix
List< Eigen::MatrixXd > FS5_matrix
precursor flux source term-5
Definition msrProblem.H:446

msrProblem::projectPPE
void projectPPE(fileName folder, label NUmodes, label NPmodes, label NFluxmodes, Eigen::VectorXi Nprecmodes, label NTmodes, Eigen::VectorXi Ndecmodes, label NCmodes)
Project using the Poisson Equation for pressure.
Definition msrProblem.C:395

msrProblem::THS3_matrix
List< Eigen::MatrixXd > THS3_matrix
temperature decay heat source term-3
Definition msrProblem.H:495

msrProblem::Dec2field
PtrList< volScalarField > Dec2field
List of pointers used to form the dec2 snapshots matrix.
Definition msrProblem.H:234

msrProblem::Prec5field
PtrList< volScalarField > Prec5field
List of pointers used to form the prec5 snapshots matrix.
Definition msrProblem.H:216

msrProblem::LD3_matrix
Eigen::MatrixXd LD3_matrix
decay heat laplacian term-3
Definition msrProblem.H:474

msrProblem::precInBool
bool precInBool
boolean variable to decide if apply prec inlet BC
Definition msrProblem.H:182

msrProblem::C_matrix
List< Eigen::MatrixXd > C_matrix
Non linear term.
Definition msrProblem.H:338

msrProblem::TXSFields
PtrList< volScalarField > TXSFields
List of pointers used to form the SP snapshosts matrix.
Definition msrProblem.H:258

msrProblem::LP5_matrix
Eigen::MatrixXd LP5_matrix
precursor laplacian term-5
Definition msrProblem.H:430

msrProblem::msrgetModesEVD
void msrgetModesEVD()
Method to compute the modes for all the fields in the MSR problem, if hombool==false the velocity mod...
Definition msrProblem.C:223

msrProblem::MD3_matrix
Eigen::MatrixXd MD3_matrix
decay heat mass term-3
Definition msrProblem.H:468

msrProblem::Amodes
PtrList< volScalarField > Amodes
List of pointers used to form the A snapshosts matrix.
Definition msrProblem.H:315

msrProblem::PS2_matrix
Eigen::MatrixXd PS2_matrix
prec_source 2
Definition msrProblem.H:375

msrProblem::Fluxmodes
PtrList< volScalarField > Fluxmodes
List of pointers used to form the flux modes.
Definition msrProblem.H:267

msrProblem::Prec2field
PtrList< volScalarField > Prec2field
List of pointers used to form the prec2 snapshots matrix.
Definition msrProblem.H:207

msrProblem::Dmodes
PtrList< volScalarField > Dmodes
List of pointers used to form the D modes.
Definition msrProblem.H:309

msrProblem::msrgetModesSVD
void msrgetModesSVD()
Method to compute the modes for all the fields in the MSR problem, if hombool==false the velocity mod...
Definition msrProblem.C:151

msrProblem::bc_prec
Eigen::MatrixXd bc_prec
matrix to store the values of precursors BC inlet conditions
Definition msrProblem.H:187

msrProblem::ST3_matrix
List< Eigen::MatrixXd > ST3_matrix
precursor stream term-3
Definition msrProblem.H:394

msrProblem::PS5_matrix
Eigen::MatrixXd PS5_matrix
prec_source 5
Definition msrProblem.H:381

reductionProblem::reductionProblem
reductionProblem()
Construct Null.
Definition reductionProblem.C:41

reductionProblem::mu
Eigen::MatrixXd mu
Row matrix of parameters.
Definition reductionProblem.H:89

reductionProblem::truthSolve
void truthSolve()
Perform a TruthSolve.
Definition reductionProblem.C:97

reductionProblem.H
Header file of the reductionProblem class.