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)

Arguments

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

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

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

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.

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.

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 IsInWake. If that is the case, JensenFactor is used to compute the penalty factor.

Examples

Run this code

# NOT RUN {
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
# }