AMApplied Mathematics2152-7385Scientific Research Publishing10.4236/am.2017.86064AM-77059ArticlesPhysics&Mathematics The Gravitational Attraction between Hemispherical Masses JonathanD. Weiss1*JSA Photonics LLC, 8019 Edith Blvd NE, Albuquerque, NM, USA* E-mail:jondweiss@aol.com08062017080682083411, May 201718, June 2017 21, June 2017© Copyright 2014 by authors and Scientific Research Publishing Inc. 2014This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/

This paper is a study of the gravitational attraction between two uniform hemispherical masses aligned such that the pair is cylindrically symmetric. Three variations are considered: flat side to flat side, curved side to curved side, and flat side to curved side. Expressions for the second and third variation are derived from the first, with the use of superposition and the well-known gravitational behavior of a spherical mass as equivalent to a point mass at its center. The study covers two masses of equal diameter and of different diameters, such that one is four times that of the other. Calculations are done for separations from zero to fifty times the radius of the larger of the two, which is effectively the asymptotic limit. It is demonstrated that at any separation, the force can be expressed as if the two hemispheres were point masses separated by a certain distance. Expressions for that distance and the location of the (fictitious) point masses within each hemisphere are presented. Unlike the case of two spherical masses, the location within their respective hemisphere is not necessarily the same for each point and both are dependent upon the separation between the two hemispheres. The calculation for the first variation is done in two ways. The first is a “brute force” multi-dimensional integral with the help of Wolfram Mathematica. The second is an axial expansion for the potential modified for off-axis locations by Legendre polynomials. With only a few terms in the expansion, the results of the second method are in extremely good agreement with those of the first. Finally, an interesting application to a split earth is presented.

Gravitational Attraction Hemispherical Solids Legendre PolynomialsSplit Earth
1. Introduction

With at least one exception, a perusal of online publications of the (Newtonian!) gravitational attraction between solid masses yields discussions of a pedagogical nature concerning spherically symmetric solids and shells. In a few cases, the force between a point mass and a disk or a hemisphere is derived for the point mass placed along the center line of the non-spherical object. The exception is a calculation of the force between two hemispheres with their flat surfaces in contact, as in the case of the earth  . In addition to considering pairs of non- spherically symmetric masses, this paper generalizes the discussion in  to the force between two hemispheres of varying separation, such that the pair forms a cylindrically-symmetric configuration. Furthermore, two additional relative orientations are analyzed, along with the one involving both flat surfaces facing each other, as in  . These involve the curved surfaces facing each other, or a curved surface facing a flat surface. As a further generalization, two sets of sizes are also considered. The first set involves both hemispheres of the same diameter, while the second involves diameters that differ from each other by a factor of four. The situation is shown in Figure 1, where 2A and 2B merge when B = A. The symmetry dictates that the net force between the two hemispheres is in the vertical direction. Calculations of the force will be done for separations between the closest two points, S, such that 0 ≤ S ≤ 50 A, essentially the asymptotic limit.

Figure 2 illustrates how the force associated with the second two arrangements can be derived from the force associated with the flat-to-flat arrangement through linear superposition of forces. Thus only the first arrangement requires a nontrivial calculation. The arrangements involving a sphere and a hemisphere are equivalent to a point source outside the hemisphere and positioned along its center line. This calculation is comparatively easy to perform.

In addition to the force calculations, it will be shown that the force can be expressed in the form of Newton’s law of gravitation for point masses, GMAMB/ (Reff)2, where Reff is an effective separation between the two point masses, each having the same mass as the associated hemisphere. Only the magnitude of the force is considered here. This situation is shown in Figure 3 for the flat-to-flat configuration. The location of each of the two fictitious point masses is within its respective hemisphere at a position αAA or αBB from the flat surface. The two α’s

refer to the fraction of the radius, A or B, where the point mass is located. If a curved surface faces down, the position will be measured from that curved surface. In all cases, the position will be calculated as a function of S and shown to approach the center of mass of the respective hemisphere at infinite separation. The separation of 50 A is fully adequate to demonstrate that point. This center-of-mass limit ((3/8) A (or B), measured from the flat surface, and (5/8) A (or B), measured from the curved surface) can be proved analytically and can be viewed as a test of the accuracy of the calculation. Their position at any separation will be calculated, as shown in the bottom of Figure 3, by replacing one of the hemispheres by its equivalent point mass at the unknown position αBB and calculating the force between it and the other hemisphere. That force between the two will be equated to the force previously calculated between the two hemispheres, thus determining αBB. Since α A A + α B B + S = R e f f , the other location is determined.

