Skip to content
Permalink
Browse files

Add power calculation to rotorDisk

  • Loading branch information...
openfoamtutorials committed May 22, 2019
1 parent 0889ff9 commit 182507ddf566b069b1d54e76bb80f63087d2d5e5
Showing with 3 additions and 0 deletions.
  1. +3 −0 src/fvOptions/sources/derived/rotorDiskSource/rotorDiskSourceTemplates.C
@@ -49,6 +49,7 @@ void Foam::fv::rotorDiskSource::calculate
scalar liftEff = 0.0;
scalar AOAmin = great;
scalar AOAmax = -great;
scalar powerEff = 0.0;

forAll(cells_, i)
{
@@ -127,6 +128,7 @@ void Foam::fv::rotorDiskSource::calculate
// Accumulate forces
dragEff += rhoRef_*localForce.y();
liftEff += rhoRef_*localForce.z();
powerEff += rhoRef_ * localForce.y() * radius * omega_;

// Transform force from local coning system into rotor cylindrical
localForce = invR_[i] & localForce;
@@ -151,6 +153,7 @@ void Foam::fv::rotorDiskSource::calculate
Info<< type() << " output:" << nl
<< " min/max(AOA) = " << radToDeg(AOAmin) << ", "
<< radToDeg(AOAmax) << nl
<< " Effective power = " << powerEff << nl
<< " Effective drag = " << dragEff << nl
<< " Effective lift = " << liftEff << endl;
}

0 comments on commit 182507d

Please sign in to comment.
You can’t perform that action at this time.