diff --git a/src/finiteVolume/fields/fvPatchFields/derived/swirlFanVelocity/swirlFanVelocityFvPatchField.C b/src/finiteVolume/fields/fvPatchFields/derived/swirlFanVelocity/swirlFanVelocityFvPatchField.C index 5939ed1fe9..acd85ee22e 100644 --- a/src/finiteVolume/fields/fvPatchFields/derived/swirlFanVelocity/swirlFanVelocityFvPatchField.C +++ b/src/finiteVolume/fields/fvPatchFields/derived/swirlFanVelocity/swirlFanVelocityFvPatchField.C @@ -79,7 +79,7 @@ void Foam::swirlFanVelocityFvPatchField::calcFanJump() deltaP /rMag[i] /fanEff_ - /rpm_*constant::mathematical::pi/30.0; + /(rpm_*constant::mathematical::pi/30.0); } } } @@ -93,22 +93,12 @@ void Foam::swirlFanVelocityFvPatchField::calcFanJump() << exit(FatalError); } 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 const vectorField tangentialVelocity(magTangU*tanDir); - this->jump_ = tangentialVelocity + UtanIn; + this->jump_ = tangentialVelocity; } } diff --git a/src/finiteVolume/fields/fvPatchFields/derived/swirlFanVelocity/swirlFanVelocityFvPatchField.H b/src/finiteVolume/fields/fvPatchFields/derived/swirlFanVelocity/swirlFanVelocityFvPatchField.H index 8d485ebf38..ed8cc6f52b 100644 --- a/src/finiteVolume/fields/fvPatchFields/derived/swirlFanVelocity/swirlFanVelocityFvPatchField.H +++ b/src/finiteVolume/fields/fvPatchFields/derived/swirlFanVelocity/swirlFanVelocityFvPatchField.H @@ -33,7 +33,7 @@ Description 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: @@ -45,7 +45,7 @@ Description Alternatively an inner and outer radii can be used instead of rEff. The 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: