/*---------------------------------------------------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | \\ / A nd | Copyright (C) 2008-2009 OpenCFD Ltd. \\/ M anipulation | ------------------------------------------------------------------------------- License This file is part of OpenFOAM. 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 2 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, write to the Free Software Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA Class Foam::incompressible::RASModels::nutkRoughWallFunctionFvPatchScalarField Description Boundary condition for turbulent (kinematic) viscosity when using wall functions for rough walls, based on turbulence kinetic energy. Manipulates the E parameter to account for roughness effects, based on KsPlus. - roughness height = sand-grain roughness (0 for smooth walls) - roughness constant = 0.5-1.0 (0.5 default) SourceFiles nutkRoughWallFunctionFvPatchScalarField.C \*---------------------------------------------------------------------------*/ #ifndef nutkRoughWallFunctionFvPatchScalarField_H #define nutkRoughWallFunctionFvPatchScalarField_H #include "nutkWallFunctionFvPatchScalarField.H" // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // namespace Foam { namespace incompressible { namespace RASModels { /*---------------------------------------------------------------------------*\ Class nutkRoughWallFunctionFvPatchScalarField Declaration \*---------------------------------------------------------------------------*/ class nutkRoughWallFunctionFvPatchScalarField : public nutkWallFunctionFvPatchScalarField { protected: // Protected data //- Roughness height scalarField Ks_; //- Roughness constant scalarField Cs_; // Protected member functions //- Compute the roughness function virtual scalar fnRough(const scalar KsPlus, const scalar Cs) const; //- Calculate the turbulence viscosity virtual tmp calcNut() const; public: //- Runtime type information TypeName("nutkRoughWallFunction"); // Constructors //- Construct from patch and internal field nutkRoughWallFunctionFvPatchScalarField ( const fvPatch&, const DimensionedField& ); //- Construct from patch, internal field and dictionary nutkRoughWallFunctionFvPatchScalarField ( const fvPatch&, const DimensionedField&, const dictionary& ); //- Construct by mapping given // nutkRoughWallFunctionFvPatchScalarField // onto a new patch nutkRoughWallFunctionFvPatchScalarField ( const nutkRoughWallFunctionFvPatchScalarField&, const fvPatch&, const DimensionedField&, const fvPatchFieldMapper& ); //- Construct as copy nutkRoughWallFunctionFvPatchScalarField ( const nutkRoughWallFunctionFvPatchScalarField& ); //- Construct and return a clone virtual tmp clone() const { return tmp ( new nutkRoughWallFunctionFvPatchScalarField(*this) ); } //- Construct as copy setting internal field reference nutkRoughWallFunctionFvPatchScalarField ( const nutkRoughWallFunctionFvPatchScalarField&, const DimensionedField& ); //- Construct and return a clone setting internal field reference virtual tmp clone ( const DimensionedField& iF ) const { return tmp ( new nutkRoughWallFunctionFvPatchScalarField(*this, iF) ); } // Member functions // Mapping functions //- Map (and resize as needed) from self given a mapping object virtual void autoMap(const fvPatchFieldMapper&); //- Reverse map the given fvPatchField onto this fvPatchField virtual void rmap ( const fvPatchScalarField&, const labelList& ); // I-O //- Write virtual void write(Ostream&) const; }; // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // } // End namespace RASModels } // End namespace incompressible } // End namespace Foam // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // #endif // ************************************************************************* //