In the early 1970s, Gerald Gregorek wrote a technical note on drag for ESTES (maker of model rockets and parts from what I gather).
The note is still very popular in the model rocket community, and it is referenced in most papers written on model rockets I've encountered. Here is one link to it (please replace with a web archive link if it ever breaks).
One equation in Gregorek's note is especially popular, and that is the equation for drag from the combined nose and body tube assembly. Here it is, with some context:
These curves can be used for all 1" models with 3:1 ogive nose cones, but what about other shapes and model rocket diameters? For other rocket configurations we would have to run more wind tunnel tests or find a mathematical way to determine the drag coefficients. Luckily, aerodynamicists have been working on this problem for many years, so short hand equations are available for our use. The one equation which we can use for the nose cone and body tube is $$C_{D_N} + C_{D_{BT}} = 1.02\ C_f \left[ 1 + \frac{1.5}{(L/d)^{3/2}} \right] > \frac{S_W}{S_{BT}} \tag{8}$$
Here, CDN and CDBT are the drag coefficients due to the nose and body tube of the rocket; Cf is the skin friction coefficient, to be found from plots over a range of Reynolds numbers; L/d is the length-to-diameter ratio; and SW and SBT are the wetted surface area and the body tube reference area, with the transverse cross-sectional area assumed in the equation.
If you've struggled to calculate drag coefficients, you can see the appeal of this equation... But there is one very big problem.
Gregorek doesn't tell you where the equation comes from. There is no rationale for it (understandable, given his audience), and there is no reference for it. All the papers that reference this equation reference only Gregorek and nobody else. Presumably the equation works for some model rockets, but why?
Can anyone say where this equation might have come from? Does anyone in the aerodynamics community recognize it? I'm especially interested in the limits to this equation---when is it valid and when does it stop working? If my rocket is 200 ft high and 12 ft wide but all else is the same, will I still get an adequate ballpark figure for my drag coefficient (assuming I don't care about Mach number or angle of attack dependences, and assuming also that the Reynolds number dependence is adequately captured in the friction coefficient calculation)?
Any pointers to the true source of Gregorek's equation (or even to an alternate equation for bigger rockets) would be appreciated!
Edit to include the bits I found from Organic Marble's hugely helpful comment:
You can arrive at an expression very close to Gregorek's using the expressions given in Ch. 6 of Hoerner's Fluid Dynamic Drag book. See pp. 6-15 thru 6-19 ("Drag of Streamline Bodies"). There, Hoerner gives the total drag of a streamlined shape (like a rocket with an aerodynamic nose) based on wetted area (this matters, see below):
The total drag, based on wetted area is consequently $$C_\mathrm{Dwet} / C_f = 1 + 1.5 (d/l)^{3/2} + 7 (d/l)^3 \tag{28}$$
About the third term, Hoerner says this:
The graph shows that the third term of this equation is practically negligible up to d/l $\approx$ 0.2.
Model rockets are usually slender, and Gregorek's examples all have l/d ratios between 10 and 20, which is to say d/l ratios between 0.1 and 0.05. In this range, and given Gregorek's amateur audience, it would have made sense to drop the third term.
Further, because Hoerner's equation is based on wetted area, and not on the frontal cross-sectional area more commonly used in drag calculations, Gregorek would have had to scale Hoerner's equation by a factor equal to the ratio of the wetted area and the cross-sectional area.
In fact, Hoerner later factors in the ratio of wetted and frontal areas, though the expressions he obtains seem less versatile than Gregorek's, since they involve certain assumptions about the (rocket's) streamline shape, which Gregorek avoids.
The only thing we haven't accounted for in Gregorek's equation at this point is his 1.02 multiplier, which suggests the true drag coefficient is 2% larger than that predicted by the simplified Hoerner equation... Maybe the 2% discrepancy comes from the third term he dropped, and the 1.02 multiplier is a way to correct for it. Maybe someone else has a better idea?