3. Aerodynamic Module Documentation

3.1. Lifting Line Theory

If \(\theta_0\) is an arbitrary span-wise location:

(1)\[\alpha (\theta_o)=\frac{2b}{\pi c(\theta_o)} \sum_1^N A_n sin(n \theta_o) + \alpha_{L=0}(\theta_o) + \sum_1^N n A_n \frac{sin(n\theta_o)}{sin(\theta_o)}\]

Each equation has \(N\) unknowns (\(A_n\)), so if there are N \(\theta_o\), we have NxN system, which in Einstein notation can be written as:

(2)\[C_{ij}A_{i}=D_{i}\]

where, \(i=0,...,N\), \(j=0,...,N\) and :

(3)\[C_{ij}= \left( \frac{2b}{\pi c(j)} + \frac{n}{sin \theta(i)} \right) sin(n \theta(i))\]
(4)\[A_i=A(i)\]
(5)\[D_i=\alpha(i)-\alpha_{L=0}(i)\]

where \(n=1,3,5,...,N-1\). Since we are considering a symmetric wing, all of the even terms would cancel each other

_images/elliptical_LLT.png

3.2. The code

Created on Mon Jul 20 17:26:19 2015

@author: Pedro Leal

aero_module.LLT_calculator(alpha_L_0_root, c_D_xfoil, N=10, b=10.0, taper=1.0, chord_root=1, alpha_root=0.0, V=1.0)[source]

Calculate the coefficients for a Wing. TODO : - Include elliptical wing

  • When alpha_L_0_root = zero, nan!
  • Include non rectangular wings
  • something else?
aero_module.Reynolds(height, V, c)[source]

Simple function to calculate Reynolds for a given height.

@author: Pedro Leal Created in Jul 17 2015

aero_module.air_properties(height, unit='feet')[source]

Function to calculate air properties for a given height (m or ft).

Sources:

Created on Thu May 15 14:59:43 2014 @author: Pedro Leal

aero_module.pressure_shell_2D(Data, chord, thickness, half_span, height, Velocity, N, txt=False)[source]

Calculate pressure field for a 2D Shell.