meilander_comp/Air.py at master · rknop/meilander_comp · GitHub

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
# -*- coding: utf-8 -*-
"""
Created on Mon Jan 27 14:17:41 2020

@author: omeil
"""

import numpy as np

class Air(object):

    def __init__(self, vo, r, mAir = (0.02897/6.022e23)):
        self.vo = np.array((vo))
        self.v = np.copy((vo))
        self.ro = np.array((r))
        self.r = np.copy((r))
        self.mAir = mAir

    def denAir(self, y):
        if (y < 11000):
            return (self.calcP1(y) / (0.2869 * (self.calcT1(y) + 273.1))) * (1 / self.mAir)

        elif (y > 25000):
            return (self.calcP3(y) / (0.2869 * (self.calcT3(y) + 273.1))) * (1 / self.mAir)

        else:
            return (self.calcP2(y) / (0.2869 * (self.calcT2(y) + 273.1))) * (1 / self.mAir)


#=============================================================================
    #found at https://www.grc.nasa.gov/WWW/K-12/airplane/atmosmet.html
    def calcT1(self, y):
        return (15.04 - (.00649 * y))

    def calcP1(self, y):
        return 101.29 * ((self.calcT1(y) + 273.1) / 288.08)

    def calcT2(self, y):
        return (-56.46)

    def calcP2(self, y):
        return 22.65 * np.exp(1.73 - .000157 * y)

    def calcT3(self, y):
        return (-131.27 + .00299 * y)

    def calcP3(self, y):
        return 2.48 * ((self.calcT3(y) + 273.1) / 288.08) ** (-11.388)