Changeset 4528 for palm


Ignore:
Timestamp:
May 11, 2020 2:14:09 PM (7 months ago)
Author:
oliver.maas
Message:

added wind_turbine_parameters-namelist parameter smearing_kernel_size

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/wind_turbine_model_mod.f90

    r4497 r4528  
    2626! -----------------
    2727! $Id$
     28! added namelist parameter smearing_kernel_size
     29!
     30! 4497 2020-04-15 10:20:51Z raasch
    2831! file re-formatted to follow the PALM coding standard
    2932!
     
    231234    REAL(wp) ::  segment_width_radial       = 0.5_wp  !< width of the segments, the rotor area is divided into
    232235                                                      !< (in radial direction, as factor of MIN(dx,dy,dz))
     236   
     237    REAL(wp) ::  smearing_kernel_size       = 2.0_wp  !< size of the smearing kernel as multiples of dx
     238
    233239    REAL(wp) ::  time_turbine_on            = 0.0_wp  !< time at which turbines are started
    234240    REAL(wp) ::  tilt_angle                 = 0.0_wp  !< vertical tilt_angle of the rotor [degree] ( positive = backwards )
     
    530536                nacelle_radius, rotor_radius, segment_length_tangential, segment_width_radial,     &
    531537                region_2_slope, speed_control, tilt_angle, time_turbine_on,                        &
    532                 tower_cd, pitch_rate,                                                              &
     538                smearing_kernel_size, tower_cd, pitch_rate,                                        &
    533539                yaw_control, yaw_speed, tip_loss_correction
    534540!                , nacelle_cd
     
    10681074!-- The kernel is set equal to twice the grid spacing which has turned out to be a reasonable
    10691075!-- value (see e.g. Troldborg et al. (2013), Wind Energy, DOI: 10.1002/we.1608):
    1070     eps_kernel = 2.0_wp * dx
     1076    eps_kernel = smearing_kernel_size * dx
    10711077!
    10721078!-- The zero point (eps_min) of the polynomial function must be the following if the integral of
Note: See TracChangeset for help on using the changeset viewer.