model.cpp

#include "model.h"
#include <cmath>

double calculatePerformance(double T,
                            double Tmin,
                            double k)
{
    if (T < Tmin)
        return 0.0;

    return 100.0 * (1.0 - std::exp(-k * (T - Tmin)));
}

double calculateEnergy(double T,
                       double massFlowRate,
                       double Cp,
                       double Tambient,
                       double efficiency)
{
    double deltaT = std::max(T - Tambient, 0.0);

    return massFlowRate * Cp * deltaT / efficiency;
}

double calculateCost(double energy,
                     double electricityPrice)
{
    return energy * electricityPrice;
}

double calculateEmissions(double energy,
                          double emissionFactor)
{
    return energy * emissionFactor;
}

ProcessResult evaluateCondition(
    double T,
    double massFlowRate,
    double electricityPrice,
    double emissionFactor)
{
    ProcessResult result;

    result.temperature = T;
    result.performance =
        calculatePerformance(T);

    result.energy =
        calculateEnergy(T, massFlowRate);

    result.cost =
        calculateCost(result.energy,
                      electricityPrice);

    result.emissions =
        calculateEmissions(result.energy,
                           emissionFactor);

    return result;
}