If this is your first visit, be sure to check out the FAQ by clicking the link above. You may have to register before you can post: click the register link above to proceed. To start viewing messages, select the forum that you want to visit from the selection below.

Go Back   Bad Astronomy and Universe Today Forum > General > Questions and Answers
Register FAQ Members List Calendar Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
  #31 (permalink)  
Old 26-April-2008, 06:30 PM
grav grav is offline
Senior Member
 
Join Date: May 2006
Posts: 2,190
Default

Well, this is astonishing. In response to Cougar's post, I changed the dimensions of the disk in the program from 100,000 light years across and 12,000 light years thick for the Milky Way to 100,000 light years across and 1,000 light years thick, which accounts only for the stars themselves and not the dust and gas which extends out further, to see what that would give for a lesser thickness. Well, the force still comes out to almost twice as great as that of a sphere and the energy 1.2 times as great. I even found for twice as many points along each axis, for eight times as many points in all, and it comes out the same. If that continues with smaller thicknesses, then the gravity of a disk plane is finite. I'm not sure what that means yet, unless the integration must somehow be performed in a much different than the ordinary way or something. It might also have something to do with the natural "clumping" that is done when considering the mass over a finite number of points in the program instead of an even spreading out of the mass throughout the disk, as would be the case when performing integrations for it. I will also recheck the program. This sort of clumping does definitely seem to be the way to go with this.
__________________
Let's put together the pieces of The Grand Puzzle . (website)

"Let's define another operator, Sz, which we won't pay any attention to."
"This transformation will automatically make zero equal zero."
"It may be true that zero equals zero -- and that is certainly an equality -- but I don't want to go into the details at this time."
Reply With Quote
  #32 (permalink)  
Old 26-April-2008, 06:54 PM
grav grav is offline
Senior Member
 
Join Date: May 2006
Posts: 2,190
Default

Quote:
Originally Posted by grav View Post
Well, this is astonishing. In response to Cougar's post, I changed the dimensions of the disk in the program from 100,000 light years across and 12,000 light years thick for the Milky Way to 100,000 light years across and 1,000 light years thick, which accounts only for the stars themselves and not the dust and gas which extends out further, to see what that would give for a lesser thickness. Well, the force still comes out to almost twice as great as that of a sphere and the energy 1.2 times as great. I even found for twice as many points along each axis, for eight times as many points in all, and it comes out the same. If that continues with smaller thicknesses, then the gravity of a disk plane is finite. I'm not sure what that means yet, unless the integration must somehow be performed in a much different than the ordinary way or something. It might also have something to do with the natural "clumping" that is done when considering the mass over a finite number of points in the program instead of an even spreading out of the mass throughout the disk, as would be the case when performing integrations for it. I will also recheck the program. This sort of clumping does definitely seem to be the way to go with this.
Whoops. I rechecked the program and it appears I changed the numbers of points found for along the z axis, but not the actual thickness of the galaxy itself, so it was still finding for the same thickness, but just using a lesser number of points along it. Sorry. So for a thickness that is 3/25 that of the diameter of the galaxy, the force at the rim is about 1.9 times greater than that of a sphere of the same mass and the energy is 1.2 times greater. For a thickness that is 1/12 as great, or 1/100 of the diameter, the force at the rim is 2.75 times as great as that of a sphere of the same mass and the energy is 1.25 times as great. I'm still going to try working with a "clumping" variable in the integrations, though.
__________________
Let's put together the pieces of The Grand Puzzle . (website)

"Let's define another operator, Sz, which we won't pay any attention to."
"This transformation will automatically make zero equal zero."
"It may be true that zero equals zero -- and that is certainly an equality -- but I don't want to go into the details at this time."
Reply With Quote
  #33 (permalink)  
Old 26-April-2008, 08:46 PM
grav grav is offline
Senior Member
 
Join Date: May 2006
Posts: 2,190
Default

