src/regionModels/thermoBaffleModels/derivedFvPatchFields/temperatureThermoBaffle/temperatureThermoBaffleFvPatchScalarField.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
   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
  13     the Free Software Foundation, either version 3 of the License, or
  14     (at your 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, see <http://www.gnu.org/licenses/>.
  23
  24 \*---------------------------------------------------------------------------*/
  25
  26 #include "temperatureThermoBaffleFvPatchScalarField.H"
  27 #include "addToRunTimeSelectionTable.H"
  28
  29 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  30
  31 namespace Foam
  32 {
  33 namespace compressible
  34 {
  35
  36
  37 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  38
  39 temperatureThermoBaffleFvPatchScalarField::
  40 temperatureThermoBaffleFvPatchScalarField
  41 (
  42     const fvPatch& p,
  43     const DimensionedField<scalar, volMesh>& iF
  44 )
  45 :
  46     turbulentTemperatureCoupledBaffleMixedFvPatchScalarField(p, iF),
  47     owner_(false),
  48     baffle_(),
  49     solidThermoType_("undefined")
  50 {}
  51
  52
  53 temperatureThermoBaffleFvPatchScalarField::
  54 temperatureThermoBaffleFvPatchScalarField
  55 (
  56     const temperatureThermoBaffleFvPatchScalarField& ptf,
  57     const fvPatch& p,
  58     const DimensionedField<scalar, volMesh>& iF,
  59     const fvPatchFieldMapper& mapper
  60 )
  61 :
  62     turbulentTemperatureCoupledBaffleMixedFvPatchScalarField
  63     (
  64         ptf,
  65         p,
  66         iF,
  67         mapper
  68     ),
  69     owner_(ptf.owner_),
  70     baffle_(ptf.baffle_),
  71     solidThermoType_(ptf.solidThermoType_)
  72 {}
  73
  74
  75 temperatureThermoBaffleFvPatchScalarField::
  76 temperatureThermoBaffleFvPatchScalarField
  77 (
  78     const fvPatch& p,
  79     const DimensionedField<scalar, volMesh>& iF,
  80     const dictionary& dict
  81 )
  82 :
  83     turbulentTemperatureCoupledBaffleMixedFvPatchScalarField(p, iF, dict),
  84     owner_(false),
  85     baffle_(),
  86     solidThermoType_()
  87 {
  88     if (!isA<directMappedPatchBase>(patch().patch()))
  89     {
  90         FatalErrorIn
  91         (
  92             "temperatureThermoBaffleFvPatchScalarField::"
  93             "temperatureThermoBaffleFvPatchScalarField\n"
  94             "(\n"
  95             "    const fvPatch& p,\n"
  96             "    const DimensionedField<scalar, volMesh>& iF,\n"
  97             "    const dictionary& dict\n"
  98             ")\n"
  99         )   << "\n    patch type '" << patch().type()
 100             << "' not type '" << directMappedPatchBase::typeName << "'"
 101             << "\n    for patch " << patch().name()
 102             << " of field " << dimensionedInternalField().name()
 103             << " in file " << dimensionedInternalField().objectPath()
 104             << exit(FatalError);
 105     }
 106
 107     const fvMesh& thisMesh = patch().boundaryMesh().mesh();
 108
 109     typedef regionModels::thermoBaffleModels::thermoBaffleModel baffle;
 110
 111     if
 112     (
 113         thisMesh.name() == polyMesh::defaultRegion
 114      && !thisMesh.foundObject<baffle>("thermoBaffle")
 115      && !owner_
 116     )
 117     {
 118         Info << "Creating thermal baffle..." << endl;
 119         baffle_.reset(baffle::New(thisMesh, dict).ptr());
 120         owner_ = true;
 121         dict.lookup("thermoType") >> solidThermoType_;
 122     }
 123 }
 124
 125
 126 temperatureThermoBaffleFvPatchScalarField::
 127 temperatureThermoBaffleFvPatchScalarField
 128 (
 129     const temperatureThermoBaffleFvPatchScalarField& ptf,
 130     const DimensionedField<scalar, volMesh>& iF
 131 )
 132 :
 133     turbulentTemperatureCoupledBaffleMixedFvPatchScalarField(ptf, iF),
 134     owner_(ptf.owner_),
 135     baffle_(ptf.baffle_),
 136     solidThermoType_(ptf.solidThermoType_)
 137 {}
 138
 139 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 140
 141
 142 void temperatureThermoBaffleFvPatchScalarField::autoMap
 143 (
 144     const fvPatchFieldMapper& m
 145 )
 146 {
 147     mixedFvPatchScalarField::autoMap(m);
 148 }
 149
 150
 151 void temperatureThermoBaffleFvPatchScalarField::rmap
 152 (
 153     const fvPatchScalarField& ptf,
 154     const labelList& addr
 155 )
 156 {
 157     mixedFvPatchScalarField::rmap(ptf, addr);
 158 }
 159
 160
 161 void temperatureThermoBaffleFvPatchScalarField::updateCoeffs()
 162 {
 163     if (this->updated())
 164     {
 165         return;
 166     }
 167
 168     const fvMesh& thisMesh = patch().boundaryMesh().mesh();
 169
 170     if
 171     (
 172         thisMesh.name() == polyMesh::defaultRegion
 173      && owner_
 174     )
 175     {
 176         baffle_->evolve();
 177     }
 178
 179     turbulentTemperatureCoupledBaffleMixedFvPatchScalarField::updateCoeffs();
 180 }
 181
 182
 183 void temperatureThermoBaffleFvPatchScalarField::write(Ostream& os) const
 184 {
 185     turbulentTemperatureCoupledBaffleMixedFvPatchScalarField::write(os);
 186
 187     const fvMesh& thisMesh = patch().boundaryMesh().mesh();
 188
 189     if (thisMesh.name() == polyMesh::defaultRegion && owner_)
 190     {
 191         os.writeKeyword("thermoBaffleModel") <<  baffle_->modelName()
 192             << token::END_STATEMENT << nl;
 193
 194         os.writeKeyword("regionName") <<  baffle_->regionMesh().name()
 195             << token::END_STATEMENT << nl;
 196
 197         os.writeKeyword("infoOutput") <<  baffle_->infoOutput()
 198             << token::END_STATEMENT << nl;
 199
 200         os.writeKeyword("active") <<  baffle_->active()
 201             << token::END_STATEMENT << nl;
 202
 203         os.writeKeyword(word(baffle_->modelName() + "coeffs"));
 204
 205         os << baffle_->coeffs() << nl;
 206
 207         os.writeKeyword("thermoType") << solidThermoType_
 208             << token::END_STATEMENT << nl;
 209
 210         os.writeKeyword(word(solidThermoType_ + "Coeffs")) << nl;
 211
 212         os << indent << '{' << nl
 213            << indent << baffle_->thermo() << nl
 214            << indent << '}' << nl;
 215     }
 216 }
 217
 218 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 219
 220 makePatchTypeField
 221 (
 222     fvPatchScalarField,
 223     temperatureThermoBaffleFvPatchScalarField
 224 );
 225
 226
 227 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 228
 229 } // End namespace compressible
 230 } // End namespace Foam
 231
 232 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 233
 234
 235 // ************************************************************************* //