Constantes.h

Go to the documentation of this file.

/*--------------------------------------------------------------------------------*\
==========================|
 \\   /\ /\   // O pen     | OpenWAM: The Open Source 1D Gas-Dynamic Code
 \\ |  X  | //  W ave     |
 \\ \/_\/ //   A ction   | CMT-Motores Termicos / Universidad Politecnica Valencia
 \\/   \//    M odel    |
 ----------------------------------------------------------------------------------
 License

 This file is part of OpenWAM.

 OpenWAM is free software: you can redistribute it and/or modify
 it under the terms of the GNU General Public License as published by
 the Free Software Foundation, either version 3 of the License, or
 (at your option) any later version.

 OpenWAM is distributed in the hope that it will be useful,
 but WITHOUT ANY WARRANTY; without even the implied warranty of
 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 GNU General Public License for more details.

 You should have received a copy of the GNU General Public License
 along with OpenWAM.  If not, see <http://www.gnu.org/licenses/>.


 \*--------------------------------------------------------------------------------*/

//---------------------------------------------------------------------------
#ifndef ConstantesH
#define ConstantesH
//#include "TBloqueMotor.h"
#ifdef __BORLANDC__
#include <vcl.h>
#endif

#include <limits>
#define _USE_MATH_DEFINES
#include <cmath>
//---------------------------------------------------------------------------

//Gas constant
namespace __R {
    const double Universal = 8314.4; // [J / (mol *K)]
    const double Air = 287.;
    const double Fuel = 55.95;
    const double Burnt = 285.4;
    const double O2 = 259.825;
    const double CO2 = 188.9207;
    const double H2O = 461.398;
    const double N2 = 296.837;
    const double Ar = 208.12;
    const double Diesel = 55.95;
    const double Gasoline = 72.42;
}
;
// [J / (kg * K)]

//Molecular weight
namespace __PM {
    const double O2 = 32.;
    const double CO2 = 44.01;
    const double H2O = 18.02;
    const double N2 = 28.01;
    const double Ar = 39.95;
    const double C = 12.01;
    const double NO2 = 46;
    const double NO = 30;
    const double CO = 28.01;
    const double UHC = 55.04;
    const double Diesel = 148.4;
    const double Gasoline = 114.8;
}
;
// [g]

namespace __cons {
    const double Pi = 3.14159265358979323846;
    const double Pi_2 = Pi / 2;
    const double Pi_4 = Pi / 4;
    const double Pi_x_2 = 2 * Pi;
    const double _1_Pi = 1 / Pi;
    const double _2_Pi = 2 / Pi;
    const double _4_Pi = 4 / Pi;
    const double SQR_4_Pi = sqrt(_4_Pi);
    const double Sigma = 5.670373e-8;
    const double ARef = 343.11;
    const double ARef2 = ARef * ARef;
    const double TRef = 292.99271;
    const double PRef = 1.0;
}
;

//const double Gamma = 1.4; ///< Air specific heat capacities ratio.
//const double ARef = 343.11; ///< Reference speed of sound. [m / s]
//const double TRef = 292.99271; ///< Reference temperature. [K]
//const double ARef2 = 117724.4721;
//const double PRef = 1.0; ///< Reference pressure. [bar]

namespace __units {

    const double _DegToRad = 2 * __cons::Pi / 360;
//GENERAL
    inline double To_kilo(double p) {
        return p * 0.001;
    }
    inline double From_kilo(double p) {
        return p * 1000;
    }
//ANGLE
    inline double DegToRad(double p) {
        return p * _DegToRad;
    }
// PRESSURE
    inline double BarToPa(double p) {
        return p * 1e5;
    }
    inline double PaToBar(double p) {
        return p * 1e-5;
    }
//ROTATIONAL SPEED
    inline double RPMToRPS(double p) {
        return p * 0.016666666666667;
    }
    inline double RPMToRad_s(double p) {
        return p * 0.104719755119660;
    }
    inline double Rad_sToRPM(double p) {
        return p * 9.549296585513721;
    }
//SPEED
    inline double m_sTokm_h(double p) {
        return p * 3.6;
    }
//TEMPERATURE
    inline double degCToK(double p) {
        return p + 273.15;
    }
    inline double KTodegC(double p) {
        return p - 273.15;
    }
}
;

//stUN __UN;

//const double Kb = 1.38054e-23; ///< Boltzmann constant. [J / K]

namespace __HFormacion {
    const double CO2 = -393510; 
    const double CO = -110530; 
    const double NO2 = 33100; 
    const double NO = 90290; 
    const double H2O = -241830; 
}
;

namespace __geom {
    inline double Circle_area(double d) {
        return d * d * __cons::Pi_4;
    }
    inline double Cylinder_volume(double d, double l) {
        return Circle_area(d) * l;
    }
    inline double Ring_area(double din, double dout) {
        return (dout * dout - din * din) * __cons::Pi_4;
    }
}
;

namespace __Gamma {

// FOR AIR AT AMBIENT CONDITIONS
    const double Cp = 1004.5;
    const double G = Cp / (Cp - __R::Air);
    const double G_1 = G - 1;
    const double G_2 = G + 1;
    const double G_3 = (G - 1) / 2;
    const double G_4 = 2 * G / (G - 1);
    const double G_5 = (G - 1) / G / 2;
    const double G_6 = 1 / (G - 1);
    const double G_7 = (3 - G) / (G + 1);
    const double G_8 = (G - 1) / G;
    const double G_9 = G / (G - 1);
    const double Cp_x2 = 2 * Cp;
    const double Cv = Cp - __R::Air;
    const double gxR = G * __R::Air;

    inline double GG(double Cp, double Cv) {
        return Cp / Cv;
    }

    inline double G1(double g) {
        return g - 1;
    }
    ;

    inline double G2(double g) {
        return g + 1;
    }
    ;

    inline double G3(double g) {
        return (g - 1) * 0.5;
    }
    ;

    inline double G4(double g) {
        return 2. * g / (g - 1);
    }
    ;

    inline double G5(double g) {
        return (g - 1) / 2. / g;
    }
    ;

    inline double G6(double g) {
        return 1 / (g - 1);
    }
    ;

    inline double G7(double g) {
        return (3 - g) / (g + 1);
    }
    ;

    inline double G8(double g) {
        return (g - 1) / g;
    }
    ;

    inline double G9(double g) {
        return g / (g - 1);
    }
}

#endif