diff --git a/test_problems/ChemEquil_ionizedGas/.cvsignore b/test_problems/ChemEquil_ionizedGas/.cvsignore new file mode 100644 index 000000000..63d34f59c --- /dev/null +++ b/test_problems/ChemEquil_ionizedGas/.cvsignore @@ -0,0 +1,7 @@ +ionizedGasEquil +Makefile +desc.txt +output.txt +csvCode.txt +ct2ctml.log +diff_test.out diff --git a/test_problems/ChemEquil_ionizedGas/Makefile.in b/test_problems/ChemEquil_ionizedGas/Makefile.in new file mode 100644 index 000000000..b8d77d8b9 --- /dev/null +++ b/test_problems/ChemEquil_ionizedGas/Makefile.in @@ -0,0 +1,110 @@ +#!/bin/sh + +############################################################################ +# +# Makefile to compile and link a C++ application to +# Cantera. +# +############################################################################# + +# addition to suffixes +.SUFFIXES : .d + +# the name of the executable program to be created +PROG_NAME = ionizedGasEquil + +# the object files to be linked together. List those generated from Fortran +# and from C/C++ separately +OBJS = ionizedGasEquil.o + +# Location of the current build. Will assume that tests are run +# in the source directory tree location +src_dir_tree = 0 + +# additional flags to be passed to the linker. If your program +# requires other external libraries, put them here +LINK_OPTIONS = @EXTRA_LINK@ + +############################################################################# + +# Check to see whether we are in the msvc++ environment +os_is_win = @OS_IS_WIN@ + +# Fortran libraries +FORT_LIBS = @FLIBS@ + +# the C++ compiler +CXX = @CXX@ + +# C++ compile flags +ifeq ($(src_dir_tree), 1) +CXX_FLAGS = -DSRCDIRTREE @CXXFLAGS@ +else +CXX_FLAGS = @CXXFLAGS@ +endif + +# Ending C++ linking libraries +LCXX_END_LIBS = @LCXX_END_LIBS@ + +# the directory where the Cantera libraries are located +CANTERA_LIBDIR=@buildlib@ + +# required Cantera libraries +CANTERA_LIBS = @LOCAL_LIBS@ -lctcxx + +# the directory where Cantera include files may be found. +ifeq ($(src_dir_tree), 1) +CANTERA_INCDIR=../../Cantera/src +else +CANTERA_INCDIR=@ctroot@/build/include/cantera +endif + +# flags passed to the C++ compiler/linker for the linking step +LCXX_FLAGS = -L$(CANTERA_LIBDIR) @LOCAL_LIB_DIRS@ @CXXFLAGS@ + +# How to compile C++ source files to object files +.@CXX_EXT@.@OBJ_EXT@: + $(CXX) -c $< -I$(CANTERA_INCDIR) $(CXX_FLAGS) + +# How to compile the dependency file +.cpp.d: + @CXX_DEPENDS@ -I$(CANTERA_INCDIR) $(CXX_FLAGS) $*.cpp > $*.d + +# List of dependency files to be created +DEPENDS=$(OBJS:.o=.d) + +# Program Name +PROGRAM = $(PROG_NAME)$(EXE_EXT) + +all: $(PROGRAM) + +$(PROGRAM): $(OBJS) $(CANTERA_LIBDIR)/libcantera.a + $(CXX) -o $(PROGRAM) $(OBJS) $(LCXX_FLAGS) $(LINK_OPTIONS) \ + $(CANTERA_LIBS) @LIBS@ $(FORT_LIBS) \ + $(LCXX_END_LIBS) + +# depends target -> forces recalculation of dependencies +depends: + $(RM) *.d .depends + @MAKE@ .depends + +.depends: $(DEPENDS) + cat *.d > .depends + +# Do the test -> For the windows vc++ environment, we have to skip checking on +# whether the program is uptodate, because we don't utilize make +# in that environment to build programs. +test: +ifeq ($(os_is_win), 1) +else + @ @MAKE@ -s $(PROGRAM) +endif + @ ./runtest + +clean: + $(RM) $(OBJS) $(PROGRAM) $(DEPENDS) .depends + ../../bin/rm_cvsignore + (if test -d SunWS_cache ; then \ + $(RM) -rf SunWS_cache ; \ + fi ) + diff --git a/test_problems/ChemEquil_ionizedGas/air_below6000K.cti b/test_problems/ChemEquil_ionizedGas/air_below6000K.cti new file mode 100755 index 000000000..b43318f0a --- /dev/null +++ b/test_problems/ChemEquil_ionizedGas/air_below6000K.cti @@ -0,0 +1,165 @@ +# +# This is the NASA thermodynamic database, which is available for download +# from http://www.galcit.caltech.edu/EDL/public/thermo.html. The original +# sources are +# +# S. Gordon and B.J. McBride, "Computer Program for Calculation of Complex +# Chemical Equilibrium Composition, Rocket Performance, Incident and +# Reflected Shocks and Chapman-Jouguet Detonations", NASA Report SP-273, 1971. +# +# B.J. McBride, S. Gordon, and M.A. Reno, "Coefficients for Calculating +# Thermodynamic and Transport Properties of Individual Species", NASA +# Report TM-4513, October 1993. +# +#----------------------------------------------------------------------------- + +ideal_gas(name = "air_below6000K", + elements = " N O E ", + species = """ N2 O2 N O NO NO+ Electron N+ O+ N2+ O2+ """, + reactions = "all", + initial_state = state(temperature = 298.15, + pressure = OneAtm, + mole_fractions = 'N2:0.78847, O2:0.21153') ) + +species(name = "Electron", + atoms = " E:1 ", + thermo = ( + NASA( [ 200.00, 1000.00], [ 2.500000000E+00, 0.000000000E+00, + 0.000000000E+00, 0.000000000E+00, 0.000000000E+00, + -7.453750000E+02, -1.172469020E+01] ), + NASA( [ 1000.00, 6000.00], [ 2.500000000E+00, 0.000000000E+00, + 0.000000000E+00, 0.000000000E+00, 0.000000000E+00, + -7.453750000E+02, -1.172469020E+01] ) + ), + note = "gas L10/92" + ) + +species(name = "N", + atoms = " N:1 ", + thermo = ( + NASA( [ 200.00, 1000.00], [ 2.500000000E+00, 0.000000000E+00, + 0.000000000E+00, 0.000000000E+00, 0.000000000E+00, + 5.610463780E+04, 4.193909320E+00] ), + NASA( [ 1000.00, 6000.00], [ 2.415942930E+00, 1.748906000E-04, + -1.190236670E-07, 3.022623870E-11, -2.036097900E-15, + 5.613377480E+04, 4.649609860E+00] ) + ), + note = "L 6/88" + ) + +species(name = "N+", + atoms = " N:1 E:-1 ", + thermo = ( + NASA( [ 200.00, 1000.00], [ 2.802694450E+00, -1.447589110E-03, + 2.771183800E-06, -2.401873520E-09, 7.808399310E-13, + 2.255752440E+05, 3.578778350E+00] ), + NASA( [ 1000.00, 6000.00], [ 2.511129670E+00, 3.464417510E-06, + -1.594269380E-08, 7.248656630E-12, -6.445014260E-16, + 2.256243400E+05, 4.927676610E+00] ) + ), + note = "L 7/88" + ) + +species(name = "NO", + atoms = " N:1 O:1 ", + thermo = ( + NASA( [ 200.00, 1000.00], [ 4.218598960E+00, -4.639881240E-03, + 1.104430490E-05, -9.340555070E-09, 2.805548740E-12, + 9.845099640E+03, 2.280610010E+00] ), + NASA( [ 1000.00, 6000.00], [ 3.260712340E+00, 1.191011350E-03, + -4.291226460E-07, 6.944814630E-11, -4.032956810E-15, + 9.921431320E+03, 6.369005180E+00] ) + ), + note = "TPIS89" + ) + +species(name = "NO+", + atoms = " N:1 O:1 E:-1 ", + thermo = ( + NASA( [ 200.00, 1000.00], [ 3.693012310E+00, -1.342291580E-03, + 2.673433950E-06, -1.026093080E-09, -6.956104920E-14, + 1.181030550E+05, 3.091266980E+00] ), + NASA( [ 1000.00, 6000.00], [ 2.945877020E+00, 1.403252600E-03, + -4.955031960E-07, 7.959489730E-11, -4.720766680E-15, + 1.182443400E+05, 6.706446410E+00] ) + ), + note = "TPIS89" + ) + +species(name = "N2", + atoms = " N:2 ", + thermo = ( + NASA( [ 200.00, 1000.00], [ 3.531005280E+00, -1.236609870E-04, + -5.029994370E-07, 2.435306120E-09, -1.408812350E-12, + -1.046976280E+03, 2.967474680E+00] ), + NASA( [ 1000.00, 6000.00], [ 2.952576260E+00, 1.396900570E-03, + -4.926316910E-07, 7.860103670E-11, -4.607553210E-15, + -9.239486450E+02, 5.871892520E+00] ) + ), + note = "TPIS78" + ) + +species(name = "N2+", + atoms = " N:2 E:-1 ", + thermo = ( + NASA( [ 200.00, 1000.00], [ 3.775407110E+00, -2.064591570E-03, + 4.757523010E-06, -3.156642280E-09, 6.705099730E-13, + 1.804811150E+05, 2.693221860E+00] ), + NASA( [ 1000.00, 6000.00], [ 3.586613630E+00, 2.530719490E-04, + 1.847782140E-07, -4.552572230E-11, 3.268180290E-15, + 1.803909940E+05, 3.095841500E+00] ) + ), + note = "TPIS89" + ) + +species(name = "O", + atoms = " O:1 ", + thermo = ( + NASA( [ 200.00, 1000.00], [ 3.168267100E+00, -3.279318840E-03, + 6.643063960E-06, -6.128066240E-09, 2.112659710E-12, + 2.912225920E+04, 2.051933460E+00] ), + NASA( [ 1000.00, 6000.00], [ 2.543636970E+00, -2.731624860E-05, + -4.190295200E-09, 4.954818450E-12, -4.795536940E-16, + 2.922601200E+04, 4.922294570E+00] ) + ), + note = "L 1/90" + ) + +species(name = "O+", + atoms = " O:1 E:-1 ", + thermo = ( + NASA( [ 200.00, 1000.00], [ 2.500000000E+00, 0.000000000E+00, + 0.000000000E+00, 0.000000000E+00, 0.000000000E+00, + 1.879352910E+05, 4.393376890E+00] ), + NASA( [ 1000.00, 6000.00], [ 2.487733170E+00, 2.176600160E-05, + -1.089558060E-08, 1.259092120E-12, 1.373167200E-16, + 1.879399650E+05, 4.461340910E+00] ) + ), + note = "L 1/90" + ) + +species(name = "O2", + atoms = " O:2 ", + thermo = ( + NASA( [ 200.00, 1000.00], [ 3.782456360E+00, -2.996734150E-03, + 9.847302000E-06, -9.681295080E-09, 3.243728360E-12, + -1.063943560E+03, 3.657675730E+00] ), + NASA( [ 1000.00, 6000.00], [ 3.660960830E+00, 6.563655230E-04, + -1.411494850E-07, 2.057976580E-11, -1.299132480E-15, + -1.215977250E+03, 3.415361840E+00] ) + ), + note = "TPIS89" + ) + +species(name = "O2+", + atoms = " O:2 E:-1 ", + thermo = ( + NASA( [ 200.00, 1000.00], [ 4.610171670E+00, -6.359519520E-03, + 1.424256240E-05, -1.209979230E-08, 3.709568780E-12, + 1.397422290E+05, -2.013268740E-01] ), + NASA( [ 1000.00, 6000.00], [ 3.316759220E+00, 1.115222440E-03, + -3.834925560E-07, 5.727846870E-11, -2.776483810E-15, + 1.398768230E+05, 5.447264760E+00] ) + ), + note = "TPIS89" + ) diff --git a/test_problems/ChemEquil_ionizedGas/air_below6000K.xml b/test_problems/ChemEquil_ionizedGas/air_below6000K.xml new file mode 100644 index 000000000..efb0bb235 --- /dev/null +++ b/test_problems/ChemEquil_ionizedGas/air_below6000K.xml @@ -0,0 +1,228 @@ + + + + + + + N O E + N2 O2 N O NO NO+ Electron N+ O+ N2+ O2+ + + + 298.14999999999998 + 101325.0 + N2:0.78847, O2:0.21153 + + + + + + + + + + + + E:1 + gas L10/92 + -1 + + + + 2.500000000E+00, 0.000000000E+00, 0.000000000E+00, 0.000000000E+00, + 0.000000000E+00, -7.453750000E+02, -1.172469020E+01 + + + + 2.500000000E+00, 0.000000000E+00, 0.000000000E+00, 0.000000000E+00, + 0.000000000E+00, -7.453750000E+02, -1.172469020E+01 + + + + + + + N:1 + L 6/88 + + + + 2.500000000E+00, 0.000000000E+00, 0.000000000E+00, 0.000000000E+00, + 0.000000000E+00, 5.610463780E+04, 4.193909320E+00 + + + + 2.415942930E+00, 1.748906000E-04, -1.190236670E-07, 3.022623870E-11, + -2.036097900E-15, 5.613377480E+04, 4.649609860E+00 + + + + + + + E:-1 N:1 + L 7/88 + 1 + + + + 2.802694450E+00, -1.447589110E-03, 2.771183800E-06, -2.401873520E-09, + 7.808399310E-13, 2.255752440E+05, 3.578778350E+00 + + + + 2.511129670E+00, 3.464417510E-06, -1.594269380E-08, 7.248656630E-12, + -6.445014260E-16, 2.256243400E+05, 4.927676610E+00 + + + + + + + O:1 N:1 + TPIS89 + + + + 4.218598960E+00, -4.639881240E-03, 1.104430490E-05, -9.340555070E-09, + 2.805548740E-12, 9.845099640E+03, 2.280610010E+00 + + + + 3.260712340E+00, 1.191011350E-03, -4.291226460E-07, 6.944814630E-11, + -4.032956810E-15, 9.921431320E+03, 6.369005180E+00 + + + + + + + E:-1 O:1 N:1 + TPIS89 + 1 + + + + 3.693012310E+00, -1.342291580E-03, 2.673433950E-06, -1.026093080E-09, + -6.956104920E-14, 1.181030550E+05, 3.091266980E+00 + + + + 2.945877020E+00, 1.403252600E-03, -4.955031960E-07, 7.959489730E-11, + -4.720766680E-15, 1.182443400E+05, 6.706446410E+00 + + + + + + + N:2 + TPIS78 + + + + 3.531005280E+00, -1.236609870E-04, -5.029994370E-07, 2.435306120E-09, + -1.408812350E-12, -1.046976280E+03, 2.967474680E+00 + + + + 2.952576260E+00, 1.396900570E-03, -4.926316910E-07, 7.860103670E-11, + -4.607553210E-15, -9.239486450E+02, 5.871892520E+00 + + + + + + + E:-1 N:2 + TPIS89 + 1 + + + + 3.775407110E+00, -2.064591570E-03, 4.757523010E-06, -3.156642280E-09, + 6.705099730E-13, 1.804811150E+05, 2.693221860E+00 + + + + 3.586613630E+00, 2.530719490E-04, 1.847782140E-07, -4.552572230E-11, + 3.268180290E-15, 1.803909940E+05, 3.095841500E+00 + + + + + + + O:1 + L 1/90 + + + + 3.168267100E+00, -3.279318840E-03, 6.643063960E-06, -6.128066240E-09, + 2.112659710E-12, 2.912225920E+04, 2.051933460E+00 + + + + 2.543636970E+00, -2.731624860E-05, -4.190295200E-09, 4.954818450E-12, + -4.795536940E-16, 2.922601200E+04, 4.922294570E+00 + + + + + + + E:-1 O:1 + L 1/90 + 1 + + + + 2.500000000E+00, 0.000000000E+00, 0.000000000E+00, 0.000000000E+00, + 0.000000000E+00, 1.879352910E+05, 4.393376890E+00 + + + + 2.487733170E+00, 2.176600160E-05, -1.089558060E-08, 1.259092120E-12, + 1.373167200E-16, 1.879399650E+05, 4.461340910E+00 + + + + + + + O:2 + TPIS89 + + + + 3.782456360E+00, -2.996734150E-03, 9.847302000E-06, -9.681295080E-09, + 3.243728360E-12, -1.063943560E+03, 3.657675730E+00 + + + + 3.660960830E+00, 6.563655230E-04, -1.411494850E-07, 2.057976580E-11, + -1.299132480E-15, -1.215977250E+03, 3.415361840E+00 + + + + + + + E:-1 O:2 + TPIS89 + 1 + + + + 4.610171670E+00, -6.359519520E-03, 1.424256240E-05, -1.209979230E-08, + 3.709568780E-12, 1.397422290E+05, -2.013268740E-01 + + + + 3.316759220E+00, 1.115222440E-03, -3.834925560E-07, 5.727846870E-11, + -2.776483810E-15, 1.398768230E+05, 5.447264760E+00 + + + + + + diff --git a/test_problems/ChemEquil_ionizedGas/ionizedGasEquil.cpp b/test_problems/ChemEquil_ionizedGas/ionizedGasEquil.cpp new file mode 100644 index 000000000..3dfa90b5c --- /dev/null +++ b/test_problems/ChemEquil_ionizedGas/ionizedGasEquil.cpp @@ -0,0 +1,78 @@ +/* + * $Author$ + * $Date$ + * $Revision$ + * + * + */ + +#ifdef SRCDIRTREE +#include "ct_defs.h" +#include "ThermoPhase.h" +#include "IdealGasMix.h" +#include "equil.h" +#else +#include "Cantera.h" +#include "IdealGasMix.h" +#include "equilibrium.h" +#endif + + +using namespace std; +using namespace Cantera; +#ifdef DEBUG_HKM_EPEQUIL +extern int debug_prnt_lvl; +#endif +int main(int argc, char **argv) { +#ifdef DEBUG_HKM_EPEQUIL + debug_prnt_lvl = 0; +#endif + try { + + IdealGasPhase* gas = new IdealGasMix("air_below6000K.xml","air_below6000K"); + + vector_fp IndVar2(6, 0.0); + IndVar2[0] = 1.5E5; + IndVar2[1] = 3.0E5; + IndVar2[2] = 9.0E5; + IndVar2[3] = 2.7E6; + IndVar2[4] = 6.7E6; + IndVar2[5] = 1.0E7; + + vector_fp IndVar1(7, 0.0); + IndVar1[0] = 1.0E-8; + IndVar1[1] = 1.0E-7; + IndVar1[2] = 1.0E-6; + IndVar1[3] = 1.0E-5; + IndVar1[4] = 1.0E-4; + IndVar1[5] = 1.0E-3; + IndVar1[6] = 1.0E-2; + int nj = 6; + int ni = 7; + + for (int j=0; jsetState_UV(IndVar2[j]+offset,1.0/IndVar1[i]); + double tkelvin = gas->temperature(); + double pres = gas->pressure(); + printf("Initial T = %g, pres = %g atm \n", tkelvin, pres/OneAtm); + beginLogGroup("topEquil", -1); + equilibrate(*gas,"UV", -1); + endLogGroup("topEquil"); + cout << report(*gas) << endl; + + tkelvin = gas->temperature(); + pres = gas->pressure(); + printf("Final T = %g, pres = %g atm\n", tkelvin, pres/OneAtm); + + } + } + delete gas; + } + + catch (CanteraError) { + showErrors(); + } + return 0; +} diff --git a/test_problems/ChemEquil_ionizedGas/output_blessed.txt b/test_problems/ChemEquil_ionizedGas/output_blessed.txt new file mode 100644 index 000000000..f370d2cb3 --- /dev/null +++ b/test_problems/ChemEquil_ionizedGas/output_blessed.txt @@ -0,0 +1,1386 @@ +Initial T = 207.284, pres = 5.8944e-09 atm + + air_below6000K: + + temperature 207.284 K + pressure 0.00059725 Pa + density 1e-08 kg/m^3 + mean mol. weight 28.8565 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy -91746.4 -2.647e+06 J + internal energy -151471 -4.371e+06 J + entropy 11980.9 3.457e+05 J/K + Gibbs function -2.57518e+06 -7.431e+07 J + heat capacity c_p 1009.36 2.913e+04 J/K + heat capacity c_v 721.23 2.081e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.78847 0.765436 -42.4943 + O2 0.21153 0.234564 -45.4397 + N 9.55809e-113 4.63943e-113 -21.2471 + O 1.18094e-56 6.54771e-57 -22.7198 + NO 1.82556e-23 1.89829e-23 -43.967 + NO+ 1.26211e-226 1.31237e-226 11.6883 + Electron 1.26211e-226 2.3837e-231 -541.807 + N+ 0 0 + O+ 0 0 + N2+ 0 0 + O2+ 0 0 + +Final T = 207.284, pres = 5.8944e-09 atm +Initial T = 207.284, pres = 5.8944e-08 atm + + air_below6000K: + + temperature 207.284 K + pressure 0.0059725 Pa + density 1e-07 kg/m^3 + mean mol. weight 28.8565 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy -91746.4 -2.647e+06 J + internal energy -151471 -4.371e+06 J + entropy 11317.4 3.266e+05 J/K + Gibbs function -2.43766e+06 -7.034e+07 J + heat capacity c_p 1009.36 2.913e+04 J/K + heat capacity c_v 721.23 2.081e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.78847 0.765436 -40.1917 + O2 0.21153 0.234564 -43.1371 + N 3.02253e-113 1.46712e-113 -20.0958 + O 3.73447e-57 2.07057e-57 -21.5685 + NO 1.82556e-23 1.89829e-23 -41.6644 + NO+ 1.66405e-227 1.73032e-227 11.9648 + Electron 1.66405e-227 3.14282e-232 -541.531 + N+ 0 0 + O+ 0 0 + N2+ 0 0 + O2+ 0 0 + +Final T = 207.284, pres = 5.8944e-08 atm +Initial T = 207.284, pres = 5.8944e-07 atm + + air_below6000K: + + temperature 207.284 K + pressure 0.059725 Pa + density 1e-06 kg/m^3 + mean mol. weight 28.8565 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy -91746.4 -2.647e+06 J + internal energy -151471 -4.371e+06 J + entropy 10654 3.074e+05 J/K + Gibbs function -2.30014e+06 -6.637e+07 J + heat capacity c_p 1009.36 2.913e+04 J/K + heat capacity c_v 721.23 2.081e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.78847 0.765436 -37.8891 + O2 0.21153 0.234564 -40.8345 + N 9.55809e-114 4.63943e-114 -18.9446 + O 1.18094e-57 6.54771e-58 -20.4172 + NO 1.82556e-23 1.89829e-23 -39.3618 + NO+ 2.2952e-228 2.3866e-228 12.2863 + Electron 2.2952e-228 4.33484e-233 -541.209 + N+ 0 0 + O+ 0 0 + N2+ 0 0 + O2+ 0 0 + +Final T = 207.284, pres = 5.8944e-07 atm +Initial T = 207.284, pres = 5.8944e-06 atm + + air_below6000K: + + temperature 207.284 K + pressure 0.59725 Pa + density 1e-05 kg/m^3 + mean mol. weight 28.8565 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy -91746.4 -2.647e+06 J + internal energy -151471 -4.371e+06 J + entropy 9990.51 2.883e+05 J/K + Gibbs function -2.16262e+06 -6.241e+07 J + heat capacity c_p 1009.36 2.913e+04 J/K + heat capacity c_v 721.23 2.081e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.78847 0.765436 -35.5865 + O2 0.21153 0.234564 -38.5319 + N 3.02253e-114 1.46712e-114 -17.7933 + O 3.73447e-58 2.07057e-58 -19.266 + NO 1.82556e-23 1.89829e-23 -37.0592 + NO+ 1.25826e-229 1.30837e-229 11.6852 + Electron 1.25826e-229 2.37642e-234 -541.81 + N+ 0 0 + O+ 0 0 + N2+ 0 0 + O2+ 0 0 + +Final T = 207.284, pres = 5.8944e-06 atm +Initial T = 207.284, pres = 5.8944e-05 atm + + air_below6000K: + + temperature 207.284 K + pressure 5.9725 Pa + density 1e-04 kg/m^3 + mean mol. weight 28.8565 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy -91746.4 -2.647e+06 J + internal energy -151471 -4.371e+06 J + entropy 9327.06 2.691e+05 J/K + Gibbs function -2.02509e+06 -5.844e+07 J + heat capacity c_p 1009.36 2.913e+04 J/K + heat capacity c_v 721.23 2.081e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.78847 0.765436 -33.2839 + O2 0.21153 0.234564 -36.2293 + N 9.55809e-115 4.63943e-115 -16.642 + O 1.18094e-58 6.54771e-59 -18.1147 + NO 1.82556e-23 1.89829e-23 -34.7566 + NO+ 2.60402e-230 2.70772e-230 12.4126 + Electron 2.60402e-230 4.9181e-235 -541.083 + N+ 0 0 + O+ 0 0 + N2+ 0 0 + O2+ 0 0 + +Final T = 207.284, pres = 5.8944e-05 atm +Initial T = 207.284, pres = 0.00058944 atm + + air_below6000K: + + temperature 207.284 K + pressure 59.725 Pa + density 0.001 kg/m^3 + mean mol. weight 28.8565 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy -91746.4 -2.647e+06 J + internal energy -151471 -4.371e+06 J + entropy 8663.62 2.5e+05 J/K + Gibbs function -1.88757e+06 -5.447e+07 J + heat capacity c_p 1009.36 2.913e+04 J/K + heat capacity c_v 721.23 2.081e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.78847 0.765436 -30.9814 + O2 0.21153 0.234564 -33.9267 + N 3.02253e-115 1.46712e-115 -15.4907 + O 3.73447e-59 2.07057e-59 -16.9634 + NO 1.82556e-23 1.89829e-23 -32.454 + NO+ 1.26648e-231 1.31691e-231 11.6917 + Electron 1.26648e-231 2.39194e-236 -541.804 + N+ 0 0 + O+ 0 0 + N2+ 0 0 + O2+ 0 0 + +Final T = 207.284, pres = 0.00058944 atm +Initial T = 207.284, pres = 0.0058944 atm + + air_below6000K: + + temperature 207.284 K + pressure 597.25 Pa + density 0.01 kg/m^3 + mean mol. weight 28.8565 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy -91746.4 -2.647e+06 J + internal energy -151471 -4.371e+06 J + entropy 8000.17 2.309e+05 J/K + Gibbs function -1.75005e+06 -5.05e+07 J + heat capacity c_p 1009.36 2.913e+04 J/K + heat capacity c_v 721.23 2.081e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.78847 0.765436 -28.6788 + O2 0.21153 0.234564 -31.6242 + N 9.55809e-116 4.63943e-116 -14.3394 + O 1.18094e-59 6.54771e-60 -15.8121 + NO 1.82556e-23 1.89829e-23 -30.1515 + NO+ 2.40463e-232 2.50039e-232 12.3329 + Electron 2.40463e-232 4.54152e-237 -541.163 + N+ 0 0 + O+ 0 0 + N2+ 0 0 + O2+ 0 0 + +Final T = 207.284, pres = 0.0058944 atm +Initial T = 414.143, pres = 1.17767e-08 atm + + air_below6000K: + + temperature 414.143 K + pressure 0.00119328 Pa + density 1e-08 kg/m^3 + mean mol. weight 28.8565 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 117856 3.401e+06 J + internal energy -1471.39 -4.246e+04 J + entropy 12482.2 3.602e+05 J/K + Gibbs function -5.05158e+06 -1.458e+08 J + heat capacity c_p 1022.57 2.951e+04 J/K + heat capacity c_v 734.438 2.119e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.78847 0.765436 -41.7103 + O2 0.21153 0.234564 -44.6571 + N 2.06855e-53 1.00406e-53 -20.8551 + O 1.8737e-25 1.03887e-25 -22.3285 + NO 5.58575e-12 5.80828e-12 -43.1837 + NO+ 4.00599e-102 4.16552e-102 11.996 + Electron 4.00599e-102 7.56595e-107 -254.374 + N+ 7.16723e-207 3.47879e-207 34.3246 + O+ 3.0296e-168 1.6797e-168 32.8512 + N2+ 6.19808e-167 6.01689e-167 13.4694 + O2+ 3.22979e-125 3.58144e-125 10.5226 + +Final T = 414.143, pres = 1.17767e-08 atm +Initial T = 414.143, pres = 1.17767e-07 atm + + air_below6000K: + + temperature 414.143 K + pressure 0.0119328 Pa + density 1e-07 kg/m^3 + mean mol. weight 28.8565 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 117856 3.401e+06 J + internal energy -1471.39 -4.246e+04 J + entropy 11818.8 3.41e+05 J/K + Gibbs function -4.77681e+06 -1.378e+08 J + heat capacity c_p 1022.57 2.951e+04 J/K + heat capacity c_v 734.438 2.119e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.78847 0.765436 -39.4077 + O2 0.21153 0.234564 -42.3545 + N 6.54134e-54 3.17512e-54 -19.7038 + O 5.92517e-26 3.28519e-26 -21.1772 + NO 5.58575e-12 5.80828e-12 -40.8811 + NO+ 5.13663e-103 5.34118e-103 12.2446 + Electron 5.13663e-103 9.70134e-108 -254.125 + N+ 2.90616e-208 1.41058e-208 33.4219 + O+ 1.22844e-169 6.81082e-170 31.9485 + N2+ 7.9474e-168 7.71508e-168 13.718 + O2+ 4.14136e-126 4.59225e-126 10.7712 + +Final T = 414.143, pres = 1.17767e-07 atm +Initial T = 414.143, pres = 1.17767e-06 atm + + air_below6000K: + + temperature 414.143 K + pressure 0.119328 Pa + density 1e-06 kg/m^3 + mean mol. weight 28.8565 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 117856 3.401e+06 J + internal energy -1471.39 -4.246e+04 J + entropy 11155.3 3.219e+05 J/K + Gibbs function -4.50205e+06 -1.299e+08 J + heat capacity c_p 1022.57 2.951e+04 J/K + heat capacity c_v 734.438 2.119e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.78847 0.765436 -37.1051 + O2 0.21153 0.234564 -40.0519 + N 2.06855e-54 1.00406e-54 -18.5525 + O 1.8737e-26 1.03887e-26 -20.026 + NO 5.58575e-12 5.80828e-12 -38.5785 + NO+ 2.82555e-104 2.93806e-104 11.6469 + Electron 2.82555e-104 5.33648e-109 -254.723 + N+ 5.05526e-210 2.45369e-210 31.6729 + O+ 2.13687e-171 1.18474e-171 30.1995 + N2+ 4.37168e-169 4.24389e-169 13.1203 + O2+ 2.27807e-127 2.52609e-127 10.1735 + +Final T = 414.143, pres = 1.17767e-06 atm +Initial T = 414.143, pres = 1.17767e-05 atm + + air_below6000K: + + temperature 414.143 K + pressure 1.19328 Pa + density 1e-05 kg/m^3 + mean mol. weight 28.8565 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 117856 3.401e+06 J + internal energy -1471.39 -4.246e+04 J + entropy 10491.9 3.028e+05 J/K + Gibbs function -4.22729e+06 -1.22e+08 J + heat capacity c_p 1022.57 2.951e+04 J/K + heat capacity c_v 734.438 2.119e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.78847 0.765436 -34.8025 + O2 0.21153 0.234564 -37.7493 + N 6.54134e-55 3.17512e-55 -17.4013 + O 5.92517e-27 3.28519e-27 -18.8747 + NO 5.58575e-12 5.80828e-12 -36.2759 + NO+ 3.93876e-105 4.09561e-105 11.9791 + Electron 3.93876e-105 7.43897e-110 -254.391 + N+ 2.22844e-211 1.08163e-211 30.8538 + O+ 9.41966e-173 5.22252e-173 29.3803 + N2+ 6.09406e-170 5.91591e-170 13.4525 + O2+ 3.17559e-128 3.52133e-128 10.5057 + +Final T = 414.143, pres = 1.17767e-05 atm +Initial T = 414.143, pres = 0.000117767 atm + + air_below6000K: + + temperature 414.143 K + pressure 11.9328 Pa + density 0.0001 kg/m^3 + mean mol. weight 28.8565 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 117856 3.401e+06 J + internal energy -1471.39 -4.246e+04 J + entropy 9828.44 2.836e+05 J/K + Gibbs function -3.95253e+06 -1.141e+08 J + heat capacity c_p 1022.57 2.951e+04 J/K + heat capacity c_v 734.438 2.119e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.78847 0.765436 -32.4999 + O2 0.21153 0.234564 -35.4467 + N 2.06855e-55 1.00406e-55 -16.25 + O 1.8737e-27 1.03887e-27 -17.7234 + NO 5.58575e-12 5.80828e-12 -33.9733 + NO+ 2.91676e-106 3.03291e-106 11.6787 + Electron 2.91676e-106 5.50875e-111 -254.691 + N+ 5.21845e-213 2.5329e-213 29.4021 + O+ 2.20585e-174 1.22298e-174 27.9287 + N2+ 4.51281e-171 4.38088e-171 13.1521 + O2+ 2.35161e-129 2.60764e-129 10.2053 + +Final T = 414.143, pres = 0.000117767 atm +Initial T = 414.143, pres = 0.00117767 atm + + air_below6000K: + + temperature 414.143 K + pressure 119.328 Pa + density 0.001 kg/m^3 + mean mol. weight 28.8565 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 117856 3.401e+06 J + internal energy -1471.39 -4.246e+04 J + entropy 9165 2.645e+05 J/K + Gibbs function -3.67776e+06 -1.061e+08 J + heat capacity c_p 1022.57 2.951e+04 J/K + heat capacity c_v 734.438 2.119e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.78847 0.765436 -30.1973 + O2 0.21153 0.234564 -33.1441 + N 6.54134e-56 3.17512e-56 -15.0987 + O 5.92517e-28 3.28519e-28 -16.5721 + NO 5.58575e-12 5.80828e-12 -31.6707 + NO+ 3.85738e-107 4.01099e-107 11.9582 + Electron 3.85738e-107 7.28527e-112 -254.411 + N+ 2.1824e-214 1.05928e-214 28.5303 + O+ 9.22504e-176 5.11462e-176 27.0569 + N2+ 5.96815e-172 5.79368e-172 13.4316 + O2+ 3.10998e-130 3.44858e-130 10.4848 + +Final T = 414.143, pres = 0.00117767 atm +Initial T = 414.143, pres = 0.0117767 atm + + air_below6000K: + + temperature 414.143 K + pressure 1193.28 Pa + density 0.01 kg/m^3 + mean mol. weight 28.8565 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 117856 3.401e+06 J + internal energy -1471.39 -4.246e+04 J + entropy 8501.55 2.453e+05 J/K + Gibbs function -3.403e+06 -9.82e+07 J + heat capacity c_p 1022.57 2.951e+04 J/K + heat capacity c_v 734.438 2.119e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.78847 0.765436 -27.8948 + O2 0.21153 0.234564 -30.8416 + N 2.06855e-56 1.00406e-56 -13.9474 + O 1.8737e-28 1.03887e-28 -15.4208 + NO 5.58575e-12 5.80828e-12 -29.3682 + NO+ 6.14722e-108 6.39201e-108 12.4242 + Electron 6.14722e-108 1.161e-112 -253.945 + N+ 1.09981e-215 5.33822e-216 27.845 + O+ 4.64894e-177 2.5775e-177 26.3716 + N2+ 9.51097e-173 9.23294e-173 13.8976 + O2+ 4.95613e-131 5.49573e-131 10.9508 + +Final T = 414.143, pres = 0.0117767 atm +Initial T = 1153.16, pres = 3.27918e-08 atm + + air_below6000K: + + temperature 1152.46 K + pressure 0.00332065 Pa + density 1e-08 kg/m^3 + mean mol. weight 28.8562 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 930594 2.685e+07 J + internal energy 598529 1.727e+07 J + entropy 13299.8 3.838e+05 J/K + Gibbs function -1.43969e+07 -4.154e+08 J + heat capacity c_p 1172.94 3.385e+04 J/K + heat capacity c_v 884.81 2.553e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.788394 0.76537 -42.7126 + O2 0.211449 0.234477 -45.755 + N 2.86312e-15 1.38975e-15 -21.3563 + O 2.23425e-05 1.23879e-05 -22.8775 + NO 0.000134494 0.000139854 -44.2338 + NO+ 9.0387e-22 9.39874e-22 11.6756 + Electron 9.0387e-22 1.70712e-26 -69.7353 + N+ 1.59898e-54 7.76114e-55 34.553 + O+ 3.0545e-41 1.69352e-41 33.0319 + N2+ 1.2126e-44 1.17717e-44 13.1967 + O2+ 3.22549e-30 3.5767e-30 10.1544 + +Final T = 1152.46, pres = 3.27723e-08 atm +Initial T = 1152.46, pres = 3.27723e-07 atm + + air_below6000K: + + temperature 1152.61 K + pressure 0.0332106 Pa + density 1e-07 kg/m^3 + mean mol. weight 28.8564 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 930634 2.685e+07 J + internal energy 598529 1.727e+07 J + entropy 12636.3 3.646e+05 J/K + Gibbs function -1.36342e+07 -3.934e+08 J + heat capacity c_p 1172.97 3.385e+04 J/K + heat capacity c_v 884.832 2.553e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.7884 0.76537 -40.4103 + O2 0.211458 0.234486 -43.4526 + N 9.11146e-16 4.42265e-16 -20.2051 + O 7.08903e-06 3.93051e-06 -21.7263 + NO 0.000134663 0.000140028 -41.9314 + NO+ 4.80823e-20 4.99972e-20 17.9386 + Electron 4.80823e-20 9.08113e-25 -63.4588 + N+ 2.72207e-53 1.32123e-53 39.6649 + O+ 5.17819e-40 2.87094e-40 38.1438 + N2+ 6.49547e-43 6.30561e-43 19.4598 + O2+ 1.72003e-28 1.90731e-28 16.4175 + +Final T = 1152.61, pres = 3.27763e-07 atm +Initial T = 1152.61, pres = 3.27763e-06 atm + + air_below6000K: + + temperature 1152.66 K + pressure 0.332118 Pa + density 1e-06 kg/m^3 + mean mol. weight 28.8565 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 930647 2.686e+07 J + internal energy 598529 1.727e+07 J + entropy 11972.9 3.455e+05 J/K + Gibbs function -1.287e+07 -3.714e+08 J + heat capacity c_p 1172.97 3.385e+04 J/K + heat capacity c_v 884.84 2.553e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.788402 0.76537 -38.1078 + O2 0.211461 0.234488 -41.1501 + N 2.88709e-16 1.40138e-16 -19.0539 + O 2.24413e-06 1.24425e-06 -20.575 + NO 0.000134716 0.000140083 -39.6289 + NO+ 4.84916e-21 5.04226e-21 17.9428 + Electron 4.84916e-21 9.1584e-26 -63.4504 + N+ 8.71412e-55 4.22961e-55 38.5179 + O+ 1.65548e-41 9.17847e-42 36.9967 + N2+ 6.5652e-44 6.37328e-44 19.464 + O2+ 1.73602e-29 1.92503e-29 16.4217 + +Final T = 1152.66, pres = 3.27775e-06 atm +Initial T = 1152.66, pres = 3.27775e-05 atm + + air_below6000K: + + temperature 1152.67 K + pressure 3.32122 Pa + density 1e-05 kg/m^3 + mean mol. weight 28.8565 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 930651 2.686e+07 J + internal energy 598529 1.727e+07 J + entropy 11309.4 3.264e+05 J/K + Gibbs function -1.21054e+07 -3.493e+08 J + heat capacity c_p 1172.97 3.385e+04 J/K + heat capacity c_v 884.842 2.553e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.788402 0.76537 -35.8052 + O2 0.211462 0.234489 -38.8475 + N 9.1356e-17 4.43436e-17 -17.9026 + O 7.09896e-07 3.936e-07 -19.4238 + NO 0.000134733 0.000140101 -37.3264 + NO+ 4.85404e-22 5.04734e-22 17.9425 + Electron 4.85404e-22 9.16762e-27 -63.4494 + N+ 2.76173e-56 1.34047e-56 37.3662 + O+ 5.24445e-43 2.90767e-43 35.8451 + N2+ 6.5764e-45 6.38415e-45 19.4636 + O2+ 1.7382e-30 1.92744e-30 16.4213 + +Final T = 1152.67, pres = 3.27779e-05 atm +Initial T = 1152.67, pres = 0.000327779 atm + + air_below6000K: + + temperature 1152.68 K + pressure 33.2124 Pa + density 0.0001 kg/m^3 + mean mol. weight 28.8565 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 930652 2.686e+07 J + internal energy 598529 1.727e+07 J + entropy 10646 3.072e+05 J/K + Gibbs function -1.13407e+07 -3.273e+08 J + heat capacity c_p 1172.97 3.385e+04 J/K + heat capacity c_v 884.843 2.553e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.788403 0.76537 -33.5026 + O2 0.211462 0.23449 -36.5449 + N 2.88951e-17 1.40255e-17 -16.7513 + O 2.24513e-07 1.24481e-07 -18.2725 + NO 0.000134738 0.000140106 -35.0238 + NO+ 4.85477e-23 5.0481e-23 17.9422 + Electron 4.85477e-23 9.169e-28 -63.4492 + N+ 8.73799e-58 4.24119e-58 36.2147 + O+ 1.6591e-44 9.1985e-45 34.6935 + N2+ 6.57884e-46 6.38652e-46 19.4634 + O2+ 1.73859e-31 1.92788e-31 16.421 + +Final T = 1152.68, pres = 0.00032778 atm +Initial T = 1152.68, pres = 0.0032778 atm + + air_below6000K: + + temperature 1152.68 K + pressure 332.124 Pa + density 0.001 kg/m^3 + mean mol. weight 28.8565 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 930653 2.686e+07 J + internal energy 598529 1.727e+07 J + entropy 9982.53 2.881e+05 J/K + Gibbs function -1.0576e+07 -3.052e+08 J + heat capacity c_p 1172.97 3.385e+04 J/K + heat capacity c_v 884.843 2.553e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.788403 0.76537 -31.2 + O2 0.211463 0.23449 -34.2424 + N 9.13803e-18 4.43553e-18 -15.6 + O 7.09996e-08 3.93655e-08 -17.1212 + NO 0.00013474 0.000140108 -32.7212 + NO+ 4.85492e-24 5.04825e-24 17.9421 + Electron 4.85492e-24 9.16928e-29 -63.4492 + N+ 2.76361e-59 1.34139e-59 35.0633 + O+ 5.2471e-46 2.90914e-46 33.5421 + N2+ 6.5795e-47 6.38717e-47 19.4633 + O2+ 1.73869e-32 1.92799e-32 16.4209 + +Final T = 1152.68, pres = 0.00327781 atm +Initial T = 1152.68, pres = 0.0327781 atm + + air_below6000K: + + temperature 1152.68 K + pressure 3321.24 Pa + density 0.01 kg/m^3 + mean mol. weight 28.8565 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 930653 2.686e+07 J + internal energy 598529 1.727e+07 J + entropy 9319.08 2.689e+05 J/K + Gibbs function -9.81128e+06 -2.831e+08 J + heat capacity c_p 1172.97 3.385e+04 J/K + heat capacity c_v 884.843 2.553e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.788403 0.76537 -28.8975 + O2 0.211463 0.23449 -31.9398 + N 2.88976e-18 1.40267e-18 -14.4487 + O 2.24523e-08 1.24486e-08 -15.9699 + NO 0.00013474 0.000140109 -30.4186 + NO+ 9.89806e-28 1.02922e-27 11.7466 + Electron 9.89806e-28 1.8694e-32 -69.6446 + N+ 1.78181e-63 8.64844e-64 27.7165 + O+ 3.38296e-50 1.87561e-50 26.1954 + N2+ 1.34144e-50 1.30222e-50 13.2678 + O2+ 3.54481e-36 3.93076e-36 10.2255 + +Final T = 1152.68, pres = 0.0327781 atm +Initial T = 3008.94, pres = 8.55633e-08 atm + + air_below6000K: + + temperature 1764.11 K + pressure 0.00544642 Pa + density 1e-08 kg/m^3 + mean mol. weight 26.9308 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 2.94317e+06 7.926e+07 J + internal energy 2.39853e+06 6.459e+07 J + entropy 14438.6 3.888e+05 J/K + Gibbs function -2.25282e+07 -6.067e+08 J + heat capacity c_p 1254.18 3.378e+04 J/K + heat capacity c_v 945.447 2.546e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.734463 0.76399 -43.5726 + O2 0.129292 0.153622 -47.1099 + N 7.0471e-08 3.6652e-08 -21.7863 + O 0.133464 0.0792902 -23.555 + NO 0.00278025 0.00309774 -45.3413 + NO+ 5.61005e-11 6.25056e-11 -0.112179 + Electron 5.61006e-11 1.13531e-15 -45.2291 + N+ 1.27882e-29 6.65089e-30 23.4428 + O+ 2.26327e-21 1.34455e-21 21.6741 + N2+ 1.38336e-25 1.43895e-25 1.65648 + O2+ 1.1543e-16 1.3715e-16 -1.88084 + +Final T = 1764.11, pres = 5.3752e-08 atm +Initial T = 1764.11, pres = 5.3752e-07 atm + + air_below6000K: + + temperature 1879.19 K + pressure 0.0576568 Pa + density 1e-07 kg/m^3 + mean mol. weight 27.099 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 2.9751e+06 8.062e+07 J + internal energy 2.39853e+06 6.5e+07 J + entropy 13729.9 3.721e+05 J/K + Gibbs function -2.28259e+07 -6.186e+08 J + heat capacity c_p 1261.81 3.419e+04 J/K + heat capacity c_v 954.994 2.588e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.738358 0.763273 -41.4151 + O2 0.135652 0.16018 -44.9212 + N 1.60763e-07 8.3094e-08 -20.7075 + O 0.121806 0.0719148 -22.4606 + NO 0.00418364 0.00463245 -43.1681 + NO+ 1.47239e-10 1.63032e-10 -1.1314 + Electron 1.4724e-10 2.9612e-15 -42.0367 + N+ 4.3077e-28 2.22644e-28 21.3292 + O+ 2.0923e-20 1.23526e-20 19.5761 + N2+ 3.12509e-24 3.23048e-24 0.621638 + O2+ 6.6046e-16 7.79863e-16 -2.88445 + +Final T = 1879.19, pres = 5.69028e-07 atm +Initial T = 1879.19, pres = 5.69028e-06 atm + + air_below6000K: + + temperature 2006.78 K + pressure 0.611351 Pa + density 1e-06 kg/m^3 + mean mol. weight 27.2925 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 3.00988e+06 8.215e+07 J + internal energy 2.39853e+06 6.546e+07 J + entropy 13025.9 3.555e+05 J/K + Gibbs function -2.31301e+07 -6.313e+08 J + heat capacity c_p 1269.31 3.464e+04 J/K + heat capacity c_v 964.672 2.633e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.742614 0.76223 -39.2679 + O2 0.142745 0.167359 -42.7409 + N 3.48946e-07 1.79082e-07 -19.634 + O 0.108396 0.0635439 -21.3704 + NO 0.00624513 0.00686606 -41.0044 + NO+ 3.67988e-10 4.04569e-10 -2.10579 + Electron 3.67992e-10 7.34837e-15 -38.8986 + N+ 1.25612e-26 6.44628e-27 19.2646 + O+ 1.73236e-19 1.01551e-19 17.5282 + N2+ 6.36212e-23 6.53005e-23 -0.369314 + O2+ 3.53454e-15 4.14396e-15 -3.84227 + +Final T = 2006.78, pres = 6.03357e-06 atm +Initial T = 2006.78, pres = 6.03357e-05 atm + + air_below6000K: + + temperature 2147.22 K + pressure 6.48861 Pa + density 1e-05 kg/m^3 + mean mol. weight 27.5144 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 3.04739e+06 8.385e+07 J + internal energy 2.39853e+06 6.599e+07 J + entropy 12327.2 3.392e+05 J/K + Gibbs function -2.34218e+07 -6.444e+08 J + heat capacity c_p 1276.53 3.512e+04 J/K + heat capacity c_v 974.348 2.681e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.747195 0.760748 -37.1302 + O2 0.150579 0.175121 -40.5685 + N 7.05856e-07 3.5933e-07 -18.5651 + O 0.0930173 0.0540888 -20.2842 + NO 0.00920811 0.010042 -38.8493 + NO+ 8.56682e-10 9.34247e-10 -3.01307 + Electron 8.56699e-10 1.69693e-14 -35.8362 + N+ 2.98345e-25 1.51872e-25 17.2712 + O+ 1.2252e-18 7.12418e-19 15.552 + N2+ 1.11649e-21 1.13672e-21 -1.29393 + O2+ 1.71647e-14 1.9962e-14 -4.73222 + +Final T = 2147.22, pres = 6.40376e-05 atm +Initial T = 2147.22, pres = 0.000640376 atm + + air_below6000K: + + temperature 2299.14 K + pressure 68.8486 Pa + density 0.0001 kg/m^3 + mean mol. weight 27.7655 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 3.08701e+06 8.571e+07 J + internal energy 2.39853e+06 6.66e+07 J + entropy 11634.5 3.23e+05 J/K + Gibbs function -2.36623e+07 -6.57e+08 J + heat capacity c_p 1283.27 3.563e+04 J/K + heat capacity c_v 983.819 2.732e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.752001 0.758717 -34.9992 + O2 0.159067 0.18332 -38.402 + N 1.28755e-06 6.49527e-07 -17.4996 + O 0.0756137 0.0435711 -19.201 + NO 0.0133167 0.0143914 -36.7006 + NO+ 1.79309e-09 1.93776e-09 -3.81689 + Electron 1.79316e-09 3.51974e-14 -32.8837 + N+ 5.23851e-24 2.64255e-24 15.3841 + O+ 6.8634e-18 3.95478e-18 13.6827 + N2+ 1.57326e-20 1.58728e-20 -2.11547 + O2+ 7.20572e-14 8.30421e-14 -5.51831 + +Final T = 2299.14, pres = 0.000679483 atm +Initial T = 2299.14, pres = 0.00679483 atm + + air_below6000K: + + temperature 2457.49 K + pressure 728.699 Pa + density 0.001 kg/m^3 + mean mol. weight 28.0401 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 3.12723e+06 8.769e+07 J + internal energy 2.39853e+06 6.725e+07 J + entropy 10948.3 3.07e+05 J/K + Gibbs function -2.37781e+07 -6.667e+08 J + heat capacity c_p 1289.28 3.615e+04 J/K + heat capacity c_v 992.754 2.784e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.756823 0.756106 -32.8686 + O2 0.167916 0.191622 -36.2366 + N 2.00897e-06 1.00353e-06 -16.4343 + O 0.056584 0.0322863 -18.1183 + NO 0.0186748 0.0199842 -34.5526 + NO+ 3.18616e-09 3.4095e-09 -4.45753 + Electron 3.1864e-09 6.19324e-14 -30.0951 + N+ 5.81834e-23 2.9063e-23 13.6608 + O+ 2.69654e-17 1.53857e-17 11.9768 + N2+ 1.58747e-19 1.58593e-19 -2.77353 + O2+ 2.41729e-13 2.75853e-13 -6.14153 + +Final T = 2457.49, pres = 0.0071917 atm +Initial T = 2457.49, pres = 0.071917 atm + + air_below6000K: + + temperature 2610.53 K + pressure 7664.71 Pa + density 0.01 kg/m^3 + mean mol. weight 28.3183 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 3.165e+06 8.963e+07 J + internal energy 2.39853e+06 6.792e+07 J + entropy 10268.9 2.908e+05 J/K + Gibbs function -2.36423e+07 -6.695e+08 J + heat capacity c_p 1294.24 3.665e+04 J/K + heat capacity c_v 1000.63 2.834e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.761291 0.753095 -30.7261 + O2 0.176463 0.199398 -34.0624 + N 2.47224e-06 1.22281e-06 -15.3631 + O 0.0372965 0.0210719 -17.0312 + NO 0.0249473 0.0264342 -32.3943 + NO+ 4.39748e-09 4.65949e-09 -4.84133 + Electron 4.39806e-09 8.46427e-14 -27.5529 + N+ 3.21521e-22 1.59024e-22 12.1898 + O+ 6.15914e-17 3.47969e-17 10.5217 + N2+ 9.61007e-19 9.50643e-19 -3.17322 + O2+ 5.73802e-13 6.48366e-13 -6.50944 + +Final T = 2610.53, pres = 0.0756448 atm +Initial T = 6437.5, pres = 1.86538e-07 atm + + air_below6000K: + + temperature 2945.65 K + pressure 0.0104079 Pa + density 1e-08 kg/m^3 + mean mol. weight 23.5318 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 7.43931e+06 1.751e+08 J + internal energy 6.39853e+06 1.506e+08 J + entropy 16362.2 3.85e+05 J/K + Gibbs function -4.0758e+07 -9.591e+08 J + heat capacity c_p 1321.41 3.11e+04 J/K + heat capacity c_v 968.084 2.278e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.630849 0.750996 -44.8654 + O2 1.39323e-06 1.89453e-06 -59.788 + N 0.0241541 0.0143771 -22.4327 + O 0.344887 0.234491 -29.894 + NO 0.000102862 0.000131163 -52.3267 + NO+ 2.91192e-06 3.71302e-06 -17.4904 + Electron 2.91277e-06 6.74603e-11 -34.8363 + N+ 7.39583e-12 4.40202e-12 12.4035 + O+ 8.35023e-10 5.67718e-10 4.94227 + N2+ 3.21595e-12 3.82836e-12 -10.0292 + O2+ 3.03228e-12 4.12326e-12 -24.9517 + +Final T = 2945.65, pres = 1.02718e-07 atm +Initial T = 2945.65, pres = 1.02718e-06 atm + + air_below6000K: + + temperature 3079.99 K + pressure 0.108465 Pa + density 1e-07 kg/m^3 + mean mol. weight 23.61 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 7.48317e+06 1.767e+08 J + internal energy 6.39853e+06 1.511e+08 J + entropy 15550 3.671e+05 J/K + Gibbs function -4.04105e+07 -9.541e+08 J + heat capacity c_p 1323.23 3.124e+04 J/K + heat capacity c_v 971.073 2.293e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.636117 0.75476 -42.6792 + O2 5.85321e-06 7.9329e-06 -56.1844 + N 0.0177433 0.0105263 -21.3396 + O 0.345878 0.234386 -28.0922 + NO 0.000248908 0.00031634 -49.4318 + NO+ 3.27998e-06 4.16849e-06 -16.9579 + Electron 3.28083e-06 7.57329e-11 -32.4739 + N+ 6.29836e-12 3.73639e-12 11.1343 + O+ 8.26462e-10 5.60036e-10 4.38167 + N2+ 4.52429e-12 5.36801e-12 -10.2053 + O2+ 9.62394e-12 1.30432e-11 -23.7105 + +Final T = 3079.99, pres = 1.07046e-06 atm +Initial T = 3079.99, pres = 1.07046e-05 atm + + air_below6000K: + + temperature 3207.45 K + pressure 1.12586 Pa + density 1e-06 kg/m^3 + mean mol. weight 23.687 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 7.52439e+06 1.782e+08 J + internal energy 6.39853e+06 1.516e+08 J + entropy 14740.4 3.492e+05 J/K + Gibbs function -3.97547e+07 -9.417e+08 J + heat capacity c_p 1324.78 3.138e+04 J/K + heat capacity c_v 973.764 2.307e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.641058 0.758149 -40.4837 + O2 2.74938e-05 3.71414e-05 -52.4587 + N 0.0116959 0.0069161 -20.2419 + O 0.346589 0.234104 -26.2294 + NO 0.000623449 0.000789772 -46.4712 + NO+ 3.3837e-06 4.28632e-06 -16.2766 + Electron 3.38439e-06 7.78696e-11 -30.1947 + N+ 3.78744e-12 2.23953e-12 9.9528 + O+ 6.56589e-10 4.43479e-10 3.9653 + N2+ 4.89296e-12 5.78656e-12 -10.2891 + O2+ 2.83815e-11 3.83399e-11 -22.2641 + +Final T = 3207.45, pres = 1.11114e-05 atm +Initial T = 3207.45, pres = 0.000111114 atm + + air_below6000K: + + temperature 3318.59 K + pressure 11.6123 Pa + density 1e-05 kg/m^3 + mean mol. weight 23.7613 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 7.55976e+06 1.796e+08 J + internal energy 6.39853e+06 1.52e+08 J + entropy 13933.4 3.311e+05 J/K + Gibbs function -3.86797e+07 -9.191e+08 J + heat capacity c_p 1326 3.151e+04 J/K + heat capacity c_v 976.084 2.319e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.645083 0.760523 -38.2726 + O2 0.000148603 0.00020012 -48.5744 + N 0.00670202 0.00395068 -19.1363 + O 0.346431 0.233265 -24.2872 + NO 0.00162924 0.00205743 -43.4235 + NO+ 3.10975e-06 3.92697e-06 -15.4005 + Electron 3.11022e-06 7.13375e-11 -28.023 + N+ 1.45279e-12 8.56351e-13 8.88668 + O+ 3.92147e-10 2.64039e-10 3.7358 + N2+ 3.76211e-12 4.43526e-12 -10.2496 + O2+ 7.586e-11 1.02157e-10 -20.5514 + +Final T = 3318.59, pres = 0.000114605 atm +Initial T = 3318.59, pres = 0.00114605 atm + + air_below6000K: + + temperature 3417.09 K + pressure 119.127 Pa + density 0.0001 kg/m^3 + mean mol. weight 23.8496 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 7.5898e+06 1.81e+08 J + internal energy 6.39853e+06 1.526e+08 J + entropy 13129.1 3.131e+05 J/K + Gibbs function -3.72736e+07 -8.89e+08 J + heat capacity c_p 1326.98 3.165e+04 J/K + heat capacity c_v 978.356 2.333e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.647745 0.760834 -36.0515 + O2 0.000876305 0.00117573 -44.5898 + N 0.00347447 0.00204054 -18.0258 + O 0.343542 0.230464 -22.2949 + NO 0.00435766 0.00548255 -40.3207 + NO+ 2.6229e-06 3.29992e-06 -14.3889 + Electron 2.62328e-06 5.9946e-11 -25.9318 + N+ 4.05713e-13 2.38264e-13 7.90602 + O+ 1.90592e-10 1.27853e-10 3.63687 + N2+ 2.27319e-12 2.67001e-12 -10.1197 + O2+ 1.87267e-10 2.5125e-10 -18.658 + +Final T = 3417.09, pres = 0.00117569 atm +Initial T = 3417.09, pres = 0.0117569 atm + + air_below6000K: + + temperature 3555.18 K + pressure 1230.49 Pa + density 0.001 kg/m^3 + mean mol. weight 24.0224 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 7.62902e+06 1.833e+08 J + internal energy 6.39853e+06 1.537e+08 J + entropy 12328.9 2.962e+05 J/K + Gibbs function -3.62024e+07 -8.697e+08 J + heat capacity c_p 1328.33 3.191e+04 J/K + heat capacity c_v 982.221 2.36e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.649921 0.757899 -33.8646 + O2 0.00420976 0.00560757 -40.8464 + N 0.00209788 0.00122321 -16.9323 + O 0.332941 0.221746 -20.4232 + NO 0.0108249 0.0135213 -37.3555 + NO+ 2.48454e-06 3.10336e-06 -13.6141 + Electron 2.48517e-06 5.63814e-11 -23.7415 + N+ 1.88247e-13 1.09757e-13 6.80913 + O+ 1.25428e-10 8.35346e-11 3.31826 + N2+ 2.02806e-12 2.36495e-12 -10.1232 + O2+ 4.96803e-10 6.6175e-10 -17.1049 + +Final T = 3555.18, pres = 0.012144 atm +Initial T = 3555.18, pres = 0.12144 atm + + air_below6000K: + + temperature 3800.16 K + pressure 12968.3 Pa + density 0.01 kg/m^3 + mean mol. weight 24.3643 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 7.69536e+06 1.875e+08 J + internal energy 6.39853e+06 1.559e+08 J + entropy 11537 2.811e+05 J/K + Gibbs function -3.61473e+07 -8.807e+08 J + heat capacity c_p 1331 3.243e+04 J/K + heat capacity c_v 989.742 2.411e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.653418 0.751284 -31.7612 + O2 0.0124851 0.0163972 -37.678 + N 0.00186299 0.00107101 -15.8806 + O 0.309477 0.203225 -18.839 + NO 0.0227506 0.0280187 -34.7196 + NO+ 3.18462e-06 3.92198e-06 -13.4287 + Electron 3.18642e-06 7.12763e-11 -21.2909 + N+ 3.11855e-13 1.79275e-13 5.41026 + O+ 1.78541e-10 1.17239e-10 2.4519 + N2+ 4.81467e-12 5.53568e-12 -10.4704 + O2+ 1.62057e-09 2.12834e-09 -16.3871 + +Final T = 3800.16, pres = 0.127987 atm +Initial T = 7052.07, pres = 2.37509e-07 atm + + air_below6000K: + + temperature 3293.99 K + pressure 0.0125025 Pa + density 1e-08 kg/m^3 + mean mol. weight 21.9058 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 1.09488e+07 2.398e+08 J + internal energy 9.69853e+06 2.125e+08 J + entropy 17407.2 3.813e+05 J/K + Gibbs function -4.63903e+07 -1.016e+09 J + heat capacity c_p 1343.53 2.943e+04 J/K + heat capacity c_v 963.974 2.112e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.51821 0.662696 -45.2973 + O2 1.57963e-07 2.30744e-07 -62.225 + N 0.160619 0.102701 -22.6487 + O 0.321096 0.23452 -31.1125 + NO 4.64565e-05 6.36352e-05 -53.7612 + NO+ 1.40272e-05 1.92139e-05 -20.4316 + Electron 1.40837e-05 3.50393e-10 -33.3296 + N+ 4.79433e-09 3.06541e-09 10.681 + O+ 5.12902e-08 3.74597e-08 2.21711 + N2+ 4.04428e-10 5.17179e-10 -11.9677 + O2+ 1.18531e-11 1.7314e-11 -28.8954 + +Final T = 3293.99, pres = 1.23391e-07 atm +Initial T = 3293.99, pres = 1.23391e-06 atm + + air_below6000K: + + temperature 3516.01 K + pressure 0.132781 Pa + density 1e-07 kg/m^3 + mean mol. weight 22.0165 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 1.10263e+07 2.428e+08 J + internal energy 9.69853e+06 2.135e+08 J + entropy 16535.4 3.641e+05 J/K + Gibbs function -4.71123e+07 -1.037e+09 J + heat capacity c_p 1346.82 2.965e+04 J/K + heat capacity c_v 969.171 2.134e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.525826 0.669053 -43.1683 + O2 5.17854e-07 7.5265e-07 -58.9387 + N 0.151375 0.0963037 -21.5842 + O 0.322655 0.234474 -29.4693 + NO 0.000104379 0.000142258 -51.0535 + NO+ 1.95395e-05 2.66298e-05 -20.2701 + Electron 1.96343e-05 4.86031e-10 -30.7834 + N+ 9.13253e-09 5.80982e-09 9.19923 + O+ 8.45593e-08 6.14472e-08 1.31404 + N2+ 1.05956e-09 1.34815e-09 -12.3849 + O2+ 4.50091e-11 6.54152e-11 -28.1553 + +Final T = 3516.01, pres = 1.31045e-06 atm +Initial T = 3516.01, pres = 1.31045e-05 atm + + air_below6000K: + + temperature 3766.66 K + pressure 1.41427 Pa + density 1e-06 kg/m^3 + mean mol. weight 22.1442 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 1.11128e+07 2.461e+08 J + internal energy 9.69853e+06 2.148e+08 J + entropy 15668.3 3.47e+05 J/K + Gibbs function -4.79042e+07 -1.061e+09 J + heat capacity c_p 1350.74 2.991e+04 J/K + heat capacity c_v 975.268 2.16e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.534543 0.676223 -41.0514 + O2 1.72905e-06 2.49852e-06 -55.6497 + N 0.140779 0.0890463 -20.5257 + O 0.324385 0.234372 -27.8249 + NO 0.000236103 0.000319927 -48.3506 + NO+ 2.7034e-05 3.66313e-05 -20.1004 + Electron 2.71895e-05 6.69171e-10 -28.2502 + N+ 1.66858e-08 1.05537e-08 7.72452 + O+ 1.35925e-07 9.82041e-08 0.425342 + N2+ 2.70465e-09 3.42145e-09 -12.8012 + O2+ 1.70141e-10 2.45852e-10 -27.3995 + +Final T = 3766.66, pres = 1.39577e-05 atm +Initial T = 3766.66, pres = 0.000139577 atm + + air_below6000K: + + temperature 4050.35 K + pressure 15.1065 Pa + density 1e-05 kg/m^3 + mean mol. weight 22.2928 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 1.12092e+07 2.499e+08 J + internal energy 9.69853e+06 2.162e+08 J + entropy 14806.5 3.301e+05 J/K + Gibbs function -4.87625e+07 -1.087e+09 J + heat capacity c_p 1355.57 3.022e+04 J/K + heat capacity c_v 982.601 2.19e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.544538 0.684274 -38.9472 + O2 5.90966e-06 8.48263e-06 -52.3543 + N 0.128599 0.0807999 -19.4736 + O 0.326244 0.234143 -26.1771 + NO 0.000538494 0.000724813 -45.6507 + NO+ 3.70145e-05 4.98205e-05 -19.9165 + Electron 3.72615e-05 9.10944e-10 -25.7343 + N+ 2.87988e-08 1.80938e-08 6.26065 + O+ 2.10966e-07 1.51403e-07 -0.442893 + N2+ 6.65571e-09 8.36349e-09 -13.213 + O2+ 6.38182e-10 9.16021e-10 -26.62 + +Final T = 4050.35, pres = 0.000149089 atm +Initial T = 4050.35, pres = 0.00149089 atm + + air_below6000K: + + temperature 4371.57 K + pressure 161.774 Pa + density 1e-04 kg/m^3 + mean mol. weight 22.4679 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 1.13163e+07 2.543e+08 J + internal energy 9.69853e+06 2.179e+08 J + entropy 13951 3.134e+05 J/K + Gibbs function -4.96714e+07 -1.116e+09 J + heat capacity c_p 1361.54 3.059e+04 J/K + heat capacity c_v 991.483 2.228e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.555954 0.693176 -36.856 + O2 2.08096e-05 2.9637e-05 -49.0466 + N 0.114619 0.0714548 -18.428 + O 0.328065 0.233616 -24.5233 + NO 0.00124059 0.00165683 -42.9513 + NO+ 4.98713e-05 6.66025e-05 -19.7098 + Electron 5.02468e-05 1.21883e-09 -23.2415 + N+ 4.59398e-08 2.86383e-08 4.81344 + O+ 3.11682e-07 2.21942e-07 -1.28183 + N2+ 1.55453e-08 1.93819e-08 -13.6146 + O2+ 2.3627e-09 3.3649e-09 -25.8051 + +Final T = 4371.57, pres = 0.00159659 atm +Initial T = 4371.57, pres = 0.0159659 atm + + air_below6000K: + + temperature 4734.02 K + pressure 1735.66 Pa + density 0.001 kg/m^3 + mean mol. weight 22.6777 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 1.14342e+07 2.593e+08 J + internal energy 9.69853e+06 2.199e+08 J + entropy 13102.7 2.971e+05 J/K + Gibbs function -5.05944e+07 -1.147e+09 J + heat capacity c_p 1368.64 3.104e+04 J/K + heat capacity c_v 1002.01 2.272e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.568792 0.702621 -34.7776 + O2 7.60334e-05 0.000107285 -45.7193 + N 0.0987443 0.0609887 -17.3888 + O 0.329365 0.232371 -22.8596 + NO 0.00289126 0.00382559 -40.2484 + NO+ 6.55575e-05 8.67412e-05 -19.468 + Electron 6.6095e-05 1.58842e-09 -20.7804 + N+ 6.57105e-08 4.05841e-08 3.39166 + O+ 4.2952e-07 3.03021e-07 -2.07919 + N2+ 3.37049e-08 4.16344e-08 -13.9971 + O2+ 8.55849e-09 1.2076e-08 -24.9388 + +Final T = 4734.02, pres = 0.0171297 atm +Initial T = 4734.02, pres = 0.171297 atm + + air_below6000K: + + temperature 5139.7 K + pressure 18630.6 Pa + density 0.01 kg/m^3 + mean mol. weight 22.9374 amu + + 1 kg 1 kmol + ----------- ------------ + enthalpy 1.15616e+07 2.652e+08 J + internal energy 9.69853e+06 2.225e+08 J + entropy 12263.1 2.813e+05 J/K + Gibbs function -5.14672e+07 -1.181e+09 J + heat capacity c_p 1376.18 3.157e+04 J/K + heat capacity c_v 1013.69 2.325e+04 J/K + + X Y Chem. Pot. / RT + ------------- ------------ ------------ + N2 0.58266 0.711602 -32.7112 + O2 0.00028938 0.000403698 -42.3671 + N 0.0812643 0.049624 -16.3556 + O 0.328809 0.229352 -21.1836 + NO 0.00681051 0.00890932 -37.5392 + NO+ 8.30639e-05 0.00010866 -19.1761 + Electron 8.37789e-05 1.99061e-09 -18.3631 + N+ 8.12574e-08 4.96179e-08 2.00744 + O+ 5.37869e-07 3.75163e-07 -2.82049 + N2+ 6.60092e-08 8.06154e-08 -14.3482 + O2+ 2.99e-08 4.17113e-08 -24.004 + +Final T = 5139.7, pres = 0.18387 atm diff --git a/test_problems/ChemEquil_ionizedGas/runtest b/test_problems/ChemEquil_ionizedGas/runtest new file mode 100755 index 000000000..be204380e --- /dev/null +++ b/test_problems/ChemEquil_ionizedGas/runtest @@ -0,0 +1,32 @@ +#!/bin/sh +# +# +temp_success="1" +/bin/rm -f output.txt outputa.txt +testName="ChemEquil_ionizedGas" + +################################################################# +# +################################################################# +CANTERA_DATA=${CANTERA_DATA:=../../data/inputs}; export CANTERA_DATA + +CANTERA_BIN=${CANTERA_BIN:=../../bin} +./ionizedGasEquil > output.txt +retnStat=$? +if [ $retnStat != "0" ] +then + temp_success="1" + echo "ERROR: $testName: ionizedGasEquil returned with bad status, $retnStat, check output" +fi + +diff -w output.txt output_blessed.txt > diff_test.out +retnStat=$? +if [ $retnStat = "0" ] +then + echo "successful diff comparison on $testName test" + temp_success="0" +else + echo "FAILED: unsuccessful diff comparison on $testName test" + echo "FAILED" > csvCode.txt +fi +