So where did I go wrong?
Remember that the average active area is computed using an integral, by averaging out the projected cross section for a given direction over all directions. The projected cross section for each direction is the surface area of an ellipse with one axis equal to , while the other axis length depends on the angle. To compute this variable axis length, I first looked at the cross section of the spheroid in the plane of the incoming direction, which is an ellipse with axes and . The projection of this ellipse perpendicular to the viewing direction is a line segment, and the length of this line segment is equal to twice the axis length we are looking for.
This is depicted in the image below:
To find the line segment in question, we first need to find out what the tangent lines to the ellipse are (blue dashed lines). These are parallel to the viewing line (solid blue line), and go through points on the ellipse (red points). The coordinates of the red points can be found using some basic analytic geometry, as I did in my old post.
The error I made was in the next step: finding the length of the projected line segment. In my old post, I took this length to be the length of the segment in between the red points (red dotted line). However, as can be seen from the figure, the red line segment is not actually perpendicular to the viewing line, so this is incorrect. In fact, it is only correct if the ellipse were a circle, which explains why the old expression is correct for the case .
Instead of using the red line segment, what we need to determine is the length of the line segment along the orange dashed line and in between the two blue dashed lines. This is a lot harder than the simple distance expression in the old post, so I had to revert to some symbolic algebra using the Python package sympy to achieve this. I won’t go into details here, as that would require me to write down too many equations. Without equations, what you need to do is (a) construct the orange line as the line perpendicular to the bottom blue line and going through the bottom red point, (b) determine the intersection point of the orange line with the top blue line, (c) compute the distance between bottom red point and this intersection point.
Assuming the coordinates of the red points are given by and , the distance expression in my old post, which read
needs to be changed to
When we substitute the expressions for and found in my old post, this leads to the following expression for the unknown axis length:
This expression is actually more compact than the old expression. As before, plugging in yields , while plugging in yields , as expected.
Using this new expression, I find the following plot for the average active cross section as a function of axis ratio:
Both curves are now indistinguishable.