FreeFOAM The Cross-Platform CFD Toolkit
saturateEvaporationModel.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration |
5  \\ / A nd | Copyright (C) 1991-2010 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 <OpenFOAM/error.H>
27 
30 
31 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
32 
33 namespace Foam
34 {
35 
36 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
37 
38 defineTypeNameAndDebug(saturateEvaporationModel, 0);
39 
41 (
42  evaporationModel,
43  saturateEvaporationModel,
44  dictionary
45 );
46 
47 
48 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
49 
50 // Construct from dictionary
52 (
53  const dictionary& dict
54 )
55 :
56  evaporationModel(dict)
57 {}
58 
59 
60 // * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
61 
63 {}
64 
65 
66 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
67 
69 {
70  return true;
71 }
72 
73 // Correlation for the Sherwood Number
75 (
76  const scalar ReynoldsNumber,
77  const scalar SchmidtNumber
78 ) const
79 {
80  return 0.0;
81 }
82 
84 (
85  const scalar diameter,
86  const scalar liquidDensity,
87  const scalar rhoFuelVapor,
88  const scalar massDiffusionCoefficient,
89  const scalar ReynoldsNumber,
90  const scalar SchmidtNumber,
91  const scalar Xs,
92  const scalar Xf,
93  const scalar m0,
94  const scalar dm,
95  const scalar dt
96 ) const
97 {
98  return max(SMALL,dt*(m0/dm - 1.0));
99 }
100 
101 
103 (
104  const scalar liquidDensity,
105  const scalar cpFuel,
106  const scalar heatOfVapour,
107  const scalar kappa,
108  const scalar Nusselt,
109  const scalar deltaTemp,
110  const scalar diameter,
111  const scalar,
112  const scalar,
113  const scalar,
114  const scalar,
115  const scalar,
116  const scalar,
117  const scalar,
118  const scalar,
119  const scalar
120 ) const
121 {
122  scalar time = GREAT;
123 
124  // the deltaTemperature is limited to not go below .5 deg K
125  // for numerical reasons.
126  // This is probably not important, since it results in an upper
127  // limit for the boiling time... which we have anyway.
128  scalar deltaT = max(0.5, deltaTemp);
129 
130  time = liquidDensity*cpFuel*sqr(diameter)/
131  (
132  6.0 * kappa * Nusselt * log(1.0 + cpFuel * deltaT/max(SMALL, heatOfVapour))
133  );
134 
135  time = max(VSMALL, time);
136 
137  return time;
138 }
139 
141 {
142  return 1;
143 }
144 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
145 
146 } // End namespace Foam
147 
148 // ************************ vim: set sw=4 sts=4 et: ************************ //