src/lagrangian/solidParticle/solidParticleCloud.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | foam-extend: Open Source CFD
   4    \\    /   O peration     | Version:     3.2
   5     \\  /    A nd           | Web:         http://www.foam-extend.org
   6      \\/     M anipulation  | For copyright notice see file Copyright
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of foam-extend.
  10
  11     foam-extend 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 3 of the License, or (at your
  14     option) any later version.
  15
  16     foam-extend is distributed in the hope that it will be useful, but
  17     WITHOUT ANY WARRANTY; without even the implied warranty of
  18     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  19     General Public License for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 \*---------------------------------------------------------------------------*/
  25
  26 #include "solidParticleCloud.H"
  27 #include "fvMesh.H"
  28 #include "volFields.H"
  29 #include "interpolationCellPoint.H"
  30
  31 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  32
  33 namespace Foam
  34 {
  35     defineParticleTypeNameAndDebug(solidParticle, 0);
  36     defineTemplateTypeNameAndDebug(Cloud<solidParticle>, 0);
  37 };
  38
  39 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  40
  41 Foam::solidParticleCloud::solidParticleCloud
  42 (
  43     const fvMesh& mesh,
  44     const word& cloudName,
  45     bool readFields
  46 )
  47 :
  48     Cloud<solidParticle>(mesh, cloudName, false),
  49     mesh_(mesh),
  50     particleProperties_
  51     (
  52         IOobject
  53         (
  54             "particleProperties",
  55             mesh_.time().constant(),
  56             mesh_,
  57             IOobject::MUST_READ,
  58             IOobject::NO_WRITE
  59         )
  60     ),
  61     rhop_(dimensionedScalar(particleProperties_.lookup("rhop")).value()),
  62     e_(dimensionedScalar(particleProperties_.lookup("e")).value()),
  63     mu_(dimensionedScalar(particleProperties_.lookup("mu")).value())
  64 {
  65     if (readFields)
  66     {
  67         solidParticle::readFields(*this);
  68     }
  69 }
  70
  71
  72 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  73
  74 void Foam::solidParticleCloud::move(const dimensionedVector& g)
  75 {
  76     const volScalarField& rho = mesh_.lookupObject<const volScalarField>("rho");
  77     const volVectorField& U = mesh_.lookupObject<const volVectorField>("U");
  78     const volScalarField& nu = mesh_.lookupObject<const volScalarField>("nu");
  79
  80     interpolationCellPoint<scalar> rhoInterp(rho);
  81     interpolationCellPoint<vector> UInterp(U);
  82     interpolationCellPoint<scalar> nuInterp(nu);
  83
  84     solidParticle::trackData td(*this, rhoInterp, UInterp, nuInterp, g.value());
  85
  86     Cloud<solidParticle>::move(td);
  87 }
  88
  89
  90 // ************************************************************************* //