From 1d9674ff45b3479a1cffeafe1be280ea3a0848e4 Mon Sep 17 00:00:00 2001 From: ignis Date: Wed, 17 Nov 2021 21:43:24 +0900 Subject: [PATCH] mass of moisture and ash contents added to N2 --- coal.py | 12 +++++++----- 1 file changed, 7 insertions(+), 5 deletions(-) diff --git a/coal.py b/coal.py index fc40600..0bde57c 100644 --- a/coal.py +++ b/coal.py @@ -15,7 +15,7 @@ fc = 0.4347 ash = 0.1055 vm_fc = 1. - moist - ash -fuelMfr = coalMfr * vm_fc +fuelMfr = coalMfr gri_species = {S.name: S for S in ct.Species.listFromFile('gri30.cti')} nasa_species = {S.name: S for S in ct.Species.listFromFile('nasa_gas.cti')} @@ -37,7 +37,7 @@ airs = [ ] for i, air in enumerate(airs): - pass # print("air", i, air.report()) + print("air {}, T = {}, mass flow rate = {}".format(i+1, air.T, air.mass)) airmix = reduce(lambda a, b: a+b, airs) print("Total Air flow rate = ", airmix.mass) @@ -45,6 +45,7 @@ print("Total Air flow rate = ", airmix.mass) """############################################################################# Dummy gaseous coal object containing elementary composition + mass of moisture and ash contents is added to N2 #############################################################################""" coal = ct.Solution(thermo='IdealGas', species=system_species) @@ -52,7 +53,7 @@ coal.TPY = 348.15, ct.one_atm, '''\ C: 81.41, H2: 5.47, O2: 10.83, - N2: 1.71, + N2: 36.9647930755, S: 0.57 ''' @@ -63,8 +64,9 @@ Heating value = Sum(Product Enthalpy of Formation) - Sum(Reactant Enthalpy of Fo #############################################################################""" # coal heating value = 5761 Kcal/kg * 4.184 kJ/kcal -coalHV = 5761 * 4184. # J/kg, with Moisture and Ash +coalHV = - 5761 * 4.184 # J/kg, with Moisture and Ash coalHVdaf = - 5761 * 4.184 / vm_fc # kJ/kg, Dry and Ash Free, Negative since exothermic +print("HV of Coal , kJ/kg = ", coalHV) print("HV of Coal(daf), kJ/kg = ", coalHVdaf) # sum product of 1kg coal enthalpy of formation - kJ/kg @@ -74,7 +76,7 @@ sum_product_hf = ( - 8919.38250 * coal.mass_fraction_dict()['S']) print("Sum(Hf_product), kJ/kg = ", sum_product_hf) -sum_coal_hf = - coalHVdaf + sum_product_hf +sum_coal_hf = - coalHV + sum_product_hf print("Sum(Hf_reactant), kJ/kg = ", sum_coal_hf)