Template:LSU CT99CommonTheme2A

From VistrailsWiki
Jump to navigation Jump to search
The printable version is no longer supported and may have rendering errors. Please update your browser bookmarks and please use the default browser print function instead.

Referencing (equivalently, Version 1 of) the above-identified integral expression for the Gravitational Potential of an Axisymmetric Mass Distribution, Trova, Huré & Hersant (2012) offer the following assessment in §6 of their paper:

"The important question we have tried to clarify concerns the possibility of converting the remaining double integral … into a line integral … this question remains open."

We also have wondered whether there is a possibility of converting the double integral in this Key Equation into a single (line) integral. This is a particularly challenging task when, as is the case with Version 1 of the expression, the integrand is couched in terms of cylindrical coordinates because the modulus of the elliptic integral is explicitly a function of both <math>~\varpi^'</math> and <math>~z^'</math>.

We have realized that if we focus, instead, on Version 2 of the expression and associate the meridional-plane coordinates of the anchor ring with the coordinates of the location where the potential is being evaluated — that is, if we set <math>~(\varpi_a, z_a) = (\varpi,z)</math> — the argument of the elliptic integral becomes,

<math>~\mu = \biggl[\frac{2}{1+\coth\eta^' }\biggr]^{1 / 2} \, ,</math>

while the integral expression for the gravitational potential becomes,



<math>~ - 2^{3/2}G \varpi^2 \iint\limits_\mathrm{config}

\biggl[\frac{ \sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} \biggl[\frac{1}{1+\coth\eta^' }\biggr]^{1 / 2} K(\mu) \rho(\eta^', \theta^') d\eta^' d\theta^' </math>



<math>~- 2^{3 / 2} G \varpi^{2} \int\limits^{\eta_\mathrm{min}}_{\eta_\mathrm{max}} \frac{K(\mu) \sinh \eta^' ~d\eta^'}{( \sinh \eta^' +\cosh \eta^' )^{1 / 2}} \int\limits_{\theta_\mathrm{min}(\eta)}^{\theta_\mathrm{max}(\eta)} \rho(\eta^', \theta^') \biggl[ \frac{d\theta^'}{(\cosh\eta^' - \cos\theta^')^{5 / 2}} \biggr] \, . </math>

Notice that, by adopting this strategy, the argument of the elliptical integral is a function only of one coordinate — the toroidal coordinate system's radial coordinate, <math>~\eta^'</math>. As result, the integral over the angular coordinate, <math>~\theta^'</math>, does not involve the elliptic integral function. Then — as is shown in an accompanying chapter titled, Attempt at Simplification — if the configuration's density is constant, the integral over the angular coordinate variable can be completed analytically. Hence, the task of evaluating the gravitational potential (both inside and outside) of a uniform-density, axisymmetric configuration having any surface shape has been reduced to a problem of carrying out a single, line integration. This provides an answer to the question posed by Trova, Huré & Hersant (2012).