As some of you may know, I've been field testing bike stuff and positioning with a power meter for 4 or 5 years now. My method of choice is Robert Chung's "Virtual Elevation", or VE protocol, sometimes known as the "Chung Method". He has a great presentation on it here: http://anonymous.coward.free.fr/wattage/cda/indirect-cda.pdf . When I first read Robert's info, I wrote up a spreadsheet that I've used since then to analyze everything from position changes to tire air pressure effects. Of course, since I wrote that spreadsheet for my own use, it's not exactly the most "user friendly" (Hey, I know what I'm supposed to do, I wrote it! ;-)...but, don't worry, everyone else is in luck since Andy Froncioni (the main tech guy behind Alphamantis and the ERO facility that recently opened at the indoor track in Carson) added a version of the same calculations (called "Aerolab") to the freeware power meter analysis software, Golden Cheetah.
So, the question with this type of testing usually comes down to just how sensitive can it really be...especially as compared to something like a wind tunnel? Admittedly, there are some limitations to this type of testing, the main one being (at present time) that the results are mostly limited to zero yaw conditions, but as we saw in one of my previous blog posts, the most common yaw angles a TT'er or triathlete encounters are usually centered around zero yaw. Using the tool to make evaluations at zero yaw still can hold a significant benefit for someone interested in improving/testing bicycle aerodynamics.
A couple years ago, Dr. Andrew Coggan published a blog post titled "A Challenge to Cycling Aerodynamicists" in which he described a field test he undertook to take up something he coined the "Tom Compton Challenge". In short, it's an effort using known geometric shapes to try to determine the "sensitivity" of the aerodynamic field testing method.
Well, last year I discovered that my preferred field testing venue for VE runs had suffered some "traffic rerouting" that had made it much less appealing for the purpose (part of that "discovery" occurred when Andy sent his test setup to me to try and the results from my first course were very mixed due to excess vehicle interference after the nearby roads had been modified). So, I started scouting around for an alternative course and luckily found one that is much closer to my home (I can ride there in just a few minutes) and that has laps that are significantly shorter than the old course (shorter laps = shorter test run time). Both of these courses are best described as a sort of "extended halfpipe", an "out and back" course having a U-shaped elevation profile that allows for turnarounds to be taken at low speeds and thus avoid braking. Since identifying the new course, and having done just a few tests on it, one thing I wanted to do was to repeat the type of testing that Andy did and attempt to characterize the potential "sensitivity" of the course using the VE method.
Using Andy's setup as a guide, I set about figuring out what sorts of objects I could use for the test. I took a quick trip to the nearby Michael's craft store and acquired some styrofoam spheres, 2", 3", and 4" in diameter.
The hole in each sphere ended up being a nice friction fit, so swapping between spheres was a very simple process.
All runs were recorded with my trusty old yellow-cap PT Pro wheel mounted on the rear, with a cover in place to turn it into a de facto disc. I prefer to use a PT for my aero field testing since it eliminates the uncertainty of variations in drivetrain resistance across the gearing, plus the PT's "coasting zero" feature allows me to have the power meter zero while soft-pedaling (i.e. turning the pedals slowly while freewheeling) down the descents of the course at least once per lap. That helps to minimize any power meter drift during the runs.
So, with the test rig sorted out, it was time to head out to the test course and do some runs! To minimize wind and traffic effects, I prefer to head out to the course early on a Saturday or Sunday morning...before the small neighborhood that the course road services begins to wake up and starts moving around. A couple weekends ago, I headed out on a Sunday morning and rode over to my test course in 5 minutes, taking a small musette bag with the spheres, a notebook, and a couple of small tools I might need. Starting at 6 am, I did the runs in the following order:
- Rod only (w/washer "endplate")
- 3" sphere
- 2" sphere
- 4" sphere
Once I returned home, it was time to download the PT files into the computer and do the VE analysis. As I described above, although it can be done in GC's Aerolab feature, I prefer to use my own home-brewed spreadsheet, mostly because I find it easier to expand the vertical scale to get a better handle on the leveling procedure, but also because I've recently added a feature that varies the on-road Crr as a function of the ambient temperature. In order to use the spreadsheet, the following inputs are required:
- Total Mass - Easy to get just by stepping on a scale with bike in hand
- Weather Conditions - this means air temp, dew point temp, and barometric pressure (to determine air density). Luckily, there's a personal weather station listed on Weather Underground literally less than a block from my test course that has updates loaded every 5 minutes. Using that station also allows for a cross-check on ambient wind conditions to make sure it stayed calm during the test runs.
- Assumed Crr - For this, I use a weighted average of the front and rear tire Crr that I've determined from roller testing. The spreadsheet then compensates for the expected Crr due to the difference in the test ambient temperature and the 20C temperature to which my Crr results are normalized.
OK then, let's get to the results. Using the procedure outlined above, my best determination for the measured CdA from the runs was as follows (in order that runs were performed):
- Rod only = .2484 m^2
- 3" sphere = .2498 m^2
- 2" sphere = .2486 m^2
- 4" sphere = .2510 m^2
Run # - Sphere Measured CdA (m^2) Calculated CdA (m^2) Difference (m^2)
2. - 3" sphere .2498 .2501 .0003
3. - 2" sphere .2486 .2488 .0002
4. - 4" sphere .2510 .2518 .0008
Another way of looking at it is the expected and measured differences from the baseline:
Run # - Sphere Meas. Diff. from Baseline(m^2) Calc. Diff. from Baseline (m^2)
2. - 3" sphere .0014 .0017
3. - 2" sphere .0002 .0004
4. - 4" sphere .0026 .0033
One last way of looking at this is from the perspective of expected change from from the previous run. Here's how that worked out:
Run # - Sphere Meas. Diff. from Previous (m^2) Calc. Diff. from Previous (m^2)
2. - 3" sphere .0014 .0017
3. - 2" sphere -.0012 -.0013
4. - 4" sphere .0024 .0029
Not bad, huh? I've always said that when using this technique I consider measurements that are within +/-.001 m^2 to be basically "tied", and the above appears to bear that assumption out as being fairly conservative. It also gives me confidence that with careful technique I should be able to easily detect CdA differences on the order of .003-.005 m^2 and greater.
Nice to see that you got this working eventually!
One question for you: what was the diameter of the rod? I gather it must have been pretty fat if a 2" diameter sphere added only 0.0004 m^2 to the CdA? (I say that because the sphere alone would have added ~2.5x that.)
Hmm, I guess I can work that out for myself: if the sphere alone should have added 0.0010 m^2, but the sphere minus a 2" section of the rod added only 0.0004 m^2, then the CdA of that part of the rod must have been 0.0006 m^2. With an assumed CdA of 1.17, the A must have been 0.0005 m^2...which divided by 0.0508 m means that the rod must have been 0.0098 m in diameter, or 9.8 mm. Is that correct?
And now that I posted that, I see that the nominal diameter of the rod was 1/2" (12.7 mm)...so considering rounding errors, etc., close enough!ReplyDelete
The dowel actually measured 12.87mm in diameter, and the spheres covered up 35mm of it.Delete
Yeah...I'm glad I finally got around to it too...if you only knew how many times I set my alarm this summer to go out and do this and then hit the snooze button when it went off ;-)
Great read Tom. Is there any way you can show what your course looks like? I'm particularly interested in the turn-around portion. I'd like to see if I can find somewhere similar in my area.ReplyDelete
Nice. One day I'll get to play with Coggan's balls.ReplyDelete
excellent post. On the last line, do you mean " I should be able to easily detect AREA differences on the order of .003-.005 m^2 and greater." instead of CdA differences.ReplyDelete
Also, as a rough estimate, if a typical cyclists represents 0.3m^2 - 0.5m^2 front area, it would seem that you can detect about a 1% change in frontal area.
Nope...CdA. CdA, also known as "drag area", is the product of the dimensionless quantity Cd and area. Since it's typically difficult to separate out the 2, they are normally just reported as the product since it doesn't reduce the utility. Oh...and thanks!Delete