• If you are citizen of an European Union member nation, you may not use this service unless you are at least 16 years old.

  • You already know Dokkio is an AI-powered assistant to organize & manage your digital files & messages. Very soon, Dokkio will support Outlook as well as One Drive. Check it out today!

View
 

Damped oscillators - the math (technical) (2013)

Page history last edited by Joe Redish 4 years, 10 months ago

Working Content > Oscillations and Waves > Harmonic oscillation > Mass on a spring > Damped oscillators

 

Prerequisites:

 

Putting a viscous damping term into Newton's 2nd law for a harmonic oscillator led to a more complex equation for the position as a function of time:

.

This relates the second derivative of x to the value of x and its first derivative (v). There are general tools to solve this in courses on differential equations, but we can get to the solution by using some physical intuition and a single rule from Calc I: the product rule. This specifies that if we take the derivative of a product of two functions, f(t) and g(t), the derivative of their product is given by the following:

 

 

Now let's apply some physical intuition. When we had no damping, the result was simply a cosine oscillation (possible shifted though that isn't shown here because we are focusing on the time dependence):

From our observations we see that the amplitude doesn't stay constant but keeps getting smaller. So let's hypothesize that it still oscillates (though possibly with a slightly different frequency) times a time-dependent amplitude:

Now we have a solution that is a product of two time dependent functions: the amplitude and the cosine term. What we now need to do is differentiate this to get the velocity and then differentiate that to get the acceleration, put the results into our F = ma equation and see what happens. This only uses the product rule, but it's messy (but straightforward if we are careful).

 

When we differentiate this once, we get the result for the velocity:

 

 

Differentiating again to get the acceleration and using the product rule on each of our two terms gives four terms, two of which are the same. Collecting these together gives:

 

Now if we put our acceleration expression and our velocity expression into our F = ma equation we get

 

This looks awful, but we notice that there are cosine terms and sine terms. If we put all the cosine terms together on one side and all the sine terms together on the other side we get:

 

This still looks intractable.  How can we possibly solve it? The easiest way is to see if we can make both of the square bracket terms vanish. Then both sides would be equal to 0 and therefore equal to each other. The right side looks simpler so let's try it first:

In the first equation the constant ω1 was in both terms so it cancels. We then solved for dA/dt and found that it is proportional to A times a constant. The function whose derivative is proportional to itself is the exponent. The result is:

 

Now, since we have decided we know A we can deal with the left side of our messy F=ma equation:

where we have just put in that we know taking the derivative of A just brings down a factor of /2. Each term has a factor of A so we can cancel that giving

 

 

Our values of γ and ω0 are given in terms of the physical parameters, but ω1 is an unknown. So we can solve for it giving:

This gives the full solution to the damped oscillator (as long as γ/2 > ω0):

 

Of course our most general solution might not start at a maximum when t = 0. So in general we have to add a phase shift giving a general solution for our damped oscillator model:

This expresses our solution in terms of the physical parameters of our system, ω0 and γ (which, you should recall, are combinations of our physical parameters k, b, and m) and two arbitrary constants -- A0 and Φ. These arbitrary constants allow us to describe a range of physical behaviors with different starting amplitudes and different starting points (initial velocities).

 

Afterthoughts: So what just happened?

Well that was strange. It didn't really look like anything you might have done in an introductory biology or chemistry class, and it didn't even look like what you might have done in a math class -- even, if we are honest about it, what you might do in a math course in ordinary differential equations where these kinds of problems are traditionally solved. We've used a  process that blends scientific modeling with mathematical structures, very much in the way we discussed it in our page Modeling with Mathematics at the beginning of this class. The idea there was that in mathematical modeling in science one starts with some knowledge of the physical world from experience, then maps it into a mathematical model where there are a lot of structures available to guide our reasoning about relationships. Then, the equations we generate are solved by doing math, but constrained and assisted by physical knowledge and interpretations. [Note that when we say "physical" here, we are referring to "the world of matter and energy" and this includes not only physics, but chemistry and biology as well.]

 

In the case of the damped oscillator, we began with physical principles -- the Newtonian framework that in the particular example we were considering led us to create a physical model (through Newton's second law) that produced a (second-order linear) differential equation: an equation that related an unknown function x(t) to its derivatives, dx/dt and d2x/dt2. Now in math classes in Ordinary Differential Equations there are algorithms taught for solving these kind of equations. They involve exponentials and complex arithmetic, and they are guaranteed to give all possible solutions. We didn't use these methods. Instead, we did something blended.

 

We had observed the phenomenon -- the damping of a physical oscillation -- and we had an idealized solution: the undamped oscillator. Physically, the more realistic case looked like it was just losing amplitude but still oscillating. So we conjectured a time-dependent amplitude and tried to see what it would imply. We wound up using a lot of math -- the product rule in calculus and a lot of equation manipulations.

 

Our blend of physics and math led us to a relationship that had a function of our unknown amplitude times a cosine equal to a different function of our unknown amplitude times a sine. Now we could have been stuck here. There is a function that when multiplied by a cosine gives a sine: the tangent (since tan(θ) = sin(θ)/cos(θ)). But that would be weird. The tangent would oscillate and go to infinity (whenever cos was 0). We don't expect an amplitude to go to infinity or look anything like this. So we rejected the possibility that our unknown multipliers transformed a cos into a sin (by not even mentioning this possibility). Instead, we tried making each side 0. This worked nicely giving us a result that made physical sense.

 

But can we trust this? Couldn't we have done something different and come up with another answer? Here, we go back to the math structure that is supporting our analysis. We can find theorems that tell us that if we have a second order linear differential equation, the general solution will have two arbitrary constants that can be set freely. The fact that we have such a solution means that we have the general solution. There is only one (though it may be expressed in a number of different looking but equivalent ways). The theorem from math tells us that our physically motivated procedure worked mathematically -- in addition to making good physical sense! Knowing that this sort of theorem holds is what gives us the license to play -- to guess physical solutions on the basis of our physical intuitions and see whether the result does the math for us. It often does! Of course we still have to check with experiment and see how well the mathematical hypothesis we started with describes real physical systems. (Sometimes very well, sometimes not.)

 

This "blending" procedure, that makes use of physical and mathematical arguments intertwined, is quite typical of the way mathematics is used to model systems in science.

 

Joe Redish 3/23/12

Comments (0)

You don't have permission to comment on this page.