42 const double pressure = outletPressure - inletPressure;
43 const double referencePower = (
flowRate + 250) * (pressure + 0.4 * (airDensity / 0.075)) / (6343 * 0.66);
44 const double kWreferencePower = referencePower * 0.7457;
45 const double efficiency1 = 0.96 * std::pow(referencePower / (referencePower + 2.2), 0.05);
46 double efficiency2 = 0.962;
47 if (kWreferencePower < 185) {
48 efficiency2 = -0.003812 * std::pow(std::log10(kWreferencePower), 4) +
49 0.025834 * std::pow(std::log10(kWreferencePower), 3) +
50 -0.072577 * std::pow(std::log10(kWreferencePower), 2) +
51 0.125559 * std::log10(kWreferencePower) + 0.850274;
54 return (kWreferencePower * (1 / efficiency1) * (1 / efficiency2)) / motorPower;
66 const double flowRate, inletPressure, outletPressure, airDensity, motorPower;