import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline
Read the posted notes for variables, dictionaries, and tuples. State if you did this.
If a dictionary has a value that is an array, we can access elements of the array as follows. (First we make the dictionary, then we access an array element.)
a = {
'CO2':np.array([0.1, 2.2, 33.3]),
'H2O':np.array([0.2, 3.8, 44.4]), # useful optional end comma
} # makes each line the same
a['CO2'][1] # --> 2.2
That is a['CO2']
is the array, and the following [1]
accesses the given element.
The heat capacity of a species $i$ is given by $$c_{p,i}(T) = R_g(a_{0,i} + a_{1,i}T + a_{2,i}T^2 + a_{3,i}T^3 + a_{4,i}T^4)$$
The following species data are given:
Species | $a_0$ | $a_1$ | $a_2$ | $a_3$ | $a_4$ |
---|---|---|---|---|---|
CO2 | 2.356773520E+00 | 8.984596770E-03 | -7.123562690E-06 | 2.459190220E-09 | -1.436995480E-13 |
H2O | 4.198640560E+00 | -2.036434100E-03 | 6.520402110E-06 | -5.487970620E-09 | 1.771978170E-12 |
O2 | 3.782456360E+00 | -2.996734160E-03 | 9.847302010E-06 | -9.681295090E-09 | 3.243728370E-12 |
N2 | 3.298677000E+00 | 1.408240400E-03 | -3.963222000E-06 | 5.641515000E-09 | -2.444854000E-12 |
CH4 | 5.149876130E+00 | -1.367097880E-02 | 4.918005990E-05 | -4.847430260E-08 | 1.666939560E-11 |
cp
that takes two arguments: a string that gives the name of the species, and the temperature.
where $x_i$ are species mole fractions.
cp_mix
that takes two arguments: a dictionary containing the names and mole fractions of species, and the temperature.
cp_mix
to evaluate the mixture heat capacity at T=900 K, for the following mixture (corresponding to a stoichiometric mixture of methane and air (CH4 + 2O2 + 7.52 N2):