28 double funCORRMass() {
29 return MassFlow * sqrt(IT_K / TREF_K) / (IP / PREF);
33 double funUnCORRMass() {
34 return MassFlowC * (IP / PREF) / sqrt(IT_K / TREF_K);
38 double funCORRRTC(
double RTC) {
39 return RTC / sqrt(IT_K / TREF_K);
43 double funUnCORRRTC() {
44 return RTC_C * sqrt(IT_K / TREF_K);
48 double funTiso(
double g,
double sig) {
49 return IT_K * pow(PR, sig * (g - 1) / g);
54 return IT_K +
pow2(MassFlow * 287 * IT_K / IP / SecIn) / 2 / Cp;
58 double funP0(
double g) {
59 return IP * pow(T0_K / IT_K, (g - 1) / g);
91 inline double K_T_H1() {
96 inline double K_H1_H2() {
101 inline double K_H2_H3() {
106 inline double K_H3_C() {
111 inline double hAk_GAS_T(
double Re_mass_tur,
double Re_shaft) {
112 return 0.306 * pow(Re_mass_tur, 0.965) * pow(Re_shaft, -0.417);
116 inline double hAk_H1_OIL(
double Re_mass_oil,
double Re_shaft,
double Pr_oil,
double mu_oil_h1) {
117 return 2.7 * pow(Re_mass_oil, 1.46) * pow(Re_shaft, -1.12) * cbrt(Pr_oil) * pow(mu_oil_h1, -0.69);
121 inline double hAk_H2_OIL(
double Re_mass_oil,
double Re_shaft,
double Pr_oil,
double mu_oil_h2) {
122 return 6.2e-3 * pow(Re_mass_oil, 3.83) * pow(Re_shaft, -2.02) * cbrt(Pr_oil) * pow(mu_oil_h2, -1.02);
126 inline double hAk_H3_OIL(
double Re_mass_oil,
double Re_shaft,
double Pr_oil,
double mu_oil_h3) {
127 return 6.6e-3 * pow(Re_mass_oil, -0.128) * pow(Re_shaft, 1.25) * cbrt(Pr_oil) * pow(mu_oil_h3, -0.876);
131 inline double hAk_H3_AIR(
double Re_mass_com) {
132 return 14.3 * pow(Re_mass_com, 0.225);
136 inline double hAk_C_AIR(
double Re_mass_com) {
137 return 2.6e-2 * pow(Re_mass_com, 0.942);
141 double Turb_Heat_Flow();
143 double Comp_Heat_Flow();
145 double Oil_Heat_Flow();
152 void InputData(
double T_AF,
double T_Humidity,
double T_MassFlow,
double T_IT_C,
double T_IP,
double T_PR,
153 double C_Humidity,
double C_MassFlow,
double C_IT_C,
double C_IP,
double C_PR,
double O_MassFlow,
double O_IT_C,
154 double O_IP,
double RTC);
158 void SolveNodeTemperatures(
double TET,
double TEC,
double TOIL);
160 void SolveHeatFlowMatix();
162 void SolveDeltaTemp();