gas_composition.h

Go to the documentation of this file.

#pragma once
#include <string>
#include <array>
 
#include "physics/gas_constants.h"
#include "physics/gas_properties.h"
using namespace gas_properties;
using namespace gas_constants;
 
namespace gas_composition {
 
struct GasFlueGasComponents {
    double mO2;
    double mN2;
    double mCO2;
    double mH2O;
    double mSO2;
};
 
class GasComposition {
  public:
    std::string   substance;
    GasProperties ch4;
    GasProperties c2h6;
    GasProperties n2;
    GasProperties h2;
    GasProperties c3h8;
    GasProperties c4h10_cnh2n;
    GasProperties h2o;
    GasProperties co;
    GasProperties co2;
    GasProperties so2;
    GasProperties o2;
    double        total_composition_weight;
    double        heating_value;
    double        heating_value_volume;
    double        specific_gravity;
    double        stoichometric_air;
    GasComposition(const std::string substance, double ch4_percent, double c2h6_percent, double n2_percent,
                   double h2_percent, double c3h8_percent, double c4h10_cnh2n_percent, double h2o_percent,
                   double co_percent, double co2_percent, double so2_percent, double o2_percent)
        : substance(substance), ch4(specificHeatCH4, 0, 0, 0, 0, 0, 0, 0, 0, 0),
          c2h6(specificHeatC2H6, 0, 0, 0, 0, 0, 0, 0, 0, 0), n2(specificHeatN2, 0, 0, 0, 0, 0, 0, 0, 0, 0),
          h2(specificHeatH2, 0, 0, 0, 0, 0, 0, 0, 0, 0), c3h8(specificHeatC3H8, 0, 0, 0, 0, 0, 0, 0, 0, 0),
          c4h10_cnh2n(specificHeatC4H10CnH2n, 0, 0, 0, 0, 0, 0, 0, 0, 0),
          h2o(specificHeatH2O, 0, 0, 0, 0, 0, 0, 0, 0, 0), co(specificHeatCO, 0, 0, 0, 0, 0, 0, 0, 0, 0),
          co2(specificHeatCO2, 0, 0, 0, 0, 0, 0, 0, 0, 0), so2(specificHeatSO2, 0, 0, 0, 0, 0, 0, 0, 0, 0),
          o2(specificHeatO2, 0, 0, 0, 0, 0, 0, 0, 0, 0) {
        double total_percent = ch4_percent + c2h6_percent + n2_percent + h2_percent + c3h8_percent +
                               c4h10_cnh2n_percent + h2o_percent + co_percent + co2_percent + so2_percent + o2_percent;
        setCH4(ch4_percent, ch4_percent / total_percent);
        setC2H6(c2h6_percent, c2h6_percent / total_percent);
        setN2(n2_percent, n2_percent / total_percent);
        setH2(h2_percent, h2_percent / total_percent);
        setC3H8(c3h8_percent, c3h8_percent / total_percent);
        setC4H10CnH2n(c4h10_cnh2n_percent, c4h10_cnh2n_percent / total_percent);
        setH2O(h2o_percent, h2o_percent / total_percent);
        setCO(co_percent, co_percent / total_percent);
        setCO2(co2_percent, co2_percent / total_percent);
        setSO2(so2_percent, so2_percent / total_percent);
        setO2(o2_percent, o2_percent / total_percent);
        setTotalCompositionWeight();
        setHeatingValue();
        setHeatingValueVolume();
        setSpecificGravity();
        setStoichometricAir();
    };
 
    double excessAirFromO2(double flue_gas_o2_percentage);
 
    double o2PercentageFromExcessAir(double excess_air);
 
    double estimateExcessAirFromO2(double flue_gas_o2_percentage);
 
    double calculateEnthalpyAtSaturation(double partial_pressure_water_vapor);
 
    double calculateSaturationTemperature(double partial_pressure_water_vapor);
 
    std::array<GasProperties*, 11> getConstituents() {
        return {&ch4, &c2h6, &n2, &h2, &c3h8, &c4h10_cnh2n, &h2o, &co, &co2, &so2, &o2};
    };
 
    double adjustedFlueGasO2ForCalcError(double excess_air, double flue_gas_o2);
 
