Skip to content
GitLab
Projects
Groups
Snippets
/
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in / Register
Toggle navigation
Menu
Open sidebar
Development
openfoam
Commits
09b91bbd
Commit
09b91bbd
authored
Jul 30, 2019
by
sergio
Committed by
Andrew Heather
Jul 30, 2019
Browse files
ENH: Small editing to alphatWallBoilingWallFunction
parent
186edcfd
Changes
2
Hide whitespace changes
Inline
Side-by-side
src/phaseSystemModels/reactingEulerFoam/derivedFvPatchFields/alphatWallBoilingWallFunction/alphatWallBoilingWallFunctionFvPatchScalarField.C
View file @
09b91bbd
...
...
@@ -1045,7 +1045,7 @@ void alphatWallBoilingWallFunctionFvPatchScalarField::updateCoeffs()
mDotL_
[
i
]
=
dmdt_
[
i
]
*
L
[
i
];
// No quenching flux
//
qq_[i] = 0.0;
qq_
[
i
]
=
0
.
0
;
this
->
operator
[](
i
)
=
(
...
...
@@ -1104,6 +1104,7 @@ void alphatWallBoilingWallFunctionFvPatchScalarField::updateCoeffs()
A1
[
i
]
=
1
.
0
;
qq_
[
i
]
=
0
.
0
;
mDotL_
[
i
]
=
0
.
0
;
dmdt_
[
i
]
=
0
.
0
;
// Turbulente thermal diffusivity for single phase.
this
->
operator
[](
i
)
=
...
...
@@ -1141,12 +1142,12 @@ void alphatWallBoilingWallFunctionFvPatchScalarField::updateCoeffs()
Info
<<
" dmdt: "
<<
gMin
((
dmdt_
))
<<
" - "
<<
gMax
((
dmdt_
))
<<
endl
;
Info
<<
" alphatlEff: "
<<
gMin
(
liquidw
*
(
*
this
+
alphaw
))
Info
<<
" alphatlEff: "
<<
gMin
(
liquidw
*
(
*
this
+
alphaw
))
<<
" - "
<<
gMax
(
liquidw
*
(
*
this
+
alphaw
))
<<
endl
;
scalar
Qeff
=
gSum
(
qEff
*
patch
().
magSf
());
Info
<<
" Effective heat transfer rate to liquid:"
<<
Qeff
<<
endl
;
Info
<<
" Effective heat transfer rate to liquid:
"
<<
Qeff
<<
endl
<<
nl
;
if
(
debug
&
2
)
{
...
...
@@ -1186,11 +1187,13 @@ void alphatWallBoilingWallFunctionFvPatchScalarField::updateCoeffs()
}
}
Info
<<
"Sub Cool faces : "
<<
nSubCool
<<
endl
;
Info
<<
"Transient faces : "
<<
nTransient
<<
endl
;
Info
<<
"Film faces : "
<<
nFilm
<<
endl
;
Info
<<
"Non Boiling faces : "
<<
nNonBoiling
<<
endl
;
Info
<<
"Total faces : "
<<
this
->
size
()
<<
endl
;
Info
<<
"Faces regime : "
<<
nl
<<
endl
;
Info
<<
" sub Cool faces : "
<<
nSubCool
<<
endl
;
Info
<<
" transient faces : "
<<
nTransient
<<
endl
;
Info
<<
" film faces : "
<<
nFilm
<<
endl
;
Info
<<
" non-Boiling faces : "
<<
nNonBoiling
<<
endl
;
Info
<<
" total faces : "
<<
this
->
size
()
<<
endl
<<
nl
;
const
scalarField
qc
(
...
...
@@ -1199,7 +1202,7 @@ void alphatWallBoilingWallFunctionFvPatchScalarField::updateCoeffs()
);
scalar
Qc
=
gSum
(
qc
*
patch
().
magSf
());
Info
<<
" Convective heat transfer:"
<<
Qc
<<
endl
;
Info
<<
" Convective heat transfer:
"
<<
Qc
<<
endl
;
const
scalarField
qFilm
(
...
...
@@ -1219,14 +1222,21 @@ void alphatWallBoilingWallFunctionFvPatchScalarField::updateCoeffs()
Info
<<
" Transient boiling heat transfer:"
<<
Qtbtot
<<
endl
;
Info
<<
"
t
DNB: "
<<
gMin
(
tDNB
)
<<
" - "
<<
gMax
(
tDNB
)
Info
<<
"
T
DNB: "
<<
gMin
(
tDNB
)
<<
" - "
<<
gMax
(
tDNB
)
<<
endl
;
scalar
QsubCool
=
gSum
const
scalarField
qSubCool
(
fLiquid
*
nSubCools
*
(
qq_
+
qe
())
*
patch
().
magSf
()
fLiquid
*
nSubCools
*
(
A1
*
alphatConv_
*
hew
.
snGrad
()
+
qe
()
+
qq
()
)
);
Info
<<
" Sub Cool boiling heat transfer:"
<<
QsubCool
scalar
QsubCool
=
gSum
(
qSubCool
*
patch
().
magSf
());
Info
<<
" Sub Cool boiling heat transfer: "
<<
QsubCool
<<
endl
;
Info
<<
" N: "
<<
gMin
(
nSubCools
*
N
)
<<
" - "
...
...
tutorials/heatTransfer/chtMultiRegionTwoPhaseEulerFoam/solidQuenching2D/system/water/fvSchemes
View file @
09b91bbd
...
...
@@ -22,58 +22,27 @@ ddtSchemes
gradSchemes
{
default Gauss linear;
/*
grad((1-alpha.gas)) leastSquares;//Gauss linear;
grad(alpha.gas) leastSquares;//Gauss linear;
grad(U.gas) leastSquares;//Gauss linear;
grad(U.liquid) leastSquares;// Gauss linear;
grad(h.gas) leastSquares;
grad(h.liquid) leastSquares;
grad(alpha.liquid) leastSquares;
grad(alpha.gas) leastSquares;
grad(rho) leastSquares;
grad(p_rgh) leastSquares;
grad(epsilon.liquid) leastSquares;
grad(k.liquid) leastSquares;
*/
}
divSchemes
{
default none;
//Gauss upwind;
default none;
"div\(phi,alpha.*\)" Gauss vanLeer;
"div\(phir,alpha.*\)" Gauss vanLeer;
"div\(alphaRhoPhi.*,U.*\)" Gauss upwind;
//limitedLinearV 1;
"div\(phi.*,U.*\)" Gauss upwind;
//limitedLinearV 1;
"div\(alphaRhoPhi.*,U.*\)" Gauss upwind;
"div\(phi.*,U.*\)" Gauss upwind;
"div\(alphaRhoPhi.*,Yi\)" Gauss upwind;
//limitedLinear 1;
"div\(alphaRhoPhi.*,(h|e|f).*\)" Gauss upwind;
//limitedLinear 1;
"div\(alphaRhoPhi.*,K.*\)" Gauss upwind;
//limitedLinear 1;
"div\(alphaPhi.*,p\)" Gauss upwind;
//limitedLinear 1;
"div\(alphaRhoPhi.*,Yi\)" Gauss upwind;
"div\(alphaRhoPhi.*,(h|e|f).*\)" Gauss upwind;
"div\(alphaRhoPhi.*,K.*\)" Gauss upwind;
"div\(alphaPhi.*,p\)" Gauss upwind;
"div\(alphaRhoPhi.*,(k|epsilon).*\)" Gauss upwind;
"div\(phim,(k|epsilon)m\)" Gauss upwind;
"div\(\(\(\(alpha.*\*thermo:rho.*\)\*nuEff.*\)\*dev2\(T\(grad\(U.*\)\)\)\)\)" Gauss linear;
/*
div(phi,U) Gauss upwind;
div(phi,K) Gauss linear;
div(phi,h) Gauss upwind;
div(phi,k) Gauss upwind;
div(phi,epsilon) Gauss upwind;
div(phi,R) Gauss upwind;
div(R) Gauss linear;
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
*/
}
laplacianSchemes
...
...
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment