CompressibleUnSteadyRhoPimple.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
    CompressibleUnSteadyRhoPimple
Description
    full-order class for a Compressible UnSteady  problem with Pimple algorithm.
SourceFiles
    CompressibleUnSteadyRhoPimple.C
\*---------------------------------------------------------------------------*/

 
#ifndef CompressibleUnSteadyRhoPimple_H
#define CompressibleUnSteadyRhoPimple_H
#include "fvCFD.H"
#include "dynamicFvMesh.H"
#include "fluidThermo.H"
#include "turbulentFluidThermoModel.H"
#include "CompressibleSteadyNS.H"
#include "CorrectPhi.H"
#include "UnsteadyProblem.H"
#include "pimpleControl.H"
#include "pressureControl.H"
#include "fvOptions.H"
#include "reductionProblem.H"
#include "ITHACAstream.H"
#include "ITHACAparameters.H"
#include "ITHACAforces.H"
#include <iostream>
 
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
/*---------------------------------------------------------------------------*\
                        Class CompressibleUnSteadyRhoPimple Declaration
\*---------------------------------------------------------------------------*/


class CompressibleUnSteadyRhoPimple: public CompressibleSteadyNS,
    public UnsteadyProblem
{
    public:
        // Constructors
        CompressibleUnSteadyRhoPimple();

        CompressibleUnSteadyRhoPimple(int argc, char* argv[]);
        ~CompressibleUnSteadyRhoPimple() {};
 
        // Functions
 
        //--------------------------------------------------------------------------
        //void truthSolve(List<scalar> mu_now);
        void truthSolve(word folder);
 
        // Momuntum equation terms
        fvVectorMatrix getUmatrix(volVectorField& U);//, Vector<double>& uresidual_v);
        fvVectorMatrix getNLTerm(volVectorField& U);
        fvVectorMatrix getViscTerm(volVectorField& U);
        volVectorField getGradP(volScalarField& p);
 
        // Energy equation terms
        fvScalarMatrix getEmatrix(volVectorField& U,
                                  volScalarField& p);//, scalar& eresidual);
        fvScalarMatrix getFluxTerm();
        volScalarField getKinEnTerm(volVectorField& U, volScalarField& p);
        fvScalarMatrix getDiffTerm();
 
        // Pressure equation terms
        fvScalarMatrix getPmatrix(fvVectorMatrix& Ueqn, volVectorField& U,
                                  volScalarField& p);
        surfaceScalarField getPhiHbyA(fvVectorMatrix& Ueqn, volVectorField& U,
                                      volScalarField& p);
        volScalarField getDivPhiHbyA(fvVectorMatrix& Ueqn, volVectorField& U,
                                     volScalarField& p);
        surfaceScalarField getRhorAUf(fvVectorMatrix& Ueqn);
        fvScalarMatrix getPoissonTerm(fvVectorMatrix& Ueqn, volScalarField& p);
 
        //--------------------------------------------------------------------------
        void changeViscosity(double mu_new);
 
        //--------------------------------------------------------------------------
        void restart();
 
        autoPtr<fluidThermo> pThermo;
        // /// Dynamic mesh field
        autoPtr<Foam::dynamicFvMesh> meshPtr;
        autoPtr<pimpleControl> _pimple;
 
        autoPtr<volScalarField> _rho;
 
        autoPtr<pressureControl> _pressureControl;
 
        autoPtr<compressible::turbulenceModel> turbulence;
 
        autoPtr<dimensionedScalar> _initialMass;
 
        autoPtr<volScalarField> _psi;
 
        autoPtr<volScalarField> _E;
        //autoPtr<surfaceVectorField> _rhoUf;

        autoPtr<volScalarField> _p0;

        autoPtr<volScalarField> _E0;

        autoPtr<volVectorField> _U0;

        autoPtr<volScalarField> _rho0;

        autoPtr<surfaceScalarField> _phi0, _phi;
        PtrList<volScalarField> Efield;

        volScalarModes Emodes;
 
        autoPtr<fvVectorMatrix> Ueqn_global;
 
        autoPtr<fvScalarMatrix> Peqn_global;
 
        autoPtr<fvScalarMatrix> Eeqn_global;
 
        autoPtr<surfaceScalarField> phiHbyA;
 
        autoPtr<surfaceScalarField> rhorAUf;
 
        autoPtr<volVectorField> HbyA;

        label middleStep;
        //label middleStepInit;

        bool middleExport;
 
        bool  correctPhi;
        bool  checkMeshCourantNo;
        bool  moveMeshOuterCorrectors;
        //const dimensionedScalar rhoMax;
        //const dimensionedScalar rhoMin;
 
 
 
};
 
#endif