2. Calculations2.1. Preliminary

First, let us be specific about the superposition of collinear forces, at least in one case. Others are similar. Assuming the form GMAMB/(Reff)2 and referring to Figure 2, panel 2A, we find that

1 [ R e f f ( c f ; S 2 ) ] 2 = 2 [ R e f f ( p f ; S 2 ) ] 2 − 1 [ R e f f ( f f ; S 2 + B ) ] 2 (1)

where, in order of their appearance, “cf” refers to curved-flat; “pf” refers to point-flat (because the sphere acts as a point source), and “ff” refers to flat-flat, as also shown in the first panel. The separation for each Reff, S2 and S2 + B, are also shown. The factor of 2 exists because the mass of the sphere is twice that of the hemisphere.

Next, we present the formulas for the magnitude of the force between a sphere and a hemisphere required in Figure 2(A), Figure 2(B), and 3 of Figure 2. In panels 2A and 2B, the flat surface of the hemisphere faces the sphere, but with A and B reversed. In panel 3, the curved side of the hemisphere faces the sphere. A simple integration yields the following for the force, Fpf, in Figure 2(A):

F p f = G M A ( 2 M B ) [ 3 A 3 ] [ ( 1 3 S 2 ) { ( 2 S 2 − A 2 ) A 2 + S 2 + 3 A S 2 + ( A ( A − S ) − 2 S 2 ) ( A + S ) } ] (2)

where we identify the multiplier of [GMA(2MB)] with (1/Reff)2, S = S2 + B, and MA and MB are the two hemispherical masses. In the case of panel 2B, A and B are reversed in the above formula. The equation for the force between the curved side facing the sphere and the sphere can be obtained by a similar integration, resulting in a similar formula. More simply, it can be obtained by subtracting from the force between two spheres the force between the sphere and the hemisphere (flat side facing the sphere), as shown in Figure 4. The force between the two spheres is simply proportional to 4 / ( S 3 + A + B ) 2 .

2.2. Force between Two Hemispheres: Flat Surface Facing Flat Surface

As was stated earlier, this situation is the only one requiring a nontrivial calculation. It was performed in two ways. The first involves a direct integration of the force between all of the points in one hemisphere and all of the points in the

other. Because integration over two volumes was required, it started out as a six-fold integration. However, azimuthal symmetry reduced it by one integral and another integral could be done analytically, thus reducing the problem to a still-unpleasant four-fold integration. The second way involved expansions of the potential created by one hemisphere along its axis for positions z ≤ A and z ≥ A. The off-axis values were then expressed in terms of the original expansion multiplied by Legendre polynomials, as in electrostatics.

3. Direct Integration Method

The direct integration method was performed with respect to Figure 5, where two hemispheres, along with the two sets of dimensions and the separation between them are shown.

The basic integral expression for the point-to-point (pp) force, Fpp, is:

F p p = G ρ A ρ B ∬ d V A d V B cos φ R 2 (3)

where VA and VB are the volume elements in their respective hemisphere; R is the separation between any two points in different hemispheres; φ is the angle between R and the vertical; and ρA(B) is the mass density of hemisphere A(B), assumed to be different for the sake of generality. By symmetry, the force is in the vertical direction. This equation can be written in an explicit form:

F p p = G ρ A ρ B ∫ 0 B d z B ∫ 0 ( B 2 − z B 2 ) 0.5 R B d R B ∫ 0 2π d φ B ∫ 0 2π d φ A ∫ 0 A d z A ∫ 0 ( A 2 − z A 2 ) 0.5 R A d R A ( z A + z B + s ) [ ( z A + z B + s ) 2 + r A 2 + r A 2 − 2 R A R B cos φ A ] 3 2 (4)

Writing the equation in the form of the product of the two point masses divided by the square of an effective separation between them scaled by the radius A(1/Reff,A)2, we obtain the following complicated expression:

( 1 R e f f , A ) 2 = 9 π ∫ 0 1 d z b ∫ 0 ( 1 − z b 2 ) 0.5 R b d R b ∫ 0 π d φ a ∫ 0 1 d z a ( z a + b z b + s a ) [ ] (5)

where

