Fix mass flow rate not dividied by the number of ovens

This commit is contained in:
Yeongdo Park 2022-12-13 20:06:47 +09:00
parent d9c6669d70
commit b06155a7cf

View file

@ -241,7 +241,9 @@ class Battery:
self.hv = hv # Base unit heat J/kg
self.normal_heat = self.heat_program.f(-1) # GJ / rev
# Energy input to battery
Q0 = self.normal_heat * 1e9 * 3 / 3600 # GJ/rev => J/s (W)
# Equivalent Fuel+Air mass flow
mdot0 = Q0 / hv # (J/s) / (J/kg) => kg/s
self.mdot0 = mdot0 # kg / s
@ -250,7 +252,7 @@ class Battery:
# ... walls_1[size-2] - chambers[size-1] - walls_0[size-1] - ovens[size-1] - walls_1[size-1] - chambers[size]
self.chambers = [
CombustionChamber(self.mdot0, self.gas,
CombustionChamber(self.mdot0/self.size, self.gas,
(self.T0, self.P0, self.X0), hA=700)
for ichamber in range(self.size+1)
]
@ -315,7 +317,7 @@ class Battery:
else:
wall_upper = None
chmbr.update_mdot(self.mdot(self.t))
chmbr.update_mdot(self.mdot(self.t)/self.size)
chmbr.update_Twall(
Twall0=(wall_lower.T_chamber if wall_lower else wall_upper.T_chamber),
Twall1=(wall_upper.T_chamber if wall_upper else wall_lower.T_chamber),
@ -425,7 +427,7 @@ class Battery:
def coke_oven_exhaust_stoichiometry (phi=1.0, return_unburned=False):
# Define the oxidizer composition, here air with 21 mol-% O2 and 79 mol-% N2
air = "O2:0.21,N2:0.79"
air = "O2:1,N2:3.762"
coke_oven_fuel = "H2:6.42, O2:0.39, N2:47.28, CH4:1.79, CO:24.25, CO2:19.72, C2H4:0.13, C2H6:0.04"
mix = ct.Solution('gri30.xml')