-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathREADME.windturbine
113 lines (100 loc) · 6.64 KB
/
README.windturbine
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
README for wind turbine drag parameterization schemes, added 20101220. JM
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!! NOTICE
!! The following paper should be cited whenever presenting results using this scheme
!! (using either the original version or any modified versions of the scheme):
!! Fitch, A. C. et al. 2012: Local and Mesoscale Impacts of Wind Farms as Parameterized in a
!! Mesoscale NWP Model. Monthly Weather Review, doi:http://dx.doi.org/10.1175/MWR-D-11-00352.1
!!
!! Anna C. Fitch, National Center for Atmospheric Research (formerly University of Bergen)
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Represents kinetic energy extracted by wind turbines and turbulence
(TKE) they produce at model levels within the rotor area. The thrust and
power coefficient curves included to calculate the momentum sink and source of TKE are
generic and an approximation to a real turbine. These coefficients should be
obtained from the turbine manufacturer for the turbines of interest and incorporated
into the code in subroutine dragcof (we could not include real curves since they
are proprietary).
References:
Fitch, A. C. et al. 2012: Local and Mesoscale Impacts of Wind Farms as Parameterized in a
Mesoscale NWP Model. Monthly Weather Review, doi:http://dx.doi.org/10.1175/MWR-D-11-00352.1
Fitch, A. C. et al. 2013: Mesoscale Influences of Wind Farms Throughout a Diurnal Cycle.
Monthly Weather Review, doi:http://dx.doi.org/10.1175/MWR-D-12-00185.1
Fitch, A. C. et al. 2013: Parameterization of Wind Farms in Climate Models.
Journal of Climate, doi:http://dx.doi.org/10.1175/JCLI-D-12-00376.1
*Compiling the code.* no different from normal
*Running the code.* There are new variables added to the phys block of
the namelist.input file. The string valued variable windturbines_spec
controls the overall operation of the scheme. If it's set to "none"
(quotes included) the wind turbine scheme is off. If it's set to
"ideal", the turbine drag scheme is in idealized mode and other namelist
variables in the block control the geometry and characteristics of the
wind farm. If it's set to the name of a file in the run directory --
for example, "windspec.in" -- then the scheme is on and the location and
characteristics of each turbine is set individually from the contents of
that file. Which scheme is active will be controlled by another variable
TBD (not yet implemented, since there's only one scheme right now).
*Idealized configuration.* If windturbines_spec is set to "ideal", the
geometry and characteristics of the farm are controlled by the following
additional variables in the phys block of the namelist.input file:
name type sample description
value
turbgridid integer 1 id of the WRF domain containing the farm
hubheight real 100. height in meters of all turbine hubs in farm
diameter real 100. diameter in meters of all turbine rotors in farm
stdthrcoef real 0.158 standing thrust coefficient of all turbines in farm
cutinspeed real 3.5 cut in speed (m/s) of all turbines in farm
cutoutspeed real 30. cut out speed (m/s) of all turbines in farm
power real 5.0 nominal power (MW) of all turbines in farm
turbpercell real 1.0 number of turbines per grid cell
ewfx integer 10 x-extent of rectangular wind farm in grid cells
ewfy integer 10 y-extent of rectangular wind farm in grid cells
pwfx integer 45 x-coordinate of grid cell in sw corner of farm
pwfy integer 45 y-coordinate of grid cell in sw corner of farm
*Real-world configuration.* The location and characteristics of each
turbine is specified individually in the file named by the variable
windturbines_spec. Each line of the file specifies one turbine.
The entries of a line are separated by spaces. The entries are listed
in order on the line and specify the following, by position on the line:
1. The WRF grid id into which the turbine is placed [integer]
2. True latitude of the wind turbine in degrees [real]
3. True longitude of the wind turbine in degrees [real]
4. Height in meters of the turbine hub [real]
5. Diameter in meters of the rotor [real]
6. Standing thrust coefficient [real]
7. Nominal power of turbine (MW) [real]
8. Cutin speed (m/s) [real]
9. Cutout speed (m/s) [real]
The location of each turbine is specified using the lat and lon elements
for the turbine's entry in the windspec file. That is, the i,j index
in the grid is computed from the true latitude (entry 2) and the true
longitude (entry 3). The coordinate in the j dimension is computed by
checking true latitude against the LAT_V field in the wrfinput initial
conditions file. The i dimension is computed by checking the true
longitude against LONG_U field in the wrfinput file.
*Pseudo-real configuration.* This mode is only for testing the real-world
specification mechanism. The contents of the windspec.in file in the
em_seabreeze2d_x directory contains a specification that is identical
to the ideal specification shown above. However, instead of true
latitude and true longitude, elements 2 and 3 of each entry specify
the i,j grid coordinates of the turbine. The code knows to treat
these as grid coordinates instead of latitude and longitude because
for idealized WRF cases, the LAT and LONG fields in wrfinput files are
set to all zero. When the model reads a wrfinput file for an idealized
case, it checks to see if the lat and lon fields of the case have data
in them or if they are just zeros. If they do not have data in them
(just have zeros) fields 2 and 3 of the windspec file are treated as i,j
coordinates are used. A sample windspec file for a pseudo-real run is
in test/em_seabreeze2d_x/windspec.in).
*Specific note for the Fitch scheme.* The Fitch scheme is based on
Blahak et al, 2010, of Wetter Jetzt GbR. It differs in that the kinetic
energy (KE) extracted is based on the thrust coefficient of the turbine (representing
the total fraction of KE extracted from the atmosphere due to the turbine), rather
than the power coefficient plus an estimated loss factor. In addition the TKE
generated is a function of wind speed. The specific implementation
comes from Anna Fitch, Alok Gupta, and Idar Barstad at Uni Bergen,
Norway. It was added to this release of WRF by Jimy Dudhia, Joe Olson
(NOAA), Julie Lundquist (U. Colorado/NREL), and John Michalakes (NREL).
It works with the MYNN PBL and modifies the QKE field representing 2xTKE
to include the TKE produced by wind turbines. QKE is advected as a part
of the scalar 4D tracer array in WRF.