This function approximates the following correlation:
where y=η ⁄ ηp (note that η refers to the hydraulic efficiency instead of total efficiency), x=log10(Eu ⁄ Eup), with the subscript p denoting the condition where the mover is operating at peak efficiency, and
Z1=(x-a) ⁄ b
Z2=(ec⋅x⋅d⋅x-a) ⁄ b
Z3=-a ⁄ b
a=-2.732094
b=2.273014
c=0.196344
d=5.267518
The approximation uses two simple polynomials stitched together by
a third one of the same order.
Care has been taken to ensure that, on the curve constructed by
if
statements, the differences of dy ⁄ dx
evaluated by different groups of coefficients at the connecting points
(i.e. at x = - 0.5 and x = + 0.5) are less than 1E-14.
This way, the derivative is still continuous to the solver even if
the solver requires a precision of 1E-10 when there are nested loops.
The correlation and the approximation have the shape as shown below (plotted by Buildings.Fluid.Movers.BaseClasses.Validation.EulerCurve).
The modified dimensionless Euler number is defined as
Eu=(Δp⋅D4) ⁄ (ρ⋅V̇2)
where Δp is the fan pressure rise in Pa, D is the fan wheel outer diameter in m, ρ is the inlet air density in kg/m3, and V̇ is the volumetric flow rate in m3/s. Note that the units in the definition do not matter to this correlation because it is the ratio of the Euler numbers that is used. Since D is constant for the same mover and ρ is approximately constant for common HVAC applications, the Euler number ratio can be simplified to
Eu ⁄ Eup=(Δp⋅V̇p2) ⁄ (Δpp⋅V̇2)
For more information regarding the correlation curve refer to EnergyPlus 9.6.0 Engineering Reference chapter 16.4 equations 16.209 through 16.218. Note that the formula is simplified here from the source document.
The svg file for the correlation equation was generated on https://viereck.ch/latex-to-svg using this script.
function correlation extends Modelica.Icons.Function; input Real x "log10(Eu/Eu_peak)"; output Real y "eta/eta_peak"; end correlation;
if-then
statements and to use
smoothMax
rather than max
.