Changeset 2152
 Timestamp:
 Feb 17, 2017 1:27:24 PM (5 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/wind_turbine_model_mod.f90
r2101 r2152 20 20 ! Current revisions: 21 21 !  22 ! 22 ! Bugfix in subroutine wtm_read_blade_tables 23 ! Addition of a tip loss model 23 24 ! 24 25 ! Former revisions: … … 104 105 LOGICAL :: speed_control = .FALSE. !< switch for use of speed controller 105 106 LOGICAL :: yaw_control = .FALSE. !< switch for use of yaw controller 107 LOGICAL :: tl_cor = .FALSE. !< switch for use of tip loss correct. 106 108 107 109 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,& 404 406 turb_cd_nacelle, turb_cd_tower, & 405 yaw_control, yaw_speed 407 yaw_control, yaw_speed, tl_cor 406 408 407 409 ! … … 1195 1197 ialpha = 1 1196 1198 1197 DO WHILE ( alpha_attack_i > alpha_attack_tab(ialpha) )1199 DO WHILE ( ( alpha_attack_i > alpha_attack_tab(ialpha) ) .AND. (ialpha <= dlen ) ) 1198 1200 ialpha = ialpha + 1 1199 1201 ENDDO … … 1369 1371 INTEGER(iwp), DIMENSION(1) :: lct=0 1370 1372 REAL(wp), DIMENSION(1) :: rad_d=0.0_wp 1373 1374 REAL(wp) :: tl_factor !< factor for tip loss correction 1371 1375 1372 1376 … … 1797 1801 ( (2.0_wp*pi) / 360.0_wp ) 1798 1802 1803 IF ( tl_cor ) THEN 1804 1805 ! Tip loss correction following Schito 1806 ! Schito applies the tip loss correction only to the lift force 1807 ! Therefore, the tip loss correction is only applied to the lift 1808 ! coefficient and not to the drag coefficient in our case 1809 ! 1810 tl_factor = ( 2.0 / pi ) * & 1811 ACOS( EXP( 1.0 * ( 3.0 * ( rr(inot)  cur_r ) / & 1812 ( 2.0 * cur_r * abs( sin( phi_rel(rseg) ) ) ) ) ) ) 1813 1814 turb_cl(rseg) = tl_factor * turb_cl(rseg) 1815 1816 END IF 1799 1817 ! 1800 1818 ! !! … … 2090 2108 ENDDO ! End of loop over turbines 2091 2109 2092 END IF 2110 END IF ! end of yaw control 2093 2111 2094 2112 IF ( speed_control ) THEN
Note: See TracChangeset
for help on using the changeset viewer.