Quote:
Originally Posted by grav
Quote:
Originally Posted by Cougar
Then there's the dark matter effect, and its distribution....
Well, dark matter was actually introduced to account for this effect, that of the rotation speed curves, but it may just simply lie within the geometry and mass distribution itself. But some relatively accurate level of the geometry and mass distribution must be known, or at least theorized according to the results of the integrations themselves. As well as this, the correct integrations must also be performed. This does not appear to be the case with some of the earlier models for galaxies, incorporating spheres for simplicity and rings used in the same way as spheres, which led to the hypothesis of dark matter to begin with. But that is what I am working on now.
By the way, to be fair, this is not saying that a galaxy is regarded in the same way as a sphere. If that were the case, the force would climb in direct proportion to the distance all the way to the rim, and the rotational speed would be constant, so that the entire galaxy should rotate together with the same period of rotation for all points. The rotational speed at the rim, figuring for the total mass and radius of the Milky Way galaxy would then be about twice as great as is observed, which would actually require some form of "negative" dark matter, or lesser matter to account for the lesser speed observed at the rim. Also, if regarded as just a thin disk of uniform density, then the rotational speed would be expected to climb very rapidly at the rim, so a leveling off of the rotational speed would also require some form of "negative" dark matter. As it is, though, the rotational speed is predicted to climb in the same way it is observed to do up to some peak and then drop back down, but instead it just levels off at the peak, as shown in this link. The difference between the predicted curve and that of a geometry incorporating a uniform density would be that of the predicted mass distribution for the galaxy. I guess my next step, then, is to find the mass distribution that would give a curve that matches that for the predicted curve for a galaxy and go from there.
__________________
Let's put together the pieces of The Grand Puzzle . (website)

"Let's define another operator, Sz, which we won't pay any attention to."
"This transformation will automatically make zero equal zero."
"It may be true that zero equals zero -- and that is certainly an equality -- but I don't want to go into the details at this time."
Reply With Quote
  #34 (permalink)  
Old 26-April-2008, 09:16 PM
grav grav is offline
Senior Member
 
Join Date: May 2006
Posts: 2,190
Default

Cougar,

You might want to check my post #30 again. I accidently erased it , so I redid it, but I know I responded differently to your second paragraph after rereading it than I did the first time, not that I remember my response to that anyway, but I think I got the other two responses the same.
__________________
Let's put together the pieces of The Grand Puzzle . (website)

"Let's define another operator, Sz, which we won't pay any attention to."
"This transformation will automatically make zero equal zero."
"It may be true that zero equals zero -- and that is certainly an equality -- but I don't want to go into the details at this time."
Reply With Quote
  #35 (permalink)  
Old 26-April-2008, 09:36 PM
Jeff Root Jeff Root is online now
Senior Member
 
Join Date: Dec 2004
Location: Minneapolis
Posts: 3,291
Default

Quote:
Originally Posted by grav View Post
I guess my next step, then, is to find the mass distribution that would
give a curve that matches that for the predicted curve for a galaxy and
go from there.
It sounds to me like you still have to get the math right before the
simulations can have any value. The fundamental problem causing
infinities must be a common one -- I'm sure it has a common solution.

-- Jeff, in Minneapolis
__________________
http://www.FreeMars.org/jeff/

"The other planets?
Well, they just happen to be there, but the point of rockets is to explore them!"
-- Kai Yeves
Reply With Quote
  #36 (permalink)  
Old 26-April-2008, 09:55 PM
grav grav is offline
Senior Member
 
Join Date: May 2006
Posts: 2,190
Default

Quote:
Originally Posted by Jeff Root View Post
It sounds to me like you still have to get the math right before the
simulations can have any value. The fundamental problem causing
infinities must be a common one -- I'm sure it has a common solution.

