Project 2

From phys660
Jump to navigationJump to search

Deterministic chaos in conservative dynamical system of two balls in one dimension with gravity

Consider a one-dimensional system that consists of two balls moving along the x-axis in response to gravity, as illustrated in the figure below:

A simple dynamical system with two degrees of freedom consists of two point masses or balls constrained to move in one dimension above a floor in a constant gravitational field. Due to assumed elastic collisions, the system exhibits discontinuities as a function of time.

Collisions between the balls and between the lower ball and the floor are completely elastic, i.e., this is a conservative Hamiltonian system.

Write a program to follow the motion of both balls, where you have to pay special attention to the collision events at which the motion changes abruptly. You can do this either

(a) Suggested for PHYS460 students: by using finite difference midpoint method (which actually yields exact, but discretized, solution for motions with constant acceleration):

or

(b) Suggested for PHYS660 students: exactly, i.e., without using finite difference equations if you can look ahead and calculate exact times and positions of the collisions.

When the collision with the floor occurs, you should reverse the velocity of the lower ball. On the other hand, when the balls collide with each other, you should use the formula for elastic collisions of two objects (derived in elementary classical mechanics) to find their velocities after the collision.

Part I for both PHYS460 and PHYS660 students

For the three cases:

construct Poincaré sections by plotting against at times when the two balls collide. Take , and use the initial conditions ; ; , and . Which, if any, of these plots indicate that the system is chaotic?

Part II for both PHYS460 and PHYS660 students

Find the position and velocities of the two balls at equal time intervals. Plot against time for these three cases.

Part III for both PHYS460 and PHYS660 students

Employ the autocorrelation function

which measures self-similarity of a signal with itself in different time intervals (here is the average value of the signal), to analyze the "signal" generated by the motion of the second ball. Note that for a discrete set of data composed of N points one should switch to "autocovariance"

which is essentially a correlation function, except that each point in is normalized by the number of points used to calculate the discrete integral.

Comment on similarities and differences in the three cases above. In general, we expect autocorellation function to be constant or oscillating for regular motions, while it decays fast to zero (typically in exponential fashion) for chaotic motions.

Part IV for PHYS660 students only

Repeart Poincaré section computation in c) for at least 11 additional sets of initial conditions. For example, you can decrease in steps of , while adjusting initial velocities so that the total energy remains the same for each set. Plot all 12 Poincaré sections on the same graph and compare with a single Poincaré section you obtained in I) using just one of these 12 sets of initial conditions.

REFERENCES