src/lagrangian/dieselSpray/parcel/parcelFunctions.C

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