src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchInjection/PatchInjection.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2009-2009 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 "PatchInjection.H"
  28 #include "DataEntry.H"
  29 #include "pdf.H"
  30
  31 // * * * * * * * * * * * * Protected Member Functions  * * * * * * * * * * * //
  32
  33 template<class CloudType>
  34 Foam::label Foam::PatchInjection<CloudType>::parcelsToInject
  35 (
  36     const scalar time0,
  37     const scalar time1
  38 ) const
  39 {
  40     if ((time0 >= 0.0) && (time0 < duration_))
  41     {
  42         return round(fraction_*(time1 - time0)*parcelsPerSecond_);
  43     }
  44     else
  45     {
  46         return 0;
  47     }
  48 }
  49
  50
  51 template<class CloudType>
  52 Foam::scalar Foam::PatchInjection<CloudType>::volumeToInject
  53 (
  54     const scalar time0,
  55     const scalar time1
  56 ) const
  57 {
  58     if ((time0 >= 0.0) && (time0 < duration_))
  59     {
  60         return fraction_*volumeFlowRate_().integrate(time0, time1);
  61     }
  62     else
  63     {
  64         return 0.0;
  65     }
  66 }
  67
  68
  69 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  70
  71 template<class CloudType>
  72 Foam::PatchInjection<CloudType>::PatchInjection
  73 (
  74     const dictionary& dict,
  75     CloudType& owner
  76 )
  77 :
  78     InjectionModel<CloudType>(dict, owner, typeName),
  79     patchName_(this->coeffDict().lookup("patchName")),
  80     duration_(readScalar(this->coeffDict().lookup("duration"))),
  81     parcelsPerSecond_
  82     (
  83         readScalar(this->coeffDict().lookup("parcelsPerSecond"))
  84     ),
  85     U0_(this->coeffDict().lookup("U0")),
  86     volumeFlowRate_
  87     (
  88         DataEntry<scalar>::New
  89         (
  90             "volumeFlowRate",
  91             this->coeffDict()
  92         )
  93     ),
  94     parcelPDF_
  95     (
  96         pdf::New
  97         (
  98             this->coeffDict().subDict("parcelPDF"),
  99             owner.rndGen()
 100         )
 101     ),
 102     cellOwners_(),
 103     fraction_(1.0)
 104 {
 105     label patchId = owner.mesh().boundaryMesh().findPatchID(patchName_);
 106
 107     if (patchId < 0)
 108     {
 109         FatalErrorIn
 110         (
 111             "PatchInjection<CloudType>::PatchInjection"
 112             "("
 113                 "const dictionary&, "
 114                 "CloudType&"
 115             ")"
 116         )   << "Requested patch " << patchName_ << " not found" << nl
 117             << "Available patches are: " << owner.mesh().boundaryMesh().names()
 118             << nl << exit(FatalError);
 119     }
 120
 121     const polyPatch& patch = owner.mesh().boundaryMesh()[patchId];
 122
 123     cellOwners_ = patch.faceCells();
 124
 125     label patchSize = cellOwners_.size();
 126     label totalPatchSize = patchSize;
 127     reduce(totalPatchSize, sumOp<scalar>());
 128     fraction_ = patchSize/totalPatchSize;
 129
 130     // Set total volume to inject
 131     this->volumeTotal_ = fraction_*volumeFlowRate_().integrate(0.0, duration_);
 132 }
 133
 134
 135 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 136
 137 template<class CloudType>
 138 Foam::PatchInjection<CloudType>::~PatchInjection()
 139 {}
 140
 141
 142 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 143
 144 template<class CloudType>
 145 bool Foam::PatchInjection<CloudType>::active() const
 146 {
 147     return true;
 148 }
 149
 150
 151 template<class CloudType>
 152 Foam::scalar Foam::PatchInjection<CloudType>::timeEnd() const
 153 {
 154     return this->SOI_ + duration_;
 155 }
 156
 157
 158 template<class CloudType>
 159 void Foam::PatchInjection<CloudType>::setPositionAndCell
 160 (
 161     const label,
 162     const label,
 163     const scalar,
 164     vector& position,
 165     label& cellOwner
 166 )
 167 {
 168     label cellI = this->owner().rndGen().integer(0, cellOwners_.size() - 1);
 169
 170     cellOwner = cellOwners_[cellI];
 171     position = this->owner().mesh().C()[cellOwner];
 172 }
 173
 174
 175 template<class CloudType>
 176 void Foam::PatchInjection<CloudType>::setProperties
 177 (
 178     const label,
 179     const label,
 180     const scalar,
 181     typename CloudType::parcelType& parcel
 182 )
 183 {
 184     // set particle velocity
 185     parcel.U() = U0_;
 186
 187     // set particle diameter
 188     parcel.d() = parcelPDF_->sample();
 189 }
 190
 191
 192 template<class CloudType>
 193 bool Foam::PatchInjection<CloudType>::fullyDescribed() const
 194 {
 195     return false;
 196 }
 197
 198
 199 template<class CloudType>
 200 bool Foam::PatchInjection<CloudType>::validInjection(const label)
 201 {
 202     return true;
 203 }
 204
 205
 206 // ************************************************************************* //