Calculating Velocity in the Millikan Oil Drop Experiment

Click For Summary
The discussion centers on calculating the velocity of an oil drop in the Millikan oil drop experiment, focusing on the differential equation governing its motion. Participants highlight the challenge of integrating the equation due to the presence of velocity on both sides, emphasizing the need to separate variables correctly. Several methods for solving the differential equation are proposed, including using definite integrals to simplify the integration process and eliminate constants. There is also debate over the correctness of solutions found online, particularly regarding the exponential terms involved. Ultimately, the conversation underscores the importance of careful algebraic manipulation and understanding the physical context of the experiment to derive the correct velocity function.
  • #31
Charles Link said:
What you have in post 12 from Wikipedia I believe is incorrect. That's why it didn't work.
In post 15, you introduce a definition for ##m ##, where ## m ## is already in the problem. You need to use some other letter.
In post 15 I didn't change m to lighten the equation. m is only the mass.
The issue is g.
the solution should be

##v(t) = \frac{mg}{6\pi n a }(1 -e)##

I don't know if you see what I mean.
 
Physics news on Phys.org
  • #32
EpselonZero said:
In post 15 I didn't change m. The issue is g.
the solution should be

##v(t) = \frac{mg}{6\pi n a }(1 -e)##

I don't know if you see what I mean.
It is not clear at all what you are computing in post 15. This differential equation is a rather simple one that can be solved in a routine manner by a couple of different methods. The method I like is given in post 28. The method in post 13 will also work.
 
  • #33
I'm using the method in post 13. It should work I guess. But I can't get the solution.
 
  • #34
EpselonZero said:
I'm using the method in post 13. It should work I guess. But I can't get the solution.
##-(du/dt)/B=u ## is the modified D.E. with the substitution ##u=A-Bv ##, so we get ## du/u =-B \, dt ##, so that
## \ln|u|=-Bt ##, so that ## u=Ce^{-Bt} ##.
## u=A-Bv ##, so that ##v=(A-u)/B ##.
Note: ## B=6 \pi \eta a /m ##.

The solution in Wikipedia is in error, (Edit: See post 41,Wikipedia solution is correct). so you may have been trying to get a solution that is incorrect.
The solution I get, by both methods is ## v(t)=\frac{A}{B}(1-e^{-Bt}) ##. See if you can get your calculations to agree with this.
 
Last edited:
  • #35
I'll retype my steps

##\int \frac{-du}{Bu} = \int dt##
##u = A-Bv, du -B dv##

##\frac{-ln |A -Bv|}{B} + C_1 = t + C_2 ##

##C_1 - C_2 = C##
##A = \frac{4 \pi a^3(p-p_1)g}{3m}, B = \frac{6 \pi na}{m}##

##-ln|\frac{ 4 \pi a^3(p-p_1)g}{3m} -\frac{6 \pi na}{m}v| + C = t\frac{6 \pi na}{m} ##

##\frac{ 4 \pi a^3(p-p_1)g}{3m} -\frac{6 \pi na}{m}v + C = e^{-t\frac{6 \pi na}{m}} ##

##m = \frac{4 \pi a^3(p-p_1)}{3}##

##\frac{ 4 \pi a^3(p-p_1)g}{3 \frac{4 \pi a^3(p-p_1)}{3}} -\frac{6 \pi na}{ \frac{4 \pi a^3(p-p_1)}{3}}v + C = e^{-t\frac{6 \pi na}{ \frac{4 \pi a^3(p-p_1)}{3}}} ##

##g -\frac{9 n}{2 a^2(p-p_1)}v + C = e^{-t\frac{6 \pi na}{ \frac{4 \pi a^3(p-p_1)}{3}}} ##

## -\frac{9 n}{2 a^2(p-p_1)}v + C = -g + e^{-t\frac{6 \pi na}{ \frac{4 \pi a^3(p-p_1)}{3}}} ##

## v(t) = \frac{2 a^2(p-p_1)}{9n} (g - e^{-t\frac{9 n}{2 a^2(p-p_1)}}) + C ##
 
  • #36
Line 7 is incorrect. ## e^{C-Bt}=C'e^{-Bt} ## is what you get on the right side.
Line 8 is also incorrect. You are redefining ## m ##.
The first 6 lines are good. Try working more carefully from there.
 
  • #37
Charles Link said:
Line 8 is also incorrect. You are redefining ## m ##.

I don't understand. m is the mass of the oil drop. I have to replace m for his mass.

The final solution should be

## v(t) = \frac{2 a^2g(p-p_1)}{9n} (1 - e^{-t\frac{9 n}{2 a^2(p)}}) ##
 
  • #38
EpselonZero said:
I don't understand. m is the mass of the oil drop. I have to replace m for his mass.
The expression with ## V_s(\rho-\rho_1)g ## has the difference in densities, which makes for the weight minus the archimedes type buoyant force. ##V_s =(4/3) \pi a^3 ## is the volume of the sphere.
## \rho_1 ## is the density of the air which is the displaced fluid in archimedes terms.

## (4/3) \pi a^3 (\rho-\rho_1) ## has the units of mass, but it is a reduced mass and is not ## m ##.
 
  • #39
Ah make sense, so ##m = \frac{4 \pi a^3}{3}p##

