Revision as of 19:53, 5 March 2010 editRHB100 (talk | contribs)Extended confirmed users2,197 edits →Gravitational potential expanded in series of Legendre polynomials: To make more concise, removed f' and f'' evaluations.← Previous edit | Revision as of 14:14, 11 March 2010 edit undoSlawekb (talk | contribs)Extended confirmed users11,467 editsm →Gravitational potential expanded in series of Legendre polynomialsNext edit → | ||
(2 intermediate revisions by the same user not shown) | |||
Line 27: | Line 27: | ||
==Gravitational potential expanded in series of Legendre polynomials== | ==Gravitational potential expanded in series of Legendre polynomials== | ||
The potential at a point '''x''' is given by | |||
For the continuous case the potential, ''V'' at point b is computed as an integral over the distributed mass, | |||
:<math>V = - \ |
:<math>V(\mathbf{x}) = - \int_{\mathbb{R}^3} \frac{G}{|\mathbf{x}-\mathbf{y}|}\ dm(\mathbf{y}).</math> | ||
where dm equals density times differential volume, where '''x''' is a vector from the ] to the differential element of mass, where '''r''' is a vector from the differential element of mass to point b i.e. <math>\mathbf{r} = \mathbf{R}-\mathbf{x} \, </math>, | |||
, and where '''R''' is a vector from the center of mass to point b as shown in figure. | |||
] | |||
The potential can be expanded in a series of ]. For simplicity, assume that an origin is chosen to coincide with the center of mass of the system. Expanding the denominator in terms of the ] gives | |||
Expanding the denominator by taking the square root of the square, carrying out the dot product, and factoring R out of the denominator we get the more useful expression, | |||
:<math>\begin{align} | |||
V &= - \int_{\mathbb{R}}^3 \frac{G}{ \sqrt{|\mathbf{x}|^2 -2 \mathbf{x} \cdot \mathbf{y} + |\mathbf{y}|^2}}\,dm(\mathbf{y})\\ | |||
:<math>- \frac{1}{R}\int G \, / \, \sqrt{1 -2 cos( \beta ) \frac{x}{R} + \left( \frac{x}{R} \right)^2}\,dm. \ </math> V can be expanded in a series of Legendre Polynomials | |||
{}&=- \frac{1}{|\mathbf{x}|}\int_{\mathbb{R}^3} G \, \left/ \, \sqrt{1 -2 \cos( \theta ) \frac{r}{|\mathbf{x}|} + \left( \frac{r}{|\mathbf{x}|} \right)^2}\right.\,dm | |||
\end{align}</math> | |||
where in the last integral we have changed to ] (''r'',θ,φ) in which ''r'' is the distance to the origin and θ is the axial angle with respect to the axis in the direction of '''x'''. | |||
The integrand is precisely the generating function for the Legendre polynomials:<ref name="AEM">C. R. Wylie, Jr. 1960,''Advanced Engineering Mathematics'' (McGraw-Hill Book Company)</ref> | |||
Consider the function, | |||
:<math> f(\frac{x}{R}) = \left(1- 2 cos(\beta)\left(\frac{x}{R}\right)+ \left(\frac{x}{R}\right) ^2 \right) ^{- \frac{1}{2}} \ .</math> | |||
: |
:<math>\left(1- 2 X Z + Z^2 \right) ^{- \frac{1}{2}} \ = \sum_{n=0}^\infty Z^n P_n(X)</math> | ||
⚫ | |||
where the ''P''<sub>''n''</sub> are the Legendre polynomials of degree ''n''. So the integral can be expanded in a series of the form | |||
:<math> \begin{align} | |||
:<math> f(\frac{x}{R})\approx P_0(cos(\beta)) + \left(\frac{x}{R}\right) \ P_1(cos(\beta)) + \left(\frac{x}{R}\right)^2 \ P_2(cos(\beta)) \ </math>. | |||
⚫ | V(\mathbf{x}) &= - \frac{G}{|\mathbf{x}|} \int \sum_{n=0}^\infty \left(\frac{r}{|\mathbf{x}|} \right)^n P_n(\cos(\theta)) \, dm\\ | ||
⚫ | {}&= - \frac{G}{|\mathbf{x}|} \int \left(1 + \left(\frac{r}{|\mathbf{x}|}\right) \cos(\theta) + \left(\frac{r}{|\mathbf{x}|}\right)^2\frac {3 \cos(\theta)^2 - 1}{2} + \cdots\right)\,dm | ||
The appearance of the first three terms of the Taylor series expansion for <math>f(\frac{x}{R})</math> suggests that the infinite series might follow the same pattern. On page 454, section 10.8 of <ref name="AEM">C. R. Wylie, Jr. 1960,''Advanced Engineering Mathematics'' (McGraw-Hill Book Company)</ref> Theorem 2 which states that <math> | |||
\end{align}</math> | |||
\left(1- 2 X Z + Z^2 \right) ^{- \frac{1}{2}} \ = \sum_{n=0}^\infty Z^n P_n(X) \ </math>, is proven. Also see ]. Applied to the notation used here Theorem 2 states that | |||
The second term in the expansion <math>\int r\cos(\theta) dm</math> is the first moment of the mass distribution about the plane through the origin that is perpendicular to the axis '''x'''; this vanishes because the origin was chosen to coincide with the center of mass. So, bringing the integral under the sign of the summation gives | |||
:<math>\left(1- 2 cos(\beta)\left(\frac{x}{R}\right)+ | |||
:<math> V(\mathbf{x}) = - \frac{GM}{|\mathbf{x}|} - \frac{G}{|\mathbf{x}|} \int \left(\frac{r}{|\mathbf{x}|}\right)^2\frac {3 \cos(\theta)^2 - 1}{2} dm + \cdots</math> | |||
We thus get the result, | |||
⚫ | |||
where the <math>P_n (cos(\beta))</math> are Legendre Polynomials of degree n and argument <math>cos( \beta)\, .</math> | |||
This is a form which is quite useful for applications as shown in section 11.7 Gravitational Potential ... of <ref name="meirovitch">Leonard Meirovitch 1970, ''Methods of Analytical Dynamics'' (McGraw-Hill Book Company)</ref> . Taking integrals of terms in the series we get: | |||
:<math>V = - \frac{G}{R} \int (1 + \left(\frac{x}{R}\right) \ cos(\beta) + \left(\frac{x}{R}\right)^2\frac {3 cos(\beta)^2 - 1}{2} \ dm + ...</math> | |||
Note that <math>\int x \ cos(\beta) dm = 0 .</math> This is true since <math>x \ cos(\beta)</math> is the component of the location with respect to the point cm of the differential element of mass along the direction of '''R''' as can be seen from the figure. But the point cm is the location of the center of mass. Therefore we get the desired result that <math>\int x \ cos(\beta) dm = 0 \ </math> and that <math>\int \left(\frac{x}{R}\right) \ cos(\beta) dm = 0 \ </math> since R is invariant over the integral. We therefore get | |||
:<math> V = - \frac{GM}{R} - \frac{G}{R} \int \left(\frac{x}{R}\right)^2\frac {3 cos(\beta)^2 - 1}{2} dm + ... </math> | |||
==Spherical symmetry== | ==Spherical symmetry== | ||
A spherically symmetric mass distribution behaves to an observer completely outside the distribution as though all of the mass were concentrated at the center, and thus effectively as a ], by the ]. On the surface of the Earth, the acceleration is given by so-called ] ''g'', approximately 9.8 m/s<sup>2</sup>, although this value varies slightly with latitude and altitude: the magnitude of the acceleration is a little larger at the poles than at the equator because the Earth is an ]. | A spherically symmetric mass distribution behaves to an observer completely outside the distribution as though all of the mass were concentrated at the center, and thus effectively as a ], by the ]. On the surface of the Earth, the acceleration is given by so-called ] ''g'', approximately 9.8 m/s<sup>2</sup>, although this value varies slightly with latitude and altitude: the magnitude of the acceleration is a little larger at the poles than at the equator because the Earth is an ]. | ||
Within a spherically symmetric mass distribution, it is possible to solve ]. Within a uniform spherical body of radius ''R'' and density σ, the gravitational force ''g'' inside the sphere varies linearly with distance ''r'' from the center, giving the gravitational potential inside the sphere, which is | Within a spherically symmetric mass distribution, it is possible to solve ]. Within a uniform spherical body of radius ''R'' and density σ, the gravitational force ''g'' inside the sphere varies linearly with distance ''r'' from the center, giving the gravitational potential inside the sphere, which is<ref>{{harvnb|Marion|Thornton|2003|loc=§5.2}}</ref> | ||
:<math>V(r) = \frac {2}{3} \pi G \sigma (r^2-3R^2),\qquad r\leq R,</math> | :<math>V(r) = \frac {2}{3} \pi G \sigma (r^2-3R^2),\qquad r\leq R,</math> | ||
Line 84: | Line 75: | ||
*] | *] | ||
== |
==Notes== | ||
<references/> | <references/> | ||
==References== | |||
* {{cite web |url=http://www.mth.uct.ac.za/omei/gr/chap5/node4.html |title=Mass in Newtonian theory |accessdate=2009-03-25 |author=Peter Dunsby |work=Tensors and Relativity: Chapter 5 Conceptual Basis of General Relativity |publisher=Department of Mathematics and Applied Mathematics University of Cape Town |date=1996-06-15 }} | * {{cite web |url=http://www.mth.uct.ac.za/omei/gr/chap5/node4.html |title=Mass in Newtonian theory |accessdate=2009-03-25 |author=Peter Dunsby |work=Tensors and Relativity: Chapter 5 Conceptual Basis of General Relativity |publisher=Department of Mathematics and Applied Mathematics University of Cape Town |date=1996-06-15 }} | ||
* {{cite web |url=http://www.eas.slu.edu/People/LZhu/teaching/eas437/gravity.ppt |title=Gravity and Earth's Density Structure |accessdate=2009-03-25 |author=Lupei Zhu Associate Professor, Ph.D. (California Institute of Technology, 1998) |publisher=Saint Louis University (Department of Earth and Atmospheric Sciences)|work=EAS-437 Earth Dynamics}} | * {{cite web |url=http://www.eas.slu.edu/People/LZhu/teaching/eas437/gravity.ppt |title=Gravity and Earth's Density Structure |accessdate=2009-03-25 |author=Lupei Zhu Associate Professor, Ph.D. (California Institute of Technology, 1998) |publisher=Saint Louis University (Department of Earth and Atmospheric Sciences)|work=EAS-437 Earth Dynamics}} | ||
* {{cite web |url=http://surveying.wb.psu.edu/sur351/geoid/grava.htm |title=The Gravity Field of the Earth|accessdate=2009-03-25 |author=Charles D. Ghilani|work=The Physics Fact Book|publisher=Penn State Surveying Engineering Program|date=2006-11-28}} | * {{cite web |url=http://surveying.wb.psu.edu/sur351/geoid/grava.htm |title=The Gravity Field of the Earth|accessdate=2009-03-25 |author=Charles D. Ghilani|work=The Physics Fact Book|publisher=Penn State Surveying Engineering Program|date=2006-11-28}} | ||
*{{Citation | last1=Thornton | first1=Stephen T. | last2=Marion | first2=Jerry B. | title=Classical Dynamics of Particles and Systems | publisher=Brooks Cole | edition=5th | isbn=978-0-534-40896-1 | year=2003}}. | |||
*{{cite book |title=Postprincipia: Gravitation for Physicists and Astronomers | *{{cite book |title=Postprincipia: Gravitation for Physicists and Astronomers | ||
|first=Peter|last=Rastall|publisher=World Scientific|year=1991|isbn=9810207786|pages=7ff.}} | |first=Peter|last=Rastall|publisher=World Scientific|year=1991|isbn=9810207786|pages=7ff.}} |
Revision as of 14:14, 11 March 2010
In celestial mechanics, the gravitational potential belonging to an individual massive object, such as the earth or the sun, is a scalar field generated by that object, see also gravitational well. It is in every way analogous to the electric potential with mass playing the role of charge; the gravitational field is analogous to the electric field, and so are their mathematics and properties.
In mathematics the gravitational potential is also known as the Newtonian potential and is fundamental in the study of potential theory.
Mathematical form
The potential V at a distance r from a point mass of mass M is
where G is the gravitational constant. V has units of energy per unit mass. V is always negative where it is defined, and as r tends to infinity, it approaches zero.
The gravitational field, and thus the acceleration of a small body in the space around the massive object, is the negative gradient of the gravitational potential. Because the potential has no angular components, its gradient is:
where r is a vector of length r pointing from the point mass towards the small body and is a unit vector pointing from the point mass towards the small body. The magnitude of the acceleration therefore follows an inverse square law:
The potential associated with a mass distribution is the superposition of the potentials of point masses. If the mass distribution is a finite collection of point masses, and if the point masses are located at the points x1, ..., xn and have masses m1, ..., mn, then the potential of the distribution at the point x is:
If the mass distribution is given as a mass measure dm on three-dimensional Euclidean space R, then the potential is the convolution of −G/|y| with dm. In good cases this equals the integral
If there is a continuous function ρ(x) representing the density of the distribution at x, so that dm = ρ(x)dx, where dx is the Euclidean volume element, then the gravitational potential is
If V is a potential function coming from a continuous mass distribution ρ(x), then ρ can be recovered using the Laplace operator Δ using the formula:
This holds pointwise whenever ρ is continuous and is zero outside of a bounded set. In general, the mass measure dm can be recovered in the same way if the Laplace operator is taken in the sense of distributions. Consequently the gravitational potential satisfies Poisson's equation. See also Green's function for the three-variable Laplace equation and Newtonian potential.
Gravitational potential expanded in series of Legendre polynomials
The potential at a point x is given by
The potential can be expanded in a series of Legendre polynomials. For simplicity, assume that an origin is chosen to coincide with the center of mass of the system. Expanding the denominator in terms of the dot product gives
where in the last integral we have changed to spherical coordinates (r,θ,φ) in which r is the distance to the origin and θ is the axial angle with respect to the axis in the direction of x.
The integrand is precisely the generating function for the Legendre polynomials:
where the Pn are the Legendre polynomials of degree n. So the integral can be expanded in a series of the form
The second term in the expansion is the first moment of the mass distribution about the plane through the origin that is perpendicular to the axis x; this vanishes because the origin was chosen to coincide with the center of mass. So, bringing the integral under the sign of the summation gives
Spherical symmetry
A spherically symmetric mass distribution behaves to an observer completely outside the distribution as though all of the mass were concentrated at the center, and thus effectively as a point mass, by the shell theorem. On the surface of the Earth, the acceleration is given by so-called standard gravity g, approximately 9.8 m/s, although this value varies slightly with latitude and altitude: the magnitude of the acceleration is a little larger at the poles than at the equator because the Earth is an oblate spheroid.
Within a spherically symmetric mass distribution, it is possible to solve Poisson's equation in spherical coordinates. Within a uniform spherical body of radius R and density σ, the gravitational force g inside the sphere varies linearly with distance r from the center, giving the gravitational potential inside the sphere, which is
which smoothly connects to the potential function for the outside of the sphere (see the figure at the top).
Potential energy
The gravitational potential (V) is the potential energy (U) per unit mass:
where m is the mass of the object. The potential energy is the negative of the work done by the gravitational field moving the body to its given position in space from infinity. If the body has a mass of 1 unit, then the potential energy to be assigned to that body is equal to the gravitational potential. So the potential can be interpreted as the negative of the work done by the gravitational field moving a unit mass in from infinity.
In some situations the equations can be simplified by assuming a field which is nearly independent of position. For instance, in daily life, in the region close to the surface of the Earth, the gravitational acceleration can be considered constant. In that case the difference in potential energy from one height to another is to a good approximation linearly related to the difference in height:
General relativity
In general relativity, the gravitational potential is replaced by the metric tensor.
See also
Notes
- Vladimirov 1984, §7.8 harvnb error: no target: CITEREFVladimirov1984 (help)
- C. R. Wylie, Jr. 1960,Advanced Engineering Mathematics (McGraw-Hill Book Company)
- Marion & Thornton 2003, §5.2 harvnb error: no target: CITEREFMarionThornton2003 (help)
References
- Peter Dunsby (1996-06-15). "Mass in Newtonian theory". Tensors and Relativity: Chapter 5 Conceptual Basis of General Relativity. Department of Mathematics and Applied Mathematics University of Cape Town. Retrieved 2009-03-25.
- Lupei Zhu Associate Professor, Ph.D. (California Institute of Technology, 1998). "Gravity and Earth's Density Structure". EAS-437 Earth Dynamics. Saint Louis University (Department of Earth and Atmospheric Sciences). Retrieved 2009-03-25.
{{cite web}}
: CS1 maint: multiple names: authors list (link) CS1 maint: numeric names: authors list (link) - Charles D. Ghilani (2006-11-28). "The Gravity Field of the Earth". The Physics Fact Book. Penn State Surveying Engineering Program. Retrieved 2009-03-25.
- Thornton, Stephen T.; Marion, Jerry B. (2003), Classical Dynamics of Particles and Systems (5th ed.), Brooks Cole, ISBN 978-0-534-40896-1.
- Rastall, Peter (1991). Postprincipia: Gravitation for Physicists and Astronomers. World Scientific. pp. 7ff. ISBN 9810207786.
- Vladimirov, V. S. (1971), Equations of mathematical physics, Translated from the Russian by Audrey Littlewood. Edited by Alan Jeffrey. Pure and Applied Mathematics, vol. 3, New York: Marcel Dekker Inc., MR0268497.