Kunz.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration | Website: https://openfoam.org
5  \\ / A nd | Copyright (C) 2011-2023 OpenFOAM Foundation
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 "Kunz.H"
28 
29 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
30 
31 namespace Foam
32 {
33 namespace compressible
34 {
35 namespace cavitationModels
36 {
39 }
40 }
41 }
42 
43 
44 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
45 
47 (
48  const dictionary& dict,
49  const compressibleTwoPhases& phases
50 )
51 :
52  cavitationModel(dict, phases),
53 
54  UInf_("UInf", dimVelocity, dict),
55  tInf_("tInf", dimTime, dict),
56  Cv_("Cv", dimless, dict),
57  Cc_("Cc", dimless, dict),
58 
59  p0_("0", dimPressure, 0)
60 {
61  correct();
62 }
63 
64 
65 // * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * * //
66 
68 Foam::compressible::cavitationModels::Kunz::mvCoeff() const
69 {
70  return Cv_*rhov()/(0.5*rhol()*sqr(UInf_)*tInf_);
71 }
72 
73 
75 Foam::compressible::cavitationModels::Kunz::mcCoeff() const
76 {
77  return Cc_*rhov()/tInf_;
78 }
79 
80 
83 {
85  phases_.mesh().lookupObject<volScalarField>("p");
86 
87  const volScalarField::Internal alphal
88  (
89  min(max(this->alphal(), scalar(0)), scalar(1))
90  );
91 
92  const volScalarField::Internal pSatv(this->pSatv());
93  const volScalarField::Internal pSatl(this->pSatl());
94 
96  (
97  mcCoeff()*sqr(alphal)
98  *max(p - pSatv, p0_)/max(p - pSatv, 0.01*pSatv),
99  - mvCoeff()*min(p - pSatl, p0_)
100  );
101 }
102 
103 
106 {
107  const volScalarField::Internal& p =
108  phases_.mesh().lookupObject<volScalarField>("p");
109 
110  const volScalarField::Internal alphav
111  (
112  min(max(this->alphav(), scalar(0)), scalar(1))
113  );
114 
115  const volScalarField::Internal alphal
116  (
117  min(max(this->alphal(), scalar(0)), scalar(1))
118  );
119 
120  const volScalarField::Internal pSatv(this->pSatv());
121  const volScalarField::Internal pSatl(this->pSatl());
122 
124  (
125  mcCoeff()*alphav*sqr(alphal)*pos0(p - pSatv)
126  /max(p - pSatv, 0.01*pSatv),
127  - mvCoeff()*alphal*neg(p - pSatl)
128  );
129 }
130 
131 
133 {}
134 
135 
137 (
138  const dictionary& dict
139 )
140 {
142  {
143  dict.lookup("UInf") >> UInf_;
144  dict.lookup("tInf") >> tInf_;
145  dict.lookup("Cv") >> Cv_;
146  dict.lookup("Cc") >> Cc_;
147 
148  return true;
149  }
150  else
151  {
152  return false;
153  }
154 }
155 
156 
157 // ************************************************************************* //
Macros for easy insertion into run-time selection tables.
Field with dimensions and associated with geometry type GeoMesh which is used to size the field and a...
Generic GeometricField class.
An ordered pair of two objects of type <T> with first() and second() elements.
Definition: Pair.H:65
virtual bool read(const dictionary &dict)=0
Read the dictionary and update.
virtual Pair< tmp< volScalarField::Internal > > mDotcvAlphal() const
Return the mass condensation and vaporisation rates as a.
Definition: Kunz.C:82
virtual void correct()
Correct the Kunz phaseChange model.
Definition: Kunz.C:132
virtual bool read(const dictionary &dict)
Read the dictionary and update.
Definition: Kunz.C:137
Kunz(const dictionary &dict, const compressibleTwoPhases &phases)
Construct for phases.
Definition: Kunz.C:47
virtual Pair< tmp< volScalarField::Internal > > mDotcvP() const
Return the mass condensation and vaporisation rates as coefficients.
Definition: Kunz.C:105
A list of keyword definitions, which are a keyword followed by any number of values (e....
Definition: dictionary.H:162
A class for managing temporary objects.
Definition: tmp.H:55
addToRunTimeSelectionTable(cavitationModel, Kunz, dictionary)
Namespace for OpenFOAM.
dimensionedScalar pos0(const dimensionedScalar &ds)
const dimensionSet dimPressure
dimensionedSymmTensor sqr(const dimensionedVector &dv)
const dimensionSet dimless
layerAndWeight min(const layerAndWeight &a, const layerAndWeight &b)
const dimensionSet dimTime
layerAndWeight max(const layerAndWeight &a, const layerAndWeight &b)
dimensionedScalar neg(const dimensionedScalar &ds)
const dimensionSet dimVelocity
dictionary dict
volScalarField & p