Volcano-watching

Calculus Level 4

The surface of a volcano is represented by the frustum of the paraboloid z = 225 x 2 y 2 z = 225 - x^2 - y^2 for 0 z 224. 0 \le z \le 224. A geological research station is located at point A = ( 1 , 0 , 224 ) A=(1,0,224) on the edge of the crater of the volcano. The main research center for the volcano is located at point B = ( 15 , 0 , 0 ) . B=(15,0,0).

The scientists studying the volcano have established a path from A A to B B that passes all the way around the volcano exactly once—so that measurements can be taken of all sides of the volcano—and has the shortest possible length just in case the volcano shows signs of erupting.

What is the angle α 0 \alpha_0 that the tangent vector t \mathbf{t} to the path makes with the horizontal tangent vector u \mathbf{u} to the surface of the volcano at point A ? A? Give your answer in degrees, correct to 1 1 decimal place.

Note: This integral will (in)definitely help.


The answer is 15.4.

This section requires Javascript.
You are seeing this because something didn't load right. We suggest you, (a) try refreshing the page, (b) enabling javascript if it is disabled on your browser and, finally, (c) loading the non-javascript version of this page . We're sorry about the hassle.

1 solution

Mark Hennings
Apr 17, 2018

It makes things simpler to perform the analysis a little more generally. Suppose that the equation of the surface of the volcano (in cylindrical polar coordinates r , θ , z r,\theta,z ) is z = R 2 r 2 ρ r R z \; = \; R^2 - r^2 \hspace{2cm} \rho \le r \le R Then the vector equation of a path from A A to B B is given as r = ( r cos θ r sin θ R 2 r 2 ) \mathbf{r} \; = \; \left(\begin{array}{c} r\cos\theta \\ r\sin\theta \\ R^2 - r^2 \end{array} \right) where r = r ( θ ) r = r(\theta) is a function such that ρ r ( θ ) R \rho \le r(\theta) \le R for all 0 θ 2 π 0 \le \theta \le 2\pi , and r ( 0 ) = ρ r(0) = \rho and r ( 2 π ) = R r(2\pi) = R . Then d r d θ = r ( θ ) ( cos θ sin θ 2 r ) + r ( sin θ cos θ 0 ) \frac{d\mathbf{r}}{d\theta} \; = \; r'(\theta) \left(\begin{array}{c} \cos\theta \\ \sin\theta \\ -2r\end{array}\right) + r\left(\begin{array}{c} -\sin\theta \\ \cos\theta \\ 0 \end{array} \right) and so d s d θ = ( 1 + 4 r 2 ) ( r ) 2 + r 2 \frac{ds}{d\theta} \; = \; \sqrt{(1 + 4r^2)(r')^2 + r^2} where s s is the arc-length. If t \mathbf{t} is the unit tangent for the path, then d s d θ t = d r d θ \frac{ds}{d\theta}\,\mathbf{t} \; = \; \frac{d\mathbf{r}}{d\theta} and so the angle α \alpha between the unit tangent of the path and the horizontal tangent to the surface is given by the equation cos α = t ( sin θ cos θ 0 ) = r ( 1 + 4 r 2 ) ( r ) 2 + r 2 \cos\alpha \; = \; \mathbf{t} \cdot \left(\begin{array}{c} -\sin\theta \\ \cos\theta \\ 0 \end{array}\right) \; = \; \frac{r}{\sqrt{(1 + 4r^2)(r')^2 + r^2}} Thus we want to find the function r = r ( θ ) r= r(\theta) which is suitably smooth and such that r ( 0 ) = ρ r(0) = \rho , r ( 2 π ) = R r(2\pi) = R and such that S [ r ] = 0 2 π ( 1 + 4 r 2 ) ( r ) 2 + r 2 d θ = 0 2 π F d θ S[r] \; = \; \int_0^{2\pi} \sqrt{(1 + 4r^2)(r')^2 + r^2}\,d\theta \; = \; \int_0^{2\pi} F\,d\theta is minimized. Since F = ( 1 + 4 r 2 ) ( r ) 2 + r 2 F \; = \; \sqrt{(1 + 4r^2)(r')^2 + r^2} does not depend explicity on θ \theta , a standard result of the Calculus of Variations is that c = F r F r = r 2 ( 1 + 4 r 2 ) ( r ) 2 + r 2 = r cos α c \; = \; F - r'\frac{\partial F}{\partial r'} \; = \; \frac{r^2}{\sqrt{(1 + 4r^2)(r')^2 + r^2}} \; = \; r\cos\alpha is constant (this relationship between r r and α \alpha is known as Clairault's Relation). At the point A A , we note that c = ρ cos α 0 c = \rho \cos\alpha_0 , where α 0 \alpha_0 is the angle that we want to determine. Thus ρ < c < ρ -\rho < c < \rho .

Rearranging this equation, we deduce that d θ d r = c r 4 r 2 + 1 r 2 c 2 \frac{d\theta}{dr} \; = \; \frac{c}{r} \sqrt{\frac{4r^2+1}{r^2 - c^2}} and a similar set of calculations to those performed in the hint question show that θ = 2 c cosh 1 ( 1 + 4 r 2 1 + 4 c 2 ) + tan 1 ( r 2 c 2 c 1 + 4 r 2 ) + k \theta \; = \; 2c \cosh^{-1}\left(\sqrt{\frac{1 + 4r^2}{1 + 4c^2}}\right) + \tan^{-1}\left(\frac{\sqrt{r^2 - c^2}}{c\sqrt{1 + 4r^2}}\right) + k for an arbitrary constant k k , and so we deduce that θ = 2 c [ cosh 1 ( 1 + 4 r 2 1 + 4 c 2 ) cosh 1 ( 1 + 4 ρ 2 1 + 4 c 2 ) ] + tan 1 ( r 2 c 2 c 1 + 4 r 2 ) tan 1 ( ρ 2 c 2 c 1 + 4 ρ 2 ) \theta \; = \; 2c\left[\cosh^{-1}\left(\sqrt{\frac{1 + 4r^2}{1 + 4c^2}}\right) - \cosh^{-1}\left(\sqrt{\frac{1 + 4\rho^2}{1 + 4c^2}}\right)\right] + \tan^{-1}\left(\frac{\sqrt{r^2 - c^2}}{c\sqrt{1 + 4r^2}}\right) - \tan^{-1}\left(\frac{\sqrt{\rho^2 - c^2}}{c\sqrt{1 + 4\rho^2}}\right) where c c is to be chosen so that A ( c ) = 2 π A(c) = 2\pi , where A ( c ) = 2 c [ cosh 1 ( 1 + 4 R 2 1 + 4 c 2 ) cosh 1 ( 1 + 4 ρ 2 1 + 4 c 2 ) ] + tan 1 ( R 2 c 2 c 1 + 4 R 2 ) tan 1 ( ρ 2 c 2 c 1 + 4 ρ 2 ) A(c) \; = \; 2c\left[\cosh^{-1}\left(\sqrt{\frac{1 + 4R^2}{1 + 4c^2}}\right) - \cosh^{-1}\left(\sqrt{\frac{1 + 4\rho^2}{1 + 4c^2}}\right)\right] + \tan^{-1}\left(\frac{\sqrt{R^2 - c^2}}{c\sqrt{1 + 4R^2}}\right) - \tan^{-1}\left(\frac{\sqrt{\rho^2 - c^2}}{c\sqrt{1 + 4\rho^2}}\right) Since A ( c ) = 2 [ cosh 1 ( 1 + 4 R 2 1 + 4 c 2 ) cosh 1 ( 1 + 4 ρ 2 1 + 4 c 2 ) ] + 1 + 4 ρ 2 ρ 2 c 2 1 + 4 R 2 R 2 c 2 A'(c) \; = \; 2\left[\cosh^{-1}\left(\sqrt{\frac{1 + 4R^2}{1 + 4c^2}}\right) - \cosh^{-1}\left(\sqrt{\frac{1 + 4\rho^2}{1 + 4c^2}}\right)\right] + \sqrt{\frac{1 + 4\rho^2}{\rho^2 - c^2}} - \sqrt{\frac{1 + 4R^2}{R^2 - c^2}} it is easy to see that A ( c ) > 0 A'(c) > 0 for all ρ < c < ρ -\rho < c < \rho , and so the desired value of c c is unique.

With ρ = 1 \rho = 1 and R = 15 R = 15 , we calculate numerically that c = 0.963915 c = 0.963915 , and hence that α 0 = cos 1 c = 15.438 8 \alpha_0 = \cos^{-1}c = 15.4388^\circ . To 1 1 decimal place, the answer is 15. 4 \boxed{15.4^\circ} .

It is worth noting that, for the optimal curve, d s d θ = r 2 c \frac{ds}{d\theta} = \frac{r^2}{c} , and hence the length of the path A B AB is s = 0 2 π r 2 c d θ = 1 15 r 4 r 2 + 1 r 2 c 2 d r = 1 2 1 225 4 u + 1 u c 2 d u = 227.99 s \; =\; \int_0^{2\pi} \tfrac{r^2}{c}\,d\theta \; = \; \int_1^{15} r\sqrt{\tfrac{4r^2+1}{r^2-c^2}}\,dr \; =\; \tfrac12\int_1^{225} \sqrt{\tfrac{4u+1}{u-c^2}}\,du \; = \; 227.99 which is shorter than the special case (which would not be handled by the Calculus of Variations) of running around the lip of the crater, then going straight down the hill; this special path has length 2 π + 1 15 1 + 4 x 2 d x = 230.95 2\pi + \int_1^{15}\sqrt{1 + 4x^2}\,dx \; =\; 230.95

If anyone wants a good reference discussing this type of thing, see Do Carmo's Differential Geometry of Curves and Surfaces (in particular, pg 258-261 of the 2nd edition)

Christopher Criscitiello - 2 years, 10 months ago

Log in to reply

Do Carmo discusses one aspect of "this type of thing". The book looks at geodesics on surfaces, which are certainly what are being studied here. However, the Calculus of Variations is a field which has wider application than geodesics...

Mark Hennings - 2 years, 10 months ago

After what seemed like way too many mistakes, I did get the answer. However, I forgot about the hint and used Wolfram to solve the integral. Wolfram gives the result in terms of simple functions of r and the F 1 F_1 Appell hypergeometric function, so I went with that. Fortunately I found an accurate function for Matlab that allowed me calculate F 1 F_1 and iterate for c. Initially I was sure circling the crater and going straight down had be the shortest. Live and learn.

Ian Leslie - 2 years, 10 months ago

Could you please check the answer for the angle α 0 = 15.38 8 \alpha_0=15.388^\circ . Given c = 0.963915 c=0.963915 , I get α 0 = 15.43887 6 \alpha_0=15.438876^\circ

Darko Simonovic - 2 years, 10 months ago

Log in to reply

Well typo-spotted! I get c = 0.9639153338593377 c = 0.9639153338593377 and α 0 = 15.43880386978524 \alpha_0 = 15.43880386978524 .

Mark Hennings - 2 years, 10 months ago

0 pending reports

×

Problem Loading...

Note Loading...

Set Loading...