Wave Propagation Part 1
Table of Contents
1. Introduction
In this series I am going to try to create a simplified model for longitudinal pressure waves propagation on fluids based on the interaction between its particles. Or main goal is to go step by step on the physics behind it, solving problems that gradually increase complexity, to have a theoretical preview of the relationship between the fluid's density and the wave's speed.
1.1. Disclaimer
I am nothing but a failed physicist physics enthusiast and the whole reason I am writing is a late attempt on trying keep the memories of some things I studied and don't really use as often as I would like to. Another reason is that I found out this cool thing called Emacs and Org mode that is making be want to go back to university again just because it is so nice to write on it!
I am not an specialist on this subject or in any subject at all and this might be all wrong. It is one of those cases that it is more about the trip then about reaching the proper destiny. If you have any real interest about this subject or any of the subjects treated on this blog and need it for an important application, you arrived on the wrong place.
Still it is not like this is all useless. You might be able to learn one thing or another and correct me, which I expect you to. Just don't expect this series to become source of any concrete scientific information or the models created here to be right. I am just as lost as you while creating this, which might be an advantage too. If I luckily end up with something interesting you have the registry of the full thought process. If I am not and this ends up being a frustrated crazy attempt, walking on circles where I just give up you will at least know exactly how not to think like when trying to solve a problem.
So I am writing this considering that someone else will read even though I am not exactly sure if it would be worth it. You are more than advised that there is no warranty that I am going to reach any conclusion at all with this or I might as well just switch to another problem that I find more interesting and never come back to this again. After all this is all just an experiment. This is nothing but a man trying to see how far he can go without having to research all around and figure out that any idea or problem he comes up with by himself was already solved by someone else before and it is public knowledge.
So I invite you to come and rediscover nature with me.
1.2. Motivation
We want to create a simplified mathematical model for wave propagation on fluids, longitudinal pressure waves to be more specific. We are specially interested on figuring out the relationship between a fluid's density and the wave propagation speed. We will start from simple problems that we can solve and extend on complexity and on number of dimensions from there.
Pressure waves propagate through the inter molecular repulsion forces, of electrostatic nature, that pull one by one in a chain effect. In real world fluids this is effect happens in some sort of 3D truss structure of interconnected particles, or a cloud, where the wave properly said starts with few particles receiving the initial pull and influencing multiple other neighbor particles and growing in a chain effect. Due to it's electrostatic we can consider classical Coulomb law with the force that decays with the inverse of the squared distance.
However the whole 3D effect would be too complex to begin with and due to that we will approach this in the most simple case. My goal with this series is to first solve this problem for the simplest case which is only the interaction of two particles. Instead of electrostatic forces let's being with a linear force. Lets begin with two masses connected by a spring that obeys Hook's law and then later see how having different potentials (how
2. Two masses connected by a spring
Consider two particles with the same mass
This is all an excuse to say that:
- The spring starts on its still state, not compressed or stretched.
- Particle 1 has an initial velocity
. - Both particles have the same mass.
This might not be much of a trivial problem as it might seem at first sight. Our goal, to say that we fully solved this problem, is determining the functions of time
Before tackling this one lets remember how the classical Harmonic Oscillator works first.
2.1. The Harmonic Oscillator
When a mass
Which is known as Hook's law, which is combined with Newton's second law there. This just means that the force applied to the object is proportional but with opposite direction of the spring's deflection. For instance if the spring is compressed, x is negative and then the force is positive, otherwise the spring is stretched, x is positive and the force is negative. So the forces varies linearly with x but it is always pointing contrary to it.
On this problem the other side of the spring is attached to a "wall". That simply means that it is an immovable object. So you can see this is an specific cause or at least an approximation of the problem we are trying to solve, where instead of two particles with the same mass we have something more like one particle with mass
Similarly to the problem described on the main section, to fully understand what is going on we need to determine
Which is the classical Differential equation of the Harmonic oscillator found on many problems including electronic circuits (RCL) besides this is a specific case of the wave equation which we might actually touch on later.
Anyway solving this (which I don't really remember the exact formal steps to do) is just a matter of asking: "What function derived twice becomes itself?", or minus itself, or… Well
What we just did there is a linear combination of the two solutions. In ODE's (Ordinary differential equations), you can add up any solutions together to have a new valid solution.
Another interesting fact is that you can add a phase difference, in this case represented by the constants
And you can now easily realize that when you take the second derivative you will be back with
One important thing to notice is that the frequency of the oscillation is already determined by the properties of the spring and the mass. In another words the boundary conditions will not affect the frequency of the movement only the amplitude and the phase. Isn't that interesting?
In case you are wondering what the heck is a boundary condition observe that when we formulated this problem we didn't talk about the particle having an initial velocity or position. Even if the spring were on its rest position, this isn't fundamental for throwing in the mechanics of this problem and solving it out, it would be still the same equation. That's where those constants come from, we can now adjust them according to whatever initial conditions we might have.
In the case of the more general problem of the two masses we are trying to solve, one boundary condition is having particle one with a
For example if
:Where
can be any integer number but we can simply use 0.for example if
and :
So the phase difference depends only of the initial position while the amplitude depends both on the position (
Using that back on the expression for A:
And
The period of the movement is the amount of time that it will take for a complete oscillation. That is whenever
And we can now say that we have fully solved the simple case of this mass + spring problem. One last thing though. You may have noticed that our general solution:
is a bit too limited. By that I mean that adding a phase to the cosine or sin part, that can be any value adjusted by the boundary conditions, you can convert sine to cosine then that phase is
Which is a sum of
2.2. Inertial Frame of Reference
So now let's try to solve the two masses problem. One argument that once came to my mind is that we could just use the particle 2 as a reference, and solve the problem from its view. It took me a while to realize that this would be completely wrong of an approach due to this thing called Inertial frame of reference. Since both particles are oscillating on this problem, particle 2 has a non zero acceleration and you can't use it as the main frame from where you make your measures from.
Rules of physics do not apply for non Inertial reference frames. How do you know if something is not a Inertial reference frame? Well if there is an inertial reference for which that one has non 0 acceleration then that is not inertial. Actually then this one will have a non 0 or varying acceleration for all other referential.
So there isn't quite much that we can use from the harmonic oscillator problem. This proves that one particle will not see the other as if it was the wall on the Harmonic Oscillator problem, even though that will be a tendency if one of the masses tends to much bigger than the other. Still that is no the problem we want to solve since it makes more sense that the particles have about the same mass for a fluid.
2.3. The ODE System
So what we are looking for is a system of differential equations that will represent both positions
Where
This doesn't change the problem at all, we are just offsetting the frame of reference in which we measure
So what we have now is:
As
Notice that:
This means, if the spring is compressed then
Which is the system we've been looking for. To solve this let's begin by subtracting the latter by the first:
Which is basically the Harmonic Oscillator equation that we already know the solution. In this case we will have:
We can also add both equations to have:
Where
Notice how they are in a phase difference of
2.4. Boundary Conditions
So we have solved the mechanics for the physics problem but now to have a better hint about the analogue wave propagation we are trying to make it would be interesting to bring up the boundary conditions since we have 4 constants to figure out on this case:
Particle 1 has an initial velocity:
The spring is initially at its rest position:
Particle 2 is at rest initially:
Since particle 1 begins with the maximum of its speed at
Also:
Similarly for
Notice that this just correcting for our
Also:
Combining this with the one for
So:
So in the end we finally have:
Where
Which fully solves the described problem!
3. Conclusion
We could have gone much far considering particles with different masses and other scenarios like "what happens if there is friction?", but for the original goal of studying longitudinal wave propagation there wont be different masses. What we should do next is try to solve this problem for 3 particles chained together, and then, try to expand for t
In any case I'm thinking on bringing up an animation of the problem we just solved maybe using something like manim. If I figure out a nifty way of hosting the gifts that would be cool.
Comments
Be the first to comment!