wflo (version 1.9)

PairPenalty: Returns the Jensen wake penalty factor for a pair of turbines.

Description

As seen from a turbine in the wind farm, computes the wake penalty factor for another turbine in that farm.

Usage

PairPenalty(x1, y1, x2, y2, Dirs, SDs)

Value

PairPenalty returns a single number between 0 and 1. If point 2 is not in the wake of point 1, the function returns 1.

Arguments

x1

must be a single value. Provide the x location of the first turbine.

y1

must be a single value. Provide the y location of the first turbine.

x2

must be a single value. Provide the x location of the second turbine.

y2

must be a single value. Provide the y location of the second turbine.

Dirs

a matrix containing average yearly wind directions. Usually, the third element of the list object FarmData will be used as this matrix.

SDs

a matrix containing average yearly wind direction standard deviations. Usually, the fourth element of the list object FarmData will be used as this matrix.

Author

Carsten Croonenbroeck

Details

First, this function uses GetAngle to compute the angle between the two points provided, as seen from point 2's point of view. It then obtains the wind direction at point 2 using GetDirInfo. After that, the distance between the two points is computed. With it, the wake cone is computed using JensenAngle to check whether point 2 is in point 1's wake using JensenTrapezoid. If that is the case, JensenFactor is used to compute the penalty factor.
Note that the penalty is the deduction to wind speed. It applies to wind power by its third power, so the user is responsible to take it to its cube himself. Profit does that automatically internally.

See Also

Use JensenFactor to see how this function operates. See FarmVars for the data object.

Examples

Run this code
Dirs <- FarmData[[3]][e$FarmVars$StartPoint:e$FarmVars$EndPoint,
e$FarmVars$StartPoint:e$FarmVars$EndPoint]
SDs <- FarmData[[4]][e$FarmVars$StartPoint:e$FarmVars$EndPoint,
e$FarmVars$StartPoint:e$FarmVars$EndPoint]
PairPenalty(0.9, 0.8, 0.6, 0.9, Dirs, SDs)
## Weak wake penalty
PairPenalty(0.1, 0.1, 0.6, 0.9, Dirs, SDs)
## No wake penalty

Run the code above in your browser using DataLab