User:Tohline/Appendix/Mathematics/ToroidalFunctions
Relationships Between Toroidal Functions
| Tiled Menu | Tables of Content | Banner Video | Tohline Home Page | |
This chapter has been put together in an effort to lay the groundwork for an evaluation of Wong's (1973) derived expression for the gravitational potential both inside and outside of a uniform-density, axisymmetric torus. The interior (i.e., <math>~\eta^' > \eta_0</math>) solution is,
<math>~U(\eta^',\theta^')\biggr|_{\mathrm{for}~\eta^' \ge \eta_0}</math> |
<math>~=</math> |
<math>~ \frac{2^{3 / 2}}{3\pi^2} \biggl(\frac{q}{a}\biggr) \frac{\sinh^3 \eta_0}{\cosh\eta_0} \biggl\{ - \frac{3\pi^2}{2^{5/ 2}} \biggl[ \frac{\sinh^2\eta^'}{(\cosh \eta^' - \cos \theta^')^2} \biggr] +~ (\cosh \eta^' - \cos \theta^')^{1 / 2} </math> |
|
|
<math>~ \times \sum\limits_{n=0}^\infty \epsilon_n \cos(n\theta^') Q_{n-1 / 2}(\cosh\eta^') B_n(\cosh\eta_0) \biggr\} \, , </math> |
Wong (1973), Eq. (2.65) |
where,
<math>~B_n(\cosh\eta_0)</math> |
<math>~\equiv</math> |
<math>~ (n+\tfrac{1}{2})P_{n+1/2} (\cosh\eta_0)Q^2_{n-1/2} (\cosh\eta_0) - (n-\tfrac{3}{2})P_{n-1/2} (\cosh\eta_0)Q^2_{n+1/2} (\cosh\eta_0) \, . </math> |
Wong (1973), Eq. (2.62) |
Our Mucking Around
Begins on p. 332 of M. Abramowitz & I. A. Stegun (1995)
Recurrence Relations
Both <math>~P_\nu^\mu</math> and <math>~Q_\nu^\mu</math> satisfy the same recurrence relations.
<math>~P_\nu^{\mu+1}(z)</math> |
<math>~=</math> |
<math>~ (z^2 - 1)^{- 1 / 2} \biggl[ (\nu - \mu)z P_\nu^\mu(z) - (\nu + \mu)P_{\nu - 1}^\mu(z) \biggr] \, ; </math> |
<math>~(z^2-1) \frac{dP_\nu^{\mu}(z)}{dz}</math> |
<math>~=</math> |
<math>~ (\nu + \mu)(\nu - \mu +1_(z^2-1)^{1 / 2} P_\nu^{\mu - 1}(z) - \mu zP_\nu^\mu(z) \, ; </math> |
<math>~(\nu - \mu + 1)P_{\nu + 1}^{\mu}(z)</math> |
<math>~=</math> |
<math>~ (2\nu+1)zP_\nu^\mu (z) -(\nu + \mu)P_{\nu-1}^\mu(z) \, ; </math> |
<math>~(z^2-1) \frac{dP_\nu^{\mu}(z)}{dz}</math> |
<math>~=</math> |
<math>~ \nu z P_\nu^{\mu }(z) - (\nu + \mu)P_{\nu-1}^\mu(z) \, . </math> |
Toroidal Functions
Relation to Elliptic Integrals
<math>~P_{-1 / 2}(z)</math> |
<math>~=</math> |
<math>~ \frac{2}{\pi} \biggl[\frac{2}{z+1}\biggr]^{1 / 2} ~K\biggl( \sqrt{ \frac{z-1}{z+1}} \biggr) \, ; </math> |
<math>~P_{-1 / 2}(\cosh\eta)</math> |
<math>~=</math> |
<math>~ \biggl[\frac{\pi}{2}~\cosh\biggl(\frac{\eta}{2}\biggr)\biggr]^{-1} ~K\biggl( \tanh \frac{\eta}{2} \biggr) \, ; </math> |
<math>~Q_{-1 / 2}(z)</math> |
<math>~=</math> |
<math>~ \biggl[\frac{2}{z+1}\biggr]^{1 / 2} ~K\biggl( \sqrt{ \frac{2}{z+1} }\biggr) \, ; </math> |
<math>~Q_{-1 / 2}(\cosh\eta)</math> |
<math>~=</math> |
<math>~2 e^{- \eta / 2} ~K(e^{-\eta} ) \, ; </math> |
<math>~P_{+1 / 2}(z)</math> |
<math>~=</math> |
<math>~ \frac{2}{\pi} \biggl[ z + \sqrt{z^2-1} \biggr]^{1 / 2} ~E\biggl( \sqrt{ \frac{2(z^2-1)^{1 / 2}}{z + (z^2-1)^{1 / 2}}} \biggr) \, ; </math> |
<math>~P_{+ 1 / 2}(\cosh\eta)</math> |
<math>~=</math> |
<math>~ \frac{2}{\pi}~e^{\eta/2} ~E\biggl( \sqrt{1-e^{-2\eta}} \biggr) \, ; </math> |
<math>~Q_{+ 1 / 2}(z)</math> |
<math>~=</math> |
<math>~z \biggl[\frac{2}{z+1}\biggr]^{1 / 2} ~K\biggl( \sqrt{ \frac{2}{z+1} }\biggr) - \biggl[2(z+1)\biggr]^{1 / 2} E\biggl( \sqrt{\frac{2}{z+1}} \biggr) \, . </math> |
When the argument, <math>~x</math>, lies in the range, <math>~-1 < x < 1</math>:
<math>~P_{-1 / 2}(x)</math> |
<math>~=</math> |
<math>~ \frac{2}{\pi} ~K\biggl( \sqrt{ \frac{1-x}{2} } \biggr) \, ; </math> |
<math>~P_{-1 / 2}(\cos\theta)</math> |
<math>~=</math> |
<math>~ \frac{2}{\pi} ~K\biggl( \sin \frac{\theta}{2}\biggr) \, ; </math> |
<math>~Q_{-1 / 2}(x)</math> |
<math>~=</math> |
<math>~ K\biggl( \sqrt{ \frac{1+x}{2} } \biggr) \, ; </math> |
<math>~P_{+1 / 2}(x)</math> |
<math>~=</math> |
<math>~ \frac{2}{\pi} \biggl[2E\biggl( \sqrt{ \frac{1-x}{2} } \biggr) - ~K\biggl( \sqrt{ \frac{1-x}{2} } \biggr) \biggr] \, ; </math> |
<math>~Q_{+ 1 / 2}(x)</math> |
<math>~=</math> |
<math>~ K\biggl( \sqrt{ \frac{1+x}{2} } \biggr) - 2E\biggl( \sqrt{ \frac{1+x}{2} } \biggr)\, ; </math> |
Piece Together
When <math>~\mu = 0</math>, and <math>~\nu = (m+ 1 / 2)</math>, the recurrence relation should be …
<math>~(m + \tfrac{3}{2})P_{m+3 / 2}(z)</math> |
<math>~=</math> |
<math>~ (2m+2)zP_{m+ 1 / 2} (z) -(m + \tfrac{1}{2} )P_{m-1 / 2} (z) </math> |
<math>~\Rightarrow ~~~(2m + 3)P_{m+3 / 2}(z)</math> |
<math>~=</math> |
<math>~ 4(m+1)zP_{m + 1 /2 } (z) - (2m + 1)P_{m-1 / 2} (z) </math> |
<math>~\Rightarrow ~~~P_{m+3 / 2}(z)</math> |
<math>~=</math> |
<math>~\biggl[ \frac{ 4(m+1)zP_{m + 1 /2 } (z) - (2m + 1)P_{m-1 / 2} (z) }{(2m + 3)} \biggr] </math> |
Drawn From Discussion of Solving the Poisson Equation
The following has been copied (May 2018) from an accompanying chapter that presents the integral representation of the Poisson equation in terms of toroidal functions.
Table 5: Green's Function in Terms of Zero Order, Half-(Odd)Integer Degree, Associated Legendre Functions of the Second Kind, <math>~Q^0_{m-1 / 2}(\chi)</math> (also referred to as Toroidal Functions) |
||||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
where: <math>~\chi \equiv \frac{\varpi^2 + (\varpi^')^2 + (z - z^')^2}{2\varpi \varpi^'}</math>
H. S. Cohl & J. E. Tohline (1999), p. 88, Eqs. (15) & (16) |
||||||||||||||||||
Note that, according to, for example, equation (8.731.5) of Gradshteyn & Ryzhik (1994), <math>~Q^0_{-m - 1 / 2}(\chi) = Q^0_{m- 1 / 2}(\chi) \, .</math> Hence, the Green's function can straightforwardly be rewritten in terms of a simpler summation over just non-negative values of the index, <math>~m</math>. |
||||||||||||||||||
Referencing equations (8.13.3) and (8.13.7), respectively, of Abramowitz & Stegun (1965), we see that for the smallest two values of the non-negative index, <math>~m</math>, the function, <math>~Q_{m- 1 / 2}(\chi)</math>, can be rewritten in terms of, the more familiar, complete elliptic integrals of the first and second kind. Specifically,
Excerpt from p. 337 of M. Abramowitz & I. A. Stegun (1995) |
||||||||||||||||||
Finally, equation (8.5.3) from Abramowitz & Stegun (1965) or equation (8.832.4) of Gradshteyn & Ryzhik (1994) — also see equation (2) of Gil, Segura & Temme (2000) — provide the recurrence relation for all other values of the index, <math>~m</math>. Specifically, for all <math>~m \ge 2</math>, <math>~Q_{m - 1 / 2}(\chi) = 4\biggl[\frac{m-1}{2m-1}\biggr] \chi Q_{m- 3 / 2}(\chi) - \biggl[ \frac{2m-3}{2m-1}\biggr] Q_{m- 5 / 2}(\chi) \, .</math> Excerpt from p. 490 of W. Guatschi (1965, Communications of the ACM, 8, 488 - 492) |
See Also
© 2014 - 2021 by Joel E. Tohline |