Changeset 2152 for palm/trunk

Timestamp:

Feb 17, 2017 1:27:24 PM (8 years ago)

Author:

lvollmer

Message:

bugfix in wtm_read_blade_tables, Addition of a tip loss model

File:

• palm/trunk/SOURCE/wind_turbine_model_mod.f90 (modified) (7 diffs)

palm/trunk/SOURCE/wind_turbine_model_mod.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Bugfix in subroutine wtm_read_blade_tables
+! Addition of a tip loss model
 ! 
 ! Former revisions:
@@ -104 +105 @@
     LOGICAL ::  speed_control = .FALSE.   !< switch for use of speed controller
     LOGICAL ::  yaw_control   = .FALSE.   !< switch for use of yaw controller
+    LOGICAL ::  tl_cor        = .FALSE.    !< switch for use of tip loss correct.
 
     REAL(wp) ::  segment_length  = 1.0_wp          !< length of the segments, the rotor area is divided into
@@ -403 +405 @@
                                   slope2, speed_control, tilt, time_turbine_on,&
                                   turb_cd_nacelle, turb_cd_tower,              &
-                                  yaw_control, yaw_speed
+                                  yaw_control, yaw_speed, tl_cor
 
 !
@@ -1195 +1197 @@
              ialpha = 1
 
-             DO WHILE ( alpha_attack_i > alpha_attack_tab(ialpha) )
+             DO WHILE ( ( alpha_attack_i > alpha_attack_tab(ialpha) ) .AND. (ialpha <= dlen ) )
                 ialpha = ialpha + 1
              ENDDO
@@ -1369 +1371 @@
        INTEGER(iwp), DIMENSION(1) :: lct=0
        REAL(wp), DIMENSION(1)     :: rad_d=0.0_wp
+       
+       REAL(wp) :: tl_factor !< factor for tip loss correction
 
 
@@ -1797 +1801 @@
                                         ( (2.0_wp*pi) / 360.0_wp )
 
+                   IF ( tl_cor )  THEN
+                   
+!--                  Tip loss correction following Schito
+!--                  Schito applies the tip loss correction only to the lift force
+!--                  Therefore, the tip loss correction is only applied to the lift 
+!--                  coefficient and not to the drag coefficient in our case 
+!--                  
+                     tl_factor = ( 2.0 / pi ) *                                &
+                          ACOS( EXP( -1.0 * ( 3.0 * ( rr(inot) - cur_r ) /     &
+                          ( 2.0 * cur_r * abs( sin( phi_rel(rseg) ) ) ) ) ) )
+                          
+                     turb_cl(rseg)  = tl_factor * turb_cl(rseg)                                  
+                                 
+                   END IF               
 !
 !--                !-----------------------------------------------------!
@@ -2090 +2108 @@
              ENDDO ! End of loop over turbines
                            
-          END IF 
+          END IF  ! end of yaw control 
           
           IF ( speed_control )  THEN