FM2Cantera/plot_tables_CH4_PV.py at master · luzyPer/FM2Cantera · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
"""
Zhen Lu 2017/07/19 <albert.lz07@gmail.com>

plot flamelet solutions, progress variable C against mixture fraction Z
"""

import numpy as np
import matplotlib.pyplot as plt
import glob

SMALL = 1.e-20

folder_name = 'tables_CH4/'
table_pre   = 'Table_'
file_end    = '.csv'

# figure and axes parameters
plot_width      =18.0 / 2.54
plot_height     =18.0 / 2.54
margin_left     =2.0 / 2.54
margin_right    =0.3 / 2.54
margin_bottom   =1.5 / 2.54
margin_top      =0.3 / 2.54
space_width     =0.0 / 2.54
space_height    =0.0 / 2.54
ftsize          =12

# use TEX for interpreter
plt.rc('text',usetex=True)
# use serif font
plt.rc('font',family='serif')
# generate the figure
fig, axes = plt.subplots(1,1,figsize=(plot_width, plot_height))

for table in glob.glob('{0}{1}*{2}'.format(folder_name,table_pre,file_end)):
    data = np.genfromtxt(table,delimiter=',',names=True)
    chist = table[len(folder_name)+len(table_pre):table.find(file_end)]
    label = r'$\chi_{st}=$'+chist

    PV = np.zeros(len(data['Z']))

    # Y_O
    PV=(data['H2O']/18.+data['CO']/28.+data['CO2']/22.)*16.

    # CO2 CO
    #PV=data['CO2']+data['CO']

    # CO2 CO H2O H2
    #PV=data['CO2']+data['CO']+data['H2O']+data['H2']

    # calculate the equivalence ratio

    # major species
    phi = (data['CO2']/44.+data['CH4']/8.+data['CO']/28.+0.5*(data['H2O']/18.+data['H2']/2.)) \
            /(data['CO2']/44.+data['O2']/32.+0.5*(data['H2O']/18.+data['CO']/28.))

    # C_O
    for i, p in enumerate(phi):
        if phi[i] < 1.0:
            PV[i] = PV[i]/16./(data['CO2'][i]/22.+data['CO'][i]/28.+data['H2O'][i]/18.+data['H2'][i]/2.+data['CH4'][i]/4.)
        else:
            PV[i] = PV[i]/16./(data['CO2'][i]/22.+data['CO'][i]/28.+data['H2O'][i]/18.+data['O2'][i]/16.)

    plt.plot(data['Z'],PV,label=label,linewidth=1.5)

axes.set_xlabel(r'$Z$',fontsize=ftsize)
axes.set_xlim(0.0,1.0)

# too many lines
#axes.legend(fontsize=ftsize,frameon=False)

plt.subplots_adjust(left   = margin_left   / plot_width,
                    bottom = margin_bottom / plot_height,
                    right  = 1.0 - margin_right / plot_width,
                    top    = 1.0 - margin_top   / plot_height,
                    wspace = space_width   / plot_width,
                    hspace = space_height  / plot_height)

plt.show()