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: Droping tangential U on the master for the swirling fan BC.

Browse Source
This commit is contained in:
sergio
2018-06-12 08:29:20 -07:00
parent c8d585b136
commit dcfaeef0ad
2 changed files with 5 additions and 15 deletions
Download Patch File Download Diff File Expand all files Collapse all files

16
src/finiteVolume/fields/fvPatchFields/derived/swirlFanVelocity/swirlFanVelocityFvPatchField.C
View File

@ -79,7 +79,7 @@ void Foam::swirlFanVelocityFvPatchField::calcFanJump()
deltaP deltaP
/rMag[i] /rMag[i]
/fanEff_ /fanEff_
/rpm_*constant::mathematical::pi/30.0; /(rpm_*constant::mathematical::pi/30.0);
} }
} }
} }
@ -93,22 +93,12 @@ void Foam::swirlFanVelocityFvPatchField::calcFanJump()
<< exit(FatalError); << exit(FatalError);
} }
magTangU = magTangU =
deltaP/rEff_/fanEff_/rpm_*constant::mathematical::pi/30.0; deltaP/rEff_/fanEff_/(rpm_*constant::mathematical::pi/30.0);
} }
// U upstream from fan
const vectorField Uup(this->patchInternalField());
// U normal to fan
const vectorField Un((Uup & patch().nf())*patch().nf());
// U tangential into fan
const vectorField UtanIn(Uup - Un);
// Calculate the tangential velocity // Calculate the tangential velocity
const vectorField tangentialVelocity(magTangU*tanDir); const vectorField tangentialVelocity(magTangU*tanDir);
this->jump_ = tangentialVelocity + UtanIn; this->jump_ = tangentialVelocity;
} }
} }

4
src/finiteVolume/fields/fvPatchFields/derived/swirlFanVelocity/swirlFanVelocityFvPatchField.H
View File

@ -33,7 +33,7 @@ Description
The U-jump is specified with a swirl component as follows: The U-jump is specified with a swirl component as follows:
Utan = deltaP/rEff/fanEff/rpm/pi/30.0; Utan = deltaP/rEff/fanEff/(rpm*pi/30.0);
where: where:
@ -45,7 +45,7 @@ Description
Alternatively an inner and outer radii can be used instead of rEff. The Alternatively an inner and outer radii can be used instead of rEff. The
Utan is as follow for r > rInner and r < rOuter Utan is as follow for r > rInner and r < rOuter
Utan = deltaP/r/fanEff/rpm/pi/30.0; Utan = deltaP/r/fanEff/(rpm/pi/30.0);
where: where: