yPlusWriter.C

Go to the documentation of this file.

/*---------------------------------------------------------------------------*\
    Copyright            : (C) 2011 Symscape
    Website              : www.symscape.com
-------------------------------------------------------------------------------
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 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/>.
 
\*---------------------------------------------------------------------------*/
 
#include "yPlusWriter.H"
#include "dictionary.H"
#include "Time.H"
#include "wallFvPatch.H"
#include "nearWallDist.H"
#include "nutWallFunctionFvPatchScalarField.H"
#include "turbulenceModel.H"
 
 
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
 
namespace Foam
{
  defineTypeNameAndDebug(yPlusWriter, 0);
  const word yPlusWriter::FIELD_NAME("yPlus");
}
 
 
// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
Foam::yPlusWriter::yPlusWriter(const dictionary & dict)
  : derivedFieldWriter(dict)
{}
 
 
// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 
 
void Foam::yPlusWriter::write(const objectRegistry& obr)
{
  const fvMesh& mesh = refCast<const fvMesh>(obr);
  
  volScalarField yPlus
    (
     IOobject
     (
      FIELD_NAME,
      obr.time().timeName(),
      mesh,
      IOobject::NO_READ,
      IOobject::NO_WRITE
      ),
     mesh,
     dimensionedScalar(FIELD_NAME, dimless, 0.0)
     );
 
  if (mesh.foundObject<turbulenceModel>(turbulenceModel::propertiesName))
    {
      const turbulenceModel & tModel = 
    mesh.lookupObject<turbulenceModel>(turbulenceModel::propertiesName);
 
      const volScalarField::GeometricBoundaryField d = nearWallDist(mesh).y();
 
      const volVectorField::GeometricBoundaryField& Ub = 
    tModel.U().boundaryField();
 
      const volScalarField::GeometricBoundaryField nutBf =
    tModel.nut()().boundaryField();
 
      const volScalarField::GeometricBoundaryField nuEffBf =
    tModel.nuEff()().boundaryField();
 
      const volScalarField::GeometricBoundaryField nuBf =
    tModel.nu()().boundaryField();
 
      const fvPatchList& patches = mesh.boundary();
 
      forAll(patches, patchi)
    {
      const fvPatch& patch = patches[patchi];
 
      if (isA<nutWallFunctionFvPatchScalarField>(nutBf[patchi]))
        {
          const nutWallFunctionFvPatchScalarField& nutPf =
        dynamic_cast<const nutWallFunctionFvPatchScalarField&>
        (
         nutBf[patchi]
         );
 
          yPlus.boundaryField()[patchi] = nutPf.yPlus();
          const scalarField& Yp = yPlus.boundaryField()[patchi];
 
          writeStatistics(Yp, nutPf.patch().name(), FIELD_NAME);
        }
      else if (isA<wallFvPatch>(patch))
        {
          yPlus.boundaryField()[patchi] =
        d[patchi]
        *sqrt
        (
         nuEffBf[patchi]
         *mag(Ub[patchi].snGrad())
         )/nuBf[patchi];
          const scalarField& Yp = yPlus.boundaryField()[patchi];
 
          writeStatistics(Yp, patch.name(), FIELD_NAME);
        }
    }
    }
 
  yPlus.write();
}
 
 
// ************************************************************************* //