[ ] = [ ( z a + b z b + s a ) 2 + ( b R b ) 2 ] 1 / 2 ( b R r b sin φ a ) 2 + ( z a + b z b + s a ) 2     + ( b R b cos φ a ) ( 1 − z a 2 ) 1 / 2 − ( ( z a + b z b + s a ) 2 + ( b R b ) 2 ) ( ( b R b sin φ a ) 2 + ( z a + b z b + s a ) 2 ) X 1 2

And where

X = ( 1 − z 2 ) − ( 2 b R b cos φ a ) ( 1 − z a 2 ) 1 / 2 + ( ( z a + b z b + s a ) 2 + ( b R b ) 2 )

Reff,A = Reff/A. All the other linear parameters in the expression are similarly scaled and use similar notation (e.g., sa ≡ S/A; b = B/A; za = ZA/A; ra = RA/A; zb = ZB/B; rb = RB/B). The integral over the azimuthal angle, φB, forthe upper hemisphere was simply replaced by 2π, and the integral over rA was doable analytically. Thus, as stated above, the 6-fold integral was reduced to a 4-fold integral. Since the 4-fold integral did not appear to be analytically evaluable, it was numerically evaluated by using the Wolfram Mathematica command “NIntegrate” as applied to multiple integration. The results of the direct integration are shown in Figure 6 for b = 1 and b = 1/4. Configuration 1 refers to the first panel in Figure 1. The basic behavior is the same for both values of b: a rapid drop off for small values of SA, followed by a transition to very small values of (1/Reff,A)2 at large values of

SA. The values for b = 1/4 are consistently larger than those for b = 1. At S = 0, for example, the two are 1.31 and 0.75, respectively. The force for b = 1 is considerably larger, however, because the product of the masses is 64 times as large. Given the value of 0.75 for (1/Reff,A)2 with b = 1, an expression for the magnitude of the force between the two hemispheres in contact with one another is (1/3) [G(πρA2)2], where ρ is the mass density of both hemispheres. This is the same result obtained in  by different methods from the one used here (see Equations (5) and 19 of that reference).

3.1. Gravitational Force between Two Halves of a Sphere

Since the gravitational force between the two halves of a solid sphere, such as the earth, is of some interest to inquiring minds  , we shall consider two applications of the result immediately above. The first concerns how the gravitational force between two contacting spheres, each of the same mass as the individual hemispheres, compares to the force between the two hemispheres. Since the ratio of the radii, sphere to hemisphere, is 1/21/3, the ratio of the two forces, spheres to hemispheres, can easily be deduced as (22/3/3) ≈ 0.53. Thus, the force between the two hemispheres is almost twice that of the two spheres. The second application is concerned with the rotation rate of the earth required to separate its two halves under only the combined action of the fictitious centrifugal force and gravity. The calculation is done with reference to Figure 7, in which is shown a sphere of radius A split along a plane perpendicular to the plane of the paper and passing through the central axis of rotation. The split sphere is rotating with angular velocity ω. The coordinates (r, φ, z) define the position of a rotating point. The total “centrifugal force”, Fc, is obtained by integrating the component of the force perpendicular to the interface plane over the hemisphere.

F c = 2 ω 2 ρ ∫ 0 A d z ∫ 0 ( A 2 − z 2 ) 1 / 2 r 2 d r ∫ − π 2 π 2 cos φ d φ = π 4 ω 2 ρ A 4 (6)

Equating Fc to [G(πρA2)2]/3, solving for ω, and expressing ω as 2πNr, where Nr is the number of rotations per second, we obtain Nr = [Gρ/(3π)]1/2. This results in ≈ 20 × 10−5 rotations/sec. required to just separate the two halves of the sphere. The earth’s rotational speed ≈ 1.16 × 10−5 rotations/sec., which ≈ (1/17) of the required speed. Of course, the separation speed would be much greater if the cohesion within the earth, which serves a kind of geological corpus callosum, were taken into account. Furthermore, tidal friction is causing the earth’s rotational speed to slowly diminish. Thus, we have nothing to worry about.

3.2. An Array of Results

