src/turbulenceModels/RAS/incompressible/derivedFvPatchFields/turbulentMixingLengthFrequencyInlet/turbulentMixingLengthFrequencyInletFvPatchScalarField.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2006-2008 OpenCFD Ltd.
   6      \\/     M anipulation  |
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of OpenFOAM.
  10
  11     OpenFOAM is free software; you can redistribute it and/or modify it
  12     under the terms of the GNU General Public License as published by the
  13     Free Software Foundation; either version 2 of the License, or (at your
  14     option) any later version.
  15
  16     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
  17     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
  18     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
  19     for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with OpenFOAM; if not, write to the Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  24
  25 \*---------------------------------------------------------------------------*/
  26
  27 #include "turbulentMixingLengthFrequencyInletFvPatchScalarField.H"
  28 #include "addToRunTimeSelectionTable.H"
  29 #include "fvPatchFieldMapper.H"
  30 #include "surfaceFields.H"
  31 #include "volFields.H"
  32 #include "RASModel.H"
  33
  34 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  35
  36 namespace Foam
  37 {
  38
  39 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  40
  41 turbulentMixingLengthFrequencyInletFvPatchScalarField::
  42 turbulentMixingLengthFrequencyInletFvPatchScalarField
  43 (
  44     const fvPatch& p,
  45     const DimensionedField<scalar, volMesh>& iF
  46 )
  47 :
  48     fixedValueFvPatchField<scalar>(p, iF),
  49     mixingLength_(0.0),
  50     kName_("undefined-k")
  51 {}
  52
  53 turbulentMixingLengthFrequencyInletFvPatchScalarField::
  54 turbulentMixingLengthFrequencyInletFvPatchScalarField
  55 (
  56     const turbulentMixingLengthFrequencyInletFvPatchScalarField& ptf,
  57     const fvPatch& p,
  58     const DimensionedField<scalar, volMesh>& iF,
  59     const fvPatchFieldMapper& mapper
  60 )
  61 :
  62     fixedValueFvPatchField<scalar>(ptf, p, iF, mapper),
  63     mixingLength_(ptf.mixingLength_),
  64     kName_(ptf.kName_)
  65 {}
  66
  67 turbulentMixingLengthFrequencyInletFvPatchScalarField::
  68 turbulentMixingLengthFrequencyInletFvPatchScalarField
  69 (
  70     const fvPatch& p,
  71     const DimensionedField<scalar, volMesh>& iF,
  72     const dictionary& dict
  73 )
  74 :
  75     fixedValueFvPatchField<scalar>(p, iF, dict),
  76     mixingLength_(readScalar(dict.lookup("mixingLength"))),
  77     kName_(dict.lookup("k"))
  78 {}
  79
  80 turbulentMixingLengthFrequencyInletFvPatchScalarField::
  81 turbulentMixingLengthFrequencyInletFvPatchScalarField
  82 (
  83     const turbulentMixingLengthFrequencyInletFvPatchScalarField& ptf
  84 )
  85 :
  86     fixedValueFvPatchField<scalar>(ptf),
  87     mixingLength_(ptf.mixingLength_),
  88     kName_(ptf.kName_)
  89 {}
  90
  91 turbulentMixingLengthFrequencyInletFvPatchScalarField::
  92 turbulentMixingLengthFrequencyInletFvPatchScalarField
  93 (
  94     const turbulentMixingLengthFrequencyInletFvPatchScalarField& ptf,
  95     const DimensionedField<scalar, volMesh>& iF
  96 )
  97 :
  98     fixedValueFvPatchField<scalar>(ptf, iF),
  99     mixingLength_(ptf.mixingLength_),
 100     kName_(ptf.kName_)
 101 {}
 102
 103
 104 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 105
 106 void turbulentMixingLengthFrequencyInletFvPatchScalarField::updateCoeffs()
 107 {
 108     if (updated())
 109     {
 110         return;
 111     }
 112
 113     // Lookup Cmu corresponding to the turbulence model selected
 114     const incompressible::RASModel& RAS =
 115         db().lookupObject<incompressible::RASModel>("RASProperties");
 116     scalar Cmu = readScalar(RAS.coeffDict().lookup("Cmu"));
 117
 118     scalar Cmu25 = pow(Cmu, 0.25);
 119
 120     const fvPatchField<scalar>& kp =
 121         patch().lookupPatchField<volScalarField, scalar>(kName_);
 122
 123     operator==(sqrt(kp)/(Cmu25*mixingLength_));
 124
 125     fixedValueFvPatchField<scalar>::updateCoeffs();
 126 }
 127
 128
 129 void turbulentMixingLengthFrequencyInletFvPatchScalarField::write
 130 (
 131     Ostream& os
 132 ) const
 133 {
 134     fvPatchField<scalar>::write(os);
 135     os.writeKeyword("mixingLength")
 136         << mixingLength_ << token::END_STATEMENT << nl;
 137     os.writeKeyword("k") << kName_ << token::END_STATEMENT << nl;
 138     writeEntry("value", os);
 139 }
 140
 141
 142 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 143
 144 makePatchTypeField
 145 (
 146     fvPatchScalarField,
 147     turbulentMixingLengthFrequencyInletFvPatchScalarField
 148 );
 149
 150
 151 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 152
 153 } // End namespace Foam
 154
 155 // ************************************************************************* //