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

fanPressureFvPatchScalarField: Corrected flow-rate

Browse Source
This commit is contained in:
Henry
2012-03-07 19:51:58 +00:00
parent a15b241d3a
commit 2c9245160f
1 changed files with 5 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
src/finiteVolume/fields/fvPatchFields/derived/fanPressure/fanPressureFvPatchScalarField.C
View File

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | \\ / O peration |
\\ / A nd | Copyright (C) 2011 OpenFOAM Foundation \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -134,17 +134,17 @@ void Foam::fanPressureFvPatchScalarField::updateCoeffs()
int dir = 2*direction_ - 1; int dir = 2*direction_ - 1;
// Average volumetric flow rate // Average volumetric flow rate
scalar aveFlowRate = 0; scalar volFlowRate = 0;
if (phi.dimensions() == dimVelocity*dimArea) if (phi.dimensions() == dimVelocity*dimArea)
{ {
aveFlowRate = dir*gSum(phip)/gSum(patch().magSf()); volFlowRate = dir*gSum(phip);
} }
else if (phi.dimensions() == dimVelocity*dimArea*dimDensity) else if (phi.dimensions() == dimVelocity*dimArea*dimDensity)
{ {
const scalarField& rhop = const scalarField& rhop =
patch().lookupPatchField<volScalarField, scalar>(rhoName()); patch().lookupPatchField<volScalarField, scalar>(rhoName());
aveFlowRate = dir*gSum(phip/rhop)/gSum(patch().magSf()); volFlowRate = dir*gSum(phip/rhop);
} }
else else
{ {
@ -157,7 +157,7 @@ void Foam::fanPressureFvPatchScalarField::updateCoeffs()
} }
// Pressure drop for this flow rate // Pressure drop for this flow rate
const scalar pdFan = fanCurve_(max(aveFlowRate, 0.0)); const scalar pdFan = fanCurve_(max(volFlowRate, 0.0));
totalPressureFvPatchScalarField::updateCoeffs totalPressureFvPatchScalarField::updateCoeffs
( (