kOmegaFieldInversionOmega.H

Go to the documentation of this file.

/*---------------------------------------------------------------------------*\
 
    DAFoam  : Discrete Adjoint with OpenFOAM
    Version : v3
 
    Description:
        kOmegaSST turbulence model with a multiplicative betaFieldInversion
        field.
 
\*---------------------------------------------------------------------------*/
 
#ifndef kOmegaFieldInversionOmega_H
#define kOmegaFieldInversionOmega_H
 
#include "RASModel.H"
#include "eddyViscosity.H"
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
namespace Foam
{
namespace RASModels
{
 
/*---------------------------------------------------------------------------*\
                      Class kOmegaFieldInversionOmega Declaration
\*---------------------------------------------------------------------------*/
 
template<class BasicTurbulenceModel>
class kOmegaFieldInversionOmega
    : public eddyViscosity<RASModel<BasicTurbulenceModel>>
{
    // Private Member Functions
 
    // Disallow default bitwise copy construct and assignment
    kOmegaFieldInversionOmega(const kOmegaFieldInversionOmega&);
    void operator=(const kOmegaFieldInversionOmega&);
 
protected:
    // Protected data
 
    // Fields
 
    volScalarField omega_;
 
    volScalarField k_;
 
    volScalarField betaFieldInversion_; 
 
    volScalarField surfaceFriction_;
 
    volScalarField surfaceFrictionData_; 
 
    volScalarField pData_; 
    
    volVectorField UData_; 
 
    volScalarField USingleComponentData_;
 
 
    //- Wall distance
    //  Note: different to wall distance in parent RASModel
    //  which is for near-wall cells only
    const volScalarField& y_;
 
    virtual void correctNut();
 
public:
    typedef typename BasicTurbulenceModel::alphaField alphaField;
    typedef typename BasicTurbulenceModel::rhoField rhoField;
    typedef typename BasicTurbulenceModel::transportModel transportModel;
 
    //- Runtime type information
    TypeName("kOmegaFieldInversionOmega");
 
    // Constructors
 
    //- Construct from components
    kOmegaFieldInversionOmega(
        const alphaField& alpha,
        const rhoField& rho,
        const volVectorField& U,
        const surfaceScalarField& alphaRhoPhi,
        const surfaceScalarField& phi,
        const transportModel& transport,
        const word& propertiesName = turbulenceModel::propertiesName,
        const word& type = typeName);
 
    //- Destructor
    virtual ~kOmegaFieldInversionOmega()
    {
    }
 
    // Member Functions
 
    //- Re-read model coefficients if they have changed
    virtual bool read();
 
    //- Return the turbulence kinetic energy
    virtual tmp<volScalarField> k() const;
 
    //- Return the turbulence kinetic energy dissipation rate
    virtual tmp<volScalarField> epsilon() const;
 
    //- Solve the turbulence equations and correct the turbulence viscosity
    virtual void correct();
};
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
} // End namespace RASModels
} // End namespace Foam
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
#ifdef NoRepository
#include "kOmegaFieldInversionOmega.C"
#endif
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
#endif
 
// ************************************************************************* //