diff --git a/src/fieldSources/basicSource/rotorDiskSource/profileModel/series/seriesProfile.C b/src/fieldSources/basicSource/rotorDiskSource/profileModel/series/seriesProfile.C index fb695078255d7f5015a3c85e9ce96a80473914ab..77dad6d37178158239de235a7583b35c5191dfcb 100644 --- a/src/fieldSources/basicSource/rotorDiskSource/profileModel/series/seriesProfile.C +++ b/src/fieldSources/basicSource/rotorDiskSource/profileModel/series/seriesProfile.C @@ -38,7 +38,24 @@ namespace Foam // * * * * * * * * * * * * Protected Member Functions * * * * * * * * * * * // -Foam::scalar Foam::seriesProfile::evaluate +Foam::scalar Foam::seriesProfile::evaluateDrag +( + const scalar& xIn, + const List<scalar>& values +) const +{ + scalar result = 0.0; + + forAll(values, i) + { + result += values[i]*cos((i + 1)*xIn); + } + + return result; +} + + +Foam::scalar Foam::seriesProfile::evaluateLift ( const scalar& xIn, const List<scalar>& values @@ -108,8 +125,8 @@ Foam::seriesProfile::seriesProfile void Foam::seriesProfile::Cdl(const scalar alpha, scalar& Cd, scalar& Cl) const { - Cd = evaluate(alpha, CdCoeffs_); - Cl = evaluate(alpha, ClCoeffs_); + Cd = evaluateDrag(alpha, CdCoeffs_); + Cl = evaluateLift(alpha, ClCoeffs_); } diff --git a/src/fieldSources/basicSource/rotorDiskSource/profileModel/series/seriesProfile.H b/src/fieldSources/basicSource/rotorDiskSource/profileModel/series/seriesProfile.H index 7c4793e28e9fbe954dd6bf72af0a22c314038665..5e2611d6c3a0503683e916f326a591333f2f2786 100644 --- a/src/fieldSources/basicSource/rotorDiskSource/profileModel/series/seriesProfile.H +++ b/src/fieldSources/basicSource/rotorDiskSource/profileModel/series/seriesProfile.H @@ -28,7 +28,7 @@ Description Series-up based profile data - drag and lift coefficients computed as sum of cosine series - Cd = sum_i(CdCoeff)*sin(i*AOA) + Cd = sum_i(CdCoeff)*cos(i*AOA) Cl = sum_i(ClCoeff)*sin(i*AOA) where: @@ -79,12 +79,21 @@ protected: // Protected Member Functions - //- Evaluate - scalar evaluate - ( - const scalar& xIn, - const List<scalar>& values - ) const; + // Evaluate + + //- Drag + scalar evaluateDrag + ( + const scalar& xIn, + const List<scalar>& values + ) const; + + //- Lift + scalar evaluateLift + ( + const scalar& xIn, + const List<scalar>& values + ) const; public: