我正在做一个关于不同复杂程度的层流预混火焰数值分析的项目。由于我是 Python 和 Cantera 的新手,我非常需要帮助。我有以下疑问:
- 我可以修改 adiabetic.py 模块吗?该模块主要给出混合物的绝热温度(从入口温度)到可以从绝热温度作为输入给出入口温度的值?
程序如下:
"""
Adiabatic flame temperature and equilibrium composition for a fuel/air mixture
as a function of equivalence ratio, including formation of solid carbon.
"""
import cantera as ct
import numpy as np
import sys
import csv
##############################################################################
# Edit these parameters to change the initial temperature, the pressure, and
# the phases in the mixture.
T = 300.0
P = 101325.0
# phases
gas = ct.Solution('gri30.xml')
carbon = ct.Solution('graphite.xml')
# the phases that will be included in the calculation, and their initial moles
mix_phases = [(gas, 1.0), (carbon, 0.0)]
# gaseous fuel species
fuel_species = 'CH4'
# equivalence ratio range
npoints = 50
phi = np.linspace(0.3, 3.5, npoints)
##############################################################################
mix = ct.Mixture(mix_phases)
# create some arrays to hold the data
tad = np.zeros(npoints)
xeq = np.zeros((mix.n_species,npoints))
for i in range(npoints):
# set the gas state
gas.set_equivalence_ratio(phi[i], fuel_species, 'O2:1.0, N2:3.76')
# create a mixture of 1 mole of gas, and 0 moles of solid carbon.
mix = ct.Mixture(mix_phases)
mix.T = T
mix.P = P
# equilibrate the mixture adiabatically at constant P
mix.equilibrate('HP', solver='gibbs', max_steps=1000)
tad[i] = mix.T
print('At phi = {0:12.4g}, Tad = {1:12.4g}'.format(phi[i], tad[i]))
xeq[:,i] = mix.species_moles
# write output CSV file for importing into Excel
csv_file = 'adiabatic.csv'
with open(csv_file, 'w') as outfile:
writer = csv.writer(outfile)
writer.writerow(['phi','T (K)'] + mix.species_names)
for i in range(npoints):
writer.writerow([phi[i], tad[i]] + list(xeq[:,i]))
print('Output written to {0}'.format(csv_file))
if '--plot' in sys.argv:
import matplotlib.pyplot as plt
plt.plot(phi, tad)
plt.xlabel('Equivalence ratio')
plt.ylabel('Adiabatic flame temperature [K]')
plt.show()
取自:http://www.cantera.org/docs/sphinx/html/cython/examples/multiphase_adiabatic.html
请问我具体的问题,因为我也不知 Prop 体要问什么。但我知道我想做什么。
最佳答案
这是一个有趣的问题。
原则上应该不会太困难。假设您已经给出了入口燃料成分、当量比和绝热火焰温度,步骤如下:
- 将温度设置为绝热火焰温度,并根据设定的 T 和 P 平衡燃料-空气混合物。
- 计算(并保存)混合物平均焓。
- 将成分设置回入口成分,并拟合温度,直到混合物平均焓与步骤 2 中的值相匹配。
关于python - 如何反转 CANTERA Python 模块中的 adiabetic.py 程序,使其输入绝热温度并给出输出入口温度?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/37851443/