mixtureFractionCalculator.hpp

#ifndef ABLATELIBRARY_MIXTUREFRACTIONCALCULATOR_HPP
#define ABLATELIBRARY_MIXTUREFRACTIONCALCULATOR_HPP
 
#include <map>
#include <memory>
#include <string>
#include <vector>
#include "eos/tChem.hpp"
#include "mathFunctions/fieldFunction.hpp"
 
namespace ablate::monitors {
 
class MixtureFractionCalculator {
   private:
    const std::shared_ptr<ablate::eos::ChemistryModel> eos;
 
    const std::vector<std::string> trackingElements;
 
    double zMixFuel, zMixOxidizer;
 
    std::vector<double> zMixCoefficients;
 
    static std::map<std::string, double> ToMassFractionMap(const std::shared_ptr<ablate::eos::EOS>& eos, const std::shared_ptr<ablate::mathFunctions::FieldFunction>& massFractions);
 
   public:
    MixtureFractionCalculator(const std::shared_ptr<ablate::eos::EOS>& eos, std::map<std::string, double> massFractionsFuel, std::map<std::string, double> massFractionsOxidizer,
                              const std::vector<std::string>& trackingElements = {});
 
    MixtureFractionCalculator(const std::shared_ptr<ablate::eos::EOS>& eos, const std::shared_ptr<ablate::mathFunctions::FieldFunction>& massFractionsFuel,
                              const std::shared_ptr<ablate::mathFunctions::FieldFunction>& massFractionsOxidizer, const std::vector<std::string>& trackingElements = {});
    template <class T>
    inline T Calculate(T* yi) const {
        double zMix = 0.;
        for (std::size_t s = 0; s < zMixCoefficients.size(); s++) {
            zMix += zMixCoefficients[s] * yi[s];
        }
        return (zMix - zMixOxidizer) / (zMixFuel - zMixOxidizer);
    }
 
    std::shared_ptr<ablate::eos::ChemistryModel> GetEos() { return eos; }
};
 
}  // namespace ablate::monitors
#endif  // ABLATELIBRARY_MIXTUREFRACTIONCALCULATOR_HPP