V(0) = 0 when C = g, right?
 
  • Like
Likes Charles Link
  • #40
I don't think ## C=g ##, but otherwise correct.
You are ok above through line 6. Try completing it with the new info.
 
  • #41
Going from line 6 to line 7 you made an algebraic error. Once you correct that, you should be able to get the correct answer. (You take ## e^x ## of both sides). ##e^{a+b}=(e^a )(e^b) ##.
and yes, your solution in post 37 (from Wikipedia) is correct, with ## \rho=3m/(4 \pi a^3) ##, and is in agreement with mine.
Note: I can see now why they wrote the solution how they did: It contains simply the radius of the drop, ##a ##, and they wanted to have ## a ## as the one parameter for the size of the drop=thereby, they eliminated the ## m ## in the exponent, and replaced it with ## (4/3) \pi a^3 \rho ##.
 
Last edited:
  • #42
I hope I don't waste your time.

I'm still stuck.

I get ##A - Bv = C'e^{-Bt}##

Then I replace A,B and m.

##\frac{(p-p_1)g}{p} - \frac{9nv}{2a^2p} = C' e^{\frac{-t9n}{2a^2p}}##

## v = \frac{2a^2p}{9n} (\frac{(p-p_1)g}{p} - C' e^{\frac{-t9n}{2a^2p}})##
 
  • Like
Likes Charles Link
  • #43
## C'=A ## because ## v=0 ## at ## t=0 ##.
I think you almost have it...
Yes, I believe you got it. The p's will cancel in numerator and denominator.
Very good. :)
 
Last edited:
  • Like
Likes happyparticle
  • #44
Wow, really? Sadly, I don't see how you can cancel the p's.

## v(t) = \frac{2a^2pg}{9n} (\frac{(p-p_1)}{p} - \frac{(p-p_1)}{p}e^{\frac{-t9n}{2a^2p}})##
 
  • #45
EpselonZero said:
Wow, really? Sadly, I don't see how you can cancel the p's.

## v(t) = \frac{2a^2pg}{9n} (\frac{(p-p_1)}{p} - \frac{(p-p_1)}{p}e^{\frac{-t9n}{2a^2p}})##
Factor out ## A ##, and you have exactly what you need, with ## v(t)=(\frac{2a^2 pg}{9n})(\frac{p-p_1}{p}) (1-e^{etc}) = \frac{2 a^2 g (p-p_1 )}{9n}(1-e^{etc}) ##.
The p's cancel. :)
 
  • Like
Likes happyparticle
  • #46
Charles Link said:
Factor out ## A ##, and you have exactly what you need, with ## v(t)=(\frac{2a^2 pg}{9n})(\frac{p-p_1}{p}) (1-e^{etc}) = \frac{2 a^2 g (p-p_1 )}{9n}(1-e^{etc}) ##.
The p's cancel. :)
You are right, after all this I didn't remember what I was looking for.

Thanks Charles and everyone else who help me. Finally! I'm going to sleep.
 
  • Like
Likes Charles Link
  • #47
EpselonZero said:
You are right, after all this I didn't remember what I was looking for.

Thanks Charles and everyone else who help me. Finally! I'm going to sleep.
Good, but there is an unresolved question.
As @hutchphd pointed out, Millikan had no need for such an equation because he only dealt with two scenarios: terminal velocity and stasis. So is the reference to Millikan in the thread title misleading?
 
  • #48
For whatever reason I don't have this in the wikipedia english page, but in french you can see here what I was looking for.

My question is titled the Millikan experiment. Ultimately, I have to find for which value of q, v=0.
 
  • #49
EpselonZero said:
For whatever reason I don't have this in the wikipedia english page, but in french you can see here what I was looking for.
The French page finds the general solution, as you have in this thread, but then takes t to infinity to find the terminal velocity. The English page avoids solving the ODE by just setting the acceleration to zero. I confess to sympathy for the English sloth.
EpselonZero said:
My question is titled the Millikan experiment. Ultimately, I have to find for which value of q, v=0.
Again, you don't need to solve the ODE for that.
So unless you were specifically instructed to find the general solution of the ODE I would say it was unnecessary - but a valuable exercise.
 
  • #50
Yeah, I had to find the general solution and as you said it's probably because it's a valuable exercise.
 
  • Like
Likes Charles Link
  • #51
Good job by @EpselonZero for staying with it until it was solved. Suggestion is to also understand the solution given in post 28. In general, this is a very standard differential equation with a solution to the homogeneous equation, (where the right side is set to zero), that comes with an arbitrary constant ## C ##, along with the particular solution. If you haven't yet had a course in differential equations, you would do well to take one as soon as you can.
 
  • #52
Alright, I'll try it.
I have done a course in differential equations, but I forgot a lot of things.
 
  • Like
Likes Charles Link
  • #53
Don't worry. Practice, practice and practice.
 
  • Like
Likes Charles Link

Similar threads

Replies
12
Views
2K
Replies
6
Views
2K
  • · Replies 11 ·
Replies
11
Views
3K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 10 ·
Replies
10
Views
6K
  • · Replies 5 ·
Replies
5
Views
2K
Replies
2
Views
3K
Replies
1
Views
3K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 2 ·
Replies
2
Views
4K