Next: B.3 Moment of inertia Up: Appendix B: Previous: B.1 Derivatives of the surface

B.2 Surface normal and surface vector quantities

The components of the unit surface normal $\hat{n}$ are


nx = $\displaystyle S^{-1} \left [ r r_{\phi} \sin{\phi} - r r_{\theta} \sin{\theta} \cos{\theta}
\cos{\phi} + r^2 \sin^2{\theta} \cos{\phi} \right ]$  
ny = $\displaystyle S^{-1} \left [-r r_{\phi} \cos{\phi} - r r_{\theta} \sin{\theta} \cos{\theta}
\sin{\phi} + r^2 \sin^2{\theta} \sin{\phi} \right ]$  
nz = $\displaystyle S^{-1} \left [r r_{\theta} \sin^2{\theta} + r^2 \cos{\theta} \sin{\theta} \right ]$ (33)

and the derivatives of the surface vector $\vec{X}$ are
$\displaystyle \frac{\partial \vec{X}}{\partial \phi}$ = $\displaystyle \left (r_{\phi}\sin{\theta}\cos{\phi} - r \sin{\theta}\sin{\phi},...
...n{\theta}\sin{\phi} + r
\sin{\theta}\cos{\phi},\: r_{\phi}\cos{\theta} \right )$ (34)
$\displaystyle \frac{\partial \vec{X}}{\partial \theta}$ = $\displaystyle \left (r \cos{\theta}\cos{\phi} + r_{\theta}\sin{\theta}\cos{\phi...
...eta}\sin{\theta}\sin{\phi},\: r_{\theta}\cos{\theta} - r \sin{\theta} \right ).$ (35)

The derivatives of the components of $\hat{n}$ with respect to $\phi$ are given by the following formula, where i = x, y,z:


\begin{displaymath}\frac{\partial n_i}{\partial \phi} = S^{-1} \left [ a_i - b_i \left (
\frac{r_{\phi}}{r} + \frac{c_i}{S^2} \right ) \right ]
\end{displaymath} (36)

where the functions ai, bi, and ci are:


a1 = $\displaystyle r^2_{\phi} \sin{\phi}
+ r r_{\phi \phi} \sin{\phi} + r r_{\phi} \...
...s{\theta} \cos{\phi} - r r_{\theta \phi} \sin{\theta} \cos{\theta}
\cos{\phi} +$  
    $\displaystyle r r_{\theta} \sin{\theta} \cos{\theta} \sin{\phi} + 2 r r_{\phi}
\sin^2{\theta} \cos{\phi} - r^2 \sin^2{\theta} \sin{\phi}$  
b1 = $\displaystyle r r_{\phi}
\sin{\phi} - r r_{\theta} \sin{\theta} \cos{\theta} \cos{\phi} +
r^2 \sin^2{\theta} \cos{\phi}$  
c1 = $\displaystyle r^2 (r_{\phi} r_{\phi \phi} + r_{\theta} r_{\theta \phi} \sin^2{\theta}
+r r_{\phi} \sin^2{\theta}$  
a2 = $\displaystyle - r^2_{\phi} \cos{\phi}
- r r_{\phi \phi} \cos{\phi} + r r_{\phi}...
...s{\theta} \sin{\phi} - r r_{\theta \phi} \sin{\theta} \cos{\theta}
\sin{\phi} -$  
    $\displaystyle r r_{\theta} \sin{\theta} \cos{\theta} \cos{\phi} + 2 r r_{\phi}
\sin^2{\theta} \sin{\phi} + r^2 \sin^2{\theta} \cos{\phi}$  
b2 = $\displaystyle - r r_{\phi}
\cos{\phi} - r r_{\theta} \sin{\theta} \cos{\theta} \sin{\phi} +
r^2 \sin^2{\theta} \sin{\phi}$  
c2 = c1  
a3 = $\displaystyle r_{\phi} r_{\theta} \sin^2{\theta}
+ r r_{\theta \phi} \sin^2{\theta} + 2 r r_{\phi} \sin{\theta}\cos{\theta}$  
b3 = $\displaystyle r r_{\theta}\sin^2{\theta} + r^2 \sin{\theta}\cos{\theta}$  
c3 = c1 (37)

The derivatives of the components of $\hat{n}$ with respect to $\theta $ are given by a formula similar to that for $\phi$:


\begin{displaymath}\frac{\partial n_i}{\partial \theta} = S^{-1} \left [ a_i - b...
...ft (
\frac{r_{\theta}}{r} + \frac{c_i}{S^2} \right ) \right ]
\end{displaymath} (38)

where the functions ai, bi, and ci are:
a1 = $\displaystyle r_{\theta} r_{\phi} \sin{\phi} + r r_{\theta \phi} \sin{\phi}
-r^...
...\cos{\theta}\cos{\phi} - r r_{\theta \theta} \sin{\theta}\cos{\theta}\cos{\phi}$  
    $\displaystyle - r r_{\theta}
\cos^2{\theta}\cos{\phi} + r r_{\theta} \sin^2{\th...
... r_{\theta} \sin^2{\theta}\cos{\phi} + 2 r^2 \sin{\theta}\cos{\theta}\cos{\phi}$  
b1 = $\displaystyle r r_{\phi}
\sin{\phi} - r r_{\theta} \sin{\theta} \cos{\theta} \cos{\phi} +
r^2 \sin^2{\theta} \cos{\phi}$  
c1 = $\displaystyle r_{\phi} r_{\theta \phi} + r_{\theta} r_{\theta \theta} \sin^2{\t...
...{\theta} + r^2_{\theta} \sin{\theta}\cos{\theta} +
r^2 \sin{\theta}\cos{\theta}$  
a2 = $\displaystyle -r_{\phi} r_{\theta} \cos{\phi} - r r_{\theta \phi} \cos{\phi}
-r...
...}\cos{\theta}\sin{\phi} -r r_{\theta \theta} \sin{\theta}\cos{\theta}\sin{\phi}$  
    $\displaystyle - r r_{\theta}
\cos^2{\theta}\sin{\phi} + r r_{\theta} \sin^2{\th...
... r_{\theta} \sin^2{\theta}\sin{\phi} + 2 r^2 \sin{\theta}\cos{\theta}\sin{\phi}$  
b2 = $\displaystyle - r r_{\phi}
\cos{\phi} - r r_{\theta} \sin{\theta} \cos{\theta} \sin{\phi} +
r^2 \sin^2{\theta} \sin{\phi}$  
c2 = c1  
a3 = $\displaystyle r^2_{\theta} \sin^2{\theta} + r r_{\theta \theta} \sin^2{\theta}
+2 r r_{\theta} \sin{\theta}\cos{\theta}$  
    $\displaystyle + 2 r r_{\theta} \sin{\theta}\cos{\theta} + r^2 \cos^2{\theta}
- r^2 \sin^2{\theta}$  
b3 = $\displaystyle r r_{\theta}\sin^2{\theta} + r^2 \sin{\theta}\cos{\theta}$  
c3 = c1 (39)


Next: B.3 Moment of inertia Up: Appendix B: Previous: B.1 Derivatives of the surface