zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

ENH: Small editing to fvSchemes and alphatWallBoilingWallFunction

Browse Source
This commit is contained in:
sergio
2019-07-30 08:49:15 -07:00
committed by Andrew Heather
parent c07ab24556
commit 962e6f6318
2 changed files with 31 additions and 52 deletions
Download Patch File Download Diff File Expand all files Collapse all files

36
src/phaseSystemModels/reactingEulerFoam/derivedFvPatchFields/alphatWallBoilingWallFunction/alphatWallBoilingWallFunctionFvPatchScalarField.C
View File

@ -1045,7 +1045,7 @@ void alphatWallBoilingWallFunctionFvPatchScalarField::updateCoeffs()
mDotL_[i] = dmdt_[i]*L[i]; mDotL_[i] = dmdt_[i]*L[i];
// No quenching flux // No quenching flux
//qq_[i] = 0.0; qq_[i] = 0.0;
this->operator[](i) = this->operator[](i) =
( (
@ -1104,6 +1104,7 @@ void alphatWallBoilingWallFunctionFvPatchScalarField::updateCoeffs()
A1[i] = 1.0; A1[i] = 1.0;
qq_[i] = 0.0; qq_[i] = 0.0;
mDotL_[i] = 0.0; mDotL_[i] = 0.0;
dmdt_[i] = 0.0;
// Turbulente thermal diffusivity for single phase. // Turbulente thermal diffusivity for single phase.
this->operator[](i) = this->operator[](i) =
@ -1145,8 +1146,8 @@ void alphatWallBoilingWallFunctionFvPatchScalarField::updateCoeffs()
<< " - " << gMax(liquidw*(*this + alphaw)) << endl; << " - " << gMax(liquidw*(*this + alphaw)) << endl;
scalar Qeff = gSum(qEff*patch().magSf()); scalar Qeff = gSum(qEff*patch().magSf());
Info<< " Effective heat transfer rate to liquid:" << Qeff Info<< " Effective heat transfer rate to liquid: " << Qeff
<< endl; << endl << nl;
if (debug & 2) if (debug & 2)
{ {
@ -1186,11 +1187,13 @@ void alphatWallBoilingWallFunctionFvPatchScalarField::updateCoeffs()
} }
} }
Info<< "Sub Cool faces : " << nSubCool << endl; Info<< "Faces regime : " << nl << endl;
Info<< "Transient faces : " << nTransient << endl;
Info<< "Film faces : " << nFilm << endl; Info<< " sub Cool faces : " << nSubCool << endl;
Info<< "Non Boiling faces : " << nNonBoiling << endl; Info<< " transient faces : " << nTransient << endl;
Info<< "Total faces : " << this->size() << endl; Info<< " film faces : " << nFilm << endl;
Info<< " non-Boiling faces : " << nNonBoiling << endl;
Info<< " total faces : " << this->size() << endl << nl;
const scalarField qc const scalarField qc
( (
@ -1199,7 +1202,7 @@ void alphatWallBoilingWallFunctionFvPatchScalarField::updateCoeffs()
); );
scalar Qc = gSum(qc*patch().magSf()); scalar Qc = gSum(qc*patch().magSf());
Info<< " Convective heat transfer:" << Qc << endl; Info<< " Convective heat transfer: " << Qc << endl;
const scalarField qFilm const scalarField qFilm
( (
@ -1219,14 +1222,21 @@ void alphatWallBoilingWallFunctionFvPatchScalarField::updateCoeffs()
Info<< " Transient boiling heat transfer:" << Qtbtot Info<< " Transient boiling heat transfer:" << Qtbtot
<< endl; << endl;
Info<< " tDNB: " << gMin(tDNB) << " - " << gMax(tDNB) Info<< " TDNB: " << gMin(tDNB) << " - " << gMax(tDNB)
<< endl; << 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; << endl;
Info<< " N: " << gMin(nSubCools*N) << " - " Info<< " N: " << gMin(nSubCools*N) << " - "

45
tutorials/heatTransfer/chtMultiRegionTwoPhaseEulerFoam/solidQuenching2D/system/water/fvSchemes
View File

@ -22,58 +22,27 @@ ddtSchemes
gradSchemes gradSchemes
{ {
default Gauss linear; 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 divSchemes
{ {
default none;//Gauss upwind; default none;
"div\(phi,alpha.*\)" Gauss vanLeer; "div\(phi,alpha.*\)" Gauss vanLeer;
"div\(phir,alpha.*\)" Gauss vanLeer; "div\(phir,alpha.*\)" Gauss vanLeer;
"div\(alphaRhoPhi.*,U.*\)" Gauss upwind;//limitedLinearV 1; "div\(alphaRhoPhi.*,U.*\)" Gauss upwind;
"div\(phi.*,U.*\)" Gauss upwind;//limitedLinearV 1; "div\(phi.*,U.*\)" Gauss upwind;
"div\(alphaRhoPhi.*,Yi\)" Gauss upwind;//limitedLinear 1; "div\(alphaRhoPhi.*,Yi\)" Gauss upwind;
"div\(alphaRhoPhi.*,(h|e|f).*\)" Gauss upwind;//limitedLinear 1; "div\(alphaRhoPhi.*,(h|e|f).*\)" Gauss upwind;
"div\(alphaRhoPhi.*,K.*\)" Gauss upwind;//limitedLinear 1; "div\(alphaRhoPhi.*,K.*\)" Gauss upwind;
"div\(alphaPhi.*,p\)" Gauss upwind;//limitedLinear 1; "div\(alphaPhi.*,p\)" Gauss upwind;
"div\(alphaRhoPhi.*,(k|epsilon).*\)" Gauss upwind; "div\(alphaRhoPhi.*,(k|epsilon).*\)" Gauss upwind;
"div\(phim,(k|epsilon)m\)" Gauss upwind; "div\(phim,(k|epsilon)m\)" Gauss upwind;
"div\(\(\(\(alpha.*\*thermo:rho.*\)\*nuEff.*\)\*dev2\(T\(grad\(U.*\)\)\)\)\)" Gauss linear; "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 laplacianSchemes