Zuber.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) 2018-2021 OpenCFD Ltd
9-------------------------------------------------------------------------------
10License
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 "Zuber.H"
31#include "phasePairKey.H"
32#include "phaseSystem.H"
33
34// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
35
36namespace Foam
37{
38namespace wallBoilingModels
39{
40namespace CHFModels
41{
44 (
46 Zuber,
48 );
49}
50}
51}
52
53
54// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
55
57(
58 const dictionary& dict
59)
60:
61 CHFModel(),
62 Cn_(dict.getOrDefault<scalar>("Cn", 0.131))
63{}
64
65
66// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
67
70(
71 const phaseModel& liquid,
72 const phaseModel& vapor,
73 const label patchi,
74 const scalarField& Tl,
75 const scalarField& Tsatw,
76 const scalarField& L
77) const
78{
79 const auto& g =
80 liquid.mesh().time().lookupObject<uniformDimensionedVectorField>("g");
81
82 const labelUList& cells = liquid.mesh().boundary()[patchi].faceCells();
83
84 const scalarField& pw = liquid.thermo().p().boundaryField()[patchi];
85
86 tmp<scalarField> trhoVapor = vapor.thermo().rhoEoS(pw ,Tsatw, cells);
87 const scalarField& rhoVapor = trhoVapor.ref();
88
89 tmp<scalarField> trhoLiq = liquid.thermo().rhoEoS(pw, Tsatw, cells);
90 const scalarField& rhoLiq = trhoLiq.ref();
91 const phasePairKey pair(liquid.name(), vapor.name());
92 const scalarField sigma
93 (
94 liquid.fluid().sigma(pair)().boundaryField()[patchi]
95 );
96
97 return
98 Cn_*rhoVapor*L*pow
99 (
100 sigma*mag(g.value())*(max(rhoLiq - rhoVapor, Zero)/sqr(rhoVapor)),
101 0.25
102 );
103}
104
105
107{
109 os.writeEntry("Cn", Cn_);
110}
111
112
113// ************************************************************************* //
Macros for easy insertion into run-time selection tables.
#define addToRunTimeSelectionTable(baseType, thisType, argNames)
Add to construction table with typeName as the key.
const uniformDimensionedVectorField & g
An Ostream is an abstract base class for all output systems (streams, files, token lists,...
Definition: Ostream.H:62
virtual tmp< scalarField > rhoEoS(const scalarField &p, const scalarField &T, const labelList &cells) const =0
Density from pressure and temperature from EoS.
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
Definition: dictionary.H:126
virtual bool write()
Write the output fields.
Generic thermophysical properties class for a liquid in which the functions and coefficients for each...
Definition: liquid.H:57
scalar sigma(scalar p, scalar T) const
Surface tension [N/m].
Definition: liquidI.H:94
Single incompressible phase derived from the phase-fraction. Used as part of the multiPhaseMixture fo...
Definition: phaseModel.H:61
const word & name() const
Definition: phaseModel.H:144
virtual const rhoThermo & thermo() const =0
Return the thermophysical model.
An ordered or unorder pair of phase names. Typically specified as follows.
Definition: phasePairKey.H:68
A class for managing temporary objects.
Definition: tmp.H:65
T & ref() const
Definition: tmpI.H:227
Base class for critical heat flux (CHF) correlation models for boiling flows.
Definition: CHFModel.H:58
A critical heat flux (CHF) correlation model based on Zuber (1958) for boiling flows.
Definition: Zuber.H:104
virtual tmp< scalarField > CHF(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: Zuber.C:70
#define defineTypeNameAndDebug(Type, DebugSwitch)
Define the typeName and debug information.
Definition: className.H:121
OBJstream os(runTime.globalPath()/outputName)
const cellShapeList & cells
Namespace for OpenFOAM.
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
Definition: hashSets.C:47
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
static constexpr const zero Zero
Global zero (0)
Definition: zero.H:131
dictionary dict
const vector L(dict.get< vector >("L"))