AirWaterCoolingUsingFlue.h

Go to the documentation of this file.

 
#ifndef TOOLS_SUITE_AIRWATERCOOLINGUSINGFLUE_H
#define TOOLS_SUITE_AIRWATERCOOLINGUSINGFLUE_H
 
#include "losses/gas_flue_gas_material.h"
 
class AirWaterCoolingUsingFlue {
  public:
    struct Output {
        Output(double excessAir, double flowFlueGas, double specHeat, double fracCondensed, double effThermal,
               double effThermalLH, double effLH, double heatRecovery, double sensibleHeatRecovery)
            : excessAir(excessAir), flowFlueGas(flowFlueGas), specHeat(specHeat), fracCondensed(fracCondensed),
              effThermal(effThermal), effThermalLH(effThermalLH), effLH(effLH), heatRecovery(heatRecovery),
              sensibleHeatRecovery(sensibleHeatRecovery) {}
 
        Output()         = default;
        double excessAir = 0, flowFlueGas = 0, specHeat = 0, fracCondensed = 0, effThermal = 0, effThermalLH = 0,
               effLH = 0, heatRecovery = 0, sensibleHeatRecovery = 0;
    };
 
    AirWaterCoolingUsingFlue() {}
 
    Output calculate(GasCompositions gasCompositions, const double heatInput, const double tempFlueGasInF,
                     const double tempFlueGasOutF, const double tempCombAirF, const double fuelTempF,
                     const double percO2, const double ambientAirTempF = 60, const double moistCombAir = 0) {
 
        const auto   res = gasCompositions.getProcessHeatProperties(tempFlueGasInF, percO2, tempCombAirF, fuelTempF,
                                                                    ambientAirTempF, moistCombAir);
        const double fracCondensed =
            (1 - (0.0000009 * pow(tempFlueGasOutF, 3.0136)) / (2.8082 - 0.1168 * percO2 * 100));
        const double effLH = (fracCondensed * 0.00935 * (1087 /*+ 0.467 * tempFlueGasOutF - tempCombAirF*/)) / 100;
        const double flowFlueGas =
            heatInput * (0.0763 * (0.6 + (res.stoichAir * (1 + res.excessAir)))) * 1000000 / res.heatValueFuel;
        const double sensibleHeatRecovery =
            res.specificHeat * (tempFlueGasInF - tempFlueGasOutF) * flowFlueGas / 1000000;
 
        return Output(res.excessAir, flowFlueGas, res.specificHeat, fracCondensed, res.availableHeat,
                      res.availableHeat + effLH, effLH, heatInput * effLH, sensibleHeatRecovery);
    }
};
 
#endif // TOOLS_SUITE_AIRWATERCOOLINGUSINGFLUE_H