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

rotorDiskSource: Added output of rotor power

Browse Source 
Patch contributed by Robert Lee
Resolves patch request https://bugs.openfoam.org/view.php?id=3262
This commit is contained in:
Henry Weller
2019-05-26 12:11:33 +01:00
parent 47d2e03ae3
commit dfe7b86e70
2 changed files with 6 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
src/fvOptions/sources/derived/rotorDiskSource/rotorDiskSource.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2011-2018 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2019 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -384,7 +384,7 @@ void Foam::fv::rotorDiskSource::createCoordinateSystem()
const scalar sumArea = gSum(area_);
const scalar diameter = Foam::sqrt(4.0*sumArea/mathematical::pi);
Info<< " Rotor gometry:" << nl
Info<< " Rotor geometry:" << nl
<< " - disk diameter = " << diameter << nl
<< " - disk area = " << sumArea << nl
<< " - origin = " << coordSys_.origin() << nl

5
src/fvOptions/sources/derived/rotorDiskSource/rotorDiskSourceTemplates.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2011-2018 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2019 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -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;
}