Oliver Sheridan-Methven | PDF - Charterhouse
Oliver Sheridan-Methven | PDF - Charterhouse
Oliver Sheridan-Methven | PDF - Charterhouse
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
<strong>Oliver</strong> <strong>Sheridan</strong>-<strong>Methven</strong> Chaos Supervisor - DL<br />
Pageites Word Count - 6460<br />
Chaos<br />
During the course of science, scientists have always looked for patterns that<br />
occur around them. The purpose of which, is to establish a pattern, and then to try to<br />
explain why this pattern might occur in nature. This ranges from finding patterns between<br />
rabbit and fox populations to the pattern between the force on an object and its<br />
acceleration, F=ma.<br />
Once a pattern has been established, physical laws and mathematical axioms<br />
are formed to create a model of the observed system. If the model produced is<br />
accurate and can produce results similar to experiments then the theory will survive; until<br />
it fails to serve as a good enough theory for predicting further results. Also, if two different<br />
theories can be explained by one unifying theory then the two will be subsumed into one<br />
new theory. Slowly phenomena in the world have been whittled down to the<br />
fundamentals explanations (axioms) we believe today, through this refining of theories by<br />
spotting further patterns in nature. However, whilst these microscopic effects can be<br />
predicted, there is still a huge and difficult science of Chaos which tries to solve the<br />
complex macroscopic systems which are operated by these microscopic laws [1,2].<br />
The basis of one physical branch of chaos is chaotic dynamics, which is the<br />
analysis and study of a dynamic (ever changing) system; which may have only a few<br />
simple equations to say how it will behave/move; such as F=ma or F=GMmr -2 . To try to<br />
solve these equations though, and to make predictions from them, may be impossible<br />
(shown later in this essay)[7]. Thus the fundamental basis of chaos, is trying to take the<br />
basic principles in nature and equations usable, and to try to make predictions from<br />
them; it would be pointless to know how everything in the world works so intricately and<br />
not be able to actually use such knowledge to solve real life problems [2]!<br />
There are many types of systems in the<br />
world (and universe) and it is difficult to say whether<br />
a system is chaotic or not. The boundary between a<br />
chaotic system and a "normal" system is tricky to<br />
define; but there are characteristics of which are<br />
often unique to chaotic systems. Normally when a<br />
system possesses more than one of these<br />
characteristics it is highly likely the system displays<br />
chaos [1].<br />
Many people would consider this cloud<br />
chamber pattern to appear chaotic.<br />
However, there is a key point to make note of with most physical and<br />
mathematical branches of chaos, which is the striving to avoid random input [1]. What<br />
this signifies is that nothing should be random and in an equation the only special values<br />
are physical constants and parameters. If the model requires some sort of random<br />
number generator then it should avoid being used; as eventually any random number<br />
generator will repeat itself or will cycle through a certain range of values [2]. Thus it will<br />
naturally input a form of predictability and a series of results which is a type of bias and<br />
predisposition (which physics and mathematics likes to avoid). This is one of the key<br />
scientific ethics in physics which is that we should study something and then see how it<br />
works and make a model to fit it [2], not to make a model which suits us and then "tamper<br />
the results" to agree with the model, which is bad science, joked about by Albert Einstein;<br />
"If the facts don't fit the theory, change the facts".<br />
Page 1 of 18
<strong>Oliver</strong> <strong>Sheridan</strong>-<strong>Methven</strong> Chaos Supervisor - DL<br />
Pageites Word Count - 6460<br />
What we might consider to be a normal system to be like then, would be one<br />
where we could reasonably easily work out, either logically or mathematically, what<br />
state the system has been in, and will be in, after some period of time. Thus given the<br />
conditions at the start we could accurately predict what it will look like at a later point in<br />
time. However, this "forecasting" is only a small branch of physics, pointed out by Richard<br />
Feynman; "Physics like to think that all you have to do is say, these are the conditions,<br />
now what happens next". However, a chaotic system can be extremely difficult or even<br />
impossible to predict well [1,2,3], such as the three body problem (explained later). Thus a<br />
property of a chaotic system is the issue with predictability; if it is difficult to predict or<br />
approximate what a later state will look like, then the system may be chaotic.<br />
This property of unpredictability can be caused by either/all of these three<br />
qualities of the system:<br />
A. Sensitive dependence<br />
B. Feedback loops<br />
C. Unsolvable equations [1,3,7,8]<br />
The last of these, C, is the simplest one; if the equations which are used to<br />
describe how the system behaves cannot be solved to give an exact solution (e.g. most<br />
nonlinear differential equations) then it may be impossible to calculate what will happen.<br />
Most often though, to resolve this problem, approximations are made to make the<br />
equations simpler; so the equations might then merge together to provide a solution.<br />
Although the result it is not 100% correct it may be a good estimate. However, if this<br />
"guessing" needs to be repeated many times to extend the prediction further into the<br />
future, slowly the answer will become increasingly more unreliable [1,9].<br />
The quality of sensitive dependence, A, is an issue not of the model's equations<br />
but rather of the data put into the model. Sensitive dependence arises with the<br />
uncertainty and error in a recorded value. If some data is input then a result is produced,<br />
however, as no data is 100% accurate [9], it is most likely a little bit off. Thus when a value<br />
is put in a range of values should be calculated, from the lowest possible value to the<br />
highest to give the spectrum of possible outcomes taking into account how wrong your<br />
measurements and data could be [9]. With many systems a little change in the data<br />
results in a little change in the result which is proportional to the error. This is the case of<br />
linear systems [3]. However, nonlinear systems are harder to work with because a little<br />
change in the data can result in a huge change in the outcome. Thus if there is a little<br />
error in the data there could be a huge variation in the results (technically quantified by<br />
Lyapunov exponents), thus if the data is not to a sufficient accuracy, the model could be<br />
unusable for making predictions [1]. James Gleick made the analogy to a pebble in a<br />
bowl and a pencil on a table. With the pebble in the bowl, if the pebble is in the bottom<br />
of the bowl then it will stay in the bottom, and if the data is wrong and the pebble is said<br />
to be an inch from the bottom then it will soon roll and settle into the bottom of the bowl<br />
and thus the system is not sensitive to it's initial conditions and it would be accurate for<br />
making predictions as to where the pebble will be. However,<br />
if you're trying to balance a pencil upright on its tip, then a<br />
tiny inaccuracy in the position of the top of the pencil can<br />
give many different results, and so this model will not be<br />
reliable for making predictions due to its sensitivity on data<br />
accuracy; which is an example of sensitive dependence [2].<br />
This is one of the reasons why it is easy to stay balanced<br />
sitting a bath but difficult to walk on a tight rope!<br />
Page 2 of 18
<strong>Oliver</strong> <strong>Sheridan</strong>-<strong>Methven</strong> Chaos Supervisor - DL<br />
Pageites Word Count - 6460<br />
The final feature of a feedback loop in a system, B, is often found hand in hand<br />
with sensitive dependence; commonly seen as two sides of the same coin. However,<br />
there are subtle differences between the two which are quite important. A feedback<br />
loop is often associated with an iterative solution to a problem, which might be forced by<br />
C. Thus if an error is present then it will still be present in the answer. The following state of<br />
the system is dictated by the previous state of the system when using iteration and so any<br />
uncertainty will be translated through the iterative feedback loop. Dependent on the<br />
system this can then be the cause of A but the fact that A is present doesn't imply that B<br />
is [1]. Thus B almost always forces A to occur, but A being present doesn't mean that B is<br />
(Figure1).<br />
Sensitive<br />
Dependence<br />
Figure1<br />
Normally Produces<br />
Sometimes Contains<br />
Feedback<br />
Loops<br />
The simplest system created which displays all of these characteristics was first<br />
discovered in biology and later analysed mathematically. This is the logistics map and<br />
was used initially to model a population of animals. What must happen is that the<br />
population must increase over time when it is small, so the first term was put in:<br />
F(X1) = λX0<br />
What this does is it increases the previous year's population, X0, by a fixed ratio,<br />
λ, to give the next years population, X1. However, this means that after enough time the<br />
population will become effectively infinite and bigger than it could possibly be. Thus a<br />
correctional term is introduced to bring down the population when it gets big and to<br />
boost it when it is small. NB the state of the population is between 0-1 for simplicity [1,2,3].<br />
F(X1) = λX0(1-X0) Ξ λ(-X0 2 + X0)<br />
Thus if this feedback loop is iterated (Figure2) then it should hopefully settle into<br />
a stable population. NB as this population represents the state of the system, this<br />
forecasting function can be used in other systems in maths, physics etc. However, when<br />
this is done over a range of values for λ then whilst it appears there may be stability in low<br />
values, past a certain point it will split into repeating solutions. This means that the system<br />
is oscillating between two states in a regular cycle which is reasonable. However, this<br />
splitting (birefringence) can occur again and again as λ increase until suddenly there is<br />
no order in the answers and there is a chaotic system! How this chaos mathematically<br />
arises is symbolically shown by the iterative process beside the logistics map in Figure2 [7].<br />
Page 3 of 18
<strong>Oliver</strong> <strong>Sheridan</strong>-<strong>Methven</strong> Chaos Supervisor - DL<br />
Pageites Word Count - 6460<br />
X 唴<br />
_<br />
F(X 1 ) = λ(-X 0<br />
2 + X 0 )<br />
Figure2<br />
λ<br />
What the logistics map symbolises for chaotic dynamics is the fact that in some of<br />
even the most simple systems; which can have just two or three factors effecting it,<br />
where one of these factors is adjustable, then it can be adjusted to produce any of three<br />
system states. First of these states is the state where nothing is happening and any initial<br />
state will soon settle into a fixed state. The second is where the system will fall into an<br />
active steady state where it will either be doing a fixed motion or will be cycling through<br />
states regularly (demonstrated in a later experiment). The final implication of the logistics<br />
map is that if the correct parameter (λ) is chosen, then the system will be unpredictable<br />
and will fall into a system of chaos that will appear to have no pattern and therefore<br />
looks random [2,7].<br />
When a system appears random the challenge is to then spot some sort of<br />
underlying pattern or trend; and to manipulate the data until such a pattern (and thus<br />
predictability) arises. Whilst many chaotic systems may not have an underlying pattern, a<br />
significant amount does. The task is to first find any strange attractors; which are points or<br />
trajectories through phase space. Phase space is just a graph of how the system is<br />
behaving as time goes by (Figure3) and where a certain point represents a<br />
characteristic, such as velocity or angle or potential energy etc. Thus one point could<br />
represent a system having a lot of kinetic energy when it is at such an angle whilst<br />
another point signifies another state the system can be in such as low kinetic energy but<br />
high temperature [1,2,3]. Thus if all paths through phase space always end up reaching a<br />
certain point then that point is a strange attractor. Similarly if a path continues orbiting a<br />
point in a circle for example, then the centre is a strange attractor. Finally, if all paths<br />
converge eventually into a certain path/trajectory (be it a straight line or a bizarrely<br />
shaped curve as in Figure3), even if initially off that path, then that trajectory is known as<br />
a strange attractor.<br />
Page 4 of 18
<strong>Oliver</strong> <strong>Sheridan</strong>-<strong>Methven</strong> Chaos Supervisor - DL<br />
Pageites Word Count - 6460<br />
Figure3<br />
Strange Attractors<br />
Figure3 shows how the measurements of a system and how it is behaving over<br />
time (top row); e.g. velocity with time, translate into phase space diagrams (bottom row)<br />
of height with gradient. Here it is evident that where the top row may be unclear in<br />
displaying how the system is behaving and what it will do with time (most prominent in<br />
the top right plot), when transformed into a phase space diagram it can become<br />
apparent what the system is doing. Thus apparent randomness, in what appears to be a<br />
chaotic system, can be translated into predictability in phase space, allowing for a<br />
chaotic system to be made predictable [1,3].<br />
The idea of there being underlying patterns in phase<br />
space which produce predictability, seems counterintuitive for<br />
chaotic systems. The fact that "randomness" can have patterns<br />
seems slightly illogical (and to an extent impossible) and thus<br />
an experiment was conducted to try to demonstrate this<br />
concept with real data collected from the lab, and not<br />
produced from theory. A very basic system encountered<br />
which was unexpectedly chaotic was that of the pendulum.<br />
For an equivalent James Gleick used the image of a child on a<br />
swing, which will be extended upon in the explanation of the<br />
experiment.<br />
If a child is on a swing and the swing is not moving<br />
then as time goes by the child will stay on the swing in the<br />
same position, which is obvious. If a child is then pushed once<br />
in the swing then they will oscillate to and from the same<br />
height, just at different angles on the swing (demonstrating the<br />
conservation of energy). Finally if friction is acknowledged then<br />
eventually the energy will dissipate from the system and the<br />
child will eventually come to a halt in the bottom of the swing<br />
position again. When these three processes are plotted in<br />
phase space (Figure4) it is much clearer to see the behaviour<br />
of the system and to spot the strange attractors in the system.<br />
NB even though the system is not yet chaotic, it has strange<br />
attractors in phase space [2].<br />
In the frictionless system the swing will carry on<br />
swinging and the phase space will form a loop (Fiugre4) which<br />
Strange Attractors<br />
Figure4<br />
Page 5 of 18
<strong>Oliver</strong> <strong>Sheridan</strong>-<strong>Methven</strong> Chaos Supervisor - DL<br />
Pageites Word Count - 6460<br />
demonstrates that after a certain period of time it will reach the same state. As it is then<br />
identical in state to how it was previously, it will follow the same course of action, and<br />
thus if a loop is formed it will continue in that loop forever. Equally, with the swing that has<br />
friction, it must eventually settle down to a fixed point which is when the swing has<br />
stopped [2] and so the phase space trajectory will always collapse into the origin<br />
(Fiugre13).<br />
However, a simple modification can be made to this system which quite<br />
unexpectedly introduces chaotic properties. This is to introduce a regular pushing force,<br />
such as a periodic torque on the swing. At first glance, this would seem just as regular as<br />
before, only just slightly more complicated mathematically. When a child is swinging,<br />
every time they reach the top of their return swing, they pull on the swing's rope (doing<br />
work on it and adding energy to the swing system) so that they stay in a regular swinging<br />
motion; where's the chaos in that? However, this is not a continuous and regular force;<br />
these are discrete and well timed packets of energy input specifically chosen for that<br />
frequency of swinging [2]. If the child were to do just the same rhythm of tugging on the<br />
rope on a longer swing, then they would pull at the wrong time, cause a jerk in the<br />
swing's motion, either speed up or slow down too quickly and would fall off. As will be<br />
shown later by the experiment, they may also swing through a full rotation and then<br />
suddenly stop at the bottom, which would certainly cause any child injury.<br />
Upon a technical analysis of what happens, is that although the pattern of<br />
forces applied may be in equilibrium, only some of the time, will this result in a system<br />
which is in equilibrium. This unexpected point was noted by James Gleick; "Even when a<br />
dampened, driven system is at equilibrium, it is not at equilibrium". This is due to the non<br />
linear relationship between the forces as time goes by [2,5]. The flow of energy in and<br />
energy out of the system has non linear equations, because the frictional torque is<br />
derived from the angular velocity, which changes with the angle of the swing as it goes<br />
from the bottom to the top of its arc. This is an example of how the system has a<br />
feedback loop and as we will see from the later experimental results, sensitive<br />
dependence.<br />
ΣT<br />
Φ<br />
L<br />
m<br />
Figure5<br />
m<br />
ω<br />
T 1 ,T 2 ,T 3<br />
The setup for the experiment was as follows<br />
(Figure5) [5]:<br />
For the experiment (in theory) there should be<br />
three sources of torque acting on the pendulum; T1, T2, T3<br />
[5]. The first is the torque from gravity which would<br />
produce simple harmonic motion (approximately for δΦ)<br />
provided it was the only force acting on the pendulum.<br />
The second is the torque due to friction which always<br />
works against the direction of the angular velocity. The<br />
final is the torque due to the driving force (child pulling<br />
on the swing's rope) adding torque in a periodic<br />
manner.<br />
Mathematically, these torques are given by:<br />
o T1 = -LmgSin(Φ)*<br />
o T2 = -kω †<br />
Commonly approximated to –LmgΦ to produce SHM<br />
Where k is a parameter value to determine the strength of<br />
the frictional torque<br />
Page 6 of 18
<strong>Oliver</strong> <strong>Sheridan</strong>-<strong>Methven</strong> Chaos Supervisor - DL<br />
Pageites Word Count - 6460<br />
o T3 = ASin(λt)<br />
This provides an external periodic torque on the system<br />
NB *the negative values denote for the first torque that it always acts to bring<br />
the pendulum back to the rest position, against the orientation of Φ, and is inwards<br />
acting. † The second negative torque signifies that the friction always acts in the<br />
orientation opposing the angular velocity of the pendulum.<br />
Thus by the laws of rotational motion where:<br />
Iα = ΣT as I = mL 2<br />
mL 2 α = T1 + T2 + T3<br />
mL 2 α = -LmgSin(Φ) –kω + ASin(λt) [5]<br />
α = -gL -1 Sin(Φ) + (-kω + ASin(λt))m -1 L -2 [7] as ω = Φ' and α = Φ''<br />
Φ'' = -gL -1 Sin(Φ) + (-kΦ' + ASin(λt))m -1 L -2<br />
The problem with these equations is that they cannot be solves explicitly to find<br />
the solution for Φ [7,8]. Thus there is no general solution for the angle with time, and<br />
therefore it is unpredictable with regards to a mathematically exact solution. This is due<br />
to the equation being a second and first order differential equation, which is nonlinear [8].<br />
This can only be approximated, to guess a solution for a value of time, by performing an<br />
iterative approximation over δt periods. This will act to guess the position quite well, but as<br />
time goes by the guess will eventually become evermore unreliable. This demonstrates<br />
how this system is mathematically chaotic. In reality though, does it actually behave<br />
chaotically? The aim of the experiment was to discover if any strange attractors appear<br />
after analysis of the data; and whether this system of unsolvable equations produces<br />
reliability. The question is whether from the chaotic equations a regular system can be<br />
found?<br />
The experiment was aimed to resemble the computer simulation provided by<br />
source [5] as closely as possible (Figure6). There were initial problems with methods of<br />
recording the position of the end of the pendulum with the corresponding time by use of<br />
video tracking due to the speed of the pendulum. This was resolved with a device with a<br />
protruding axle which was attached to the pivot of the pendulum. The device has a<br />
relatively much lower rotational friction compared to the driving force device (motor) of<br />
the pendulum. Thus whilst this setup was not ideal, as the added friction was relatively<br />
minor effect, it could be neglected. The driving torque was provided by an electrical<br />
motor powered by a signal generator giving out an AC produced by a varying voltage<br />
which followed a sine wave. This produced a varying torque which exactly matched T3.<br />
NB this was verified by the cathode ray oscilloscope during the experiment.<br />
Page 7 of 18
<strong>Oliver</strong> <strong>Sheridan</strong>-<strong>Methven</strong> Chaos Supervisor - DL<br />
Pageites Word Count - 6460<br />
Figure6<br />
The CRO (above) verified the voltage<br />
acting on the pendulum was a sine<br />
wave and that the signal generator<br />
(below) caused an AC through the<br />
motor.<br />
The motor here provided the driving<br />
torque to act on the pendulum. The<br />
nuts and bolts at the end<br />
approximate a pendulum as closely<br />
as possible.<br />
The blue device records the angle of<br />
the pendulum with time and has a<br />
very low frictional torque. The axes<br />
were joined together by the tube<br />
device at the top of the pendulum.<br />
NB – Whilst the essential principles for the pendulum were from [5], the experiment was designed entirely by myself.<br />
However, the experiment’s setup was constructed with the aid of the <strong>Charterhouse</strong>’s Physics technicians. Yet, all data was<br />
collected and analysed by myself.<br />
M<br />
CRO<br />
A<br />
C<br />
m<br />
Output data of angle v s<br />
time<br />
Φ vs t<br />
The pendulum was approximated as closely as it could by a light metal rigid<br />
rod. A mass of two large metal washers was bolted onto the end of the pendulum so that<br />
the bulk of the mass was at the end of the rod (approximating a pendulum). When the<br />
signal generator was switched on the data was then recorded into "data studio", and<br />
then moved to Microsoft Office Excel (Figure7).<br />
Figure7<br />
Page 8 of 18
<strong>Oliver</strong> <strong>Sheridan</strong>-<strong>Methven</strong> Chaos Supervisor - DL<br />
Pageites Word Count - 6460<br />
The parameters of the experiment were:<br />
o The amplitude of the voltage from the signal generator and thus the<br />
maximum torque which could be driven into the pendulum.<br />
o The length and mass of the pendulum.<br />
o The frequency of the signal generator and thus how frequently the driving<br />
torque changes direction.<br />
Since the aim was to determine the nature of how birefringence from stable to<br />
chaotic might occur within the system due to the variation of just one of the parameters,<br />
two of the parameters were forces into set values. The only parameter which could be<br />
accurately changed through a relatively large range of values was the signal frequency,<br />
thus the others were kept at constant values during the collection of data. The signal<br />
amplitude was maintained at 5.6V + 0.1V and the length of the pendulum was always<br />
0.14m, which did not change during experiments, however, the error in this length was<br />
+0.003m.<br />
One of the preliminary experiments was to determine the nature of friction in<br />
the setup. Whether the frictional torque was a first, second, or even zero order<br />
polynomial function of the angular velocity. To discover this, a large blob of putty was<br />
formed into a ball around the axis, before all the data was collected, and was spun up<br />
to a very fast angular velocity (c200 rad/s), by a separate power supply. Then the power<br />
was disconnected and the observed angular deceleration was recorded (Figure8).<br />
Figure8<br />
Page 9 of 18
<strong>Oliver</strong> <strong>Sheridan</strong>-<strong>Methven</strong> Chaos Supervisor - DL<br />
Pageites Word Count - 6460<br />
The significance of the straight line in Figure8 during the deceleration period<br />
shows that the frictional torque acting on the system acted on a zero order polynomial<br />
function of the angular velocity. What this meant was that there was an invariant<br />
(constant) frictional torque acting on the pendulum whenever it was in movement.<br />
Once this preliminary was conducted, the experimental data could be<br />
collected (Figure7). This was done over a set of values for the signal frequency from<br />
1.01Hz – 1.31Hz +0.01Hz. A typical spread sheet of data would contain a table of the<br />
initial conditions, the time, angle, angular velocity and a corrected angle. The corrected<br />
angle was a manipulation of the recorded angle; to transform it such that if the<br />
pendulum has rotated through 360 0 (2π c ) although the angles are recorded differently,<br />
then this new corrected angle would take the same value when the pendulum was in<br />
the same position. The reason for this should become clear after the analysis has been<br />
performed and the phase space trajectories plotted.<br />
Then a data set of the<br />
recorded angle against time was<br />
taken (Figure9) which could show<br />
how much the pendulum had<br />
rotated during the course of the<br />
data collection. What this shows is<br />
that during the blue run the<br />
pendulum had no preferred<br />
rotation, whilst in the purple run,<br />
the pendulum tended to<br />
repeatedly rotate through a full<br />
anticlockwise turn. This<br />
demonstrates how although all the<br />
conditions can appear the same to a<br />
high degree of accuracy, the results<br />
can soon diverge (Figure10), due to<br />
the sensitive feedback system which<br />
this proves to be.<br />
Figure9<br />
Figure10<br />
This sensitive dependence in<br />
the experiment shows how tiny<br />
discrepancies in the accuracy can<br />
result in differing outcomes and<br />
behaviours. This is best displayed with<br />
the corrected angles plot with time<br />
(Figure11).<br />
Demonstration of how identical states can quickly diverge<br />
into different final states. This was a weather simulation by<br />
Lorenz using only a few simple equations but the feedback<br />
loops caused this sensitive dependence.<br />
In Figure11 the top graph traces how the repeat run, the purple run, during the<br />
first five seconds almost identically follows the behaviour of the first run, the blue run.<br />
What this demonstrates is how the conditions in the start were (to a very high accuracy)<br />
the same. However, the middle graph follows how the two runs behaved in the<br />
proceeding five seconds, and how the two behaviours were starting to separate, shown<br />
by the waves becoming ever increasingly out of phase. Very soon after this difference<br />
had set in, it is clear to see how these two now slightly differing runs quickly become<br />
incoherent with respect to the other. During the fifteenth to the twentieth second<br />
(bottom graph), the two runs were in two completely different states, the blue run is<br />
following the motion of a normal pendulum, swinging back and forth, whilst the purple<br />
Page 10 of 18
<strong>Oliver</strong> <strong>Sheridan</strong>-<strong>Methven</strong> Chaos Supervisor - DL<br />
Pageites Word Count - 6460<br />
run is swinging completely through a full rotation in a regular cycle. Thus although even<br />
the two runs had the same initial conditions, after a period of time (due to their sensitive<br />
dependence and feedback loops found in their equations) had settled into two different<br />
states/rhythms.<br />
0-5s<br />
5-10s<br />
Figure11<br />
15-20s<br />
There are ultimately three ways that the pendulum could behave (Figure12)<br />
depending on the frequency chose. If the frequency is too high then it will only cause a<br />
rapid and small vibration of the pendulum and overall this effect will cancel itself out,<br />
thus it will act approximately to a normal pendulum with friction and will come to a halt.<br />
This would cause an implicit trajectory through the phase space. This would also happen<br />
with a frequency that is too low, and then once the pendulum has settled into a<br />
stationary state the driving force would just tilt the pendulum slowly from side to side.<br />
Overall this would still produce an implicit trajectory. The second possibility, which is the<br />
most familiar to people, is to produce a symplectic trajectory. This would be where the<br />
frequency chosen would produce a regular oscillation, just as a child does on a swing to<br />
maintain reaching a constant height. In phase space this would form regular orbits and<br />
loops without touching the constraints (borders) of the phase space. The final possibility<br />
would be to produce an explicit<br />
trajectory where the pendulum would<br />
Figure12<br />
"randomly" gain speed and lose speed<br />
at certain points in the rotation. If a child<br />
were to do this on a swing then they<br />
would get higher and higher and<br />
possibly suddenly slow down at one<br />
point in the swing, then slow down to a<br />
stop again, giving a chaotic<br />
appearance when observed. Also they<br />
have the possibility of swinging over the<br />
top of the swing and then back round<br />
Explicit Implicit Symplectic<br />
Page 11 of 18
<strong>Oliver</strong> <strong>Sheridan</strong>-<strong>Methven</strong> Chaos Supervisor - DL<br />
Pageites Word Count - 6460<br />
again and then returning to a normal pendulum action. In phase space this would be<br />
when the trajectory travels off the boundary in the phase space and then would return<br />
on the other side of the phase space. This is where the angle of the child changes from<br />
positive to negative once he crosses over the top of the swing. The purpose of correcting<br />
the angles was technically so that when a point went off of one boundary/edge, the<br />
corrected angle would cause the point to appear on the opposite side of the phase<br />
space.<br />
The transition from one of these trajectories into another type would be a form<br />
of birefringence of the system. Thus a further and more detailed aim was to see whether<br />
adjusting the frequency of the experiment's signal could produce these results. The first<br />
was to choose a very high frequency; which as explained earlier would have no overall<br />
effect but vibration. Thus the frequency was set to an excessively high frequency and<br />
was manually raised to a perpendicular position (90 0 ). Then it was released and the<br />
motion plotted and the<br />
trajectory through phase space<br />
traced (Figure13).<br />
What Figure13 displays<br />
is how the pendulum went from<br />
perpendicular and slowly swung<br />
through the centre four times,<br />
each time slowing down and<br />
then becoming stationary.<br />
Through phase space it traced<br />
the implicit trajectory, under the<br />
influence of the strange attractor<br />
at the origin. Therefore by<br />
adjusting the frequency the<br />
trajectory could be predicted.<br />
Thus at this frequency the<br />
chaotic system was easily<br />
approximated to a predictable<br />
one. Then the aim was then to<br />
construct the other two types of<br />
trajectories; symplectic and<br />
explicit, by causing a<br />
birefringence of the states<br />
through the variation of the<br />
signal frequency.<br />
Figure13<br />
The proceeding experiments were the runs where the frequency of the signal<br />
significantly affected the behaviour of the pendulum (i.e. not implicit trajectories);<br />
causing what appeared to be erratic behaviour. The first frequency (Figure14) was<br />
1.01Hz (there was a constant error in all frequencies of +0.01Hz) which initial caused an<br />
explicit trajectory (top left), behaving chaotically, until at around twenty seconds, it<br />
settled into an implicit trajectory (bottom left), following a regular pendulum motion.<br />
Overall the transition from chaotic into regular is best displayed by the angle time graph<br />
(top right) and gives the two overall trajectories (symplectic and explicit) through phase<br />
space (bottom right).<br />
Page 12 of 18
<strong>Oliver</strong> <strong>Sheridan</strong>-<strong>Methven</strong> Chaos Supervisor - DL<br />
Pageites Word Count - 6460<br />
Explicit<br />
Symplectic<br />
Explicit<br />
Figure14<br />
Symplectic<br />
The significance of the pendulum having the capability of making the sudden<br />
transition from an explicit trajectory to a symplectic trajectory is that it is the<br />
manifestation of the feedback loop and the sensitive dependence. Where there was a<br />
feedback loop, the system was easily repeating the previous trajectory again and again.<br />
However, where it was a sensitive system, when the system diverged from this trajectory<br />
slightly, it was attracted by the strange attractor significantly enough to make the<br />
transition of trajectories into the symplectic trajectory. This demonstrates how tiny<br />
variations have the potential to cause large differences in states.<br />
Once this sensitivity within the system was established, the intension is to<br />
increase the frequency to find a birefringence of trajectories or a region of stability. At<br />
the frequencies around 1.14Hz (Figure15) the phase spaces began to have the majority<br />
of the trajectory following the explicit path and the minority having a symplectic form.<br />
Figure15<br />
Symplectic<br />
Explicit<br />
Page 13 of 18
<strong>Oliver</strong> <strong>Sheridan</strong>-<strong>Methven</strong> Chaos Supervisor - DL<br />
Pageites Word Count - 6460<br />
Similar results were obtained with increasing frequencies until 1.31Hz was used<br />
(Figure16). When this frequency was applied an extremely regular explicit trajectory was<br />
formed. The pendulum swung very periodically through a full rotation in the same<br />
orientation. This was performed in a period four motion of the pendulum, due to the<br />
unique frequency that was able to resonate the pendulum, such that it produced a very<br />
regular oscillation pattern.<br />
Figure16<br />
The resulting phase space diagram in Figure16 represents how from the previous<br />
phase spaces, which appeared chaotic and without well defined trajectories (Fig14&15),<br />
a well selected frequency could produce a stable state. This is demonstrated by the<br />
extremely well defined trajectory which has great clarity (bottom graph in Figure16),<br />
forming a closed loop (joined at the edges of the phase space). However, how does this<br />
regularity arise from the chaotic birefringence<br />
and from amid the surrounding chaos? This can Figure17<br />
be said to be a unique solution/frequency<br />
which could be the one specially chosen to<br />
repeat this circuit of rotation; the<br />
correspondingly lower frequency is the one<br />
commonly used by a child on a swing [2].<br />
However, from the logistics map which<br />
parallels this system, only the chaos arose. It<br />
didn't appear that there was any stability past<br />
the first onset of chaos. However, on further<br />
inspection there are streaks of stability<br />
(Figure17), emerging as clear lines amongst the<br />
black. These can correspond to frequencies<br />
which induce stability [1,2,3,7], such as 1.31Hz.<br />
These clear streaks represent the unique<br />
frequencies which can resonate the pendulum<br />
into a regular rhythm, such as 1.31Hz.<br />
Page 14 of 18
<strong>Oliver</strong> <strong>Sheridan</strong>-<strong>Methven</strong> Chaos Supervisor - DL<br />
Pageites Word Count - 6460<br />
Figure18<br />
Computer Simulation by [5]<br />
Experimental Results<br />
So upon comparison of the<br />
experimental results to the computer<br />
simulation (Figure18) it is evident that the<br />
simulations can demonstrate many of the<br />
aspects of that occur in reality within the<br />
system. However, it would be difficult<br />
(almost impossible) to determine, from<br />
the four simple equations which govern<br />
the torque and the angular acceleration,<br />
the nature of the phase space<br />
trajectories and the strange attractors<br />
present. Secondly, it has proved to be a<br />
high quality manifestation of the logistics<br />
map and its qualities of birefringence and<br />
isolated regions of stability (1.31Hz) in such<br />
a relatively simply governed system<br />
(T1,T2,T3).<br />
Ultimately, this experiment<br />
demonstrates the capabilities that<br />
strange attractors have: that a plot in<br />
phase space can reduce the system's<br />
apparent randomness and<br />
unpredictability into order and patterns. From these patterns predictability can be<br />
induced which is the ultimate aim in experimental science; combining microscopic<br />
equations into macroscopic predictions, and the knowledge of what will happen next [2].<br />
Throughout the history of physics and natural philosophy though, there has<br />
always been one topic that used to take precedence over all others. During the days<br />
when religion and science were intermingled, the stars and planets were of divine<br />
significance. The movement of the heavenly spheres (planets) and the study of<br />
astronomy were the professions of the educated few (natural philosophers,<br />
mathematicians and some members of the church). Beginning with the theory of a<br />
geocentric universe, this was changed by Nicolaus Copernicus, proposing the<br />
heliocentric model for calculations in his book, On the Revolution of Celestial Spheres.<br />
Following this Galileo Galilei more fully supported this idea until Johannes Kepler<br />
published in Harmony of the World the collected data which supported, to the point of<br />
proof, the heliocentric model. Finally the issue appeared to be settled definitively in Isaac<br />
Newton's Mathematical Principles of Natural Philosophy. This had proved to be the most<br />
problematic issue: predicting where the divine bodies will appear at a given time [4].<br />
The key to Newton's success in resolving the issue definitively was to take<br />
Galileo's concept of inertia and extend upon it with mathematical description. These are<br />
known as Newton's three laws of motion which describe how things move. The first two<br />
laws can be described just by the second law; that when a force acts on a body its<br />
velocity will change. This is mathematically described by:<br />
F = dPdt -1 =ma<br />
The third law says for every impulse (F◦t) there is an equal and oppositely<br />
directed impulse. This law was not crucial though in his solution for a heliocentric model.<br />
Page 15 of 18
<strong>Oliver</strong> <strong>Sheridan</strong>-<strong>Methven</strong> Chaos Supervisor - DL<br />
Pageites Word Count - 6460<br />
The ultimate key to the success of the Newtonian model was his first and second law and<br />
the knowledge of Kepler's laws of planetary motion. Newton used Kepler's third law to<br />
deduce his law of universal gravitation:<br />
F = -GMmr -2 the negative sign denotes the force is inwards acting [9].<br />
From these laws there have been many derivations of how the planets move,<br />
which vary from the geometric proofs demonstrated in Principia and Feynman's Lost<br />
Lecture to algebraic proofs, using complicated calculus, such as in Mechanics and<br />
Vectors by Terry Heard [4,7,10]. The results from all of these are that the planets move in<br />
ellipses around the sun, the eccentricity of which can vary immensely; earth travels in a<br />
nearly circular orbit whilst Halley's Comet travels in a very eccentric orbit. Not only does<br />
theory predict this, but simulations support this (Figure19) and observations confirm this to<br />
a reasonable degree of accuracy.<br />
Figure19<br />
Here are two computer simulations using only the equations F = ma and F = -GMmr -2 which generate very accurate<br />
orbits (left) for earth and various comets (right).<br />
The key to the success in the predictions of these equations and simulations was<br />
the nature in which they were analysed. During all of these solutions there are only two<br />
bodies/planets, never more than two. The importance of this is that at all points in the<br />
orbit there is only once source of the gravitation force acting on any one body, which is<br />
the one produced by the other body. The significance is that during the geometrical<br />
analysis that all changes in velocity (acceleration) are in towards the other body directly.<br />
During the algebraic analysis, this means that the acceleration's component vector<br />
which is perpendicular to the axis connecting the two planets (a2) is equal to zero<br />
(Figure20). This is the fundamental principle<br />
from which all other derivations about the<br />
a 2<br />
planet's motion are made possible [4,7,10].The Figure20<br />
problem when three or more bodies of<br />
considerable mass are present is that the<br />
a 2 = 0<br />
principle that only centripetal acceleration<br />
exists (only a1) is frequently incorrect (i.e. a2 ≠<br />
m<br />
0). Consequently then the geometric<br />
derivations and algebraic demonstrations<br />
a 1<br />
which relied on only a2 existing are reduced<br />
to false models [7].<br />
M<br />
Page 16 of 18
<strong>Oliver</strong> <strong>Sheridan</strong>-<strong>Methven</strong> Chaos Supervisor - DL<br />
Pageites Word Count - 6460<br />
The issue with a three or more body problem is that a new set of nonlinear<br />
equations are created, and a chaotic system is born (Figure20). Mathematically the<br />
equations are all unsolvable and can only be guessed by iterative functions (e.g. "the<br />
disturbing function"). However, normally with the huge masses, distances, velocities and<br />
forces involved, these models soon become unreliable. This unreliability is what induces<br />
the sensitive dependence within the system whilst the iterative process causes the<br />
feedback loops. This chaotic system is demonstrated by the simulation in Figure21,<br />
provided by source [6].<br />
Figure21<br />
No Strange Attractors Exist<br />
The significance of three body problems having no strange attractors over the<br />
vast range of situations is that the ability of astronomy to produce reliable predictions<br />
becomes effectively impossible. When three bodies are present they could settle into<br />
some dynamic orbit within each other (which is irregular); they could collide with each<br />
other catastrophically! There is even the possibility that they call all collapse into each<br />
other and fling out one of the planets with a speed faster than the relative escape<br />
velocity. This makes efforts to forecast how our solar system formed and how it will<br />
continue, fundamentally ineffective. The exact course of our solar system, and even<br />
universe, can never be determined as a whole. No grand underlying pattern can<br />
possibly emerge and thus the accurate nature of the solar system, its future and past, is<br />
indeterminable [7].<br />
Over the course of science, the purpose of investigation and experimentation<br />
has been to discover underlying patterns. From these trends, physical laws are produced<br />
in ever increasing accuracy. Our ability to forecast the most intricate events has grown in<br />
strength. However, whilst our knowledge of the microscopic world has grown, that of the<br />
macroscopic has reduced. With ever increasing knowledge of how systems behave on<br />
their own, it emerges that the giant mesh of systems that exist within the world form an<br />
ever decreasingly predictably structure. Chaos has proven that this is the case for even<br />
the most basic of systems such as the chaotic pendulum; extending through to<br />
historically the most philosophically important questions about the planets and their<br />
motions. This new and emerging science of chaos and chaotic dynamics has forced us<br />
to question the knowledge which we have traditionally believed to govern how things<br />
behave, finely summarised by James Gleick; "Knowledge of the fundamental equations<br />
no longer [seems] to be the right kind of knowledge at all".<br />
Page 17 of 18
<strong>Oliver</strong> <strong>Sheridan</strong>-<strong>Methven</strong> Chaos Supervisor - DL<br />
Pageites Word Count - 6460<br />
Sources:<br />
[1] Chaos: A Very Short Introduction, Leonard Smith<br />
[2] Chaos, James Gleick<br />
[3] Introducing Chaos, Ziauddin Sardar & Iwona Abrams<br />
[4] Feynman's Lost Lecture, David & Judith Goodstein<br />
[5] http://myphysicslab.com/pendulum2.html - Chaotic Pendulum<br />
14/04/2010<br />
[6] http://myphysicslab.com/molecule3.html - Three Body Problem<br />
14/04/2010<br />
[7] Discussions with Ivan Hoffman DeVisme (IDV)<br />
[8] Discussions with Edward Reid (EJR)<br />
[9] Discussions with Dean Johnson (DRHJ)<br />
[10] Mechanics and Vectors, Terry Heard<br />
Page 18 of 18