-- Jeff, in Minneapolis
The math is right, although I am still trying to find out the predicted or observed mass distribution of galaxies, yes. The infinity is for the edge of a flat disk plane of uniform density, not one with some thickness, but as far as I know, the exact solution for a thick disk is not known either, only approximated, so that is why I am just using the program unless I can manage to work it out through the integrations. But even then, the program does the same thing as the integration, just point by point, and I can vary the mass distribution at will instead of reperforming the integrations accordingly, and it's much quicker to boot.
__________________
Let's put together the pieces of The Grand Puzzle . (website)

"Let's define another operator, Sz, which we won't pay any attention to."
"This transformation will automatically make zero equal zero."
"It may be true that zero equals zero -- and that is certainly an equality -- but I don't want to go into the details at this time."
Reply With Quote
  #37 (permalink)  
Old 27-April-2008, 12:53 AM
TomT TomT is online now
Senior Member
 
Join Date: Jan 2005
Posts: 616
Default

Quote:
Originally Posted by grav View Post
Here is a link I found a little while ago that summarizes some methods. It looks like I have already tried most of the ones mentioned. Is this pretty much all that has been explored of the subject so far?
Hi grav,
I think you are on the right track in referring to the above paper. Nicholson's calculation was discussed extensively in this forum here:

Galaxy Rotation Curve Calculation.


Nicholson's solution avoided the dreaded elliptical integral solution by using a numerical interation method in a computer program he used while working on the Space Shuttle Program. Dont stop reading the thread after the first parts, as it takes a turn later, and Nicholson joins the discussion.

TomT
Reply With Quote
  #38 (permalink)  
Old 27-April-2008, 02:39 AM
grav grav is offline
Senior Member
 
Join Date: May 2006
Posts: 2,190
Default

Thanks, Tom. Good to hear from you again. I remember you helped out with this when I first went through it. Hh and publius have since given me a two week crash course for finding integrations and I have learned a lot from them. I have been trying to use it for this, but you say it's not necessary, at least not with a computer program for it, or by using iterations, eh? I have read through some of Nicholson's papers before. Maybe this time I can follow it better. I'll read through it again and try to use any calculations he provides to see what it comes up with.
__________________
Let's put together the pieces of The Grand Puzzle . (website)

"Let's define another operator, Sz, which we won't pay any attention to."
"This transformation will automatically make zero equal zero."
"It may be true that zero equals zero -- and that is certainly an equality -- but I don't want to go into the details at this time."
Reply With Quote
  #39 (permalink)  
Old 27-April-2008, 06:40 AM
Jeff Root Jeff Root is online now
Senior Member
 
Join Date: Dec 2004
Location: Minneapolis
Posts: 3,291
Default

Quote:
Originally Posted by grav View Post
The math is right,
I don't think it can be right if it produces infinities, even if only for
the unrealistic case of zero thickness. I suspect that the added
calculations for points out of the plane are partially masking the
error which is apparent at zero thickness.

Quote:
Originally Posted by grav View Post
The infinity is for the edge of a flat disk plane of uniform density, not
one with some thickness,
If it happens at the edge, it should happen everywhere. I suspect
that in the interior, the error in one direction is masked by the error
in the opposite direction.

Can you describe how the infinities come about in your program?
Does your program actually calculate the attraction between a mass
at one point and the same mass at the same point? Or does it
calculate attraction between masses at two points that are very
close together? Or what?

-- Jeff, in Minneapolis
__________________
http://www.FreeMars.org/jeff/

"The other planets?
Well, they just happen to be there, but the point of rockets is to explore them!"
-- Kai Yeves
Reply With Quote
  #40 (permalink)  
Old 27-April-2008, 08:00 AM
mugaliens's Avatar
mugaliens mugaliens is offline
Senior Member
 
Join Date: Dec 2005
Location: Dortmund
Posts: 5,961
Default

Quote:
Originally Posted by grav View Post
They would still both lead to an infinite acceleration of gravity and therefore an infinite orbital speed for the point at the center of the mass that is right at the disk edge for a straight forward calculation of the disk plane. The calculation just blows up there.
I strongly disagree, for one simple reason - you have two stars, one on the edge of the circular field, the other positioned such that half the mass is inside it's orbit, the other half outside it's orbit. That mass is finite, and so is the gravitational attraction towards the center. Nothing goes to zero. Nothing approaches infinity.

