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: PDRsetFields adjustments

Browse Source 
- missing minThreshold in activePressureForceBaffleVelocity

- update names of cyclics

- increase tolerance for zero-thickness baffles

- velocity field "U" instead of "UBet"

- viscosity field "nut" instead of "mut"
This commit is contained in:
Jonathan Cranford
2020-04-06 15:27:03 +02:00
committed by Mark Olesen
parent b6a953cfc9
commit 8cf346a9c5
2 changed files with 12 additions and 11 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
applications/utilities/preProcessing/PDRsetFields/PDRarraysAnalyse.C
View File

@ -673,15 +673,15 @@ void Foam::PDRarrays::addBlockage
p.patchType = obs.blowoff_type;
p.blowoffPress = obs.blowoff_press;
p.blowoffTime = obs.blowoff_time;
if (obs.span.x() < 1e-5)
if (obs.span.x() < 0.01)
{
indir = 1;
}
else if (obs.span.y() < 1e-5)
else if (obs.span.y() < 0.01)
{
indir = 2;
}
else if (obs.span.z() < 1e-5)
else if (obs.span.z() < 0.01)
{
indir = 3;
}

17
applications/utilities/preProcessing/PDRsetFields/PDRarraysCalc.C
View File

@ -934,9 +934,9 @@ void Foam::PDRarrays::calculateAndWrite
);
write_uniformField
(
"mut", 0, MUT_WALL_FN,
"nut", 0, NUT_WALL_FN,
meshIndexing, patches,
dimDynamicViscosity, casepath
dimViscosity, casepath
);
// combustFlag is 1 in rectangular region, 0 or 1 elsewhere
// (although user could set it to another value)
@ -1312,9 +1312,9 @@ void write_uniformField
// outer
{
os.beginBlock("outer");
if (fieldName == "alphat" || fieldName == "mut")
if (fieldName == "alphat" || fieldName == "nut")
{
// Different b.c. for alphat & mut
// Different b.c. for alphat & nut
os.writeEntry("type", "calculated");
}
else
@ -1346,10 +1346,10 @@ void write_pU_fields
{
// Velocity field
{
OFstream os(casepath / pars.timeName / "Ubet");
OFstream os(casepath / pars.timeName / "U");
os.precision(outputPrecision);
make_header(os, "", volVectorField::typeName, "Ubet");
make_header(os, "", volVectorField::typeName, "U");
os.writeEntry("dimensions", dimVelocity);
@ -1406,6 +1406,7 @@ void write_pU_fields
os.writeEntry("cyclicPatch", word(patchName + "Cyclic_half0"));
os.writeEntry("openFraction", 0); // closed
os.writeEntry("openingTime", p.blowoffTime);
os.writeEntry("minThresholdValue", p.blowoffPress);
os.writeEntry("maxOpenFractionDelta", 0.1);
os.writeEntry("forceBased", "false");
os.writeEntry("opening", "true");
@ -1865,8 +1866,8 @@ void write_blocked_face_list
if (p.patchType > 0) // Panel
{
os.beginBlock(setName);
os.writeEntry("wallPatchName", word(patchName + "Wall"));
os.writeEntry("cyclicMasterPatchName", patchName);
os.writeEntry("wallPatch", word(patchName + "Wall"));
os.writeEntry("cyclicMasterPatch", word(patchName + "Cyclic_half0"));
os.endBlock();
}
}