Kutadeladze.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 | www.openfoam.com
6  \\/ M anipulation |
7 -------------------------------------------------------------------------------
8  Copyright (C) 2021 OpenCFD Ltd
9 -------------------------------------------------------------------------------
10 License
11  This file is part of OpenFOAM.
12 
13  OpenFOAM is free software: you can redistribute it and/or modify it
14  under the terms of the GNU General Public License as published by
15  the Free Software Foundation, either version 3 of the License, or
16  (at your option) any later version.
17 
18  OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
19  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
20  FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
21  for more details.
22 
23  You should have received a copy of the GNU General Public License
24  along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
25 
26 \*---------------------------------------------------------------------------*/
27 
28 #include "Kutadeladze.H"
31 #include "phasePairKey.H"
32 #include "phaseSystem.H"
33 
34 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
35 
36 namespace Foam
37 {
38 namespace wallBoilingModels
39 {
40 namespace nucleateFluxModels
41 {
42  defineTypeNameAndDebug(Kutadeladze, 0);
44  (
45  nucleateFluxModel,
46  Kutadeladze,
47  dictionary
48  );
49 }
50 }
51 }
52 
53 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
54 
55 Foam::wallBoilingModels::nucleateFluxModels::Kutadeladze::Kutadeladze
56 (
57  const dictionary& dict
58 )
59 :
61  Cn_(dict.getOrDefault<scalar>("Cn", 5.66e-10)),
62  an_(dict.getOrDefault<scalar>("an", 2.5)),
63  bn_(dict.getOrDefault<scalar>("bn", 1)),
64  n_(dict.getOrDefault<scalar>("n", 1))
65 {}
66 
67 
68 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
69 
72 (
73  const phaseModel& liquid,
74  const phaseModel& vapor,
75  const label patchi,
76  const scalarField& Tl,
77  const scalarField& Tsatw,
78  const scalarField& L
79 ) const
80 {
81  const auto& p = liquid.mesh().lookupObject<volScalarField>("p");
82  const scalarField& pb = p.boundaryField()[patchi];
83 
84  const labelUList& cells = liquid.mesh().boundary()[patchi].faceCells();
85 
86  tmp<scalarField> trhoVapor = vapor.thermo().rhoEoS(pb, Tsatw, cells);
87  const scalarField& rhoVapor = trhoVapor.ref();
88 
89  tmp<scalarField> trhoLiq = liquid.thermo().rhoEoS(pb, Tsatw, cells);
90  const scalarField& rhoLiq = trhoLiq.ref();
91 
92  const phasePairKey pair(liquid.name(), vapor.name());
93  const scalarField sigma
94  (
95  liquid.fluid().sigma(pair)().boundaryField()[patchi]
96  );
97 
98  const fvPatchScalarField& Tw =
99  liquid.thermo().T().boundaryField()[patchi];
100 
101  const scalarField kappaLiquid(liquid.kappa(patchi));
102 
103  tmp<scalarField> tCpliq = liquid.thermo().Cp(pb, Tsatw, cells);
104  const scalarField& Cpliquid = tCpliq();
105 
106  const scalarField muLiquid(liquid.mu(patchi));
107 
108  const scalarField deltaTsub
109  (
110  pow(max((Tw-Tsatw), scalar(0)), an_)
111  );
112 
113  return
114  Cn_*kappaLiquid*pow(Cpliquid, 1.5)*pow(rhoLiq,1.28)*pow(pb,1.75)
115  *deltaTsub
116  /
117  (pow(muLiquid, 0.625)*pow(sigma,0.9)*pow(L,1.5)*pow(rhoVapor,1.5));
118 }
119 
120 
122 (
123  Ostream& os
124 ) const
125 {
127  os.writeEntry("Cn", Cn_);
128  os.writeEntry("an", an_);
129  os.writeEntry("bn", bn_);
130  os.writeEntry("n", n_);
131 }
132 
133 
134 // ************************************************************************* //
Foam::fvPatchField< scalar >
Foam::phaseModel
Single incompressible phase derived from the phase-fraction. Used as part of the multiPhaseMixture fo...
Definition: phaseModel.H:54
L
const vector L(dict.get< vector >("L"))
p
volScalarField & p
Definition: createFieldRefs.H:8
Foam::wallBoilingModels::nucleateFluxModels::Kutadeladze::qNucleate
virtual tmp< scalarField > qNucleate(const phaseModel &liquid, const phaseModel &vapor, const label patchi, const scalarField &Tl, const scalarField &Tsatw, const scalarField &L) const
Calculate and return the nucleation-site density.
Definition: Kutadeladze.C:72
Foam::liquid::Cp
scalar Cp(scalar p, scalar T) const
Liquid heat capacity [J/(kg K)].
Definition: liquidI.H:46
Foam::wallBoilingModels::nucleateFluxModel
Base class for nucleation flux models.
Definition: nucleateFluxModel.H:54
Foam::tmp
A class for managing temporary objects.
Definition: PtrList.H:61
Foam::wallBoilingModels::nucleateFluxModels::defineTypeNameAndDebug
defineTypeNameAndDebug(exponential, 0)
Foam::liquid
Generic thermophysical properties class for a liquid in which the functions and coefficients for each...
Definition: liquid.H:54
Foam::phaseModel::thermo
virtual const rhoThermo & thermo() const =0
Access const to phase thermo.
Foam::tmp::ref
T & ref() const
Definition: tmpI.H:227
Foam::Field< scalar >
Foam::phasePairKey
An ordered or unorder pair of phase names. Typically specified as follows.
Definition: phasePairKey.H:65
Foam::pow
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
Definition: dimensionedScalar.C:75
Foam::max
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
Definition: hashSets.C:47
dict
dictionary dict
Definition: searchingEngine.H:14
Kutadeladze.H
Foam::liquid::kappa
scalar kappa(scalar p, scalar T) const
Liquid thermal conductivity [W/(m K)].
Definition: liquidI.H:82
Foam::dictionary
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
Definition: dictionary.H:123
os
OBJstream os(runTime.globalPath()/outputName)
addToRunTimeSelectionTable.H
Macros for easy insertion into run-time selection tables.
Foam
Namespace for OpenFOAM.
Definition: atmBoundaryLayer.C:33
Foam::wallBoilingModels::nucleateFluxModels::Kutadeladze::write
virtual void write(Ostream &os) const
Write.
Definition: Kutadeladze.C:122
Foam::phaseModel::name
const word & name() const
Definition: phaseModel.H:140
uniformDimensionedFields.H
Foam::UList< label >
Foam::wallBoilingModels::nucleateFluxModels::addToRunTimeSelectionTable
addToRunTimeSelectionTable(nucleateFluxModel, exponential, dictionary)
Foam::liquid::sigma
scalar sigma(scalar p, scalar T) const
Surface tension [N/m].
Definition: liquidI.H:94
Foam::vtk::write
void write(vtk::formatter &fmt, const Type &val, const label n=1)
Component-wise write of a value (N times)
Definition: foamVtkOutputTemplates.C:36
cells
const cellShapeList & cells
Definition: gmvOutputHeader.H:3
Foam::Ostream
An Ostream is an abstract base class for all output systems (streams, files, token lists,...
Definition: Ostream.H:56
sigma
dimensionedScalar sigma("sigma", dimMass/sqr(dimTime), transportProperties)
Foam::GeometricField< scalar, fvPatchField, volMesh >
Foam::basicThermo::rhoEoS
virtual tmp< scalarField > rhoEoS(const scalarField &p, const scalarField &T, const labelList &cells) const =0
Density from pressure and temperature from EoS.
Foam::liquid::mu
scalar mu(scalar p, scalar T) const
Liquid viscosity [Pa s].
Definition: liquidI.H:70