62 #define _USE_MATH_DEFINES
68 const double Universal = 8314.4;
69 const double Air = 287.;
70 const double Fuel = 55.95;
71 const double Burnt = 285.4;
72 const double O2 = 259.825;
73 const double CO2 = 188.9207;
74 const double H2O = 461.398;
75 const double N2 = 296.837;
76 const double Ar = 208.12;
77 const double Diesel = 55.95;
78 const double Gasoline = 72.42;
85 const double O2 = 32.;
86 const double CO2 = 44.01;
87 const double H2O = 18.02;
88 const double N2 = 28.01;
89 const double Ar = 39.95;
90 const double C = 12.01;
91 const double NO2 = 46;
93 const double CO = 28.01;
94 const double UHC = 55.04;
95 const double Diesel = 148.4;
96 const double Gasoline = 114.8;
102 const double Pi = 3.14159265358979323846;
103 const double Pi_2 = Pi / 2;
104 const double Pi_4 = Pi / 4;
105 const double Pi_x_2 = 2 * Pi;
106 const double _1_Pi = 1 / Pi;
107 const double _2_Pi = 2 / Pi;
108 const double _4_Pi = 4 / Pi;
109 const double SQR_4_Pi = sqrt(_4_Pi);
110 const double Sigma = 5.670373e-8;
111 const double ARef = 343.11;
112 const double ARef2 = ARef * ARef;
113 const double TRef = 292.99271;
114 const double PRef = 1.0;
126 const double _DegToRad = 2 * __cons::Pi / 360;
128 inline double To_kilo(
double p) {
131 inline double From_kilo(
double p) {
135 inline double DegToRad(
double p) {
136 return p * _DegToRad;
139 inline double BarToPa(
double p) {
142 inline double PaToBar(
double p) {
146 inline double RPMToRPS(
double p) {
147 return p * 0.016666666666667;
149 inline double RPMToRad_s(
double p) {
150 return p * 0.104719755119660;
152 inline double Rad_sToRPM(
double p) {
153 return p * 9.549296585513721;
156 inline double m_sTokm_h(
double p) {
160 inline double degCToK(
double p) {
163 inline double KTodegC(
double p) {
173 namespace __HFormacion {
174 const double CO2 = -393510;
175 const double CO = -110530;
176 const double NO2 = 33100;
177 const double NO = 90290;
178 const double H2O = -241830;
183 inline double Circle_area(
double d) {
184 return d * d * __cons::Pi_4;
186 inline double Cylinder_volume(
double d,
double l) {
187 return Circle_area(d) * l;
189 inline double Ring_area(
double din,
double dout) {
190 return (dout * dout - din * din) * __cons::Pi_4;
198 const double Cp = 1004.5;
199 const double G = Cp / (Cp - __R::Air);
200 const double G_1 = G - 1;
201 const double G_2 = G + 1;
202 const double G_3 = (G - 1) / 2;
203 const double G_4 = 2 * G / (G - 1);
204 const double G_5 = (G - 1) / G / 2;
205 const double G_6 = 1 / (G - 1);
206 const double G_7 = (3 - G) / (G + 1);
207 const double G_8 = (G - 1) / G;
208 const double G_9 = G / (G - 1);
209 const double Cp_x2 = 2 * Cp;
210 const double Cv = Cp - __R::Air;
211 const double gxR = G * __R::Air;
213 inline double GG(
double Cp,
double Cv) {
217 inline double G1(
double g) {
222 inline double G2(
double g) {
227 inline double G3(
double g) {
228 return (g - 1) * 0.5;
232 inline double G4(
double g) {
233 return 2. * g / (g - 1);
237 inline double G5(
double g) {
238 return (g - 1) / 2. / g;
242 inline double G6(
double g) {
247 inline double G7(
double g) {
248 return (3 - g) / (g + 1);
252 inline double G8(
double g) {
257 inline double G9(
double g) {