This module aims to provide an open source and comprehensive approach to combustion problems and further isentropic expansion using the NASA polynomials To do it, this module is separated into 3 parts: the thermodynamic properties, the Algorithms and Other functionalities.
The original script is developed/maintained by C. Fuster, M. Soria, A. Miró, et all.
NASA polynomials are downloaded from Burcat mirrors link check his documentation.
There is no more documentation in this package that the comments on the functions, but you can check the master thesis developed with the Matlab code for more info. Request it to author of this package or download it from UPCcommons when the institution publishes it.
Thermodynamic properties are accessible using HGSsingle() or HGSprop() (see their own documentation using help())
help(fun)
HGSsingle
*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*
res = hgs_single(species, prop, T, P)
*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*
hgs_single returns the property of a species
*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*
Inputs:
-----------------------------------------------------------------------------
species --> String or numbers of species
prop --> Property requested (see below)
T --> [K] Temperature
P --> [bar] Pressure
*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*
Outputs:
-----------------------------------------------------------------------------
res --> Property result
mm [g/mol]
cp [kJ/(mol*K)]
cv [kJ/(mol*K)]
h [kJ/mol]
s [kJ/(mol*K)]
g [kJ/mol]
*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*
* Matlab HGS 2.0
* By Caleb Fuster, Manel Soria and Arnau Miró
* ESEIAAT UPC
HGSprop
*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*
var = hgs_prop(species, n, T, P, *args)
*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*
hgs_prop returns the properties of the mixture of gasses
*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*
Inputs:
-----------------------------------------------------------------------------
species --> String or numbers of species
prop --> Property requested (see below)
T --> [K] Temperature
P --> [bar] Pressure
*args --> Expected return: 'Mm' 'Cp' 'Cv' 'H' 'S' 'G' 'Rg' 'gamma' 'a'
If it is empty, all the properties will be
return
*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*
Outputs:
-----------------------------------------------------------------------------
var --> Property result
mm [g/mol]
cp [kJ/K]
cv [kJ/K]
h [kJ]
s [kJ/K]
g [kJ]
Rg [kJ/(kg*K)]
gamma
a [m/s]
*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*
* Matlab HGS 2.0
* By Caleb Fuster, Manel Soria and Arnau Miró
* ESEIAAT UPC
This module includes an equilibrium algorithm (HGSeq()), combustion algorithm (HGStp()) and isentropic expansion algorithm (HGSisentropic()) (see their own documentation using help())
help(fun)
HGSeq
*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*
species, n, Gmin = hgs_eq(species, n0, T, P, **kwargs)
*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*
hgs_eq calculates the species mols equilibrium at a certain temperature
and pressure
*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*
Inputs:
-----------------------------------------------------------------------------
species --> String or numbers of species
n0 --> [mol] Initial mixture
T --> [K] Temperature. Could be a single value or an array.
P --> [bar] Pressure
**kwargs --> opti_eq= Options for the minimize Scipy function
*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*
Outputs:
-----------------------------------------------------------------------------
species --> Species
n --> [mol] Final mixture
Gmin --> [kJ] Minimum Gibbs free energy
*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*
* Matlab HGS 2.0
* By Caleb Fuster, Manel Soria and Arnau Miró
* ESEIAAT UPC
HGStp
*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*
Tp, n, species, flag = hgs_tp(species, n0, tipo, V0, P, **kwargs)
*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*
hgs_tp calculates the reaction temperature considering dissociation,
and the products composition in equilibrium
*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*
Inputs:
-----------------------------------------------------------------------------
species --> String or numbers of species
n0 --> [mols] Number of mols of each species
tipo --> Entry type that defines the state of the input.
It can be 'T' or 'H'
V0 --> Entry that should be for type:'T' V0=T [K] input temperature
'H' V0=H [kJ] input enthalpy
P --> [bar] Mixture pressure
**kwargs --> opti_eq= Options for the minimize Scipy function
opt_sec= Dictionary with the options for the secant method.
"xmin": [K] Temperature minimum for the solver;
"xmax" [K] Temperature maximum for the solver;
"maxiter" Max iterations for the solver;
"epsx" Diferential T where the solver reachs the solution;
"epsy" Diferential S where the solver reachs the solution;
"fchange" T difference where secant method is
changed by bisection method;
"tipo" Select between: 'Frozen' for frozen flow
'Shifting' for shifting flow
"info" Detailed info == 1; No info == 0.
"dTp" Improve the velocity with the approximation of
parabola. +- dTp
opt_sec = {"xmin": 300, "xmax": 4000, "maxiter": 200,
"epsx": 0.1, "epsy": 1, "tipo": "Shifting",
"fchange": 5, "info": 0, "dTp": 100}
*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*
Outputs:
-----------------------------------------------------------------------------
Tp --> [K] Final temperature
n --> [mol] Final mixture
species --> String or numbers of species
flag --> Solver error detection:
1 Solver has reached the solution
-1 Solver failed. Maximum iterations
-2 Solver failed. Initial sign change not found
*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*
* Matlab HGS 2.0
* By Caleb Fuster, Manel Soria and Arnau Miró
* ESEIAAT UPC
HGSisentropic
*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*
Tp, n, species, v2, M2, flag = hgs_isentropic(species, n0, T0, P0, P1, **kwargs)
*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*
hgs_tp calculates the outlet variables for an isentropic expansion
*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*
Inputs:
-------------------------------------------------------------------------------
species --> String or numbers of species
n0 --> [mols] Number of mols of each species
T0 --> [K] Initial temperature
P0 --> [bar] Inlet pressure
P1 --> [bar] Exit pressure
**kwargs --> opti_eq= Options for the minimize Scipy function
opt_sec= Dictionary with the options for the secant method.
"xmin": [K] Temperature minimum for the solver;
"xmax" [K] Temperature maximum for the solver;
"maxiter" Max iterations for the solver;
"epsx" Diferential T where the solver reachs the solution;
"epsy" Diferential S where the solver reachs the solution;
"fchange" T difference where secant method is
changed by bisection method;
"tipo" Select between: 'Frozen' for frozen flow
'Shifting' for shifting flow
"info" Detailed info == 1; No info == 0.
"dTp" Improve the velocity with the approximation of
parabola. +- dTp
opt_sec = {"xmin": 300, "xmax": 4000, "maxiter": 200,
"epsx": 0.1, "epsy": 1, "tipo": "Shifting",
"fchange": 5, "info": 0, "dTp": 100}
*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*
Outputs:
-----------------------------------------------------------------------------
Tp --> [K] Exit temperature
n --> [mols] Species resultant mols
species --> String or numbers of species
v2 --> [m/s] Velocity of the mixture
M2 --> [Mach] Mach of the mixture
flag --> Solver error detection:
1 Solver has reached the solution
-1 Solver failed. Maximum iterations
-2 Solver failed. Initial sign change not found
*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*+*
* Matlab HGS 2.0
* By Caleb Fuster, Manel Soria and Arnau Miró
* ESEIAAT UPC
This module includes other functionalities to facilitate the user life. You can:
Update the database:
HGSresetDataBase.mCreate/Erase your own mixtures:
HGSaddMixture("NAME", ["species1", "species2",...], [20, 80,...])
HGSsubtMixture("NAME")And search and print info from a species:
HGSfind("NAME")
HGSprintInfo("NAME")