In this section, we present a graphical array of results, in addition to the one already presented in Figure 6. The graphical array will be of (1/Reff,A)2 vs. SA for the various configurations in Figure 1 and of Alpha (the spelled out version of the Greek letters in Figure 2) for those same configurations. Figure 8 is a more detailed view of the behavior of (1/Reff,A)2 for SA ≤ 3, b = 1, and for its three distinguishable configurations. The dotted line is added to show the linear behavior in the immediate vicinity of SA = 0. The ordering in size is consistent with the amount of mass in one hemisphere that is close to the other hemisphere. That order is flat-to-flat, flat-to-curved, curved-to-curved. Figure 9 is similar to Figure 8, but with b = 0.25. For these graphs, (1/Reff,A)2 is larger than in Figure 8 because one of the hemispheres is so much smaller. The ordering in size in Figure 9, in general, is also consistent with the amount of mass in one hemisphere that is close to the other. Curiously, however, there is a cross over between con-

figurations 2A and 2B at SA = 0.25, which is maintained down to SA = 0. Figure 10 and Figure 11 are the same as Figure 8 and Figure 9, except that the entire range of values are shown on a log-log scale. The dashed line in Figure 10 shows how the behavior at large values of SA varies as 1 / S A 2 . Figure 12 and Figure 13 illustrate the behavior of Alpha for various configurations and the two values of b (1 and 1/4). Alpha 1 (or Alpha 1/4) is the value of Alpha for a hemispherical radius = 1 (or 1/4). As expected from the comment in the introduction, this parameter approaches (3/8) or (5/8) for large values of SA, depending on whether

the flat or curved side faces the other hemisphere. That condition is denoted by “flat” or “curved” in Figure 12. The curves in Figure 12 and Figure 13 can be interpreted as indicating that the force between the two hemispheres decays away more slowly than 1/S2 and only approaches that behavior at very large vales of S (as shown in Figure 10 and Figure 11). As Alpha diminishes, the increasing separation between the two fictitious point masses with S is partially offset by their move toward each other within their respective hemisphere.

3.3. The Use of Power Expansions and Legendre Polynomials

Without the use of high-speed computers and sophisticated computational software, such as Wolfram Mathematica or, perhaps, MATLAB, it would be extremely difficult to perform the four-fold integration that makes up Equation (5). This is particularly true because of the coupling among the variables of integration. However, the use of a power expansion of the on-axis gravitational potential of hemisphere A (or B), followed by a term-by-term multiplication by Legendre polynomials to produce off-axis values of the potential simplifies the calculation. This is commonly done when dealing with electrostatics problems containing azimuthal symmetry, and it is discussed in many textbooks  . This method can be applied to the current problem because, as in electrostatics, it requires a solution of Laplace’s equation. The only possible drawback is that the result is an infinite series, which may require a largenumber of terms to produce the desired accuracy. Fortunately, this will not be true of the number of terms required for an accurate solution to the current problem.

This method of solution for configuration 1 is presented in conjunction with Figure 14, where a hemisphere of radius A, assumed to be the source of the potential, Φ(R, θ), is shown with the hemisphere of radius B, assumed to be in the gravitational field of hemisphere A.

Once the potential, Φ(r, θ), is determined, its negative z-derivative (the z- component of its gravitational field) is calculated and integrated over hemisphere B to determine the total force between the two hemispheres. In the case of Figure 14, when R ≤ A, the on-axis potential can be written as a power series in Rl, whereas for R ≥ A, it can be written as a power series in R−(l+1). Thus, the potential at angle θ is of the form:

Φ ( R , θ ) = ∑ l = 0 ∞     C l R l P l ( cos θ ) ;     R ≤ A (7a)

Φ ( R , θ ) = ∑ l = 0 ∞     D l R − ( l + 1 ) P l ( cos θ ) ;     R ≥ A (7b)

Cl and Dl are parameters to be determined and P l ( cos θ ) is the Legendre Polynomial of order l, all of which are unity at θ = 0. In order to proceed, we must first determine the on-axis potential. After a simple integration it can be shown to be:

Φ ( R , 0 ) = − G ( 2 π ρ a ) [ ( A + R ) 3 3 R − ( A 2 + R 2 ) 3 2 3 R − A R − A 2 2 ] (8)

Ignoring constants in the above expression, which contribute nothing to the force, results in the following expansions for R ≤ A and R ≥ A, respectively.

Φ ( R , 0 ) = − G ( 2 π 3 ρ a A 3 ) ( 1 A ) [ − ( 3 2 ) ( R A ) + ( R A ) 2 − ( 3 8 ) ( R A ) 3 + ( 1 16 ) ( R A ) 5 − ( 3 128 ) ( R A ) 7 + ⋯ ] (9a)

