From b06155a7cfa3eea5cd7fd29a2067b37755823ed8 Mon Sep 17 00:00:00 2001 From: Yeongdo Park Date: Tue, 13 Dec 2022 20:06:47 +0900 Subject: [PATCH] Fix mass flow rate not dividied by the number of ovens --- Battery.py | 8 +++++--- 1 file changed, 5 insertions(+), 3 deletions(-) diff --git a/Battery.py b/Battery.py index 174495e..e7b5812 100644 --- a/Battery.py +++ b/Battery.py @@ -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')