We're dealing with finite entities, here, and the model is patterned after, and close to, our own Milky way. Clumping doesn't enter into the picture, nor do the spiral arms, because it's at such a distance, and there are so many clumps and arms, that while these might cause pertubations in our orbit, those are minimal unless we get fairly close to a clump. If that happens, then, yes - we might get tossed outwards (slingshot) or decayed inwards (orbital braking).
__________________
I am Mugs, of the Alien clan of Usa, Nordamerica, a Terran, of Sol.

Mine: "Perception isn't reality. It's merely an abstraction thereof, and quite often not a very good one at that."

Heinlein's: "Staying young requires the unceasing cultivation of the ability to unlearn old falsehoods." "Freedom begins when you tell Ms. Grundy to go fly a kite."
Reply With Quote
  #41 (permalink)  
Old 27-April-2008, 08:09 AM
mugaliens's Avatar
mugaliens mugaliens is offline
Senior Member
 
Join Date: Dec 2005
Location: Dortmund
Posts: 5,961
Default

Quote:
Originally Posted by grav View Post
The thing about the gravity at the edge of a disk plane, though, is that since it comes out to 4 G (M/A) [atanh(1) - 1], the gravity is just right on the "edge" of infinity.
You're really missing something, here, grav...

Instead of stars, let's use steel ball bearings, 1" in diameter, arranged in a circular plane. Let's use the same numbers as in the first problem I stated, and set things in motion such that all initial orbits are circular given the initial conditions. Let's assume we're in the center of that supermassive void they discovered so the local gravitational gradient from even the closest object is zero.

Forget the math you're using for just one second and answer this question: Is gravity "just right on the 'edge' of infinity?"

I'll answer it for you: Hardly.

In fact, it'll be so low that the outermost ball bearing will have an orbital velocity measured in micrometers per second, if not hour.

If we blow this back up to the model I proposed, the proportions increase along with it, but nothing, especially the gravitational tug at the edge of the disk, get's anywhere close to infinity, any more than it does for our own Milky way.

Again, the model is patterned after our Milky Way, and for all intents and purposes, when you get beyond about 1% of the Milky Way's radius, things can be approximated by using a thin disk model.

Thus, I conclude that since the only substantial difference between my model and the Milky Way is local clumping, I think you're tripping over your calculations, somewhere...[/quote]

Thanks, Jeff Root, and Cougar, for posts 27 and 28.
__________________
I am Mugs, of the Alien clan of Usa, Nordamerica, a Terran, of Sol.

Mine: "Perception isn't reality. It's merely an abstraction thereof, and quite often not a very good one at that."

Heinlein's: "Staying young requires the unceasing cultivation of the ability to unlearn old falsehoods." "Freedom begins when you tell Ms. Grundy to go fly a kite."
Reply With Quote
  #42 (permalink)  
Old 27-April-2008, 07:27 PM
grav grav is offline
Senior Member
 
Join Date: May 2006
Posts: 2,190
Default

To answer these last three posts at once, hopefully, the integration for a disk plane does indeed come out to a result of infinite acceleration at the very edge for a point particle. A larger sphere at the edge, half in and half out should experience the same thing, at least for the point at its center, or considering the sphere a point mass at the center. In other words, if I were to perform separate integrations for every point in the sphere toward the disk, the point in the center would still blow up regardless of the rest of the points, and other points in the sphere shouldn't make much difference adding or subtracting with their finite contributions, then, I wouldn't think.

If we took a bunch of spheres and arranged them in a disk, then they would have thickness, and the acceleration toward the disk would be finite.

The computer program doesn't blow up at the edge, since it is calculating over clumped points within the disk instead of a mass spread evenly throughout the disk, as with an integration. However, as I make the thickness of the disk smaller and smaller, the gravity at the edge climbs higher and higher.