and

Φ ( R , 0 ) = − G ( 2 π 3 ρ a A 3 ) ( 1 A ) [ ( A R ) + ( 3 8 ) ( A R ) 2 − ( 1 16 ) ( A R ) 4 + ( 3 128 ) ( A R ) 6 − ( 3 256 ) ( A R ) 8 + ⋯ ] (9b)

From these two equations, the constants Cl and Dl can be directly determined. To obtain Φ(R, θ), we multiply each term in the above two expansions by the appropriate Legendre polynomial, express the expansion in terms of the hemispheric coordinates z, r, φ and, to obtain the force, integrate its negative z-deriv- ative over z from s to ( s + ( B 2 – r 2 ) 1 / 2 ) and over r from 0 to B. The azimuthal integration over φ is simply replaced by 2π. In addition, we note that cos θ = z / ( z 2 + r 2 ) 1 / 2 . When all of this and a few other manipulations are completed, we identify (1/Reff,A)2 for R ≤ A as:

1 ( R e f f , A ) 2 = ( 3 b 3 ) ∫ 0 b r a d r a ∫ S a S a + ( b 2 − r a 2 ) 1 / 2 d z a ∂ ∂ z a [ 3 2 ( z a 2 + r a 2 ) 1 2 − ( z a 2 + r a 2 ) P 2 ( z a ( z a 2 + r a 2 ) 1 2 ) + ⋯ ] ;       R ≤ A (10a)

We note the use of scaled coordinates and the fact that the za-integral is particularly simple because it is the integral of a derivative. The integral over ra was done numerically. For R ≥ A, the results are similar:

1 ( R e f f , A ) 2 = ( 3 b 3 ) ∫ 0 b r a d r a ∫ S a S a + ( b 2 − r a 2 ) 1 / 2 d z a ∂ ∂ z a [ − ( z a 2 + r a 2 ) − 1 2 − ( 3 8 ) ( z a 2 + r a 2 ) − 1 P 1 ( z a ( z a 2 + r a 2 ) 1 2 ) + ⋯ ] ;     R ≥ A (10b)

For simplicity, only the first two terms of each expansion are shown in Equations (10a) and (10b). Figure 15 and Figure 16 show the results of this method for b = 1 and b = 0.25, respectively, along with those using Wolfram Mathematica. With only a very few terms of the expansion, the two sets are virtually indistinguishable on the graphs.

4. Summary

This paper consists of a study of the gravitational attraction between two inline hemispheres over a range of separations from zero to 50 times the radius of the larger of the two hemispheres. Two relative radii were chosen, 1:1 and 0.25:1 as sufficiently representative, along with three relative orientations: flat side facing flat side, curved side facing curved side, and flat side facing curved side. Initially, the first was calculated through a “brute force”, four-fold numerical integration with the help of Wolfram Mathematica. The second and third orientations were derived from the first by a superposition of the first with a simpler, easily calculated, intermediate configuration. The results were expressed in terms of the canonical form of two point masses separated by an effective distance. That “effective distance” was calculated, as was the location of the “point masses” within the hemispheres. An application of the 1:1, flat-to-flat configuration was made to a split, rotating earth. Finally, the flat-to-flat configuration was analyzed using a power series expansion for the on-axis potential and Legendre polynomials. With the use of very few terms in the expansion, the results agreed extremely well with the four-fold integration.

Based on the methods of calculation presented here, future work could involve other pairs of masses that form a cylindrically symmetric configuration. Elongating or compressing one or both hemispheres along their individual axes immediately comes to mind.

Acknowledgements

The author is grateful to Joseph S. Accetta, CEO of JSA Photonics, and William F. Filter, retired Sandia National Laboratories physicist, for reviewing the manuscript and making valuable suggestions for its improvement.

Cite this paper

Weiss, J.D. (2017) The Gravitational Attraction between Hemispherical Masses. Applied Mathematics, 8, 820-834. https://doi.org/10.4236/am.2017.86064

ReferencesGray, A. (1906) Solutions of Physical Problems. Proceedings of the Royal Philosophical Society of Glasgow, 37, 192.If the earth is cut into two halves, what will be the gravitational force between them? https://www.quora.com/If-the-earth-is-cut-into-two-halves-what-will-be-the-gravitational-force-between-themJackson, J.D. (1999) Classical Electrodynamics. 3rd Edition, John Wiley & Sons, New York.