HuaXu.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 "HuaXu.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 HuaXu,
48 );
49}
50}
51}
52
53
54// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
55
57(
58 const dictionary& dict
59)
60:
62 Kburn_(dict.getOrDefault<scalar>("Kburn", 1.5))
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 scalarField alphaLiq(liquid.alpha(patchi));
83
84 const labelUList& cells = liquid.mesh().boundary()[patchi].faceCells();
85
86 const scalarField& pw = liquid.thermo().p().boundaryField()[patchi];
87
88 tmp<scalarField> trhoVapor = vapor.thermo().rhoEoS(pw, Tsatw, cells);
89 const scalarField& rhoVapor = trhoVapor.ref();
90
91 tmp<scalarField> trhoLiq = liquid.thermo().rhoEoS(pw, Tsatw, cells);
92 const scalarField& rhoLiq = trhoLiq.ref();
93
94
95 tmp<volScalarField> tCp = liquid.thermo().Cp();
96 const volScalarField& Cp = tCp();
97 const fvPatchScalarField& Cpw = Cp.boundaryField()[patchi];
98
99 const phasePairKey pair(liquid.name(), vapor.name());
100 const scalarField sigma
101 (
102 liquid.fluid().sigma(pair)().boundaryField()[patchi]
103 );
104
105 const scalarField Pe
106 (
107 pow(sigma, 0.75)
108 /
109 (
110 alphaLiq
111 * pow025(mag(g.value())*(rhoLiq-rhoVapor))
112 * sqrt(rhoVapor)
113 )
114 );
115
116 const scalarField Ja
117 (
118 rhoLiq*Cpw*max(Tsatw - Tl, scalar(0))/(rhoVapor*L)
119 );
120
121 return Kburn_*(scalar(1) + 0.345*Ja/pow025(Pe));
122}
123
124
126(
127 Ostream& os
128) const
129{
131 os.writeEntry("Kburn", Kburn_);
132}
133
134
135// ************************************************************************* //
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
scalar Cp(scalar p, scalar T) const
Liquid heat capacity [J/(kg K)].
Definition: liquidI.H:46
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
A critical heat flux (CHF) sub-cooling correlation model based on Hua-Xu (2000) for boiling flows.
Definition: HuaXu.H:105
virtual tmp< scalarField > CHFSubCool(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: HuaXu.C:70
#define defineTypeNameAndDebug(Type, DebugSwitch)
Define the typeName and debug information.
Definition: className.H:121
const tmp< volScalarField > & tCp
Definition: EEqn.H:4
const volScalarField & Cp
Definition: EEqn.H:7
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
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensionedScalar sqrt(const dimensionedScalar &ds)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
dimensionedScalar pow025(const dimensionedScalar &ds)
dictionary dict
const vector L(dict.get< vector >("L"))