I know it doesn't seem to make much sense that the gravity at the edge of a finite mass disk should be infinite, but I have two theories about this. The first is with the thickness. With zero thickness, the particle does not lie against a surface as it would with a sphere or thick disk. Since the force varies with the square of the distance, we might need two dimensions to cancel out the infinities. That can be tested by finding for the potential energy, which only varies with the inverse of the distance, should become finite against the one dimension of the disk's edge. The other is with the clumping. I'm thinking stars, even gas and dust, are not spread evenly throughout a geometry, so the minimum possible distance between points must be taken into account.

I am wondering about the precise relationship between the resulting force and energy in all of this.
__________________
Let's put together the pieces of The Grand Puzzle . (website)

"Let's define another operator, Sz, which we won't pay any attention to."
"This transformation will automatically make zero equal zero."
"It may be true that zero equals zero -- and that is certainly an equality -- but I don't want to go into the details at this time."
Reply With Quote
  #43 (permalink)  
Old 27-April-2008, 08:38 PM
Cougar's Avatar
Cougar Cougar is offline
Senior Member
 
Join Date: Aug 2003
Location: The Wild West
Posts: 3,806
Default

I'm not sure what "integration" you're talking about. All you need is the distance from your test particle to each object in the galaxy and the approximate mass of each object. I think you could simplify it to a single central black hole 5% to 10% of the total mass, and the rest solar-sized masses. Let the computer add up the resultant accelerations....

This "infinite gravity" is obviously an error.
__________________
Everyone is entitled to his own opinion, but not his own facts.
Reply With Quote
  #44 (permalink)  
Old 27-April-2008, 09:31 PM
grav grav is offline
Senior Member
 
Join Date: May 2006
Posts: 2,190
Default

Quote:
Originally Posted by Cougar View Post
I'm not sure what "integration" you're talking about. All you need is the distance from your test particle to each object in the galaxy and the approximate mass of each object. I think you could simplify it to a single central black hole 5% to 10% of the total mass, and the rest solar-sized masses. Let the computer add up the resultant accelerations....

This "infinite gravity" is obviously an error.
For a central black hole, we would still perform the same integration for the disk, and then just add the gravity of a point mass for the central body, which wouldn't add much considering the mass of the rest of the galaxy. In any case, here's the integration for a particle on the rim of a flat disk plane from post #14.

Quote:
If we integrate for the acceleration of gravity of a disk plane on a particle, we have Int G (M/A) dx dy (R-x) / d^3, where M/A is the mass per area which is uniform, R is the distance of the particle from the center, d is the distance from the particle to each point in the disk, and (x,y) are the coordinates of a point using the center of the disk as the origin. From this, we get

Int G (M/A) dx dy (R-x) / sqrt[(R-x)^2 +y^2]^3

= Int G (M/A) dx (R-x) * [2 * sqrt(r^2 - x^2) / (R-x)^2 / sqrt[(R-x)^2 + (r^2 - x^2)] ] for y = -sqrt(r^2 - x^2) to sqrt(r^2 - x^2)

but the last integration for dx does lead to complicated elliptic integrals, three different kinds, in fact, so let's avoid that by finding a solution for a specific R, such as R=r, for acceleration of gravity acting on a particle right at the edge of the disk. Then we get

Int 2 G (M/A) dx sqrt(r^2 - x^2) / (r-x) / sqrt[(r-x)^2 - r^2 -x^2]

= Int 2 G (M/A) dx sqrt(r+x) sqrt(r-x) / (r-x) / sqrt[2r * (r-x)]

= Int 2 G (M/A) dx sqrt(r+x) / (r-x) / sqrt(2r)

= 2 G (M/A) / sqrt(2r) * [ 2 sqrt(2r) atanh[sqrt(r+x) / sqrt(2r)] - 2 sqrt(r+x)]

from x = -r to r, so

= 4 G (M/A) [ atanh(1) - 1]

