OpenEngine/parameters.py at main · WiktorBoltrukiewicz/OpenEngine · 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
"""
parameters.py — Obliczanie parametrow fizycznych strony gazowej w danym punkcie dyszy.

Wersja uproszczona (bez chlodzenia regeneracyjnego).
Oblicza tylko parametry potrzebne do eksportu siatki dla Ansysa:
  - temperatura adiabatyczna scianki (T_aw)
  - liczba Macha (M)
"""

import numpy as np


def compute_gas_parameters(YSol, params, i):
    """
    Oblicza parametry strony gazowej w i-tym punkcie siatki.

    Parameters
    ----------
    YSol : np.ndarray, shape (n, 3) — rozwiazanie ODE [N, P, T]
    params : dict — slownik ze wszystkimi parametrami silnika
    i : int — indeks punktu (0-indexed)

    Returns
    -------
    T_aw : float — temperatura adiabatyczna scianki [K]
    M : float — liczba Macha [-]
    """
    gamma = params['gamma_arr'][i]
    Prcg = params['Prcg_arr'][i]

    N = YSol[i, 0]
    Ts = YSol[i, 2]

    # --- Temperatura adiabatyczna scianki ---
    T_aw = Ts * (
        (1 + Prcg ** 0.33 * ((gamma - 1) / 2) * N)
        / (1 + ((gamma - 1) / 2) * N)
    )

    M = np.sqrt(N)

    return T_aw, M