DASpalartAllmarasFv3.H

Go to the documentation of this file.

/*---------------------------------------------------------------------------*\
 
    DAFoam  : Discrete Adjoint with OpenFOAM
    Version : v4
 
    Description:
    Child class for the SpalartAllmarasFv3 model used in OpenFOAM 2.4 and before
 
    This file is modified from OpenFOAM's source code
    src/TurbulenceModels/turbulenceModels/RAS/SpalartAllmaras/SpalartAllmaras.H
 
    OpenFOAM: The Open Source CFD Toolbox
 
    Copyright (C): 2011-2016 OpenFOAM Foundation
 
    OpenFOAM License:
 
        OpenFOAM is free software: you can redistribute it and/or modify it
        under the terms of the GNU General Public License as published by
        the Free Software Foundation, either version 3 of the License, or
        (at your option) any later version.
    
        OpenFOAM 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 General Public License
        for more details.
    
        You should have received a copy of the GNU General Public License
        along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
 
\*---------------------------------------------------------------------------*/
 
#ifndef DASpalartAllmarasFv3_H
#define DASpalartAllmarasFv3_H
 
#include "DATurbulenceModel.H"
#include "addToRunTimeSelectionTable.H"
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
namespace Foam
{
 
/*---------------------------------------------------------------------------*\
                       Class DASpalartAllmarasFv3 Declaration
\*---------------------------------------------------------------------------*/
 
class DASpalartAllmarasFv3
    : public DATurbulenceModel
{
 
protected:
 
    dimensionedScalar sigmaNut_;
    dimensionedScalar kappa_;
    dimensionedScalar Cb1_;
    dimensionedScalar Cb2_;
    dimensionedScalar Cw1_;
    dimensionedScalar Cw2_;
    dimensionedScalar Cw3_;
    dimensionedScalar Cv1_;
    dimensionedScalar Cv2_;
 
 
    tmp<volScalarField> chi() const;
 
    tmp<volScalarField> fv1(const volScalarField& chi) const;
 
    tmp<volScalarField> fv2(
        const volScalarField& chi,
        const volScalarField& fv1) const;
 
    tmp<volScalarField> fv3(
        const volScalarField& chi,
        const volScalarField& fv1) const;
 
    tmp<volScalarField> fw(const volScalarField& Stilda) const;
 
 
 
    volScalarField& nuTilda_;
    volScalarField nuTildaRes_;
 
    volScalarField pseudoNuTilda_;
    fvScalarMatrix pseudoNuTildaEqn_;
 
    const volScalarField& y_;
 
    volScalarField betaFINuTilda_;
 
    label solveTurbState_ = 0;
 
    // fvSolutions parameters
    scalar relaxNuTildaEqn_ = 1.0;
    dictionary solverDictNuTilda_;
 
public:
    TypeName("SpalartAllmarasFv3");
    // Constructors
 
    //- Construct from components
    DASpalartAllmarasFv3(
        const word modelType,
        const fvMesh& mesh,
        const DAOption& daOption);
 
    //- Destructor
    virtual ~DASpalartAllmarasFv3()
    {
    }
 
    // Member functions
 
    tmp<volScalarField> DnuTildaEff() const;
 
    virtual void correctModelStates(wordList& modelStates) const;
 
    virtual void correctNut();
 
    virtual void correctBoundaryConditions();
 
    virtual void updateIntermediateVariables();
 
    virtual void correctStateResidualModelCon(List<List<word>>& stateCon) const;
 
    virtual void addModelResidualCon(HashTable<List<List<word>>>& allCon) const;
 
    virtual void calcResiduals(const dictionary& options);
 
    virtual void correct(label printToScreen);
 
    virtual void invTranProdNuTildaEqn(
        const volScalarField& mySource,
        volScalarField& pseudoNuTilda);
 
    virtual void constructPseudoNuTildaEqn();
 
    virtual void rhsSolvePseudoNuTildaEqn(const volScalarField& nuTildaSource);
 
    virtual void calcLduResidualTurb(volScalarField& nuTildaRes);
 
    virtual void getFvMatrixFields(
        const word varName,
        scalarField& diag,
        scalarField& upper,
        scalarField& lower);
 
    virtual void getTurbProdOverDestruct(volScalarField& PoD) const;
 
    virtual void getTurbConvOverProd(volScalarField& CoP) const;
};
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
} // End namespace Foam
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
#endif
 
// ************************************************************************* //