Here, since x = -r to r, it reduces to just zero for the whole thing for -r, but for +r, the part in atanh becomes sqrt(r+x)/sqrt(2r) = 1, and atanh(q) = (1/2)[ln(1+q) - ln(1-q)], so if q=1, then - ln(1-q) = -ln(0) = - (-infinity) = infinity. We have nothing else that might cancel out the infinity here, so that is the result for a disk plane. The gravity is infinite at the edge of the disk plane.

Now, I am thinking that finite masses cannot result in infinite force, so the only thing I see that might cause this is due to considering any and all point masses that lie directly in front of the particle, even within zero distance, like singularities. I am thinking that a better method would be to allow some distance between the particle and the closest possible point mass. For the gravity of something like a planet or the sun, this could possibly be as small as the distance between atoms, so wouldn't make much difference in the integration for those, at least that of spheres. With galaxies, however, especially being disks, I'm thinking an average or minimum possible distance between stars must be accounted for, at the very least being that of the diameter of a typical star.
__________________
Let's put together the pieces of The Grand Puzzle . (website)

"Let's define another operator, Sz, which we won't pay any attention to."
"This transformation will automatically make zero equal zero."
"It may be true that zero equals zero -- and that is certainly an equality -- but I don't want to go into the details at this time."
Reply With Quote
  #45 (permalink)  
Old 28-April-2008, 03:13 AM
Jeff Root Jeff Root is online now
Senior Member
 
Join Date: Dec 2004
Location: Minneapolis
Posts: 3,291
Default

grav,

I'm going to reply to a couple of things you just said in a separate post in a
few minutes, but right now I want to suggest a test I just thought of.

Can you have your program calculate the force on a plane shape other than
a disk? Like, the force on a point at the middle of an edge of a square? Or
on the edge of a U-shape? Like this:

Code:
  _____       _____
  |    |     |    |
  |    |     |    |
  |    |__.__|    |
  |               |
  |_______________|
The period indicates the location of the test point.

From your descriptions, I expect that your program will give a result of
infinite force toward the bottom of the U, and the force from the upper
part of the U, while calculated, will have no effect on the result. If it
does have an effect, that could tell you something about the problem.

You might want to use an L-shape instead of a U. The lopsidedness
could provide more diagnostic data.

-- Jeff, in Minneapolis
__________________
http://www.FreeMars.org/jeff/

"The other planets?
Well, they just happen to be there, but the point of rockets is to explore them!"
-- Kai Yeves
Reply With Quote
  #46 (permalink)  
Old 28-April-2008, 03:54 AM
Jeff Root Jeff Root is online now
Senior Member
 
Join Date: Dec 2004
Location: Minneapolis
Posts: 3,291
Default

Quote:
Originally Posted by grav View Post
The computer program doesn't blow up at the edge, since it is
calculating over clumped points within the disk instead of a mass
spread evenly throughout the disk, as with an integration. However,
as I make the thickness of the disk smaller and smaller, the gravity
at the edge climbs higher and higher.
That should make the problem much easier to resolve than if it only
showed up when it reaches zero thickness. This way, you will get
a definite indication that you have found the problem and fixed it,
rather than applying some arbitrary restriction which makes the
infinities go away but doesn't necessarily give correct results.

Quote:
Originally Posted by grav View Post
With zero thickness, the particle does not lie against a surface as it
would with a sphere or thick disk.
I don't understand the geometry you are describing when you say
"...lie against a surface...".

Quote:
Originally Posted by grav View Post
Since the force varies with the square of the distance, we might need
two dimensions to cancel out the infinities.
I think you shouldn't be producing infinities in the first place. Cancelling
them out sounds like a hack, and a particularly dangerous one.

-- Jeff, in Minneapolis
__________________
http://www.FreeMars.org/jeff/

"The other planets?
Well, they just happen to be there, but the point of rockets is to explore them!"
-- Kai Yeves
Reply With Quote
  #47 (permalink)  
Old 28-April-2008, 04:03 AM
Jeff Root Jeff Root is online now
Senior Member