Untitled

def SpecificHeat_Cp(T: float) -> float:  # Функция поиска удельной теплоемкости
    return Cp * (T/T0[0]) ** fi  # Удельная теплоемкость

def DynamicViscosity_mu(T: float) -> float:  # Функция поиска динамической вязкости
    return mu * (T/T0[0]) ** n  # Динамическая вязкость

def CoeffThermalConductivity_Lambda(T: float) -> float:  # Функция поиска коэффициента теплопроводности
    return Lambda * (T/T0[1]) ** z  # Коэффициент теплопроводности

def findDefiningTemperature (lam_or_turb: float,T: float, M: float) -> float: # Функция поика температуры восстановления и определяющей температуры
    rl = 0.83
    rt = 0.88
    r = rl
    if (lam_or_turb == 1):
        r = rt
    Tr = T * (1 + r * (k - 1) * 0.5 * M ** 2)
    T_z = 0.5 * (Tst + T) + 0.22 * (Tr - T)
    T_z_previous = 0
    while abs(T_z - T_z_previous) > (T_z * 0.001):
        T_z_previous = T_z
        Cpi = SpecificHeat_Cp(T_z)
        mui = DynamicViscosity_mu(T_z)
        Lambdai = CoeffThermalConductivity_Lambda(T_z)
        Prl = (Cpi * mui) / Lambdai
        r = np.sqrt(Prl)
        if (lam_or_turb == 0):
            r = Prl ** (1/3)
    return Tr,T_z

# Для граней 1-4 температура восстановления и определяющая температура
l = 0
t = 1
Trl1, Tl1_z = findDefiningTemperature(l,T1,M1)
Trl2, Tl2_z = findDefiningTemperature(l,T2,M2)
Trl3, Tl3_z = findDefiningTemperature(l,T3,M3)
Trl4, Tl4_z = findDefiningTemperature(l,T4,M4)
Trt1, Tt1_z = findDefiningTemperature(t,T1,M1)
Trt2, Tt2_z = findDefiningTemperature(t,T2,M2)
Trt3, Tt3_z = findDefiningTemperature(t,T3,M3)
Trt4, Tt4_z = findDefiningTemperature(t,T4,M4)