Here’s my model so far …
Let h and h1 be the height of the windway exit and entrance in mm,
w and w1 be the width of the windway exit and entrance in mm,
L be the length of the windway (the windway, not the whole whistle) in mm,
rho (ρ) be the density of air in kg/m^3,
Q be the air flow in L/min,
P be the mouth pressure in mm H2O,
v be the air speed leaving the windway.
If we know the flow rate, we can calculate the air speed in m/s as:
v = Q / (0.060 * h * w)
(The 0.060 is to convert from L/min to mL/sec.)
If we know the mouth pressure, we can calculate the air speed in m/s as:
v = sqrt(2 * 9.807 * P / (rho * K))
(The 9.807 is to convert from mm H20 to Pascals.)
Here, K is the flow loss coefficient, a dimensionless constant that is specific to a whistle. It is the sum of three parts:
K = Kb + Kdw + Ke
Kb comes from Bernoulli’s principle. If the cross-sectional area of the mouth is A0 mm^2, and the area of the windway entrance is A1 = h1 * w1,
Kb = 1 - (A1/A0)^2
Kdw comes from the Darcy-Weisbach equation:
Kdw = fd * T * L/Dh
Here, fd is the Darcy friction factor, which is close enough to 0.04 for pretty much everything we’re interested in.
Dh is the hydraulic diameter of the windway entrance, Dh = 2h1w1/(h1+w1).
T is a taper factor. Assuming most of the taper happens in the height, and the taper is more-or-less linear:
T = (1 + h/h1)/2 * (1 + h^2/h1^2)/2
Ke is an empirical adjustment. As yet, I don’t have a handle on this value. For the calibrators and the Feadog it is around 0.35; for the old Gen, Killarney, and Mellow D, it is around zero. When Terry put the flowmeter after a calibrator instead of before, Ke was around zero, so Ke may be nothing more than compensation for some idiosyncrasy of the flowmeter.
For the discharge coefficient, Cd = 1/sqrt(K)
For the “resistance”, resistance = sqrt(P)/Q = 0.2657 * h * w / sqrt(rho * K) = 0.2657 * h * w * Cd / sqrt(rho)