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

flowRateInletVelocity extrapolated: Removed reverse flow and correct only the normal component

Browse Source 
Improved stability and convergence.
This commit is contained in:
Henry Weller
2017-07-17 12:18:47 +01:00
parent ba50e560aa
commit 4fa0192e1d
2 changed files with 24 additions and 12 deletions
Download Patch File Download Diff File Expand all files Collapse all files

30
src/finiteVolume/fields/fvPatchFields/derived/flowRateInletVelocity/flowRateInletVelocityFvPatchVectorField.C
View File

@ -156,31 +156,43 @@ void Foam::flowRateInletVelocityFvPatchVectorField::updateValues
{
const scalar t = db().time().timeOutputValue();
tmp<vectorField> n = patch().nf();
const vectorField n(patch().nf());
if (extrapolateProfile_)
{
vectorField newValues(this->patchInternalField());
vectorField Up(this->patchInternalField());
scalar flowRate = flowRate_->value(t);
scalar estimatedFlowRate = -gSum(rho*(this->patch().Sf() & newValues));
// Patch normal extrapolated velocity
scalarField nUp(n & Up);
// Remove the normal component of the extrapolate patch velocity
Up -= nUp*n;
// Remove any reverse flow
nUp = min(nUp, 0.0);
const scalar flowRate = flowRate_->value(t);
const scalar estimatedFlowRate = -gSum(rho*(this->patch().magSf()*nUp));
if (estimatedFlowRate/flowRate > 0.5)
{
newValues *= (mag(flowRate)/mag(estimatedFlowRate));
nUp *= (mag(flowRate)/mag(estimatedFlowRate));
}
else
{
newValues -=
((flowRate - estimatedFlowRate)/gSum(rho*patch().magSf()))*n;
nUp -= ((flowRate - estimatedFlowRate)/gSum(rho*patch().magSf()));
}
this->operator==(newValues);
// Add the corrected normal component of velocity to the patch velocity
Up += nUp*n;
// Correct the patch velocity
this->operator==(Up);
}
else
{
const scalar avgU = -flowRate_->value(t)/gSum(rho*patch().magSf());
operator==(n*avgU);
operator==(avgU*n);
}
}

6
src/finiteVolume/fields/fvPatchFields/derived/flowRateOutletVelocity/flowRateOutletVelocityFvPatchVectorField.C
View File

@ -160,10 +160,10 @@ void Foam::flowRateOutletVelocityFvPatchVectorField::updateValues
Up -= nUp*n;
// Remove any reverse flow
nUp = min(nUp, 0.0);
nUp = max(nUp, 0.0);
const scalar flowRate = flowRate_->value(t);
const scalar estimatedFlowRate = -gSum(rho*(this->patch().magSf()*nUp));
const scalar estimatedFlowRate = gSum(rho*(this->patch().magSf()*nUp));
if (estimatedFlowRate/flowRate > 0.5)
{
@ -171,7 +171,7 @@ void Foam::flowRateOutletVelocityFvPatchVectorField::updateValues
}
else
{
nUp -= ((flowRate - estimatedFlowRate)/gSum(rho*patch().magSf()));
nUp += ((flowRate - estimatedFlowRate)/gSum(rho*patch().magSf()));
}
// Add the corrected normal component of velocity to the patch velocity