  private:
    void setCH4(double composition_percent, double composition_by_volume) {
        ch4 = GasProperties(specificHeatCH4, kCh4MolecularWeight, kCh4SpecificWeight, composition_percent,
                    composition_by_volume, kCh4O2Generated, kCh4HeatingValue, kCh4HeatingValueVolume,
                    kCh4H2oGenerated, kCh4Co2Generated);
    }
 
    void setC2H6(double composition_percent, double composition_by_volume) {
        c2h6 = GasProperties(specificHeatC2H6, kC2h6MolecularWeight, kC2h6SpecificWeight, composition_percent,
                     composition_by_volume, kC2h6O2Generated, kC2h6HeatingValue, kC2h6HeatingValueVolume,
                     kC2h6H2oGenerated, kC2h6Co2Generated);
    }
 
    void setN2(double composition_percent, double composition_by_volume) {
        n2 = GasProperties(specificHeatN2, kN2MolecularWeight, kN2SpecificWeight, composition_percent,
                   composition_by_volume, 0, 0, 0, 0, 0);
    }
 
    void setH2(double composition_percent, double composition_by_volume) {
        h2 = GasProperties(specificHeatH2, kH2MolecularWeight, kH2SpecificWeight, composition_percent,
                   composition_by_volume, kH2O2Generated, kH2HeatingValue, kH2HeatingValueVolume,
                   kH2H2oGenerated, 0);
    }
 
    void setC3H8(double composition_percent, double composition_by_volume) {
        c3h8 = GasProperties(specificHeatC3H8, kC3h8MolecularWeight, kC3h8SpecificWeight, composition_percent,
                     composition_by_volume, kC3h8O2Generated, kC3h8HeatingValue, kC3h8HeatingValueVolume,
                     kC3h8H2oGenerated, kC3h8Co2Generated);
    }
 
    void setC4H10CnH2n(double composition_percent, double composition_by_volume) {
        c4h10_cnh2n = GasProperties(specificHeatC4H10CnH2n, kC4h10Cnh2nMolecularWeight, kC4h10Cnh2nSpecificWeight,
                        composition_percent, composition_by_volume, kC4h10Cnh2nO2Generated,
                        kC4h10Cnh2nHeatingValue, kC4h10Cnh2nHeatingValueVolume,
                        kC4h10Cnh2nH2oGenerated, kC4h10Cnh2nCo2Generated);
    }
 
    void setH2O(double composition_percent, double composition_by_volume) {
        h2o = GasProperties(specificHeatH2O, kH2oMolecularWeight, kH2oSpecificWeight, composition_percent,
                    composition_by_volume, 0, 0, 0, kH2oH2oGenerated, 0);
    }
 
    void setCO(double composition_percent, double composition_by_volume) {
        co = GasProperties(specificHeatCO, kCoMolecularWeight, kCoSpecificWeight, composition_percent,
                   composition_by_volume, kCoO2Generated, kCoHeatingValue, kCoHeatingValueVolume, 0,
                   kCoCo2Generated);
    }
 
    void setCO2(double composition_percent, double composition_by_volume) {
        co2 = GasProperties(specificHeatCO2, kCo2MolecularWeight, kCo2SpecificWeight, composition_percent,
                    composition_by_volume, 0, 0, 0, 0, kCo2Co2Generated);
    }
 
    void setSO2(double composition_percent, double composition_by_volume) {
        so2 = GasProperties(specificHeatSO2, kSo2MolecularWeight, kSo2SpecificWeight, composition_percent,
                    composition_by_volume, 0, 0, 0, 0, 0);
    }
 
    void setO2(double composition_percent, double composition_by_volume) {
        o2 = GasProperties(specificHeatO2, kO2MolecularWeight, kO2SpecificWeight, composition_percent,
                   composition_by_volume, kO2O2Generated, 0, 0, 0, 0);
    }
 
    void                 setTotalCompositionWeight();
    void                 setHeatingValue();
    void                 setHeatingValueVolume();
    void                 setSpecificGravity();
    void                 setStoichometricAir();
    GasFlueGasComponents getMassFlueGasComponents(double excess_air);
};
}; // namespace gas_composition