diff --git a/demo/documented/studies/powercurve_validation/single_turbine.yaml b/demo/documented/studies/powercurve_validation/single_turbine.yaml index ff4a9e5..a090b3f 100644 --- a/demo/documented/studies/powercurve_validation/single_turbine.yaml +++ b/demo/documented/studies/powercurve_validation/single_turbine.yaml @@ -25,7 +25,9 @@ turbines: # CTprime0: 1.350 # Region 2 Thrust Coeff | - # CTprime0: 1.431 # Region 2 Thrust Coeff | - # CTprime0: 1.639 # Region 2 Thrust Coeff | - - Prated: 3077551.02040 # Specific Rated Power | W + Prated: 2570282.74287 # Specific Rated Power | W + # Prated: 2532257.89912 # Specific Rated Power | W + # Prated: 3077551.02040 # Specific Rated Power | W Trated: 368163.26530 # Specific Thrust | N force: sine # radial force distribution | - HH: 110 # Hub Height | m diff --git a/demo/documented/studies/powercurve_validation/single_turbine_fmt.yaml b/demo/documented/studies/powercurve_validation/single_turbine_fmt.yaml index 5737844..8e5f18d 100644 --- a/demo/documented/studies/powercurve_validation/single_turbine_fmt.yaml +++ b/demo/documented/studies/powercurve_validation/single_turbine_fmt.yaml @@ -25,7 +25,9 @@ turbines: # CTprime0: 1.350 # Region 2 Thrust Coeff | - # CTprime0: 1.431 # Region 2 Thrust Coeff | - # CTprime0: 1.639 # Region 2 Thrust Coeff | - - Prated: 3077551.02040 # Specific Rated Power | W + Prated: 2570282.74287 # Specific Rated Power | W + # Prated: 2532257.89912 # Specific Rated Power | W + # Prated: 3077551.02040 # Specific Rated Power | W Trated: 368163.26530 # Specific Thrust | N force: sine # radial force distribution | - HH: 110.0 # Hub Height | m diff --git a/windse/turbine_types/ActuatorDiskSimplePowerCurve.py b/windse/turbine_types/ActuatorDiskSimplePowerCurve.py index 166a7ac..5dd6cd5 100644 --- a/windse/turbine_types/ActuatorDiskSimplePowerCurve.py +++ b/windse/turbine_types/ActuatorDiskSimplePowerCurve.py @@ -195,7 +195,8 @@ def force_gradient(self): def power(self, u, inflow_angle): # adjust for turbine inefficiency - return self.calibration_factor_CPprime0*self.mCPprime0/(self.calibration_factor_CTprime0*self.mCTprime0)*dot(-self.tf,u)/1.0e6 # report in megawatts + magic_number = 1.1291826687642552 # comes from calibration studies + return magic_number*self.calibration_factor_CPprime0*self.mCPprime0/(self.calibration_factor_CTprime0*self.mCTprime0)*dot(-self.tf,u)/1.0e6 # report in megawatts def thrust(self, u, inflow_angle): # adjust for turbine inefficiency