src/lagrangian/dieselSpray/parcel/parcelFunctions.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 "parcel.H"
  27 #include "mathematicalConstants.H"
  28
  29 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  30
  31 namespace Foam
  32 {
  33
  34 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  35
  36 // The diameter based Reynolds number
  37 scalar parcel::Re
  38 (
  39     const vector& U,
  40     const scalar nu
  41 ) const
  42 {
  43     return mag(Urel(U))*d_/nu;
  44 }
  45
  46 // The diameter based Reynolds number
  47 scalar parcel::Re
  48 (
  49     const scalar rho,
  50     const vector& U,
  51     const scalar mu
  52 ) const
  53 {
  54
  55     return rho*mag(Urel(U))*d_/mu;
  56 }
  57
  58 // The diameter based Weber number
  59 scalar parcel::We
  60 (
  61     const vector& U,
  62     const scalar rho,
  63     const scalar sigma
  64 ) const
  65 {
  66     return 0.5*rho*pow(mag(Urel(U)),2)*d_/sigma;
  67 }
  68
  69
  70 scalar parcel::Sc
  71 (
  72     const scalar mu,
  73     const scalar rho,
  74     const scalar massDiffusion
  75 ) const
  76 {
  77     return mu/(rho*massDiffusion);
  78 }
  79
  80
  81 scalar parcel::Sc
  82 (
  83     const scalar nu,
  84     const scalar massDiffusion
  85 ) const
  86 {
  87     return nu/massDiffusion;
  88 }
  89
  90
  91 scalar parcel::Pr
  92 (
  93     const scalar cp,
  94     const scalar mu,
  95     const scalar kappa
  96 ) const
  97 {
  98     return cp*mu/kappa;
  99 }
 100
 101
 102 scalar parcel::N(const scalar rho) const
 103 {
 104     return 6.0*m_/(rho*pow(d_, 3.0)*mathematicalConstant::pi);
 105 }
 106
 107
 108 scalar parcel::Vd() const
 109 {
 110     return pow(d_, 3.0)*mathematicalConstant::pi/6.0;
 111 }
 112
 113
 114 scalar parcel::V(const scalar rho) const
 115 {
 116     return m_/rho;
 117 }
 118
 119
 120 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 121
 122 } // End namespace Foam
 123
 124 // ************************************************************************* //