As an appendix promised for the article How fast helium balloon floats…, and also appendix for the follow up article Hint for solving quadratic drag force on a damping oscillator, I’m going to show the phase plot or phase diagram of the oscillation of quadratically damped floating balloon.
Originally, I want to make it interactive, but I realized the computation is complicated, it might make your browser hangs. I am using GeoGebra instead to make some screenshots of it.
Our original equation of motions is this:
In the previous article we call as the damping coefficient for quadratic drag force. Because we are very clear now that this article only consider quadratic case, let’s just call the coefficient .
The equation of motion above have interesting property. It will converge to a stable position, and and , due to losing energy. That means, the center of the coordinate system is the position where the object converges and stop moving.
Any oscillation with linear restoring force and quadratic damping drag force , can be reduced into a generic equation above. If the equation have constant force , we shall find where A is the amplitude value that can be used to do a coordinate shift from the original reference frame to a reference frame where position 0 is where the object will be at rest.
So, it’s suffice for us to find the behaviour of above equation, because any similar motions can be reduced to our equation above.
The term in the drag force , the absolute value of , is just a convenient way for us to say that the direction of the drag force is always opposing the direction of the velocity. With the magnitude of drag force proportional to .
So, our equation above can be broken down into positive (positive x) direction and negative (negative x) direction.
We solved the phase space equation in the previous article. We also defined more variables that will make our life simpler.
The phase space relations conveniently becomes:
Some interesting insight from this. When we normalized the phase space variables and , it seems that the phase is deterministic. The constant and is entirely defined by the substitution of and . That means, for any such system with this behavior, we can transform the initial condition to fit such phase plot. Then find the solution and behavior, then transform back from this standard plot.
Now, let’s define our standard initial condition.
First, we define that our standard coordinate system is a position when the motion has converges.
Second, we define that our initial condition start when the object is displaced from , with maximum displacement. We can pick either in the left side or right side of origin. But I prefer to start from the left side. By starting from the left side, we will use the positive equation of motion (where velocity direction is positive).
If the distance between initial position and is (the maximum amplitude). Then the initial position . The initial acceleration becomes the maximum acceleration possible, and we can easily calculate it using Newton’s First Law in the initial condition. Or, shortly
After defining the initial condition above, we have completely mapped a phase from our standard phase space and into the actual phase space of the object and .
To illustrate, we are going to do a transformation . this will map into . From the initial condition, we got the constant
The constant for our standard phase is determined by it’s original position.
If we use Lambert W function to find at , we should find two values of . One at when the object began to move and another one at , which is when the object stop and change it’s velocity direction.
By the way, we got the above expression using Lambert W substitution: into , with .
Most of numerical function only implements the principal branch. But the equation above has two solution, which means it will use two branch of function, the principal and branch. Luckily for us, the function returns the other value, and not again. This is probably because the absolute value of has to be bigger than due to the damping.
With just this information we are able to plot the phase. The equation is:
It looks like this (plotted using GeoGebra):
As you can see that the maximum of is slanted a little bit left. We can calculate this position using the acceleration relation above, setting .
The second iteration of the phase is when the velocity is in the negative direction. So, we use the second set of the equation:
The initial value is the value that we find before: . From there, we calculate
Then, we can calculate .
The plot looks like this:
By this point, you probably have found some patterns.
Our equation basically tells us the same shape of phase, but it is rotated, and scaled, due to the Lambert W function behaviour.
Let’s do a small recap. Our most basic relation tells us this:
It will have two zeroes, will be negative and it’s absolute value is bigger than the positive
The peak is only determined by , which is .
is determined only by it’s initial position.
Or the other way around is also valid. It’s starting initial position is determined by .
The equation when is negative is actually the same equation but with a transform and . After that, the equation just alternates it’s sign with different value of .
We can reformulate which equation/constants to use based on which iteration number.
Note, that if we define , there is only one value possible, which is 1. We can imagine that the original amplitude is so far away at that it technically doesn’t cross axis. Let’s use this as a base.
We got the next , by flipping the equation, so the amplitude now is . I’m sure you got the pattern if you continue through multiple iteration. We got this relation:
Choice of index is completely arbitrary. With this kind of choice, that means, when is even, velocity is positive. When is odd, velocity is negative.
From this recursive relation, you can make a simple program that calculate all the necessary points of interests.
It’s also worth noting that subsequent values of is only determined by the initial values in the range of . If were to happen in the range it is probably not the first iteration, and you can determine the the original that is capable of generating that value.
Simply speaking, if you make a phase plot with then you will get the same phase plot if you draw the same plot, but starting at that is generated from . So, hypothetically speaking, if we have a canonical phase plot and the other phase plot can be transformed into this canonical phase plot, that means it is possible to define a standard function that will cover the whole dynamics.
Other insights from the above recursive relations is the maxima/minima of each cycle, from the value of . We can infer that eventually will converge to 1. In which the will converge to 0.
This is a plot of graph with with 5 iterations:
Same plot with but now with , also with 5 iterations as you can see, same plot is constructed:
We can go a little bit further deep. There is an interesting relation that we don’t have if the plot phase is a simple harmonic cycles. If we choose two arbitrary index of a different parity , with the same value of . We have an interesting relationship between positive and negative . Obviously in the different half cycle, it will corresponds with a different value that maps to the same . With simple harmonics, this is not interesting because the value will be the same, but with this plot phase, we will have extra identity. Let’s suppose that is the position where velocity is positive, meanwhile is the position where velocity is negative.
We just introduced a new variable and hyperbolic angle . For every same between consecutive half cycles, this relation emerges. The definition of the new variables:
We can do lots of things with this relation. For example, one of the things I immediately recognize is the term shows up twice in the identity above and the definition of . That means we can use the identity to approximate something using recursion.
We know that in the final iteration, the value of , because the peak has to be in the middle. We can then inspect every where . By setting either or to 1, then, quantifies how far the current cycle is with the final cycles. This is set unambigously because either Odd or even will corresponds to a positive value of . For example if then it must be due to the recursion definition above. If then it must be , causing .
If we compare with , then which is the maximum value. When , . We can immediately know which cycle is earlier or later. We can even compare different to get relative cycle between the two. The idea is very powerful, that we can think about as the absolute cycle lengths. Then and can be divided to get a relative cycles. We can even know, if then that means it is the very first initial cycle.
Now, we are going to try to interpret the hyperbolic angle . We can immediately recognize that it is possible to set the right hand side to zero. But the cycle length can not be zero. So the zero must comes from . The reason why it happens, maybe because of these:
First, whatever that is. But in order for that to be valid. It also must satisfy
Because we assumed that and from a different parity of , that means this can only happen if and it is in the final cycle. We actually proved otherwise. The then cannot be zero if we are inspecting two cycles from a different parity.
If we define, the absolute cycle length of cycle n as :
Then our previous is actually a square root/geometric average between two absolute cycles:
While can be defined as half the difference between two absolute cycles, and can be interpreted as hyperbolic angle difference:
The quantity can becomes a natural choice to understand which cycle is currently early or late if we compare them. We can retrieve the cycle’s hyperbolic angle from any absolute cycle length , either by calculating it directly with the natural logarithm, or by using Like what we understand before. But now, we can conveniently give it a name. It arises intuitively that when is small, that means the cycles is at it’s end and very similar.
From these results and understanding, we can then try to generalize these identities.
It follows that it turns out this identity is valid for any time that we observe two different cycles with the same .
However, remember that the sign alternates whenever if it’s not the same parity. Currently, the above identity holds if odd and even.
It is a direct result, that if we compare two different in the same even cycle (positive velocity) with :
In the odd cycle (negative velocity)
From both equation, we can’t determine if we set , which is make sense because the only time it happens is when .
If both are zeroes of the equation, then is the full range/length of of that cycle. We can call this the span . Meanwhile the quantity and turns out to be always positive due to the slanted shape of the phase plot. It is the average of absolute value of left and right amplitude difference. Let’s call this .
Regardless of if it is an odd or even cycle, is actually a good quantity to compare against difference cycles. If is small or near zero, that means it is closer to final cycle. If we remember our definition of absolute cycle length . Then the value is a relative cycle length between the same span of cycle .
With this property, we can easily scale any phase plot with a given span into another span. If we treat span as a continous variables (we don’t have any reason now, but for the sake of arguments, let’s pretend), we can take partial derivative. A span is invariant under origin translation. That means, in any given cycle, if we know the span and the motion constant , we can determine where the canonical origin should be from the value.
Next, we are going to observe the quantity .
The is an equivalent of but it is in the form of hyperbolic angle. Since we have :
It’s value may signify some importance with the hyperbolic geometry. The range of values of , depends on it’s span and amplitude difference in a given cycle . The minimum value , the maximum value . It’s very similar with the phase of circular function like, . I’m very suspicious. The midpoint of is a very special value where .
From, that we have an interesting hyperbolic identity, where any hyperbolic angle can be described with relative to the midpoint. Then, we have a quantity to for us to actually which phase the oscillation is currently in with arbitrary value.
The new quantity is the relative hyperbolic angle from the midpoint angle. From this form, we can intuitively understand that, the corresponds to how big the value of is compared with the current . If is negative, that mean is less than is magnitude, (also the other way around by inference).
You can check the graph of and . It is very convenient that can be used to show the phase angle of the cycle. Check that the zero value of does match with the maximum of .
Now, because we had a new definition of , we can add a new definition, analoguous with the relation and . Since is a normalized value, then we shall define a normalized absolute cycles , with the definition:
Here’s the graph of both and .
Because of how and cramped together, I only included the first three cycles. Exponential line above corresponds to . Exponential line below corresponds to . Each color of the line corresponds to the phase cycle plot directly below.
From what I observed in the graph. I understand that:
- The closer to cross , that means the closer the cycle to the final state. Cycle with is defined with everywhere (the furthest from final state)
- The closer to cross , that means the closer the cycle to the final state. Cycle with is defined have undefined normalized cycle length. (because division by zero). However we can intuitively says that if the value is bigger in the same phase, that means it is furthest from the final state.
From the phase and cycle plot, we summarized the following proposed definition and identities to understand the cycles.
Recursive relation, starting from
- The span is the length from to in a given cycle n.
- The amplitude is the maximum displacement from the origin point to a point where the velocity is zero in a given cycle.
- The absolute amplitude difference is the difference between the left amplitude and right amplitude of the given cycle, relative to the origin point.
The absolute cycle length quantifies how far the current cycle into the final cycle state .
Absolute cycle length preserves ordering. Whenever , if compared in the same value of X. Maximum value is .
Relative cycle length is a quantity to measure the ordering of the cycle , if compared in the same value of X.
Cycle average is a quantity that measures the average of two cycles.
Hyperbolic phase angle quantifies the phase of motion in the current cycle n. It can be used for ordering in the same cycle (but with some care about the parity).
The mid of the half cycle phase corresponds to when the acceleration is zero.
The minima and maxima of the phase in a given cycle is completely determined from the current span and amplitude difference .
The phase angle difference is the difference between two phases.
The normalized phase angle is a translated phase where the origin is . The normalized phase angle has a special property that the sign can also determined the sign of the equation of motions (combination of sign of acceleration, and velocity).
The normalized phase angle difference is the same phase angle difference, .
From the normalized phase angle, we can construct a normalized cycle length . The contour of preserves ordering between cycles.
For any and , with , with the same squared velocity , it follows:
Above equation assumes both odd parity cycle . If one of the arguments is even cycle, switch the sign of . If both is even cycle, switch both sign.
If expanded, odd b and a:
odd b, even a:
even b, even a:
even b, odd a:
If it is the same cycle, relation reduced to:
If both are zeroes in the same cycle, relation is simplified into:
Some defition have connecting relation:
The relation below happens in the same cycle n:
The relation below connects between , , and in the same cycle n: