zhyang-dev/openfoam
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

BUG: Fixed parallel operation of activePressureForceBaffle

Browse Source
This commit is contained in:
sergio
2016-06-02 09:44:40 -07:00
parent 91d4bab09c
commit 9b0cd898a6
1 changed files with 17 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

21
src/finiteVolume/fields/fvPatchFields/derived/activePressureForceBaffleVelocity/activePressureForceBaffleVelocityFvPatchVectorField.C
View File

@ -279,8 +279,6 @@ void Foam::activePressureForceBaffleVelocityFvPatchVectorField::updateCoeffs()
{
valueDiff -=p[nbrFaceCells[facei]]*mag(initCyclicSf_[facei]);
}
Info<< "Force difference = " << valueDiff << endl;
}
else //pressure based
{
@ -293,8 +291,20 @@ void Foam::activePressureForceBaffleVelocityFvPatchVectorField::updateCoeffs()
{
valueDiff -= p[nbrFaceCells[facei]];
}
}
Info<< "Pressure difference = " << valueDiff << endl;
reduce(valueDiff, sumOp<scalar>());
if (Pstream::master())
{
if (fBased_)
{
Info<< "Force difference = " << valueDiff << endl;
}
else
{
Info<< "Pressure difference = " << valueDiff << endl;
}
}
if (mag(valueDiff) > mag(minThresholdValue_) || baffleActivated_)
@ -322,7 +332,10 @@ void Foam::activePressureForceBaffleVelocityFvPatchVectorField::updateCoeffs()
openFraction_ = max(min(1 - 1e-6, openFraction_), 1e-6);
}
Info<< "Open fraction = " << openFraction_ << endl;
if (Pstream::master())
{
Info<< "Open fraction = " << openFraction_ << endl;
}
scalar areaFraction = 0.0;