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):

Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle v_{1,i+1} = v_{1,i} - g \Delta t }

Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle x_{1,i+1} = x_{1,i} + \frac{v_{1,i+1} + v_{1,i}}{2} \Delta t }

Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle v_{2,i+1} = v_{2,i} - g \Delta t }

Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle x_{2,i+1} = x_{2,i} + \frac{v_{2,i+1} + v_{2,i}}{2} \Delta t }

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:

  • Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle m_2 = 0.5 m_1 }
  • Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle m_2 = m_1 }
  • Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle m_2 = 9 m_1 }

construct Poincaré sections by plotting Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle v_2 } against Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle x_2 } at times when the two balls collide. Take Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle m_1 = 1 } , and use the initial conditions Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle x_1 = 1 } ; Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle v_1 = 0 } ; Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle x_2 = 3 } , and Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle v_2 = 0 } . 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 Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle x_2 } against time Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle t } for these three cases.

Part III for both PHYS460 and PHYS660 students

Employ the autocorrelation function

Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle C(\tau) = \int_0^\infty [x(t) - \bar{x}][x(t+\tau)-\bar{x}] dt }

which measures self-similarity of a signal with itself in different time intervals (here Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \bar{x} } is the average value of the signal), to analyze the "signal" Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle x_2(t)} 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"

Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle R(\tau) = \frac{1}{N-\tau} \sum_{n=1}^{N-\tau} x_n x_{n+\tau} }

which is essentially a correlation function, except that each point in Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle R(\tau) } 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 Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle x_2 = 3 } in steps of Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle 0.1 } , 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

Retrieved from "https://wiki.physics.udel.edu/wiki_phys660/index.php?title=Project_2&oldid=558"