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

BUG: Fixes issue #1606

Browse Source 
Proper counting of regime on faces of the patch
This commit is contained in:
sergio
2020-02-28 11:56:06 -08:00
parent 7e8a215506
commit 18e53c3c06
1 changed files with 9 additions and 11 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@ -6,7 +6,7 @@
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2015-2018 OpenFOAM Foundation
Copyright (C) 2018 OpenCFD Ltd
Copyright (C) 2018-2020 OpenCFD Ltd
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@ -1005,11 +1005,6 @@ void alphatWallBoilingWallFunctionFvPatchScalarField::updateCoeffs()
if (debug & 2)
{
scalar nSubCool(0);
scalar nTransient(0);
scalar nFilm(0);
scalar nNonBoiling(0);
scalarField nSubCools(this->size(), 0);
scalarField nTransients(this->size(), 0);
scalarField nFilms(this->size(), 0);
@ -1021,26 +1016,27 @@ void alphatWallBoilingWallFunctionFvPatchScalarField::updateCoeffs()
switch (regimeTypes[i])
{
case regimeType::subcool:
nSubCool++;
nSubCools[i] = 1;
break;
case regimeType::transient:
nTransient++;
nTransients[i] = 1;
break;
case regimeType::film:
nFilm++;
nFilms[i] = 1;
break;
case regimeType::nonBoiling:
nNonBoiling++;
nNonBoilings[i] = 1;
break;
}
}
scalar nSubCool(gSum(nSubCools));
scalar nTransient(gSum(nTransients));
scalar nFilm(gSum(nFilms));
scalar nNonBoiling(gSum(nNonBoilings));
Info<< "Faces regime : " << nl << endl;
@ -1048,7 +1044,9 @@ void alphatWallBoilingWallFunctionFvPatchScalarField::updateCoeffs()
Info<< " transient faces : " << nTransient << endl;
Info<< " film faces : " << nFilm << endl;
Info<< " non-Boiling faces : " << nNonBoiling << endl;
Info<< " total faces : " << this->size() << endl << nl;
Info<< " total faces : "
<< nSubCool + nTransient + nFilm + nNonBoiling
<< endl << nl;
const scalarField qc
(