src/lagrangian/basic/particle/particle.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2004-2011 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
  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 "particle.H"
  27 #include "transform.H"
  28
  29 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  30
  31 Foam::label Foam::particle::particleCount_ = 0;
  32
  33 const Foam::scalar Foam::particle::trackingCorrectionTol = 1e-5;
  34
  35 const Foam::scalar Foam::particle::lambdaDistanceToleranceCoeff = 1e3*SMALL;
  36
  37 namespace Foam
  38 {
  39     defineTypeNameAndDebug(particle, 0);
  40 }
  41
  42
  43 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  44
  45 Foam::particle::particle
  46 (
  47     const polyMesh& mesh,
  48     const vector& position,
  49     const label cellI,
  50     const label tetFaceI,
  51     const label tetPtI
  52 )
  53 :
  54     mesh_(mesh),
  55     position_(position),
  56     cellI_(cellI),
  57     faceI_(-1),
  58     stepFraction_(0.0),
  59     tetFaceI_(tetFaceI),
  60     tetPtI_(tetPtI),
  61     origProc_(Pstream::myProcNo()),
  62     origId_(getNewParticleID())
  63 {}
  64
  65
  66 Foam::particle::particle
  67 (
  68     const polyMesh& mesh,
  69     const vector& position,
  70     const label cellI,
  71     bool doCellFacePt
  72 )
  73 :
  74     mesh_(mesh),
  75     position_(position),
  76     cellI_(cellI),
  77     faceI_(-1),
  78     stepFraction_(0.0),
  79     tetFaceI_(-1),
  80     tetPtI_(-1),
  81     origProc_(Pstream::myProcNo()),
  82     origId_(getNewParticleID())
  83 {
  84     if (doCellFacePt)
  85     {
  86         initCellFacePt();
  87     }
  88 }
  89
  90
  91 Foam::particle::particle(const particle& p)
  92 :
  93     mesh_(p.mesh_),
  94     position_(p.position_),
  95     cellI_(p.cellI_),
  96     faceI_(p.faceI_),
  97     stepFraction_(p.stepFraction_),
  98     tetFaceI_(p.tetFaceI_),
  99     tetPtI_(p.tetPtI_),
 100     origProc_(p.origProc_),
 101     origId_(p.origId_)
 102 {}
 103
 104
 105 Foam::particle::particle(const particle& p, const polyMesh& mesh)
 106 :
 107     mesh_(mesh),
 108     position_(p.position_),
 109     cellI_(p.cellI_),
 110     faceI_(p.faceI_),
 111     stepFraction_(p.stepFraction_),
 112     tetFaceI_(p.tetFaceI_),
 113     tetPtI_(p.tetPtI_),
 114     origProc_(p.origProc_),
 115     origId_(p.origId_)
 116 {}
 117
 118
 119 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 120
 121 void Foam::particle::transformProperties(const tensor&)
 122 {}
 123
 124
 125 void Foam::particle::transformProperties(const vector&)
 126 {}
 127
 128
 129 Foam::scalar Foam::particle::wallImpactDistance(const vector&) const
 130 {
 131     Info<< "particle::wallImpactDistance" << endl;
 132
 133     return 0.0;
 134 }
 135
 136
 137 // * * * * * * * * * * * * * * Friend Operators * * * * * * * * * * * * * * //
 138
 139 bool Foam::operator==(const particle& pA, const particle& pB)
 140 {
 141     return (pA.origProc() == pB.origProc() && pA.origId() == pB.origId());
 142 }
 143
 144
 145 bool Foam::operator!=(const particle& pA, const particle& pB)
 146 {
 147     return !(pA == pB);
 148 }
 149
 150
 151 // ************************************************************************* //