DecibelsMethod.h

Go to the documentation of this file.

#pragma once
 
class DecibelsMethod {
  public:
    struct Output {
        Output(const double leakRateEstimate, const double annualConsumption)
            : leakRateEstimate(leakRateEstimate), annualConsumption(annualConsumption) {}
 
        const double leakRateEstimate, annualConsumption;
    };
 
    DecibelsMethod(const double operatingTime, const double linePressure, const double decibels,
                   const double decibelRatingA, const double pressureA, const double firstFlowA,
                   const double secondFlowA, const double decibelRatingB, const double pressureB,
                   const double firstFlowB, const double secondFlowB)
        : operatingTime(operatingTime), linePressure(linePressure), decibels(decibels), decibelRatingA(decibelRatingA),
          pressureA(pressureA), firstFlowA(firstFlowA), secondFlowA(secondFlowA), decibelRatingB(decibelRatingB),
          pressureB(pressureB), firstFlowB(firstFlowB), secondFlowB(secondFlowB) {}
 
    Output calculate() {
        /*
        double operatingTime;
        double linePressure; // X
        double decibels; // Y
        double decibelRatingA; // Y1
        double pressureA; // X1
        double firstFlowA; // Q11
        double secondFlowA; // Q21
        double decibelRatingB; // Y2
        double pressureB; // X2
        double firstFlowB; // Q12
        double secondFlowB; // Q22
        */
 
        const double denominator = (pressureB - pressureA) * (decibelRatingB - decibelRatingA);
        const double leakRateEstimate =
            ((pressureB - linePressure) * (decibelRatingB - decibels)) / denominator * firstFlowA +
            ((linePressure - pressureA) * (decibelRatingB - decibels)) / denominator * secondFlowA +
            ((pressureB - linePressure) * (decibels - decibelRatingA)) / denominator * firstFlowB +
            ((linePressure - pressureA) * (decibels - decibelRatingA)) / denominator * secondFlowB;
        const double           annualConsumption = (leakRateEstimate * operatingTime * 60) / 1000;
        DecibelsMethod::Output output(leakRateEstimate, annualConsumption);
 
        return output;
    }
 
  private:
    double operatingTime;
    double linePressure;   // X
    double decibels;       // Y
    double decibelRatingA; // Y1
    double pressureA;      // X1
    double firstFlowA;     // Q11
    double secondFlowA;    // Q21
    double decibelRatingB; // Y2
    double pressureB;      // X2
    double firstFlowB;     // Q12
    double secondFlowB;    // Q22
};