src/lagrangian/distributionModels/RosinRammler/RosinRammler.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 "RosinRammler.H"
  27 #include "addToRunTimeSelectionTable.H"
  28
  29 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  30
  31 namespace Foam
  32 {
  33     namespace distributionModels
  34     {
  35         defineTypeNameAndDebug(RosinRammler, 0);
  36         addToRunTimeSelectionTable(distributionModel, RosinRammler, dictionary);
  37     }
  38 }
  39
  40 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  41
  42 Foam::distributionModels::RosinRammler::RosinRammler
  43 (
  44     const dictionary& dict,
  45     cachedRandom& rndGen
  46 )
  47 :
  48     distributionModel(typeName, dict, rndGen),
  49     minValue_(readScalar(distributionModelDict_.lookup("minValue"))),
  50     maxValue_(readScalar(distributionModelDict_.lookup("maxValue"))),
  51     d_(readScalar(distributionModelDict_.lookup("d"))),
  52     n_(readScalar(distributionModelDict_.lookup("n")))
  53 {
  54     check();
  55 }
  56
  57
  58 Foam::distributionModels::RosinRammler::RosinRammler(const RosinRammler& p)
  59 :
  60     distributionModel(p),
  61     minValue_(p.minValue_),
  62     maxValue_(p.maxValue_),
  63     d_(p.d_),
  64     n_(p.n_)
  65 {}
  66
  67
  68 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
  69
  70 Foam::distributionModels::RosinRammler::~RosinRammler()
  71 {}
  72
  73
  74 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  75
  76 Foam::scalar Foam::distributionModels::RosinRammler::sample() const
  77 {
  78     scalar K = 1.0 - exp(-pow((maxValue_ - minValue_)/d_, n_));
  79     scalar y = rndGen_.sample01<scalar>();
  80     scalar x = minValue_ + d_*::pow(-log(1.0 - y*K), 1.0/n_);
  81     return x;
  82 }
  83
  84
  85 Foam::scalar Foam::distributionModels::RosinRammler::minValue() const
  86 {
  87     return minValue_;
  88 }
  89
  90
  91 Foam::scalar Foam::distributionModels::RosinRammler::maxValue() const
  92 {
  93     return maxValue_;
  94 }
  95
  96
  97 // ************************************************************************* //