4.2.1.9. Tail fin Aerodynamics Theory

4.2.1.9.1. Notations

Tail fin aerodynamic reference point

The tail fin aerodynamic reference point, \(\boldsymbol{x}_\text{ref}\), is the point where the aerodynamic loads are calculated on the tail fin. The structural solver computes the instantenous position, velocity, acceleration, of the reference point at each time step. The initial position of the reference point with respect to the tower top is a user input. Typical choices are the leading edge/apex of the fin or a point close to the center of pressure at zero angle of attack. The other aerodynamic inputs (e.g. aerodynamic moment coefficient) need to be consistent with the choice of the reference point.

Tail fin coordinate system

The inertial and tail fin coordinate systems are illustrated in Fig. 4.9. The transformation matrix from the inertial coordinate system to the tail fin coordinate system is given by \(\boldsymbol{R}_\text{tf,i}\).

../../../_images/TailFinCoord.png

Fig. 4.9 Coordinate systems and velocity vectors used for the tail fin aerodynamics

The reference orientation (when the structure is un-deflected), the transformation matrix is:

(4.1)\[\boldsymbol{R}_\text{tf,i} = \operatorname{EulerConstruct}(\theta_\text{bank}, \theta_\text{tilt}, \theta_\text{skew})\]

For a common application with a vertical fin, the three angles are zero. TODO: The order of the angles might be different in the current implementation (3-2-1) instead of (1-2-3) above)

Velocities

The following velocity vectors (3D vectors in global coordinates) are defined (see Fig. 4.9):

  • \(\boldsymbol{V}_\text{wind}\): Undisturbed Wind speed vector at the reference point

  • \(\boldsymbol{V}_\text{dist}\): Disturbed wind speed vector at the reference point (the disturbed wind contains the influence of the tower on the flow). AeroDyn has internal methods to compute \(\boldsymbol{V}_\text{dist}\) from \(\boldsymbol{V}_\text{wind}\). For now, we use “wind” but in the future we might use “dist”. In the theory below we would simply replace all the “wind” by “dist”.

  • \(\boldsymbol{V}_\text{elast}\): Structural translational velocity vector at the reference point

  • \(\boldsymbol{V}_\text{ind}\): Induced velocity from the wake at the reference point (assumed to be zero for now)

  • \(\boldsymbol{\omega}\): Structural rotational velocity of the fin

All velocities (except for \(\boldsymbol{V}_\text{ind}\) and \(\boldsymbol{V}_\text{dist}\) which are computed internally by AeroDyn) are provided as input to the aerodynamic solver. The relative wind experienced by the airfoil is given by:

(4.2)\[\boldsymbol{V}_\text{rel} = \boldsymbol{V}_\text{wind} -\boldsymbol{V}_\text{elast} +\boldsymbol{V}_\text{ind}\]

Angle of attack

The angle of attack is defined in the \(x_\text{tf}-y_\text{tf}\) plane of the tail fin coordinate systems as illustrated in Fig. 4.10.

../../../_images/TailFinAirfoilCoord.png

Fig. 4.10 Tail fin airfoil coordinate system and definition of angle of attack in the x-y plane

We write \(V_{\text{rel},\perp}\) the projection of \(\boldsymbol{V}_\text{rel}\) in this plane. The angle of attack is given by the components of this vector:

(4.3)\[\alpha = \arctan\frac{V_{\text{rel},y_\text{tf}}}{V_{\text{rel},x_\text{tf}}}\]

In this implementation, the function atan2 is used to compute the angle of attack.

Loads

If the dimensionless coefficients are known, they can be projected in the \(x_\text{tf}-y_\text{tf}\) plane as follows:

(4.4)\[C_{x_\text{tf}}(\alpha) = -C_l(\alpha) \sin\alpha + C_d(\alpha)\cos\alpha ,\quad C_{y_\text{tf}}(\alpha) = C_l(\alpha) \cos\alpha + C_d(\alpha)\sin\alpha\]

and the loads are therefore given by:

(4.5)\[f_{x_\text{tf}} = \frac{1}{2}\rho V_{\text{rel},\perp}^2 A \,C_{x_\text{tf}}(\alpha) ,\quad f_{y_\text{tf}} = \frac{1}{2}\rho V_{\text{rel},\perp}^2 A \,C_{y_\text{tf}}(\alpha) ,\quad m_{z_\text{tf}} = \frac{1}{2}\rho V_{\text{rel},\perp}^2 Ac \, C_m(\alpha)\]

Once the loads are known in the tail fin coordinate systems, they are transferred to the inertial system as follows:

(4.6)\[\begin{split}\left.\boldsymbol{f}\right|_{i} = \boldsymbol{R}_\text{tf,i}^t \left.\boldsymbol{f}\right|_\text{tf} = \boldsymbol{R}_\text{tf,i}^t \begin{bmatrix} f_{x_\text{tf}}\\ f_{y_\text{tf}}\\ 0\\ \end{bmatrix} ,\qquad \left.\boldsymbol{m}\right|_{i} = \boldsymbol{R}_\text{tf,i}^t \left.\boldsymbol{m}\right|_\text{tf} = \boldsymbol{R}_\text{tf,i}^t \begin{bmatrix} 0\\ 0\\ m_{z_\text{tf}}\\ \end{bmatrix}\end{split}\]

Induced velocity

The induced velocity from the wake at the reference point will affect the relative wind and therefore the angle of attack of the tail fin. Different models are implemented to compute this induced velocity. As a first approximation, this velocity may be set to zero (corresponding to the input TFinIndMod=0):

(4.7)\[\boldsymbol{V}_\text{ind}=0\]

The rotor-averaged induced velocity can also be used as an estimate (TFinIndMod=1). It is computed as the mean induced velocity over all the blade and aerodynamic nodes

(4.8)\[\boldsymbol{V}_\text{ind}=\frac{1}{n_B n_r}\sum_{i_b=1..n_B} \sum_{i_r=1..n_r} \boldsymbol{V}_{\text{ind},\text{blade}}[i_b, i_r]\]

Where \(\boldsymbol{V}_{\text{ind},\text{blade}}[i_b, i_r]\) is the induced velocity vector for blade \(i_b\) and at the radial node \(i_r\). NOTE: This averaging corresponds to what is done for the disk-average of the inflow in AeroDyn. In the future, we can use something weighted by the radius, or using precomputed coefficients, as done by Envision.

More advanced models could set the induced velocity to zero when outside of the wake boundary, or include a tower-shadow-like wake model. Such option is not yet available.

4.2.1.9.2. Polar-based model

In the polar-based model, the user provides the aerodynamic coefficients \(C_l, C_d, C_m\), as tabulated data, functions of the angle of attack. The aerodynamic moment is assumed to be provided at the reference point. A common practice is to use the center of pressure at zero angle of attack for polar data, so the user might want to chose such a point as the reference point of the fin. The tabulated data are provided as part of the list of airfoils given with AFNames in the AeroDyn input file. The user only needs to indicate the index TFinAFIndex within the list AFNames to indicate which polar to use for the tail fin.

4.2.1.9.3. Unsteady slender body model

The unsteady slender body (USB) model is documented in [ad-HW22].

The theory will be implemented and documented in a future release.