2009 Lectures
2009 Lectures
ill be
l um w its’. August-November, 2009
r ricu 0 ‘un
0 1 cu in1
d
PH1 covere Dr. Deshmukh, P.C.
)
(un
Dr. Kadhane, U.R.
Dr. Harish Kumar, N.
Dr. Vijayan. C.
Department of Physics
Indian Institute of
Technology Madras
Chennai 600036
[email protected]
[email protected]
[email protected]
[email protected]
PCD-09
A very hearty welcome to
‘PH101: Physics I’ course!
PCD-09
Alternative way to arrive at Newton’s III law,
following a very different approach
- the principle of translational invariance
in homogeneous space.
PCD-09
In contemporary physics,
SYMMETRY is placed ahead of LAWS OF
NATURE.
PCD-09
w ill
lu m
ricu
cu r
1 0 1 0
P H i n 1
d
be covere
n )
(u
s ’.
n it
‘ u
Today’s lecture:
Brief overview of the course contents of the
Ten Units
PCD-09
Course material is available at:
http://www.physics.iitm.ac.in/~labs/amp/
ATOMIC & MOLECULAR PHYSICS
Department of Physics,
Indian Institute of Technology - Madras
"One must always
do what one really
Dr. P. C. Deshmukh cannot."
Professor
Department of Physics, IIT Madras
- Niels Bohr
(1885-1962)
Teaching Courses: PH101 pdf files
uploaded (lectures, resources etc.)
Fragmentary Tale of The Atom
•
PCD-09
Big Bang
‘SCIENCE’ is ~ between
12 and 14
relatively young! billion
years ago
Just a few hundred
years old, perhaps
The Solar
a few thousand System ~
years or so, 4.5 billion
years old
but not older!
Humans ~
few million
years
Time
PCD-09 http://www.ccsf.edu/Departments/History_of_Time_and_Life/PDFs/BigBang24x24.pdf
supernova: 11 Billion Years Old!
PCD-09
Why do obj
ects ‘fall’ ?
PCD-09
Why do objects ‘fall’ ?
PCD-09
Albert Einstein: “We owe a lot to Indians, who
taught us how to count, without which no worthwhile
scientific discovery could have been made.”
http://www.physics.iitm.ac.in/~labs/amp/kerala-astronomy.pdf
PCD-09
Unit 1 : Equations of Motion.
Principle of Causality and Newton’s
I & II Laws.
∫
t1
L ( q , &
q , t )dt = extremum Principle of Variation
PCD-09
Central problem in ‘Mechanics’: How is the
‘mechanical state’ of a system described,
and how does this ‘state’ evolve with time?
‘position’ and ‘velocity’: both needed
to specify the mechanical state of a system?
Causality
&
Determinism
I Law II Law
ur r
F = ma Effect is proportional to the Cause.
Linear Response. Principle of causality.
PCD-09
Unit 1 :
ur r
F = ma Linear Response Formalism
Cause-Effect Relationship
t2
∫
t1
L ( q , &
q , t ) dt = extremum Principle of Variation
Lagrange’s equations:
One 2nd Order equation
Hamilton’s equations: ,
Two first order equations
PCD-09
“When we ask advice, we are usually
looking for an accomplice. ”
Joseph-Louis Lagrange
(1736-1813)
⎛k⎞
x = −⎜ ⎟ x
&&
⎝m⎠
&& 1
Q = −( )Q
LC
PCD-09
Resonances Enrico Caruso
1873 - 1921
Enrico Caruso was said to be able to shatter
a crystal goblet by singing a note of just the
right frequency at full voice.
In 2005, the Discovery Channel television
show MythBusters recruited rock singer
and vocal coach Jamie Vendera to hit some crystal
ware with his best shot.
http://www.sciam.com/article.cfm?id=fact-or-fiction-opera-singer-can-shatter-glass&sc=rss
PCD-09
Tacoma Narrows Bridge (Washington state), 1.9 km
length, one of the largest suspended bridges built at the
time.
PCD-09
http://www.math.harvard.edu/archive/21b_fall_03/tacoma/index.html
Wave motion in one
dimension.
Wave equation and
travelling wave
solutions.
PCD-09
Unit 3 :
Y
Plane polar
and
cylindrical polar
coordinate
systems.
eˆϕ
eˆρ
eˆy
ϕ eˆρ = cos ϕ eˆx + sin ϕ eˆy
O X eˆϕ = − sin ϕ eˆx + cos ϕ eˆy
eˆx
⎛ eˆ ρ ⎞ ⎛ co s ϕ sin ϕ 0 ⎞ ⎛ eˆ ⎞
⎜ eˆ ⎟ = ⎜ − sin ϕ co s ϕ
⎟ ⎜ ˆx
0 ⎟ ey ⎟
⎜⎜ ϕ ⎟⎟ ⎜ ⎜⎜ ⎟⎟
⎝ eˆ z ⎠ ⎜⎝ 0 0 ⎟
1 ⎠ ⎝ eˆ z ⎠
PCD-09
Unit 4: Kepler Problem.
Laplace-Runge-Lenz vector, ‘Dynamical’ symmetry.
Conservation principle ↔ Symmetry relation.
Kepler Problem.
Laplace-Runge-Lenz vector, ‘Dynamical’ symmetry.
Conservation principle ↔ Symmetry relation.
b
a
Rosette motion
For this potential and
associated field:
no precession of orbit.
PCD-09
ur
×L
Laplace Runge Lenz Vector is ê ρ ur
p
constant for a strict 1/r potential.
ur ur ur
A = p × L − mkeˆρ ur
r ur A
ur ⎡ v× L ⎤ r
A = mk ⎢ − eˆρ ⎥ = mka
⎣ k ⎦ uur
r ur dA
r v× L For = 0,
a= − eˆρ r
dt
k dp k
r = − 2 eˆρ
LRL vector is often defined as a dt ρ
DYNAM ICAL
which is also, obviously, a constant.
r SYM M ETRY
One can show that a = e, Eccentricity
vector
SYMMETRY CONSERVATION
PRINCIPLE
Homogeneity of time Energy
x
y = r sin θ sin ϕ
x = r sin θ cos ϕ Y
ϕ
0≤r<∞
0 ≤θ ≤π x = r sin θ cos ϕ
ρ = r sin θ
0 ≤ ϕ < 2π y = r sin θ sin ϕ
X
z = r cos θ
PCD-09
Transformations of the unit vectors
PCD-09
Unit 6: Inertial and non-inertial reference frames.
Moving coordinate systems. Pseudo forces.
Inertial and non-inertial reference frames.
Deterministic cause-effect relations in inertial frame,
and their modifications in a non-inertial frame.
PCD-09
ur
ω Western
North Pacific
Coriolis
1792 - 1843
o n th e N or th er n
te o f fr e e fall
An object in a
s t a
r d th e E a s t, i n the
t s d e f l e c t ed towa ar d We st!
r e g e lected t o w
hemisphe er e w o ul d g et d e f
er n h e m i s p h
South
g t hr ough 1 0 0
0 a n o b j e c t f a llin
a l a ti t ud e o f 60 .
At thr ou g h ~ 1 c m
is de f l e c t e d
meters
Western
South Pacific
PCD-09 http://agora.ex.nii.ac.jp/digital-typhoon/start_som_wnp.html.en
Unit 7: Galilean & Lorentz transformations.
Special Theory of Relativity.
1
γ=
v2
1- 2
c r
v
γ → 1 as v → 0.
y dψ r
= u • ∇ψ
ˆ
ds
h uˆ = lim
uur uur
δ r dr
=
δ s→0 δ s ds
uur
h: dependent δ s = δ r , tiny
variable
increament
x & y:
independent uur
variables ds = dr , differential
ur ur r ur r uur
∫∫∫
volume
dτ ⎡⎣∇ • A(r ) ⎤⎦ = ∫∫
surface
A(r ) • da
region enclosing
that
region
PCD-09
Unit 10: Fluid Flow, Bernoulli’s Principle
Ishant Sharma
Swing Bowling
Equation of motion for fluid flow.
Definition of curl, vorticity,
irrotational flow and circulation.
Steady flow. Bernoulli’s principle.
Introduction to applications of Gauss’ law
and Stokes’ Theorem in Electrodynamics.
∂ρ (r , t )
∇ • J (r , t ) + =0
∂t
PCD-09
Curl of the vector
r r r r r uur
∫ A ( r ) • dl = ∫∫ ( ∇ × A ) • dS
William Thomson,
1st Baron Kelvin
(1824-1907)
This theorem is named after George
Gabriel Stokes (1819–1903), although
George Gabriel the first known statement of the theorem
Stokes is by William Thomson (Lord Kelvin) and
(1819–1903)
appears in a letter of his to Stokes in
PCD-09
July 1850.
2
p(r ) v
Ψ= +φ + = constant for a given streamline
ρ 2
To read about Bernoulli family, refer to:
http://www.york.ac.uk/depts/maths/histstat/people/bernoulli_tree.htm
http://library.thinkquest.org/22584/temh3007.htm
ur r uur
∫∫ B(r ) • dS = 0 There is no
Magnetic Monopole
ur r uur ∂φB
∫c E ( r ) • dl = −
∂t
Faraday-Lenz-Maxwell
ur r uur Ampere-Maxwell
∂φE
∫c B(r ) • dl = μ0ε 0 ∂t + μ0 I enclosed Next
Fluid Dynamics, Electrodynamics, class:
Maxwell’s equations:
Another course!! Unit 1
PCD-09
Unit 1: Equations of Motion
Lagrangian/Hamiltonian formulation.
Application to SHO.
1
PCD-09
ATOMIC & MOLECULAR PHYSICS
Department of Physics, Indian Institute of Technology - Madras
Chennai, India - 600 036 Email : [email protected]
Course material is available at:
http://www.physics.iitm.ac.in/~labs/amp/ "One must
Dr. P. C. Deshmukh always do
what one
Professor
really
Department of Physics, IIT Madras cannot."
Affilialiation
•
- Niels Bohr
Publications (1885-1962
Teaching Courses: PH101 pdf files uploaded (courses
lectures, resources etc.)
Sponsored Projects
•
Collaboration
•
[email protected]
2
PCD-09
Learning goals:
3
PCD-09
Learning goals:
4
PCD-09
Big Bang
‘SCIENCE’ is ~ between
12 and 14
relatively young! billion
years ago
Just a few hundred
years old, perhaps
The Solar
a few thousand System ~
years or so, 4.5 billion
years old
but not older!
Humans ~
few million
years
Time
5
PCD-09 http://www.ccsf.edu/Departments/History_of_Time_and_Life/PDFs/BigBang24x24.pdf
supernova: 11 Billion Years Old!
6
PCD-09
Why do obj
ects ‘fall’ ?
7
PCD-09
“Gravity is not
responsible for people
falling in love”
- Albert Einstein
8
PCD-09
“We owe a lot to the Indians, who taught us how
to count, without which no worthwhile scientific
discovery could have been made.” - Einstein
Scientific outlook in early days……….
ARYABHATTA: first to propound the theory that
earth has a spherical shape. (5th Century)
the famous Indian tradition…..
BRAHMAGUPTA: Estimated, with significant
accuracy, the circumference of the earth.
(7th Century)
9
PCD-09
This paper is available at:
http://www.physics.iitm.ac.in/~labs/amp/kerala-astronomy.pdf
Nicolus Copurnicus
1473-1543
10
Central problem in ‘Mechanics’: How is the
‘mechanical state’ of a system described, and how
does this ‘state’ evolve with time? Formulations due to
Galileo/Newton, Lagrange and Hamilton.
Why ‘position’ and ‘velocity’ are both needed to specify the
mechanical state of a system?
q
What is meant by ‘Equation of Motion’ ?
Galileo Galilee
1564 - 1642 13
PCD-09
What is ‘equilibrium’?
What causes departure from ‘equilibrium’?
Galileo Galilei Isaac Newton
1564 - 1642 (1642-1727)
Causality
&
Determinism
I Law II Law
ur r
F = ma Effect is proportional to the Cause.
Linear Response. Principle of causality.
14
PCD-09
ur r
F = ma Effect is proportional to the Cause.
Linear Response. Principle of causality.
15
From Carl Sagan’s ‘Cosmos’:
“Like Kepler, he (Newton) was not immune to the
superstitions of his day… in 1663, at his age of twenty, he
purchased a book on astrology, …. he read it until he
came to an illustration which he could not understand,
because he was ignorant of trigonometry. So he
purchased a book on trigonometry but soon found himself
unable to follow the geometrical arguments, So he found a
copy of Euclid’s Elements of Geometry, and began to read.
Two years later he invented the differential calculus.”
16
PCD-09
f ( x)
Differential of a function.
Derivative of a function
Slope of the curve.
x
It is a quantitative measure of how sensitively the
( x ) responds to changes in the
function ff(x)
independent variable x .
δf
f ( x0 )
⎡ df ⎤ δx
δ f =⎢ ⎥ δx
⎣ dx ⎦ x0
x
Tangent to the curve. x0
18
Dimensions of the derivative of the function: [f][x]-1
PCD-09
One often has to deal with functions of many variables,
such as the height above a flat horizontal surface of a
handkerchief that is stretched out in a warped surface,
…… y
h One must
then define
h: dependent
variable
‘Partial
derivatives’
x & y: of a function
independent
of more than
variables
one variable.
x
Or, the temperature on a flat surface that has various
heat sources spread out on that surface – such as a
tiny hot filament here, a hotter there, a tiny ice cube
here, a tiny beaker of liquid helium there, etc. 19
PCD-09
δx δx
h( x0 + , y0 ) − h( x0 − , y0 )
⎡ ∂h ⎤ 2 2 ⎡δ h ⎤
⎢ ∂x ⎥ = lim = lim ⎢ ⎥
⎣ ⎦ ( x0 , y0 ) δ x →0 δx δ x →0 δ x
⎣ ⎦ y0
y
h
h: dependent
variable Partial
x & y: derivatives of a
independent function of
variables more than one
variable.
x
δy δy
⎡ ∂h ⎤ h( x0 , y0 + ) − h( x0 , y0 − ) ⎡δ h ⎤
= lim 2 2 = lim ⎢ ⎥
⎢ ∂y ⎥ δ y →0 δy δ y →0 δ y
⎣ ⎦ ( x0 , y0 ) ⎣ ⎦ x0 20
PCD-09
Time-Derivatives of
δ t = t2 − t1 position
r &
r (t1 ) uur uur
δr = r δr velocity.
YI r r v = lim velocity
r (t2 ) − r (t1 ) δ t →0 δ t
OI XI r r
r r r (t2 ) − r (t1 )
ZI r (t2 ) v = lim
δ t →0 δt
‘Equation of Motion’
Rigorous relationship uur
between r δv acceleration
a = lim
position, δ t →0 δ t
r r
velocity r v(t2 ) − v(t1 )
and acceleration. a = lim
δ t →0 δt 21
PCD-09
I Law: Law of Inertia
What is ‘equilibrium’?
Relative to whom?
- frame of reference
Equilibrium means
‘state of rest’,
or of uniform
motion along a Galileo Galilee
straight line. 1564 - 1642
Isaac
Newton Linear Response.
(1642-1727) Principle of causality/determinism.
23
PCD-09
We have introduced the first two laws
of Newton as fundamental principles.
As is often done, we can introduce
Newton’s III law simply by stating
that
‘Action and Reaction are Equal and
Opposite’.
In Newtonian scheme of
mechanics, this is introduced
as a ‘fundamental’ principle
–i.e., as a law of nature.
24
PCD-09
Newton’s III law :
‘Action and Reaction are Equal and Opposite’.
ur ur
F 12 = − F 21
In Newtonian scheme of uur uur
d p1 d p2
mechanics, this is introduced =−
dt dt
as a ‘fundamental’ principle d ur ur
–i.e., as a law of nature. dt
( )
p1 + p 2 = 0
26
(1879 – 1955) (1882 – 1935) (1902 – 1995)
PCD-09
Consider an isolated N-particle System
Total internal force on k th particle
ur N ur
due to all other particles : F k = ∑ F ik
i =1, i ≠ k
Space : Homogeneous
Consider simulteneous displacement of the entire
uur
N-particle system through δ s
without affecting relative positions.
PASSIVE VIEW:
uur
Displacement of the entire SPACE through - δ s.
Virtual Displacement in homogeneous space
(translational invariance)
27
PCD-09
PASSIVE VIEW:
uur
Displacement of the entire SPACE through - δ s.
Virtual Displacement in homogeneous space
(translational invariance)
⎧ N ur ⎫ uur ⎧ N N ur ⎫ uur
0 = δ W = ⎨∑ F k ⎬ • δ s = ⎨∑ ∑ F ik ⎬ • δ s
⎩ k =1 ⎭ ⎩ k =1 i =1,i ≠ k ⎭
28
PCD-09
Conservation of
Space : Homogeneous. Internal forces linear momentum,
do *NO* work in this displacement! and thus Newton’s
III law, is a
⎧ N ur ⎫ uur ⎧ N N ur ⎫ uur property of
0 = δ W = ⎨∑ F k ⎬ • δ s = ⎨∑ ∑ F ik ⎬ • δ s. translational
⎩ k =1 ⎭ ⎩ k =1 i =1,i ≠ k ⎭
invariance.
For the above relation to hold good for arbitrary displacements,
we must have :
r ⎧ N ur ⎫ N dp
uur
N uur
ur ur
dP r
d dP
0 = ⎨∑ F k ⎬ = ∑ k = ∑ pk =
⎩ k =1 ⎭ k =1 dt dt k =1 dt
=0
ur dtuur
r dP
0=
uur
dt d p2 d p1
=−
Symmetry Conservation Law
urdt ur dt
F 12 = − F 21 29
PCD-09
Are the conservation principles consequences of the laws of
nature? Or, are the laws of nature the consequences of the
symmetry principles that govern them?
30
Emilly Noether: Symmetry Conservation Laws
31
The connection between SYMMETRY and
CONSERVATION PRINCIPLES brought out in
the previous example,
becomes even more transparent in an
alternative scheme of MECHANICS.
Common knowledge:
when a ray of light meets the edge of a medium,
it usually does not travel along the direction of
incidence
- gets reflected and refracted.
t1 t0 t2
t
38
PCD-09
t2
(q2,p2)
.
⎧ ∂L ∂L ⎫ ⎧ ∂L ∂L d ⎫
t2 t2
0 = ∫ ⎨ δ q + δ q& ⎬ dt = ∫ ⎨ δ q + (δ q )⎬ dt
t1 ⎩
∂q ∂q& ⎭ t1 ⎩
∂q ∂q& dt ⎭
42
differential and
integral of a
product of two
functions.
d ⎧ df ⎫ ⎧ dg ⎫
{ f ( x) g ( x)} = ⎨ ⎬ g ( x) + f ( x) ⎨ ⎬
dx ⎩ dx ⎭ ⎩ dx ⎭
∴
⎧ dg ⎫ d ⎧ df ⎫
f ( x) ⎨ ⎬ = { f ( x) g ( x)} - ⎨ ⎬ g ( x)
⎩ dx ⎭ dx ⎩ dx ⎭
43
PCD-09
⎧ dg ⎫ d ⎧ df ⎫
f ( x) ⎨ ⎬ = { f ( x) g ( x)} - ⎨ ⎬ g ( x)
⎩ dx ⎭ dx ⎩ dx ⎭
x1 ⎩ dx ⎭ x1 ⎩
dx ⎭
x2 x
⎧ dg ⎫ ⎧ df ⎫
2
∫
x1
f ( x) ⎨ ⎬ dx = f ( x2 ) g ( x2 ) − f ( x1 ) g ( x1 ) − ∫ ⎨ ⎬ g ( x) dx
⎩ dx ⎭ x1 ⎩
dx ⎭
44
PCD-09
⎧ ∂L ∂L ⎫ ⎧ ∂L ∂L d ⎫
t2 t2
0 = δ S = ∫ ⎨ δ q + δ q ⎬ dt = ∫ ⎨ δ q +
& (δ q )⎬ dt
t1 ⎩
∂q ∂q& ⎭ t1 ⎩
∂q ∂q& dt ⎭
⎧ ∂L ⎫
t2
⎧ ∂L d (δ q ) ⎫
t2
i.e. 0 = δ S = ∫ ⎨ δ q ⎬ dt + ∫ ⎨ ⎬ dt
t1 ⎩
∂q ⎭ t1 ⎩
∂q& dt ⎭
x1 ⎩ dx ⎭ x1 ⎩
dx ⎭
⎧ d ⎛ ∂L ⎞ ⎫
t2
⎧ ∂L ⎫ ⎡ ∂L ⎤
t2 t2
0 = ∫ ⎨ δ q ⎬ dt + ⎢ δ q ⎥ − ∫ ⎨ ⎜ ⎟ δ q ⎬ dt
t1 ⎩
∂q ⎭ ⎣ ∂q& ⎦ t1 t1 ⎩ dt ⎝ ∂q& ⎠ ⎭ 45
⎧ ∂L ∂L ⎫ ⎧ ∂L ∂L d ⎫
t2 t2
0 = δ S = ∫ ⎨ δ q + δ q ⎬ dt = ∫ ⎨ δ q +
& (δ q )⎬ dt
t1 ⎩
∂q ∂q& ⎭ t1 ⎩
∂q ∂q& dt ⎭
⎧ ∂L ⎫
t2
⎧ ∂L d (δ q ) ⎫
t2
i.e. 0 = δ S = ∫ ⎨ δ q ⎬ dt + ∫ ⎨ ⎬ dt
t1 ⎩
∂q ⎭ t1 ⎩
∂q& dt ⎭
x1 ⎩ dx ⎭ x1 ⎩
dx ⎭
⎧ d ⎛ ∂L ⎞ ⎫
t2
⎧ ∂L ⎫ ⎡ ∂L ⎤
t2 t2
0 = ∫ ⎨ δ q ⎬ dt + ⎢ δ q ⎥ − ∫ ⎨ ⎜ ⎟ δ q ⎬ dt
t1 ⎩
∂q ⎭ ⎣ ∂q& ⎦ t1 t1 ⎩ dt ⎝ ∂q& ⎠ ⎭ 46
δ q(t ) at time t p
@t
t1 < t < t2 (q1,p1)
@t
at time t1
(q2,p2)
at time t2
t 1 < t < t2
δ q ( t1 ) =δ q ( t2 ) = 0
47
PCD-09
⎧ d ⎛ ∂L ⎞ ⎫
t2
⎧ ∂L ⎫ ⎡ ∂L ⎤
t2 t2
0 = ∫ ⎨ δ q ⎬ dt + ⎢ δ q ⎥ − ∫ ⎨ ⎜ ⎟ δ q ⎬ dt
t1 ⎩
∂q ⎭ ⎣ ∂q& ⎦ t1 t1 ⎩ dt ⎝ ∂q& ⎠ ⎭
⎧ ∂L d ⎛ ∂L ⎞ ⎫
t2
⎡ ∂L ⎤
t2
i.e. 0 = ⎢ δ q⎥ + ∫ ⎨ − ⎜ ⎟ ⎬ δ qdt
⎣ ∂q& ⎦ t1 t1 ⎩ ∂q dt ⎝ ∂q& ⎠ ⎭
Arbitrary variation
Zero !
between the end points.
δ q ( t1 ) =δ q ( t2 ) = 0
∂L d ⎛ ∂L ⎞
Hence, − ⎜ &⎟=0
∂q dt ⎝ ∂q ⎠
Lagrange ' s Equation 48
PCD-09
∂L d ⎛ ∂L ⎞
− ⎜ & ⎟ = 0 Lagrange ' s Equation
∂q dt ⎝ ∂q ⎠
∂L
=−
∂V
= F , the force
L(q, q& , t ) = f1 (q& 2 ) + f 2 (q)
∂q ∂q
∂L m 2
= mq& = p, the momentum L(q, q& , t ) = q& − V (q)
∂q& 2
ur
dp dP r =T - V
i.e., = F : in 3D: =F
dt dt Interpretation of L as T-V
Newton ' s II Law gives equivalent
∂L correspondence with
p= , generalized momentum Newtonian formulation.
∂q& 49
PCD-09
When the Lagrangian is not an
L = L(q, q& , t ) ∂L
⎛ ∂L ⎞ ⎛ d ∂L ⎞ explicit function of time, =0
0 = ⎜ ⎟−⎜ ∂t
⎟
∂
⎝ ⎠ ⎝
q dt ∂q& ⎠ and
dL ⎛ ∂L ⎞ ⎛ ∂L ⎞ ∂L d ⎡ ⎧ ∂L ⎫ ⎤
= ⎜ ⎟q + ⎜ ⎟q +
& && ⎢ ⎨ ⎬ q& - L ⎥ = 0
dt ⎝ ∂q ⎠ ⎝ ∂q& ⎠ ∂t dt ⎣ ⎩ ∂q& ⎭ ⎦
dL ⎛ d ⎧ ∂L ⎫ ⎞ ⎛ ∂L ⎞ ∂L ⎡ ⎧ ∂L ⎫ ⎤
= ⎜ ⎨ ⎬⎟ q& + ⎜ & ⎟ q&& + ⎢ ⎨ ⎬ q& - L ⎥ is CONSTANT:
dt ⎝ dt ⎩ ∂q& ⎭ ⎠ ⎝ ∂q ⎠ ∂t ⎣ ⎩ ∂q& ⎭ ⎦
ENERGY
d ⎛ ⎧ ∂L ⎫ ⎞ ∂L
= ⎜ ⎨ ⎬ q& ⎟ + Hamiltonian's Principal Function
dt ⎝ ⎩ ∂q& ⎭ ⎠ ∂t ⎡ ⎧ ∂L ⎫ ⎤
H = ⎢ ⎨ ⎬ q& - L ⎥ = pq& - L
d ⎛ ⎧ ∂L ⎫ ⎞ ∂L ⎣ ⎩ ∂q& ⎭ ⎦
⎜ ⎨ ⎬ q& − L ⎟ = − When there are N degrees
dt ⎝ ⎩ ∂q& ⎭ ⎠ ∂t of freedom,
N
H = ∑ pi q&i − L
i =1 50
Time is homogeneous:
Lagrangian of a closed system
does not depend explicitly on time.
⎡ ∂L ⎤
⎢q& ∂q& - L ⎥ is a CONSTANT.
⎣ ⎦
Conservation of Energy is
thus connected with the
Hamiltonian
symmetry principle
“ENERGY”
regarding invariance with
Summation over i: respect to temporal
degrees of freedom translations.
N
⎛ ∂L ⎞
N
H = ∑ pi q&i − L = ∑ ⎜ ⎟ q&i − L
i =1 i =1 ⎝ ∂q
&i ⎠
PCD-09
Hamiltonian / Hamilton’s Principal Function 51
In an inertial frame,
Time is homogeneous; Space is homogenous and isotropic
∂L ∂ ⎡ m 2 ⎤
Now, p = = ⎢ q − V (q) ⎥ = mq&
&
∂q& ∂q& ⎣ 2 ⎦
& = q& m = 2T
qp 2
H = 2T - L = 2T - (T - V ) = T + V
Total Energy
53
PCD-09
Hamiltonian (Hamilton’s Principal Function) of a system
When there are many degrees of freedom, one takes the sum
∂L ∂L
dH = ∑ pk dq&k + ∑ q&k dpk −∑ dqk −∑ dq&k
k k k ∂qk k ∂q
&k
∂L
= ∑ q&k dpk −∑ dqk
k k ∂qk
= ∑ q&k dpk −∑ p& k dqk
k k
We have used:
∂L d ⎛ ∂L ⎞
− ⎜ & ⎟ = 0 Lagrange ' s Equation
∂q dt ⎝ ∂q ⎠
∂L d ⎛ ∂L ⎞ d
= ⎜ &⎟= p = p&
∂q dt ⎝ ∂q ⎠ dt
54
PCD-09
Hamiltonian (Hamilton’s Principal Function) of a system
When there are many degrees of freedom, one takes the sum
∂L ∂L
dH = ∑ pk dq&k + ∑ q&k dpk −∑ dqk −∑ dq&k
k k k ∂qk k ∂q
&k
∂L
= ∑ q&k dpk −∑ dqk
k k ∂qk
= ∑ q&k dpk −∑ p& k dqk
k k
But , H = H ( pk , qk )
∂H ∂H
Hamilton’s
so dH = ∑k ∂p k ∑k ∂q dqk
dp +
k k
equations of ∂H ∂H
Hence ∀ k : = q&k and = − p& k
motion ∂pk ∂qk
55
PCD-09
xi mu m ?
u m – M a ?
E xtre m im u m
Min
P o in t ?
Saddle ari ty… …
Station
cep ts l!
d c on 1 l e v e
a n ce P H 1 0
Adv e y on d
w e l l b
– 56
REFERENCES
• ‘Action’ : dimensions
‘angular momentum’ :
h: Max Planck :
fundamental quantity
in Quantum Mechanics
Illustrations: use of Lagrange’s / Hamilton’s equations
to solve simple problems in Mechanics
59
PCD-09
Manifestation of simple phenomena
in different unrelated situations
Dynamics of
spring–mass systems,
¾ pendulum,
¾ oscillatory electromagnetic circuits,
¾ bio rhythms,
¾ share market fluctuations …
1581:
Observations on the
swaying chandeliers
at the Pisa cathedral.
q: Generalized Coordinate
⎛ ∂L ⎞
p=⎜ ⎟
p: Generalized Momentum
⎝ ∂q& ⎠
L = L(q, q& , t )
H = H ( q, p, t )
62
PCD-09
L = L(q, q& , t ) Lagrangian
2nd order
∂L d ⎛ ∂L ⎞
− ⎜ &⎟=0 Lagrange ' s Equation differential
∂q dt ⎝ ∂q ⎠ equation
H = ∑ q&k pk − L
k
TWO
H = H ( q, p, t ) Hamiltonian
1st order
∂H ∂H differential
∀ k: = qk and
& = − p& k
∂pk ∂qk equations
Hamilton ' s Equations
63
PCD-09
Mass-Spring Simple Harmonic Oscillator
L = L(q, q& , t ) Lagrangian
m 2 k 2
L = T − V = q& − q
2 2
∂L d ⎛ ∂L ⎞
− ⎜ &⎟=0 Lagrange ' s Equation
∂q dt ⎝ ∂q ⎠
k d ⎛m ⎞
i.e. − 2q − ⎜ 2q& ⎟ = 0
2 dt ⎝ 2 ⎠
2nd order
− kq − mq&& = 0
Newton’s differential
mq&& = −kq
Lagranges’ equation
Equation of Motion for a simple harmonic oscillator
64
PCD-09
k
mq&& = −kq q&& = − q
m
Linear relation between restoring force and displacement
for spring-mass system: k
&&x =− x
m
Hooke’s law, after Robert Hooke (1635-1703), (a
contemporary of Newton), who empirically discovered this
relation for several elastic materials in 1678.
http://www-groups.dcs.st-and.ac.uk/~history/PictDisplay/Hooke.html
65
PCD-09
Mass-Spring Simple Harmonic Oscillator
Do we get same result using HAMILTONIAN approach ?
L = L(q, q& , t ) Lagrangian Note! Begin
H = H ( q, p, t ) Hamiltonian Always with
m 2 k 2 the
Lagrangian : L = T − V = q& − q LAGRANGIAN.
t u m 2 2
en
mom ⎛ ∂L ⎞
p = ⎜ ⎟ = mq&
⎝ ∂q& ⎠ L = L(q, q& , t )
m 2 k 2
H = pq& − L = mq& − q& + q H = H (q, p, t )
2
2 2
m 2 k 2
H = q& + q VERY
2 2
p2 k 2 IMPORTANT !
H= + q
2m 2
for a simple harmonic oscillator 66
PCD-09
m 2 k 2
H = pq& − L = mq& − q& + q
2
2 2
m k
H = q& 2 + q 2
2 2 ∂H 2 p p
p2 k 2 q& = = =
H= + q ∂p 2m m
2m 2
L = L(q, q& , t ) ∂H k
and p& = − = − 2q
H = H ( q, p, t ) ∂qk 2
67
PCD-09
L = L(q, q& , t )
H = H ( q, p, t ) ∂H 2 p p
q& = = =
∂p 2m m
VERY
IMPORTANT ! ∂H k
and p& = − = − 2q
∂qk 2
(i.e. f = −kq)
Hamilton ' s Equations
p& = − kq TWO first order equations
68
PCD-09
Be careful about how you write the Lagrangian and
the Hamiltonian for the Harmonic oscillator!
m 2 k 2 L = L(q, q& , t )
Lagrangian : L = T − V = q& − q
2 2
⎛ ∂L ⎞
p=⎜ ⎟ = mq&
⎝ ∂q& ⎠ H = H ( q, p, t )
p2 k 2
H= + q VERY
2m 2 IMPORTANT !
Generalized Momentum is interpreted
only as p = ⎛⎜ ∂L ⎞⎟ , and not a product of mass with velocity
⎝ ∂q& ⎠ 69
PCD-09
S: support
θ V = mgh = mg l(1 − cos θ )
l cos θ l : length
h = l − l cos θ
mg = l(1 − cos θ )
E: equilibrium
L = L(r , θ , r&,θ&)
Remember this! L = T −V
ALWAYS, the first 1 2 &2
thing to do is to set-up
L = m(r + r θ ) − mg l(1 − cos θ )
& 2
2
the Lagrangian in terms 1 2 &2
of the generalized L = ml θ − mg l(1 − cos θ )
coordinates and the 2
1 2 &2
generalized velocities. L = ml θ − mg l + mg l cos θ 70
PCD-09
2
1 2 &2
L = ml θ − mg l + mg l cos θ
2
r = l : fixed length
Subsequently, we can
find the generalized
momentum for each
degree of freedom. ∂L
= pr =0
∂r&
∂L 2 &
∂L d ⎛ ∂L ⎞ = pθ = ml θ
− ⎜ &⎟=0
∂q dt ⎝ ∂q ⎠ ∂θ&
71
PCD-09
S: support Application to a simple pendulum
l cos θ θ l : length
L = L(r ,θ , r&,θ&) = T − V
1 2 &2
L = ml θ − mg l + mg l cos θ
2
h = l − l cos θ
mg = l(1 − cos θ )
E: equilibrium
∂ L d ⎛ ∂L ⎞
∂L ∂L − ⎜ &⎟=0
= pr =0 = 0. ∂q dt ⎝ ∂q ⎠
∂r& ∂r
∂L
= pθ = ml 2θ& ∂L = − mgl sin θ ≈ − mglθ d
∂θ&
∂θ − mglθ − (ml 2θ&) = 0
dt
(1) q&& = −α q ml 2θ&& = − mglθ
g
(2) Solution: q = Aeiω0t + Be − iω0t ω0 = g
Substitute (2) in (1) ⇒ ω0 = α l &&
θ =− θ
l 72
PCD-09 Does the r equation give you anything useful?
For the simple pendulum oscillating in the
gravitational field where the acceleration
due to gravity is g, we must, and do, get
the same answer regardless of which
approach we employ:
(1) Newtonian g
θ&& = − θ
(2) Lagrangian l
(3) Hamiltonian
Note! We haven’t used ‘force’,
g
ω0 = ‘tension in the string’ etc. in the
Lagrangian and Hamiltonian
l approach!
73
PCD-09
Hamiltonian Approach Application to a simple pendulum
2 &
⎡ ∂L ⎤ p = ml θ
H = ⎢ ∑ q&i - L ⎥ = ⎡⎣ ∑ q&i pi - L ⎤⎦ θ
2 &&
⎣ ∂q&i ⎦ p&θ = ml θ
∂H H = θ pθ − ml θ + mg l − mg l cos θ
& 1 2 &2
p&θ = − 2
∂θ ∂H = −mgl ⎡ ∂ (cos θ )⎤ ∂(cos θ )
∂θ ⎢ ∂ (cos θ ) ⎥ ∂θ
⎣ ⎦
= − mgl (− sin θ ) ≈ + mglθ
∂H
≈ mglθ
∂θ
∂H
&pθ = ml 2θ&& = − = − mglθ g
∂θ
g ω 0 =
θ&& = − θ l
l
74
PCD-09
Lagrangian and Hamiltonian Mechanics has
very many applications. All problems in
‘classical mechanics’ can be addressed using
these techniques.
75
PCD-09
However, the accuracy of the experimental
determination of (q,p) is limited by the
unavoidable coupling between the target and
the probe (Heisenberg’s principle of
uncertainty).
Flip-Flop Square
Wave oscillator
1
PCD-09
Learning Goals
unstable
stable
unstable
stable neutral
stable
neutral
a c e d b
4
PCD-09
Is a point mid-way between two
equal positive point charges a
point of equilibrium for a unit
point positive test charge?
+ +
Can it be unambiguously
classified as a point or ‘stable’ /
‘unstable’ equilibrium?
5
PCD-09
meaning of small oscillations U(x) x
‘Zero’, at equilibrium
∂U
U ( x) = U ( x0 ) + ( x − x0 ) +
∂x x0
1 ∂ 2U U(x)
+ ( x − x0 ) +
2
2! ∂x 2 x0
1 ∂ 3U
+ ( x − x0 )3 + ...
3! ∂x 3 x0
x
Potential for a
Approximations, close to x0 Linear harmonic oscillator
1 ∂ 2U dU
U ( x) ≈ U ( x0 )+
1
( x − x0 ) = kx 2
2 F =− = − kx
2! ∂x 2 2 dx
x0
by choosing U ( x0 ) = 0 k
x=− x
&&
m
6
PCD-09
Uniform circular motion
and SHM
Intrinsic www.answers.com
natural frequency
ω = 2πν
Shadow of the red dot
www.physics.uoguelph.ca
would constitute SHM 7
PCD-09
Unlike what happens in
a resistor, the current
and voltage in an
inductance L and in a
capacitor C do not peak
dV
I is proportional to , together.
dt
not to V, as in the case of a resistor.
VL Voltage lags the current in a
capacitor by 900,
I
but
VC
leads the current in an inductor by
8
PCD-09
the same amount.
Qmac
Vmax =
C
dI d 2Q &&
VL = − L = − L 2 = − LQ
dt dt −VL + VC = 0
& d d dV d Q Q
2
I = Q = Q = (CV ) = C +L 2 + = 0
dt dt dt
dt C
&& 1
I is proportional to
dV
,
Q = −( )Q
dt LC
not to V, as in the case of a resistor.
VL
I Voltage lags the current in a capacitor by 900,
VC but leads the current in an inductor by the same amount.
9
PCD-09
Electro-mechanical analogues
Qmin
10
PCD-09
k
x=− x
&&
m Electro-mechanical analogues:
&& 1
Q = −( )Q Inductance mass, inertia
LC Capacitance 1/k, compliance
(1) q&& = −α q
iω0t − iω0t
(2) Most general solution: q = Ae + Be
Substitute (2) in (1) ⇒ ω0 = α
k Question:
ω0 =
m Could we have
associated L with
1/k and C with m?
1
ω0 =
LC 11
PCD-09
g
θ&& = − θ (1) sim ple pendulum
l
k
x=− x
&& (2) spring-m ass system
m
Q&& = − ( 1 ) Q electrical LC circuit oscillator
(3)
LC
(1) q&& = − α q
(2) M ost general solution: q = Ae iω 0 t + Be − iω 0 t
(1) q&& = −α q
(2) Most general solution: q = Ae iω0 t + Be − iω 0 t
13
PCD-09
The most general solution is
q = Ae iω0 t +B e − iω0 t where ω 0 = α
the frequency is governed by α ;
A and B are determined by initial conditions
14
PCD-09
M ean kinetic energy
k
x=− x
&& spring-mass system KE =
1
mx& 2
m 2
x = A sin(ω0t + ϕ )
k M ean potential energy
ω0 =
T
m ∫
t =0
f ( t ) dt
PE =
1 2
kx
M ean: f (t ) = T 2
∫
t =0
dt
1
<cos 2 (ω t)dt> = = <sin 2 (ω t)dt>
2
1 1 1 1
KE = m x& 2 = m ( Aω 0 ) 2 < >= A 2 mω 0 2
2 2 2 4
1 2 1 1 1
PE = kx = kA 2 < >= kA 2
2 2 2 4
k
Note that < PE >=< KE > , since ω 0 2
= 15
PCD-09 m
&& 1
Q = −( )Q
LC
(1) q&& = −α q
(2) Most general solution: q = Aeiω0t + Be − iω0t
Substitute (2) in (1) ⇒ ω0 = α
16
PCD-09
Longitudinal oscillations
1 2 Longitudinal Oscillations
Frictionless support
PCD-09
Reference: Berkeley’s Mechanics are the same? 17
ao: relaxed length of the springs Tension exerted by each string
a: instantaneous stretched length AT EQUILIBRIUM
T0 = k (a − a0 )
T0=k(a-a0)
Transverse oscillations
in
‘vertical’
T=k(l-a0)
l and in
x
θ ‘horizontal’
a a planes
View in the plane of vibration
The total restoring force along − x
is − 2T sin θ
x
1 2 mx&& = −2T sin θ = −2k (l − a0 )
l
18
PCD-09
ao: relaxed length of the springs
a: instantaneous stretched length
T0=k(a-a0)
phase: θ = ω t ± kx ± Δ
On a surface of constant phase: dθ =0 i.e. ω dt ± kdx = 0
± kdx=-ω dt i.e. kdx= m ω dt
dx ω
when =m
dt k
dx
<0
dt
← a wave travelling to the left
when
dx
>0
dt
→ a wave travelling to the right
21
PCD-09
q (t ) = q0 cos {ωt ± (kx + Δ )} = q 0 cos{ω t ± kx ± Δ}
phase: θ = ω t ± kx ± Δ
when phase: θ = ω t + kx ± Δ,
on a surface of constant phase,
dx
<0
dt
← a wave travelling to the left
when phase: θ = ω t − kx ± Δ,
on a surface of constant phase,
dx
>0
dt
→ a wave travelling to the right
22
PCD-09
f(x-vt) represents a pulse traveling to the right
dx dx
> 0, i.e. as a positive quantity
dt dt
→ a wave travelling to the right
24
The total energy E is constant: conservative forces
k c
where ω0 = 2
γ=
m 2m
If EM & Gravitational forces are conservative,
and all forces are made up of fundamental forces,
Then,
why is friction dissipative?
Just what is ‘lost’, and why?
26
PCD-09
All net interactions in nature are superpositions
of fundamental interactions, namely the nuclear
(‘strong’ interaction), electro-weak
(electromagnetic/nuclear ‘weak’ interaction),
and gravity.
The electrodynamics/electroweak/gravitational
interactions are conservative.
Eq.[#]
k c
where ω0 = γ=
2
m 2m
We seek a solution in the form: x = Ae qt Eq.[*]
Why seek this form? and inquire what conditions would result on q if
Eq.[*] is to be admitted as a solution of Eq.[#]
q2 = − γ − γ 2 − ωο 2
CASE 1:
OVERDAMPED OSCILLATOR
When γ > ω0 ,
γ 2 − ωο 2 is a real number whose value is < γ ,
so both q1 and q 2 become 'real' and essentially 'negative'
31
PCD-09
x(t ) = A1e q1 t + A2 e q2 t k
q1,2 = − γ ± γ 2 − ωο 2 , ω0 =2
m
When γ > ω0 ,
c
γ − ωο
2 2
is a real number whose value is < γ , γ=
so both q1 and q 2 become 'real' and essentially 'negative'
2m
x(t ) = A1e q1 t
+ A2 e q2 t
hence
&x(t ) = q1 A1e q1 t + q2 A2 e q2 t q1 x(t = 0) − x& (t = 0) = A2
q1 − q2
Hence, q1 x(t = 0) − x& (t = 0)
x(t = 0) − = A1
x(t = 0) = A1 + A2 q1 − q2
x& (t = 0) = q1 A1 + q2 A2 ‘Overshoot’ is not possible.
q1 x(t = 0) = q1 A1 + q1 A2 Oscillations being completely
killed, this oscillator is called
x& (t = 0) = q1 A1 + q2 A2 ‘OVERDAMPED’. 32
PCD-09
k
ω0 =2
CASE 2
m
UNDERDAMPED OSCILLATOR
When γ < ω0 ,
γ=
c γ − ωο is an imaginary number
2 2
2m
q1 = − γ + i ω0 2 − γ 2 = −γ + iω where ω = ω0 2 − γ 2
q2 = − γ - i ω0 2 − γ 2 = −γ − iω i.e., ω < ω0
by an amount determined by γ
[ − γ −iω ] t
[ − γ + iω ]t
x = A1e + A2 e
x=e −γ t
{A e1
+ iω t
+ A2 e
− iωt
}
x = e −γ t {( A1 + A2 ) cos(ωt ) + i ( A1 − A2 ) sin(ωt ) }
33
PCD-09
UNDERDAMPED OSCILLATOR k
When γ < ω0 , γ 2 − ωο 2 is an imaginary number
ω0 =
2
m
c
γ=
[ − γ −iω ] t
[ − γ + iω ]t
x = A1e + A2 e
2m
x=e −γ t
{A e1
+ iω t
+ A2 e
− iωt
}
x = e −γ t {( A1 + A2 ) cos(ωt ) + i ( A1 − A2 ) sin(ωt ) }
34
PCD-09
CASE 2 UNDERDAMPED OSCILLATOR k
ω0 =
2
} γ=
2m
x = e −γ t {( A1 + A2 ) cos(ωt ) + i ( A1 − A2 ) sin(ωt ) }
Introduce two new parameters B & θ instead of A1 and A2.
A1 + A2 = B sin θ
i ( A1 − A2 ) = B cos θ iBe + iθ iBe − iθ
A1 = − , A2 = +
2 2
x(t ) = Be −γ t
{sin θ cos(ωt ) + cos θ sin(ωt ) }
x(t ) = Be −γ t sin(ωt + θ )
35
PCD-09
UNDERDAMPED OSCILLATOR
γ < ω0 , γ 2 − ωο 2 :imaginary ω 2 = k
0
m
x(t ) = Be −γ t {sin θ cos(ωt ) + cos θ sin(ωt ) } c
γ=
x(t ) = Be −γ t sin(ωt + θ ) 2m
• Frequency ω < ω0
ω = ω0 2 − γ 2
i.e., ω < ω0 by an amount determined by γ 36
PCD-09
k
ω0 =2
UNDERDAMPED OSCILLATOR m
c
γ=
2m
When γ < ω0 , γ 2 − ωο 2 is an imaginary number
x(t ) = Be −γ t
sin(ωt + θ )
x(t ) = Be −γ t
sin(ωt + θ )
Unlike the ‘overdamped oscillator’ (no oscillations),
we do have oscillations that are ‘damped’, not
‘killed’; hence called UNDERDAMPED
OSCILLATIONS
40
PCD-09
Case 3: ‘CRITICAL DAMPING’ q1,2 = − γ ± γ 2 − ωο 2
Overdamped Oscillator
When γ > ω0 , i.e. c 2 > 4mk ,
γ 2 − ωο 2 is a real number whose value is < γ ,
so both q1 and q 2 become 'real' and essentially 'negative'
42
PCD-09
Underdamped Oscillator x = A1e q1 t
+ A2 e q2 t
Overdamped: oscillations
are `killed’
Underdamped:
oscillations are damped
critically damped:
overshoots equilibrium
in finite time
45
PCD-09
Forced oscillations
Restoring force, damping force and driving force
46
PCD-09
Forced oscillations
Restoring force, damping force and driving force
F = mx&& = − kx − cx& + Fdr
i.e.,
c k Fdr
x + x& + x =
&&
m m m
For a simple pendulum with damping ,
&& c & g Fdr
θ+ θ+ θ=
ml l ml
&& R & 1 Vdr
for an LCR oscillator, Q + Q + Q= or
L LC L
&& & 1
LQ + RQ + Q = Vdr
C 47
PCD-09
c k Fdr
F = mx&& = − kx − cx& + Fdr or &&
x + x& + x =
m m m
Actual form of the solution depends on the functional form of Fdr
iθ iΩt
) iΩt
Fdr = Fο e e =Fe
Substituting for x& and &&
x: )
) where F = Fο eiθ
[−Ω + i(2γΩ) + ω ] x(t ) = (F / m)eiΩt
2 2
0
) iΩt )
i.e.[(ω0 −Ω ) + i2γΩ ] A e
2 2
= (F / m)eiΩt
50
PCD-09
General case, including damping: iθ iΩt
) iΩt
Fdr = Fο e e =Fe
F0 )
x + 2γ x& + ω02 x =
&& exp i (Ωt + θ ) where F = Fο eiθ
m
) iΩt
x(t ) = Ae
)
with A = A0 ei (θ −φ )
) iΩt )
[(ω − Ω ) + i 2γΩ ] A e
2
0
2
= ( F / m)eiΩt
)
)
A=
F /m { ;
}
{
(ω0 − Ω ) + i 2γΩ
2 2
}
i.e. A0 e i (θ −φ )
=
{ }
F0 eiθ / m
{( 0 ) + i 2γΩ
ω 2
− Ω 2
}
)
as F = F0 eiθ A0 = A0 (Ω)
51
PCD-09
cancel
A0 e i (θ −φ )
=
{ F0 eiθ / m } ; A0 e −iφ
=
{F0 / m}
(ω − Ω ) + i 2γΩ
2
0
2
( 0 ) + i 2γΩ
ω 2
− Ω 2
e −iφ
=
{F0 /(mA0 )}
( 0 ) + i 2γΩ
ω 2
− Ω 2
{0
F ( mA0 }( 0
) ω − Ω ) {F / ( mA )} 2γΩ
2 2
/
cos φ = and sin φ =
0 0
(ω − Ω
2
0
2 2
) + 4γ 2 Ω 2 ( 0 −Ω
ω 2
)
2 2
+ 4γ 2 Ω 2
52
PCD-09
{0
F ( mA0 }( 0
) ω − Ω ) {F / ( mA )} 2γΩ
2 2
/
cos φ = and sin φ =
0 0
(ω − Ω
2
0 )
2 2
+ 4γ 2 Ω 2 ( 0 −Ω
ω 2
)
2 2
+ 4γ 2 Ω 2
⎧ 2γΩ ⎫ 2γΩ
and φ = tan ⎨ 2 2⎬
; tan φ = 2
−1
ω
⎩ 0 − Ω ⎭ ω 0 − Ω 2
Fascinating
Nature of the solution depends on γ applications in
and on the proximity of Ω to ω0 . mechanical, electrical
and many other
physical systems. 54
PCD-09
) iΩt
x(t ) = Ae ( F0 / m)
x(t ) = ei (Ωt +θ −φ ) .
) i (θ −φ ) (ω02 − Ω 2 ) 2 + 4γ 2 Ω 2
with A = A0 e
F0
Aο (Ω) =
m (ω − Ω
2
0 )
2 2
+ 4γ 2 Ω 2
A0 = A0 (Ω)
Two frequencies
are of interest ω0 Intrinsic, natural frequency.
Ω External, under our control!
56
Reference: Fowles ‘Analytical Mechanics’; Our notation is slightly different!
PCD-09
Condition for Resonance F0
Aο (Ω) =
when is
dA0
= 0? m (ω − Ω
2
0 )
2 2
+ 4γ 2 Ω 2
dΩ
ω0 Intrinsic, natural frequency.
Ω External, under our control!
1 F0
dA0
− { 2(ω 2
0 − Ω 2
)(-2Ω )+ 8 γ 2
Ω}
when is = 2m = 0?
dΩ {(ω02 − Ω2 )2 + 4γ 2Ω2 }
3/ 2
F0
Aο (Ω) MAXIMUM = m
( ω − ( ω − 2γ )) + 4γ 2 (ω02 − 2γ 2 )
2
2 2 2
0 0
F0 i.e.
A0 (Ω) MAXIMUM = m F0
= 2γ A0 (Ω) MAXIMUM ω0 2 − γ 2
2γ ω0 − γ
2 2 m
59
PCD-09
F0
Using: = 2γ A0 (Ω) MAXIMUM ω0 2 − γ 2
m
When γ << ω0 : Ω r ≈ ω0
A0 (Ω)
Less Damping
Ωr More Damping
Ω 61
PCD-09
Frequency of the Driving Force
F0 i (Ωt +θ )
x + 2γ x& + ω x = e
Thus the solution for && 2
0 becomes
m
( F0 / m)
x(t ) = ei (Ωt +θ −φ ) .
(ω02 − Ω 2 ) 2 + 4γ 2 Ω 2
g 1
ω0 = and ωο = , for LC -circuit
l LC
for simple pendulum
64
PCD-09
A0 (Ω) MAXIMUM γ
Aο (Ω) ≈
Energy is proportional to
( ω0 − Ω ) + γ 2
2
2γ “RESONANCE WIDTH”
ω ω0
Define: Q= ≈ (for the case of weak damping)
2γ 2γ
Quality Factor 65
PCD-09
The Tacoma Narrows Bridge in Washington state, was
with 1.9 km length one of the largest suspended bridges
built at the time. The bridge connecting the Tacoma
Narrows channel collapsed in a dramatic way on Thursday
November 7, 1940. Winds at about 50-70 km/hr produced
an oscillation which eventually broke the construction.
Forced/Driven
Damped
Oscillator
http://www.math.harvard.edu/archive/21b_fall_03/tacoma/index.html 66
PCD-09
Google: MRI picture
67
PCD-09
Solutions of the
oscillator problem play
a fundamental, crucial
role in DSP, information
transmission, etc.
68
A wave packet,
or a wave Wave / Pulse propagation
pulse, is made
up by
superposing a
y = f ( x)
y = f ( x − vt )
large number
of sinusoidal Unit 2:
waves.
Oscillators
Each Resonances
component Waves
wave has a
wavelength X
λ Wave pulse at t=0
vt Wave pulse at t > 0
and the Peak is at x=0
component Peak is at x=vt
travels at its
own phase Shape: same
velocity Medium: Non-dispersive
v = νλ .
69
PCD-09
Pulse shapes ---- Fourier Analysis
Fourier :
70
PCD-09
⎡ ⎛x⎞ ⎛ x ⎞⎤
Plot the function: f ( x) = 2 ⎢ H ⎜ ⎟ − H ⎜ − 1⎟ ⎥ − 1
⎣ ⎝L⎠ ⎝ L ⎠⎦
when x ∈ [ 0, 2 L ]
H (x) =Unit
0 for
2: x<0 Heaviside step function
=Oscillators
1 for x>0 “Unit step function”
Resonances
Waves
2L
x
x 0
X=0 L
-1
⎡ ⎛x⎞ ⎛ x ⎞⎤
Square Wave: f ( x) = 2 ⎢ H ⎜ ⎟ − H ⎜ − 1⎟ ⎥ − 1
⎣ ⎝L⎠ ⎝ L ⎠ ⎦ 71
PCD-09
∞
4 1
f (θ ) = ∑
π n =1,3,5,.... n
sin ( nθ )
72
8 ⎡ 1 1 1 1 ⎤
f (θ ) = sin θ − sin(3θ ) + sin(5θ ) − sin(7θ ) + sin(9θ ) − ......
π 2 ⎢⎣ 32 52 72 92 ⎥⎦
Applications:
DSP
ψ ( x, t ) = f ( x − vt )
73
We worked with the function f = f ( x)
⎡ ⎛x⎞ ⎛ x ⎞⎤
Square Wave: f ( x) = 2 ⎢ H ⎜ ⎟ − H ⎜ − 1⎟ ⎥ − 1
⎣ ⎝L⎠ ⎝ L ⎠⎦
ω is a function of k , given as ω ( k ), vφ = vφ ( k ),
the functional form is different for different systems
Actually, it is the MEDIUM that is non-dispersive.
Then, we get
ψ ( z , t ) = Amod ( z , t ) cos(ωavet − k ave z )
where Amod ( z , t ) = 2 A cos(ωmod t − kmod z )
1 1
ωmod = (ω1 − ω2 ) ; kmod = ( k1 − k2 )
2 2
Ref: Berkeley/
1 1
Vol.3/ also, ωave = (ω1 + ω2 ) ; kave = ( k1 + k2 )
Page270 2 2 78
PCD-09
ψ ( z , t ) = Amod ( z , t ) cos(ωavet − kave z ) At what speed does
where Amod ( z , t ) = 2 A cos(ωmod t − kmod z ) the modulation
propagate?
1 1
ωmod = (ω1 − ω2 ) ; kmod = ( k1 − k2 )
2 2
1 1
also, ωave = (ω1 + ω2 ) ; kave = ( k1 + k2 )
2 2
To follow a given modulation wave crest of the
modulation amplitude Amod ( z , t ),
we need to maintain a constant value of (ωmod t − k mod z )
i.e., in time dt , z must increase by dz in such a way that
d (ωmod t − k mod z ) = (ωmod dt − k mod dz ) = 0
To satisfy this, the modulation must propogate at:
dz ωmod ω1 − ω2 δω d ω
= v mod = = = ≈ = v g = 'group velocity'
PCD-09
dt k mod k1 − k 2 δ k dk 79
If all the components of the wave-packet
travel at the same speed, the ‘shape’ of
the wave-packet propagates without
distortion.
80
x
A Refraction
θ1 a Why does the light
‘1’ B
B’ ray go along the path
A B C,
‘2’ θ2 b and not along
C A B’ C
d
Time taken for light
( ) ( )
1/ 2 1/ 2
a 2
+x 2
b 2
+ (d − x) 2
to travel the path
t= + A B C:
v1 v2
dt 1 1 2
( ) ( )
−1/ 2 1 1 2 −1/ 2
= a + x2 (2 x) − b + (d − x) 2 (2(d − x))
dx v1 2 v2 2
dt 1 x 1 (d − x) sin θ1 sin θ 2
= − 0= −
(
dx v1 a 2 + x 2 ) ( )
1/ 2 1/ 2
v 2 b 2 + (d − x ) 2 v1 v2
sin θ1 v1
Zero (Fermat’s principle) = = nRefractive Index
sin θ 2 v 2 81
Refractive index, n:
Ratio of phase velocity of light in
vacuum to that in the medium
c νλvac λvac
n= = =
vϕ νλϕ λϕ
Thuban
1
PCD-09
Unit 3
Plane polar and cylindrical
polar coordinate systems.
Motivation: Description of
planetary motion in solar
system, and other dynamics
wherein the central interaction
has a center of symmetry, or
at least an axis of symmetry. 2
PCD-09
Unit 3: Learning goals
Develop confidence in vector methods and coordinate
transformations and learn to use an appropriate
coordinate system to simplify analysis.
Vector methods are fundamental to physics and not
merely a matter of using powerful mathematical tools.
Physical quantities are tensors of various ranks and
depend on how their components transform under the
rotation of a coordinate frame of reference.
We shall get well acquainted with the algebra of the
vectors, which are tensors of rank 1, and learn to
distinguish between ‘polar’ (such as position, velocity,
acceleration, force, etc.) and ‘axial’ (such as magnetic
field, angular momentum, angular velocity, etc.)
vectors. 3
PCD-09
Originally published in
1884, this wonderful
fantasy--written three
decades before Einstein's
general theory of relativity
and a century before string
theory--explores the world
of many dimensions by
imagining a world of only
two.
EDWIN A. ABBOTT
When its protagonist--
narrator, A Square, is
visited by a sphere, his
mind is stretched, and
ultimately his placid world is
turned inside out.
4
http://findarticles.com/p/articles/mi_m1134/is_6_117/ai_n27907784
PCD-09
Ptolemy’s The sun and the planets
coordinate were considered to move on
system a small circle (called
accounted for ‘epicycle’) whose center
the retrograde would move on a large circle
apparent (called ‘deferent’).
motion of all
planets known
at his time.
http://csep10.phys.utk.edu/astr161/lect/retrograde/aristotle.html
5
PCD-09
Claudius Ptolemaeus (AD100-170) (called Ptolemy).
June:
Galileo is sentenced to prison for an indefinite term. Seven of
ten cardinals presiding at his trial sign the sentencing order. Galileo signs a
formal recantation. Galileo is allowed to serve his term under house-arrest in
the home of the archbishop of Siena.
11
PCD-09
z
Definition of a vector: “magnitude” and “direction”
ur ur
A • uˆ = A cos θ
ur THE
LEANING
A TOWER OF
PISA, ITALY
PCD-08
uˆ 14
scalar or inner or dot product
r r r r
A • B = AB cos θ = B • A (Commutative)
r r r r r r r
A • ( B + C ) = A • B + A • C ( Distributive)
r r r r r r
m( A • B ) = (mA) • B = A • (mB ) = mAB cos θ ( Associtive)
¾ Invariant under ROTATION of basis, i.e. coordinate
rotation
¾ ……. This is how a ‘scalar’ *MUST* be defined !
eˆi • eˆ j = δ ij r r
A • B = Aiδ ij B j
= 1, when i = j
= 0, when i ≠ j = Ai Bi
“God created the integers, all else is the work of man.”
– Leopold Kroneker (1823-1891) 16
PCD-09 http://www-groups.dcs.st-and.ac.uk/~history/Biographies/Kronecker.html
vector or cross product
r r r r
A × B = AB sin θ nˆ = − B × A ( Not Commutative)
r
n̂ B r θ
B sin
directed area
θ A
r r r r r r r
A × ( B + C ) = A × B + A × C ( Distributive)
( ) ( )
A × B × C ≠ A × B × C ( Not Associtive)
( ) ( ) (
A × B × C = A • C B − A • BC = A • C B − A • B C )
( ) ( ) (
A × B × C = A • C B − A • BC = B A • C − C A • B )
(A × B)× C = (C • A)B − (C • B )A 17
PCD-09
Y
How do vectors
Y/
transform under
rotation of a coordinate VY
VY/ X/
system?
r VX/
V = V x eˆ x + V y eˆ y X
VX
ê z
eˆ y '
θ
ê y
θ
ê x
eˆ x '
Positive rotation about z axis
θ
ê x
eˆ x ' ⎡ eˆ ' ⎤ ⎡ cos θ sin θ 0 ⎤ ⎡ eˆx ⎤
⎢ x⎥ ⎢ ⎥ ⎢ ⎥
⎢ eˆy' ⎥ = ⎢ − sin θ cos θ 0 ⎥ ⎢eˆy ⎥
⎢ ⎥ ⎢ 0 1 ⎥⎦ ⎢⎣ eˆz ⎥⎦
⎢⎣ eˆz ' ⎥⎦ ⎣ 0
Vz/ =Vx ⎡⎣eˆz/ ⋅ eˆx ⎤⎦ +Vy ⎡⎣eˆz/ ⋅ eˆy ⎤⎦ +Vz ⎡⎣eˆz/ ⋅ eˆz ⎤⎦
Father of experimental
Physics
Grandpa of Engineering?
22
PCD-09
Why is an image in a mirror
inverted left to right
but not top to bottom?
23
PCD-09
Discrete Symmetries OBSERVE
Parity REFLECTION
Charge Conjugation ⎛ − x ⎞ ⎛ x ⎞
(matter – antimatter) ⎜ ⎟ ⎜ ⎟
Time Reversal ⎜ y ⎟ ← y
⎜ ⎟
⎜ z ⎟ ⎜ z ⎟
⎝ ⎠ ⎝ ⎠
z M
I
x’ R
R
x O
y
R
T
⎡ ⎛ R11 R12 R13 ⎞ ⎛ x1 ⎞⎤ ⎛ R11 R12 R13 ⎞ ⎛ y1 ⎞
⎢⎜ ⎟ ⎜ ⎟⎥ ⎜ ⎟ ⎜ ⎟
(R X ) T (R Y ) = ⎢ ⎜ R 21 R 23 R 23 ⎟ ⎜ x2 ⎟⎥ ⎜ R 2 1 R 2 3 R 2 3 ⎟ ⎜ y2 ⎟
⎜ ⎟ ⎜ ⎟⎥ ⎜ R ⎟ ⎜ ⎟
⎣⎢ ⎝ R 3 1 R 32 R 33 ⎠ ⎝ x3 ⎠⎦ ⎝ 31 R 32 R 33 ⎠ ⎝ y3 ⎠
T T
⎛ x1 ⎞ ⎛ R11 R12 R13 ⎞ ⎛ R11 R12 R13 ⎞ ⎛ y1 ⎞
⎜ ⎟ ⎜ ⎟ ⎜ ⎟⎜ ⎟
= ⎜ x2 ⎟ ⎜ R 21 R 23 R 23 ⎟ ⎜ R 2 1 R 2 3 R 2 3 ⎟⎜ y2 ⎟
⎜ x ⎟ ⎜ R ⎟ ⎜ R R 3 3 ⎟⎠ ⎜⎝ y 3 ⎟⎠
⎝ 3 ⎠ ⎝ 31 R 32 R 33 ⎠ ⎝ 31 R 32
T
⎛ R11 R12 R13 ⎞ ⎛ R11 R12 R13 ⎞ ⎛ y1 ⎞
⎜ ⎟ ⎜ ⎟ ⎜ ⎟
= ( x1 x2 x 3 )⎜ R 21 R 23 R 23 ⎟ ⎜ R 21 R 23 R 23 ⎟ ⎜ y 2 ⎟
⎜ R 3 3 ⎠⎟ ⎜⎝ R 3 1 ⎟ ⎜ y ⎟
⎝ R 31 R 32 R 32 R 33 ⎠ ⎝ 3 ⎠
T
⎛ R11 R12 R13 ⎞ ⎛ R11 R12 R13 ⎞
⎜ ⎟ ⎜ ⎟ Orthogonal Matrix
⎜ R21 R23 R23 ⎟ ⎜ R21 R23 R23 ⎟ = 1 Inverse=Transpose
⎜R R33 ⎟⎠ ⎜R R33 ⎟⎠
⎝ 31 R32 ⎝ 31 R32 27
Matrix representation of the transformation operator
For a rotation
TRANSFORMATION ⎛ −1 0 0 ⎞ For
⎜ ⎟ inversion,
(SQUARE) MATRIX RT R = RRT = 1 ⎜ 0 −1 0 ⎟ det R = −1
⎜ 0 0 −1⎟
and det R = 1 ⎝ ⎠ 28
Rotation & Reflection
r’= R r are fundamentally different
symmetries!
For rotation, For parity,
det |R| = +1 det |R| = −1
29
ORTHOGONAL
MATRICES
⎛ R11 R1 2 R1 3 ⎞T ⎛ R1 1 R1 2 R1 3 ⎞
⎜ ⎟ ⎜ ⎟
⎜ R 21 R 23 R 23 ⎟ ⎜ R 21 R 23 R 2 3 ⎟ = I 3× 3
⎜R R 3 3 ⎟⎠ ⎜⎝ R 3 1 R 3 3 ⎟⎠
⎝ 31 R 32 R 32
D e te r m in a n t o f th e m a tr ix is ± 1
R1 1 R1 2 R1 3
R 21 R 23 R 23 = ± 1
R 31 R 32 R 33
+ R O T A T IO N S
- P A R IT Y
30
⎛ x1 ⎞ ⎛ y1 ⎞ Set of all rotations
r ⎜ ⎟ r ⎜ ⎟
= ⎜ x2 ⎟ , y =
x
⎜ x ⎟
⎜ y
⎜ y ⎟
2 ⎟, forms a group,
⎝ 3 ⎠ ⎝ 3 ⎠
r
x T
= ( x1 x2 x 3 ), y
r T
= ( y1 y2 y3 ) O(3)
r r r r Invariance of scalar O: “orthogonal.”
Rx • Ry = x • y
(R X )T (R Y ) = X T
Y under rotations The subset with
T
⎡ ⎛ R11
⎢⎜
R12 R13 ⎞
⎟
⎛ x1 ⎞ ⎤
⎜ ⎟⎥
⎛ R11
⎜
R12 R13 ⎞
⎟
⎛ y1 ⎞
⎜ ⎟
determinant +1 is
• ⎜ R 21
⎢ ⎜ R 21
⎜
R 23 R 23 ⎟ ⎜ x 2 ⎟⎥ R 23 R 23 ⎟ ⎜ y 2
⎜ y ⎟
⎟ also a group:
⎣⎢ ⎝ R 3 1 R 32 R 3 3 ⎠⎟ ⎜ x ⎟⎥
⎝ 3 ⎠⎦
⎜
⎝ R 31 R 32 R 3 3 ⎠⎟ ⎝ 3 ⎠
SO(3)
⎛ y1 ⎞
⎜ ⎟ This part: S: “special.”
X T Y = ( x1 x2 x 3 )⎜ y 2 ⎟
⎜ y ⎟
⎝ 3 ⎠
OPTIONAL. It is usual to refer
⎛ R11 R12 R13 ⎞T ⎛ R11 R12 R13 ⎞ to the elements of
⎜ ⎟ ⎜ ⎟
⎜ R 21 R 23 R 23 ⎟ ⎜ R 21 R 23 R 2 3 ⎟ = I 3× 3 this subgroup as
⎜ R R 3 3 ⎟⎠ ⎜ R R 3 3 ⎟⎠
⎝ 31 R 32 ⎝ 31 R 32
R11 R12 R13
proper rotations;
R 21 R 23 R 23 = ± 1 Applications in one excludes
R 31 R 32 R 33
somewhat advanced those with
topics! negative
+ R O T A T IO N S Sophus Lie 31
- P A R IT Y determinant.
Polar vectors and pseudo- or axial-vectors
under inversion,
r r
r → −r ,
ur ur
p → −p
r r ur
but if l = r × p,
r r
then, under inversion, l → −l
Axial vector (pseudo vector) is a vector which does not transform
like a position vector under reflection. Its components are
governed by a different transformation law with respect to
rotation of the coordinate system. 32
r r r r r ur
l = r× p
C = A× B
angular momentum
ur
R R O R p
MI
r
r r
l right − hand −cross − product r r ur
r ur ⊗ l = r× p
= r image × p image 33
Proper- or Polar-vectors (odd parity)
Polar vector
3D
inversion
Pseudo vector
3D inversion
34
Reflection ≠ Rotation
d
[ Product of two functions ]
dt
To get acceleration, d d
we have to do that
twice! dt dt 40
(eˆρ , eˆϕ ) are not constant vectors,
Y
whereas (eˆx , eˆy ) are.
eˆρ = eˆρ ( ρ , ϕ )
eˆϕ1
“Unit Circle” eˆρ 1 eˆϕ = eˆϕ ( ρ , ϕ )
1
ϕ X
O
∂eˆρ
= − sin ϕ ( cos ϕ eρ − sin ϕ eϕ ) + cos ϕ ( sin ϕ eρ + cos ϕ eϕ )
) ) ) )
∂ϕ
)
= eϕ
42
PCD-09
Geometrical Consider (eˆρ , eˆϕ )
determination Y
at two neighboring
∂ê ρ points, infinitesimally
of
∂ϕ −eˆρ 1 eˆρ 2 close to each other.
Y
−eˆρ 1 eˆρ 2
2
eˆρ = cos ϕ eˆx + sin ϕ eˆy
“Unit Circle” eˆϕ = − sin ϕ eˆx + cos ϕ eˆy
δϕ 1 eˆρ 1
ϕ X
O
∂eˆρ ∂eˆρ
= 0, = eˆϕ ,
∂ρ ∂ϕ
∂eˆϕ ∂eˆϕ
= 0, = −eˆρ
PCD-09
∂ρ ∂ϕ 44
If ξ = ξ (u ) and u = φ ( x), chain rule
dξ
then will be a measure of the d ξ ⎛ d ξ ⎞ ⎛ du ⎞
dx =⎜ ⎟⎜ ⎟
sensitivity of ξ to changes in x : dx ⎝ du ⎠ ⎝ dx ⎠
If ξ = ξ (u , v)
where u = u ( x), v = v( x),
dξ ⎛ ∂ξ ⎞ ⎛ du ⎞ ⎛ ∂ξ ⎞ ⎛ dv ⎞
the rate at which ξ will = ⎜ ⎟⎜ ⎟ + ⎜ ⎟⎜ ⎟
dx ⎝ ∂u ⎠ ⎝ dx ⎠ ⎝ ∂v ⎠ ⎝ dx ⎠
change with respect to x
will be given by :
If ξ = ξ (u , v, x) where u = u ( x), v = v( x), the rate at
which ξ will change with respect to x will be given by :
dξ ⎛ ∂ξ ⎞⎛ du ⎞ ⎛ ∂ξ ⎞ ⎛ dv ⎞ ⎛ ∂ξ ⎞
= ⎜ ⎟⎜ ⎟ + ⎜ ⎟ ⎜ ⎟ + ⎜ ⎟ . &
Eg. L(q (t ), q (t ), t )
dx ⎝ ∂u ⎠⎝ dx ⎠ ⎝ ∂v ⎠ ⎝ dx ⎠ ⎝ ∂x ⎠ 45
PCD-09
Example: Velocity of a molecule of water
in a river – depends on ‘where’ and
‘when’ both!
dv d
= v ( r (t ) , t )
dt dt
Time Dependence: both
implicit and explicit
v ( x ( t ), y ( t ), z ( t ), t )
dv d
=
dt dt
d v ∂ v dx ∂ v dy
= + +
dt ∂ x dt ∂ y dt
∂ v dz ∂ v
+
∂ z dt ∂t 46
Elemental area Position vector & Velocity
in plane polar coordinates in plane polar coordinates
dA = ρ d ρ dϕ ur
ρ = ρ eˆρ
Y
uuur
ρ dϕ d ρ = (d ρ )eˆρ + ρ deˆρ
ur
dρ r r d ρ d ( ρ eˆρ )
v = ρ& = =
dϕ dt dt
ρ dρ deˆρ
= eˆρ + ρ
ϕ dt dt
∂eˆρ ∂eˆρ
R 2π 2
R
X = 0, = eˆϕ ,
∂ρ ∂ϕ
∫ρ =0 ϕ∫=0 ρ ρ ϕ = π = π 2
d d 2 R
2 ∂eˆϕ ∂eˆϕ
= 0, = −eˆρ
∂ρ ∂ϕ 47
PCD-09
∂eˆρ ∂eˆρ Motion of a particle
= 0, = eˆϕ , in plane polar coordinates
∂ρ ∂ϕ
∂eˆϕ ∂eˆϕ deˆρ ∂eˆρ
= 0, = −eˆρ = ϕ& = eˆϕϕ&
∂ρ ∂ϕ dt ∂ϕ
and chain rule
Time-dependence deˆϕ ∂eˆϕ
of unit vectors = ϕ& = −eˆρϕ&
dt ∂ϕ
ur
r r d ρ d ( ρ eˆρ ) d ρ deˆ
&
v=ρ = = = eˆρ + ρ ρ
dt dt dt dt
= ρ& eˆρ + ρϕ& eˆϕ
48
Radial velocity and Azimuthal velocity
PCD-09
r
v = ρ& eˆ ρ + ρϕ& eˆ ϕ instantaneous velocity
deˆρ ∂eˆρ
= ϕ& = eˆϕϕ&
dt ∂ϕ
and
deˆϕ ∂eˆϕ
= ϕ& = −eˆρϕ&
dt ∂ϕ
acceleration
r
r dv deˆρ deˆϕ
a = = ρ&& eˆ ρ + ρ& + ρϕ
& & eˆϕ + ρϕ&& eˆϕ + ρϕ&
dt dt dt
r
⇒ a =(ρ&& − ρϕ& ) eˆ ρ + (2ρϕ
2
& & + ρϕ&&)eˆ ϕ
49
PCD-09
Position vector is d rr = d ρ eˆ ρ + ρ d ϕ eˆ ϕ + dz eˆ z
r
r = ρ eˆ ρ + z eˆ z Z right handed triad
ρdρdφêz
of unit vectors
{eˆ ρ ,eˆ ϕ ,eˆ z }
dz
ρ
ρdφdzêρ
r
dρdz êϕ z
φ
Y
dφ ρdφ
X
dρ V = ρdρdφdz
50
PCD-09
General Reference on Vector analysis :
Kepler Problem.
Laplace-Runge-Lenz vector, ‘Dynamical’ symmetry.
Conservation principle ↔ Symmetry relation.
SYMMETRY CONSERVATION
LAWS
4
PCD-09
Space is homogenous and isotropic
∂L d ⎛ ∂L ⎞ ∂L
since − ⎜ &⎟ = 0, this means i.e. = p is conserved.
∂q dt ⎝ ∂q ⎠ ∂q&
i.e., is independent of time, is a constant of motion
Nr (i) →
Displacement is arbitrary, hence, ∑ Fk = 0;
k =1
→ d N ur (i) d →
i.e. 0 = ∑ pk = P.
dt k = 1 dt
r
Thus P is conserved in the absence of external forces.
The conservation of momentum was secured in ‘Unit 1’ on
the basis of ‘translational invariance in homegenous
space’, and not on the basis of Newton’s III law.
PCD-09
ur
dP r
=0
dtuur uur
d p2 d p1
=−
urdt ur dt
F 12 = − F 21
Given the symmetry related to
translation in homogenous
space, could you have
discovered Newton’s 3rd law?
7
Kepler-Bohr problem
ρ
ϕ
∂L d ⎛ ∂L ⎞ We explore
− ⎜ &⎟=0 the
∂q dt ⎝ ∂q ⎠
connection
1
L = T − V = m( ρ& 2 + ρ 2ϕ& 2 ) − V (ρ )
between
2 symmetry
∂L and
= 0, ϕ is cyclic conservation
∂ϕ laws further
∂L
pϕ = = constant = mρ 2ϕ& = angular momentum
∂ϕ&
8
PCD-09
How did Kepler deduce that planetary orbits are
ellipses around the sun ?
Johannes Kepler
1571- 1630
9
PCD-09
Johannes Kepler Galileo Galilei Isaac Newton
1571- 1630 1564 - 1642 (1642-1727)
b
a
Rosette motion
For this potential and
associated field:
no precession of orbit.
ur ur ur this figure.
A = p × L − m keˆ ρ
17
PCD-09
ur ur ur
A = p × L − mkeˆρ
ur ur ur ur
dA dp u
r u
r dL deˆρ r dL r
= × L + p× − mk
dt dt dt dt τ= = 0.
dt
Central force: No torque!
ur ur
d A d p ur deˆρ
= × L − mk
dt dt dt
deˆρ ∂eˆρ
But = ϕ& = eˆϕϕ&
dt ∂ϕ
ur ur
d A d p ur
Hence = × L − mkeˆϕϕ&
dt dt
18
PCD-09
ur ur
ur ur ur d A d p ur
A = p × L − mkeˆρ ⇒ = × L − mkeˆϕϕ&
dt dt
ur
dp
What is the form of the force: ?
dt
If we know it, we can plug it in the above
expression
- but for us to do be able to do that,
we must know it!
We must know the form of the interaction!
19
PCD-09
ur ur ur
A = p × L − mkeˆρ Force:
r
? dp k
ur ur = − ˆ
eρ
d A d p ur dt ρ 2
= × L − mkeϕϕ
ˆ &
dt
ur dt Newton told us!
dA ⎛ k ⎞ u
r
= ⎜ − 2 eˆρ ⎟ × L − mkϕ& eˆϕ
dt ⎝ ρ ⎠
ur
dA ⎛ k ⎞
= ⎜ − 2 eˆρ ⎟ × ( mρ 2ϕ& eˆz ) − mkϕ& eˆϕ
dt ⎝ ρ ⎠ ur
= ( − keˆρ ) × ( mϕ& eˆz ) − mkϕ& eˆϕ
dA
right handed triad
of unit vectors: {eˆ ρ ,eˆ ϕ ,eˆ z }
dt
ur
r
−eˆρ × eˆz = eˆϕ = ( mkϕ& eˆϕ ) − mkϕ& eˆϕ = 0
dA
20
PCD-09
dt
ur
dA r
Conclusion: =0
dt
ρ
ϕ
ur ur ur
perigee A = p × L − m keˆ ρ apogee
ur ur
A. L = 0 LRL vector is in the plane of the orbit
ur
Find its direction at perigee and at apogee. A
Direction of the LRL vector is: focus to perigee.
ur
:Must remain constant
A
– no matter where the planet is!
This is precisely what FIXES the orbit!
21
PCD-09
r
dp k
Force: = − 2 eˆρ Newton told us!
dt ρ
A = p × L − m keˆ ρ
ur
For L (angular momentum vector)
to be conserved,
ur uur
any central force would do.
A For
dA
= 0,
[Geometrical Symmetry] dt
1 one requires
L = T − V = m ( ρ& 2 + ρ 2ϕ& 2 ) − V ( ρ ) r
2 dp k
= − 2 eˆ ρ
k dt ρ
V (ρ ) = −
ρ D Y N AM IC AL
Reference: Goldstein’s ‘Classical Mechanics’, Section 9, Chapter 3.
SY M M ET R Y
23
PCD-09
Symmetry & Conservation Principles!
Eugene
Paul
Wigner
Emmily Amalie (1902-1995)
Noether
24
(1882 – 1935)
Properties of an ellipse: ρ y ρ'
locus of all points ϕ
such that ρ + ρ ' = 2 a x
ρ ' = 2a − ρ 2λ
( λ − ) + = ρ
2 2 2
2 x y '
Pythagoras Theorem
ρ ,ϕ
relation a b
for
ellipse 2λ
b = a −λ
2 2 2
25
ρ ,ϕ ρ + ρ ' = 2a ρ y ρ'
ϕ
relation ρ ' = 2 a − ρ
x
for 2λ
ellipse
(2λ − x ) + y = ρ '
2 2 2
( 2λ − ρ cos ϕ ) + y = ( 2 a − ρ )
2 2 2
( 2λ − ρ cos ϕ ) + ( ρ sin ϕ ) = ( 2 a − ρ )
2 2 2
a2 − λ 2 b2
ρ= =
a − λ cos ϕ a − λ cos ϕ
⎛λ⎞
1 − ⎜ ⎟ cos ϕ
1 a − λ cos ϕ ⎝ a⎠ 1 − e cos ϕ
= = =
ρ b 2
⎛b ⎞
2
⎛ b2 ⎞
⎜ ⎟ ⎜ ⎟
⎝ ⎠
a ⎝ a ⎠
27
ur ur ur
A = p × L − m keˆ ρ T he ( ρ , ϕ )
ur ur ur relation
A. L = 0 l =| L | describes the
r ur familiar ellipse.
Examine: r • A = ρ Acosϕ
r ur r ur ur LRL vector
r • A = r • p × L − mk ρ provides an
r ur ur alternative way to
= r × p • L − mk ρ
determine the
= l2 − mk ρ ‘orbit’ equation for
1 mk ⎡ A ⎤ the Kepler
ρ
= 2
l ⎢⎣1 + m k cos ϕ ⎥⎦ problem.
28
PCD-09
ur ur ur
A = p × L − mkeˆρ 1 1 − e cos ϕ “usual”
=
r ur ρ ⎛ b 2
⎞ ellipse
ur ⎡ v× L ⎤ r ⎜ ⎟
A = mk ⎢ − eˆρ ⎥ = mka ⎝ a ⎠ “LRL”
⎣ k ⎦ ellipse
r ur
r v× L 1 mk ⎡ A ⎤
where a = − eˆρ = 2 ⎢1 + cos ϕ ⎥
k ρ l ⎣ mk ⎦
r
a =e , The above ρ , ϕ relation
orbit's eccentricity. expresses the ellipse with the
eccentricity identified as defined
above. Details : bit laborious!
We thus get equation to the orbit using
the LRL vector properties. 29
ur ur ur Conservation law
A = p × L − m keˆ ρ associated with
Laplace R unge Lenz V ector : ‘dynamical/accidental’
1 symmetry.
constant for a strict potential.
r Remember Noether?
Wilhelm
Carl David Tolmé Lenz
Runge 1888 -1957
Pierre-Simon Laplace
1856 - 1927 Symmetry of the H atom: ‘old’
1749 - 1827
quantum theory. En ~ n-2
30
PCD-09
Further reading, beyond PH101
• Symmetry & Conservation laws play an
important role in understanding the very frontiers
of Physics.
• The implications go as far as testing the
‘standard’ model of physics, and exploring if
there is any physics beyond the standards
model.
31
Continuous Symmetries - Translation, Rotation
Discrete Symmetries
- P : Parity
- C : Charge Conjugation
- T : Time Reversal
32
Electroweak Unification
The PHOTON, along with the W and Z provides unification of
weak and Electromagnetic (Electrodynamic) interactions.
Mass of W & Z particles: ~80 and 90 Gev, respectively;
Photon is mass-less.
33
Lorentz symmetry: associated with another fundamental
symmetry of particle physics, the PCT symmetry.
http://www.relativitycalculator.com/laws_of_conservation.shtml
37
PCD-09
“In the judgment of the most competent living
mathematicians, Fräulein Noether was the most significant
creative mathematical genius thus far produced since the
higher education of women began.”
Learning goals:
x
y = r sin θ sin ϕ
x = r sin θ cos ϕ Y
ϕ
x = r sin θ cos ϕ
ρ = r sin θ
y = r sin θ sin ϕ
X
z = r cos θ
PCD-09 2
Transformations of the unit vectors
Y
ϕ
Recognize the distance
between the old position
ρ = r sin θ and the new position to be rdθ
PCD-09 4
Z
ρ eˆr ρ dϕ = r sin θ dϕ
êϕ The point
dθ êθ
is displaced
to a new point
r on the sphere r = constant
and on the surface of the
PCD-09 5
dV =
Volume spanned Z
by the
( dr )(rdθ )(r sin θ dϕ )
three = r si n θ dr dθ dϕ
2
displacements
through dr
dr
rdθ r
r sin θ dϕ r
θ dθ
Y
φ dϕ
X
PCD-09 6
uur
dS r = rdθ eˆθ × r sin θ dϕ eˆϕ = r 2 sin θ dθ dϕ eˆr
uur
Z dS ϕ = dreˆr × rdθ eˆθ = rdrdθ eˆϕ
Determination of r
infinitesimal vector dr dSθ = r sin θ d ϕ eˆϕ × dreˆr
surface elements
made up of r = r sin θ drd ϕ eˆθ
incremental sides r
θ dθ Tiny volume contained
Y between the three pairs
of surfaces
φ dϕ
dV =
( dr )( rdθ )(r sin θ dϕ )
X
= r 2si n θ dr dθ dϕ
PCD-09 7
Note that (eˆρ , eˆϕ )
Do you remember how Y we got,
are not constant vectors,
in Unit_3,
the CHANGES δ ê whereas (eˆx , eˆy ) are.
in the unit vectors for
the plane polar coordinates? Consider (eˆρ , eˆϕ )
Use the same technique! at two neighboring
δϕ points, infinitesimally
eˆρ = cos ϕ eˆx + sin ϕ eˆy
eˆϕ = − sin ϕ eˆx + cos ϕ eˆy ϕ close to each other.
X
O
∂eˆρ ∂eˆρ
eˆρ 2 − eˆρ 1 = δϕ eˆϕ eˆϕ 2 − eˆϕ 1 = −δϕ eˆρ = 0, = eˆϕ ,
∂ρ ∂ϕ
∂eˆρ ∂eˆϕ ∂eˆϕ ∂eˆϕ
= eˆϕ = −eˆρ = 0, = −eˆρ
∂ϕ ∂ϕ ∂ρ ∂ϕ
PCD-09 8
Z
ρ eˆr 1st position
eˆr 2 − eˆr1 = δθ eˆθ The point
−eˆr 2nd position
is displaced
dθ êθ
In the limit to a new point
δθ → 0, r on the sphere
of same radius
∂eˆr θ r = constant
= eˆθ
∂θ and in the same plane
ϕ = constant
Y
ϕ
Distance between the 1st
position and 2nd position
is rdθ
X
PCD-09 9
Partial derivatives of the unit vectors with
respect to the coordinates:
For example:
∂eˆϕ
= − cos ϕ eˆx − sin ϕ eˆy
∂ϕ
= − cos ϕ ( sin θ cos ϕ eˆr + cos θ cos ϕ eˆθ − sin ϕ eˆϕ )
− sin ϕ ( sin θ sin ϕ eˆr + cos θ sin ϕ eˆθ + cos ϕ eˆϕ )
= − sin θ eˆr − cos θ eˆθ
But, How?
Pictures from
What is relativity? Landau and Rumer; Dover 1993
PCD-09 13
z
êϕ
r êθ
θ
rcosθ
rs y
inθ
φ r sinθ cosϕ
r sinθ sinϕ
x
¾ the surface of constant r sphere of radius r
¾ the surface of constant θ surface of a cone
¾ the surface of constant ϕ plane perpendicular to xy
¾ constant r and θ ‘mouth’ of the cone
¾ constant r and φ arc of a circle
PCD-09
¾ constant θ and φ line segment 14
Motion in in spherical polar:
Velocity and acceleration
Infinitesimal displacement
r
Position vector r = reˆr
r
dr = dreˆr + rdeˆr
r
dr = dreˆr + rd θ eˆθ + r sin θ d ϕ eˆϕ
r ∂ eˆr ∂ eˆr
⇒ dr = dreˆr + r δθ + r δϕ
∂θ ∂ϕ
PCD-09 15
Motion in in spherical polar:
Velocity and acceleration
r
dr = dreˆr + rd θ eˆθ + r sin θ d ϕ eˆϕ
r
& ˆr + rθ&eˆθ + r sin θϕ& eˆϕ
v = re
r
r dv &
a= = ( r - rθ - r sin θϕ ) eˆr
&& 2 2
& 2
dt
+ (2 r&θ − r sin θ cos θϕ + rθ&&) eˆθ
& & 2
PCD-09 16
Unit 6: Inertial and non-inertial reference frames.
Moving coordinate systems. Pseudo forces.
Inertial and non-inertial reference frames.
Deterministic cause-effect relations in inertial frame,
and their modifications in a non-inertial frame.
ur r
• The same force F = ma explains the
linear response (effect/acceleration is
linearly proportional to the cause/interaction)
relationship in all inertial frames.
PCD-09 7
Galilean relativity
Time t is the same in the
red frame
and
in the green, double-primed frame.
PCD-09 8
Time t is the same Galilean relativity
in the red frame r ur ur
and in the green, r ⎛ d ⎞ ⎛ d ⎞ ⎛ d ⎞ ⎛ d ⎞
double-primed ⎜ ⎟ ⎜ ⎟ r (t ) = ⎜ ⎟ ⎜ ⎟ r ''(t ) + f
r (t ) ⎝ dt ⎠ ⎝ dt ⎠ ⎝ dt ⎠ ⎝ dt ⎠
frame. Y r r r ur
r ''(t ) a = a "+ f
F
O X uuuur r 1 ur 2
Z OO '' = vt + f t
2 Y’’
r r r 1 ur 2 F”
r (t ) = r ''(t ) + vt + f t
2 X’’
O’’
⎛ d ⎞r ⎛ d ⎞ ur r 1 ur
⎜ ⎟ r (t ) = ⎜ ⎟ r ''(t ) + v + f ( 2t ) Z’’
⎝ dt ⎠ ⎝ dt ⎠ 2
⎛ d ⎞r ⎛ d ⎞ ur r ur ur uur r ur
⎜ ⎟ r (t ) = ⎜ ⎟ r ''(t ) + v + f t F " = ma " = m(a − f )
⎝ dt ⎠ ⎝ dt ⎠ ur r ur uur ur
F " = m(a − f ) = F − F pseudo
‘Real’/’Physical’
force/interaction
PCD-09 9
Galilean relativity
r ur uur r ur
r (t ) r F " = ma " = m(a − f )
r ''(t ) ur r ur ur ur
Y
F " = m(a − f ) = F − F pseudo
F
O X uuuur r 1 ur 2
OO '' = vt + f t
2 Y’’
Z F”
X’’
Z’’ O’’
Laws of Mechanics are not the same in a FRAME OF
REFERENCE that is accelerated with respect to an
inertial frame. The force/interaction which explained the
acceleration in the inertial frame does not account for the
acceleration in the accelerated frame of reference.
Same cause-effect relationship does not explain dynamics.
PCD-09 10
uur Galilean relativity
ma " =
Time t is the same in the
red frame and in the green, double-primed frame.
ur ur ur
F " = F real / physical − F pseudo
Note! An observer in frame F” may not
be aware of which frame is in –
whether it is an inertial or an
accelerated one with respect to an
inertial one!
Real effects of pseudo-forces!
P. Chaitanya Das, G. Srinivasa Murthy, Gopal Pandurangan
and P.C. Deshmukh
Resonance, Vol. 9, Number 6, 74-85 (2004)
http://www.ias.ac.in/resonance/June2004/pdf/June2004Classroom1.pdf
PCD-09 11
Galilean relativity: Time t is the same in all frames of references.
RELATED ISSUES:
Weightlessness
http://www.gozerog.com/
PCD-09 16
Not Just Another Old Flame !
PCD-09 18
Galilean relativity
Time t is the same in the
red frame
and
in the green, double-primed frame.
PCD-09 20
We have considered earlier the
translational motion of the new frame
relative to the inertial frame FI; so for
the current discussion on rotating
frame of reference, we shall ignore
that.
PCD-09 21
ZI
Galilean relativity Time t is the same in the
C red frame
ZR
and
in the purple, rotating, frame.
YR
F
YI
OI
OR
First, we consider the position
FR
XI vector of a particle at rest in the
XR
rotating frame.
PCD-09
24
Galilean relativity
Time t is the same in
the red frame
and in the purple,
rotating, frame.
The time-derivative is
NOTE! different in the red frame
and
in the purple, rotating, frame.
⎛d ⎞ ⎛d ⎞
⎜ ⎟ ≠⎜ ⎟
⎝ dt ⎠ I ⎝ dt ⎠ R
PCD-09 25
1. Position vectors in the two frames are always the same!
r r
rI = rR
2. Of course, their respective Cartesian components are
different!
Operator ⎛d ⎞ ⎛d ⎞ ur
Equivalence: ⎜ ⎟ =⎜ ⎟ + ω ×
⎝ dt ⎠ I ⎝ dt ⎠ R
PCD-09 29
⎛d ⎞ ⎛d ⎞ ur ⎛d ⎞ r ⎛d ⎞ r ur r
⎜ ⎟ =⎜ ⎟ + ω × ⎜ ⎟ r =⎜ ⎟ r + ω × r
⎝ dt ⎠ I ⎝ dt ⎠ R ⎝ dt ⎠ I ⎝ dt ⎠ R
Operating twice:
⎛ d ⎞ ⎛ d ⎞ r ⎛ d ⎞ ⎧⎛ d ⎞ r ur
r⎫
⎜ ⎟ ⎜ ⎟ r = ⎜ ⎟ ⎨⎜ ⎟ r + ω × r⎬ +
⎝ dt ⎠ I ⎝ dt ⎠ I ⎝ dt ⎠ R ⎩⎝ dt ⎠ R ⎭
ur
⎧⎛ d ⎞ r ur r ⎫
+ ω ×
⎨⎜ ⎟ r + ω × r ⎬
⎩⎝ dt ⎠ R ⎭
⎛ d 2 ⎞ r ⎛ d 2 ⎞ r ⎛ d ⎞ ur r ur ⎛d ⎞ r
⎝ dt ⎠ I ⎝ dt ⎠ R ⎝ dt ⎠ R
(
⎜ 2 ⎟ r = ⎜ 2 ⎟ r +⎜ ⎟ ω × r +ω × ⎜ ⎟ r )⎝ dt ⎠ R
ur ur r
(
+ ω × ω × r )
Multiplying by mass ‘m’, we shall get quantities
PCD-09
that have dimensions of ‘force’.
30
ur
⎛d ⎞ r2
⎛d ⎞ r 2
⎛ dω ⎞ r ur ⎛d ⎞ r
m ⎜ 2 ⎟ r = m ⎜ 2 ⎟ r − m ⎜⎜ ⎟⎟ × r − 2mω × ⎜ ⎟ r +
⎝ dt ⎠ R ⎝ dt ⎠ I ⎝ dt ⎠ R ⎝ dt ⎠ R
ur ur r
- mω × ω × r ( )
r r r r ⎛d ⎞ r r r r
FR = FI − Fω& − 2mω × ⎜ ⎟ r - mω × (ω × r )
⎝ dt ⎠ R
‘Coriolis force’ ‘Centrifugal force’
‘Leap second’ term Gaspard Gustave
de Coriolis term
1792 - 1843
PCD-09 31
Why are Leap Seconds Used?
Leap seconds are added to keep the atomic clocks
synchronized with the Earth's rotation. The Earth rotates at a
slower pace over time while the atomic clocks do not slow
down. The International Earth Rotation and Reference
System Service (IERS) (IERS) decides when to introduce a
leap second in UTC (Coordinated Universal Time).
On one average day the difference between atomic clocks
and Earth's rotation is around 0.002 seconds, or around 1
second every 1.5 years.
The time to do one rotation differs from day to day and from
year to year. The Earth was slower than atomic clocks by:
0.16 seconds in 2005; 0.30 seconds in 2006; 0.31 seconds in
2007; and 0.32 seconds in 2008. It was only 0.02 seconds
slower in 2001.
http://www.timeanddate.com/time/leapseconds.html ; 14th October, 2009
32
The atomic clocks are occasionally instructed to add an
extra second, known as the leap second, to synchronize
the atomic clocks with the Earth's observed rotation.
Real
Effects of
Pseudo
A Forces
Please read the paper at the following internet weblink:
PCD-09
http://www.physics.iitm.ac.in/~labs/amp/homepage/dasandmurthy.htm
34
Force = Rate of change of momentum
What interaction is making the instantaneous
momentum change?
R = 10m
Very fundamental question!!
ω = 0.95 rad/sec
v = (-10.0325,2.755) m/s The answer
P. Chaitanya Das, G. determines
Srinivasa Murthy,
Gopal Pandurangan
our notion of
and P.C. Deshmukh a
‘fundamental
'The real effects of
pseudo-forces',
Resonance, Vol. 9,
Number 6, 74- interaction’
85(2004)
(http://www.ias.ac.in/res
onance/June2004/pdf/
June2004Classroom1.
PCD-09
pdf) 35
The plane of oscillation of
the Foucault pendulum is
seen to rotate due to the
Coriolis effect. The plane
rotates through one full
rotation in 24 hours at
poles, and in ~33.94
“ Foo-Koh” hours at a latitude of 450
(Latitude of Paris is
~490).
ur ur ur ur ⎛ d ⎞ r ur ur r
F R = F I − F ω& − 2mω × ⎜ ⎟ r - mω × ω × r
⎝ dt ⎠ R
( )
‘Leap second’ term ‘Coriolis force’ ‘Centrifugal force’
PCD-09 36
To estimate the REAL effects of the PSEUDO-FORCES
(cf: centrifugal and Coriolis), we use in the following
figure a Cartesian coordinate system with reference to a
point on the earth’s surface.
δ Depends on
the latitude.
PCD-09 40
r ⎛d ⎞ r r
r r
(
r
)
FCoriolis = −2mω × ⎜ ⎟ r ω = ω • eˆy eˆy + (ω • eˆz ) eˆz
⎝ dt ⎠ R
By plugging in the magnitude and direction of the
⎛d ⎞ r
velocity ⎜ ⎟ r
⎝ dt ⎠ R
of an object in the above
expression, one can easily determine the effect of
the Coriolis force. The effects are quite dramatic!
How would you estimate the Coriolis deflection of an
object ‘falling’ at a point on earth’s surface? Note that
it would have a velocity along −eˆz .
eˆy = eˆNorth , pointing toward North.
eˆz = eˆvertical ; eˆx is toward EAST (−eˆy ) × (−eˆz ) = eˆx (East)
PCD-09 41
ur
ω Western
North Pacific
Coriolis
1792 - 1843
ll o n th e N o r thern
a s t a te o f fr ee f a , i n the
An o b je c t in ward th e E a st
t s d e f l e c t ed t o r d West!
he r e g e t ed t o w a
hemisp e w o ul d g et d eflec
th er n h e m i s pher
So u
c t fra l ling
r t hr⎛od ⎞ r ugh 100
r At ar latitude ofr 60 0 a n ob j e
FCoriolis = −2mω × ⎜ ⎟ r
( )
ω = ω • eˆy eˆy +ef(lωec•teeˆzd) eˆtzhrough ~1 cm. ⎝ dt ⎠ R
meters is d
Western
South Pacific
PCD-09 http://agora.ex.nii.ac.jp/digital-typhoon/start_som_wnp.html.en 42
http://daphne.palomar.edu/pdeen/Animations/34_Coriolis.swf
r r ⎛d ⎞ r
FCoriolis = −2mω × ⎜ ⎟ r
⎝ dt ⎠ R
http://www.physics.iitm.ac.in/~labs/amp/homepage/dasandmurthy.htm
1. P. Chaitanya Das, G. Srinivasa Murthy, Gopal Pandurangan
and P.C. Deshmukh
'The real effects of pseudo-forces', Resonance, Vol. 9, Number 6, 74-85(2004)
(http://www.ias.ac.in/resonance/June2004/pdf/June2004Classroom1.pdf)
PCD-09 44
Unit 7: Galilean & Lorentz transformations.
Special Theory of Relativity
Learning goals:
PCD-09 1
Unit 7: Galilean & Lorentz Transformations.
Special Theory of Relativity.
Smoking is
Hendrik Antoon Lorentz
1853-1928 injurious to
Galileo Galilei
1564 - 1642 health!
Albert Einstein
PCD-09 1879-1955 2
Loop : ×
Stationary B
Lorentz
force (a) Clockwise Current
(b) Counterclockwise Current
predicts: √ (c) No Current
×
Strength of B decreased. B I
Nothing is moving,
but still, current seen!!! Decreasing B↓
Einstein:
dB Special Theory of Relativity
I ∝
dt What matters is relative motion
‘Relativity’ applied to EM fields.
PCD-09 4
Einstein: Special Theory of Relativity
(1905: http://lorentz.phl.jhu.edu/AnnusMirabilis/)
PCD-09 6
Why did the chicken cross the road?
PCD-09 7
ur ur r ur
F = q[ E + v x B]
uur
v is obviously not the same in different frames of references moving
with respect to each other, even if their relative velocity is constant.
Nevertheless, laws of physics must be the same in all inertial frames
of references. The consistency in the laws of physics in all inertial
frames of references is then achieved by Einstein's bold proposal:
1
speed of light, c= constant, in all inertial
μ 0ε 0
frames, contrary to Galilean principles,
ur ur uur uur
and x, y, z , t ↔ x',y',z',t', and E , B ↔ E ', B ' .
P. Chaitanya Das, G. Srinivasa Murty, K. Satish Kumar, T A. Venkatesh and P.C. Deshmukh
'Motion of Charged Particles in Electromagnetic Fields and Special Theory of
Relativity', Resonance, Vol. 9, Number 7, 77-85 (2004)
PCD-09 http://www.ias.ac.in/resonance/July2004/pdf/July2004Classroom3.pdf 8
1. Maxwell’s equations are correct in all inertial frames of references.
1
2. Maxwell’s formulation predicts that EM waves travel at the speed c= .
μ 0ε 0
3. HENCE, light (EM waves) travels at the constant speed c= 1
in all inertial frames of references. μ 0ε 0
PCD-09 10
Just what does it mean to say that “Light (EM waves)
travels at the constant speed in all inertial frames of
references” ?
ê y ê y'
r c=
1
μ 0ε 0
L v
R
~~ ~~
ê x ê x'
OI The blue frame moves
ê z ê z' toward the right at a constant
Inertial
frame velocity fc where 0<f<1.
The laser guns L & R fire light toward the blue frame. Both the
beams of light, from L & from R, are seen to arrive in the blue
frame at the speed c and not c±v. COUNTER-INTUITIVE !
PCD-09 11
Light (EM waves) travels at the constant speed in all inertial frames of references.
Experimental proof: A. A. Michelson and E. W. Morley, "On the Relative Motion of the
Earth and the Luminiferous Ether," American Journal of Science, 34, 333-345 (1887).
L F
R
r
v
S
Light sensors Stationary observer
L F
R
r
v
F would give out GREEN LIGHT if R
senses light before L, or it would give
Stationary S out RED LIGHT if both L & R sense light
observer simultaneously (and also give out red
light if L senses light before R).
Now, both the firecrackers explode simultaneously as seen by S just
as M crosses S. Let us now catalogue S’s observations:
(1) S detects both the flashes simultaneously.
(2) S is half-way between the firecrackers when the two explosions
occurred.
(3) Light from both explosions travels at equal speed toward S.
S must expect F to flash RED LIGHT. Contradiction can be
removed only by
However, recognizing that the
M would ‘expect’ F to flash out GREEN LIGHT. bursts are not
simultaneous in M’s frame !
PCD-09 15
Moving observer M
L F
R
r
v
F would give out GREEN LIGHT if R
Stationary S senses light before L, or it would give
observer out RED LIGHT if both L & R sense light
simultaneously (and also give out red
light if L senses light before R).
M would expect to see GREEN LIGHT, and he must conclude that the
two bursts were not simultaneous. Instead, the right burst occurred
before the left one!
STR upshot!
We must revise our notions of ‘simultaneity’.
Events that are ‘simultaneous’ in one frame of reference S are not
so in another frame of reference S’ that is moving relative to S.
PCD-09 16
STR
Time Dilation
Length Contraction
Lorentz Transformations
PCD-09 17
We shall first describe what is known as a ‘LIGHT CLOCK’
‘TIME KEEPING COUNTER’
⎣⎝ ⎠ ⎦ ⎣⎝ ⎠ ⎦
Time Dilation
2h / c Δt '
⇒ Δt = =
1− v / c
2 2
1− β 2 Conclusions do not depend on
where β =v/c. the use of the ‘Light Clock’. Any
clock would give the same
result.
PCD-09 20
2h / c Δτ
Δt = =
clock frame, S ′ 1 − v2 / c2 1− β 2
r where β =v/c
v and Δτ = Δt ' = PROPER TIME
Δt > Δτ
Time Dilation
“TIMEINTERVAL " Moving clocks go slow;
BETWEEN TWO time interval between
TICKS IS THE
SHORTEST IN two ticks is longer when
THE CLOCK’s measured in a frame in
REST FRAME.” which the clock is moving "
PCD-09 22
y We consider two stars in deeeeeep space,
an orange-star and a blue-star.
r
v
x
S
In frame S, the stars are at rest at a length L apart from
each other. The rocket moves from the star O to the star
B. An observer in frame S carries out measurements of
lengths and time intervals.
y’
x’
S’: Rocket frame
PCD-09 24
y
The rocket moves
r from the star O to r
star B at velocity
v in frame S. v
x
S In frame S, the objects are at rest at a length L apart from
each other.
This LENGTH L is therefore the “PROPER LENGTH”, l.
PCD-09 26
In the ROCKET-FRAME-S', it is the two stars that move to the left at
L'
speed v= , where Δ t' is the PROPER TIME (Δ τ ) measured in
Δ t'
S ' for the blue star to travel the LENGTH L '.
L l L' L' L'
v= = = = =
Δ t Δ t Δ t' Δ τ
( 1-β 2 ) ( Δt )
⇒ L ' = L 1-β 2 = l 1-β 2 ≤ l
y’
This is LORENTZ ( LENGTH )
CONTRACTION
x’
S’: Rocket frame
PCD-09 27
1902 Nobel Prize in Physics
Hendrik Antoon Lorentz
1853-1928
"in recognition of the extraordinary service they
rendered by their researches into the influence of
magnetism upon radiation phenomena"
Lorentz contraction!
Pieter Zeeman
1865-1943
http://www.bun.kyoto-u.ac.jp/~suchii/lorentz.tr.html
PCD-09 28
LORENTZ transformations (x,y,z,t) to (x’,y’,z’,t’)
Requirements:
Transformation equations must agree
with Galilean transformations when
v<<<c.
Transform both space and time
coordinates.
Ensure that speed of light is same in all
inertial frames of references.
PCD-09 29
Y Origins O and O’ of Y’ r
the two frames S v
X and S’ coincide at X’
Z
S t=0 and t’=0. Z’
S’
z'=z z=z' 1
=
1− β 2
⎛ vx ⎞ ⎛ vx' ⎞
t'= γ ⎜ t- 2 ⎟ t= γ ⎜ t'+ 2 ⎟ Note : γ → 1 as v → 0.
⎝ c ⎠ ⎝ c ⎠
Lorents transformations transform the space-time
PCD-09
coordinates of ONE EVENT. 30
References
Main Reference:
‘Physics for Scientists and Engineers’, II Edition, by
Randall D. Knight (Pearson, Addison-Wesley, 2007)
http://www.ias.ac.in/resonance/July2004/pdf/July2004Classroom3.pdf
‘twin paradox’”
http://mentock.home.mindspring.com/twins.htm
28/09/09 21:50
PCD-09 31
Seeta and Geeta are identical twins. Twin Paradox
Geeta stays at home,
3
and Seeta travels in a rocket at a speed 5 c for 4 yrs
measured in the rocket-clock (proper time).
Geeta’s home-based clock measures the corresponding
time interval as Δt = Δτ ; β = v/c = 3 .
1− β 2 5
( )
2
3 c
v2 5 32 4
1− 2 = 1− = 1− 2 =
c c2 5 5
5 5
Δt = Δτ = (4 yrs ) = 5 yrs “TIMEINTERVAL
4 4 BETWEEN TWO
Δτ = Δt ' = PROPER TIME TICKS IS THE
SHORTEST IN
Δt > Δτ (Time Dilation).
THE CLOCK’s
REST FRAME.”
PCD-09 32
Seeta now turns around, and returns at same
speed, taking another 4 years
(measured in rocket clock) to return,
during which Geeta’s clock advances by
another 5 years.
G
S
PCD-09 35
Sum of the velocities as per Galilean relativity:
3 3 6
c+ c= c 〉 c!
5 5 5
Impossible for Seeta to get that speed > c
However, this is not how relative velocity is
added! One must use Lorentz, not Galilean
relativity.
PCD-09 36
r uuuur
Z r (t ) ur O O ' = ( u x eˆ x ) t
Y r '( t ')
x ' = γ ( x − u xt )
O X
S y'= y
Z’
z'= z
Y’
ux x
O’ X’ t ' = γ (t − 2 )
c
r S’
r d r dx dy dz 1
v= = eˆ x + eˆ y + eˆ z γ =
2
dt dt dt dt 1− 2
u
ur c
ur d r '
v' =
dt '
dx' dy' dz'
= eˆ x + eˆ y + eˆ z
dt' dt' dt'
PCD-09 37
r uuuur
Z r (t ) ur O O ' = ( u x eˆ x ) t
Y r '( t ')
x ' = γ ( x − u xt )
O X
ux x
S Z’ t ' = γ (t − 2 )
c
r Y’ 1
r dr γ =
v= O’ u2
dt X’ 1− 2
S’ c
ur
ur d r ' dx' dy' dz'
v' = = eˆ x + eˆ y + eˆ z
dt ' dt' dt' dt'
dx' d ( γ ( x − u x t )) d ( x − u xt )
= = dx' v x − u x
dt' d ( γ ( t − u x x )) d ( t − u x x ) ⇒ =
dt' 1 − u x v x
c2 c2
c2
PCD-09 38
dx' v x − u x
=
dt' 1 − u x v x
c2 If the frame of reference
S’ is moving in the
negative x direction, we
shall get: dx' v x + u x
=
dt' ux vx
1+
c2
vx + ux
vx ' =
ux vx
1+ 2
c
vy ' = vy
vz ' = vz Twins!
PCD-09 39
G vx + ux
vx ' =
ux vx
S 1+ 2
c
G vy ' = vy
S vz ' = vz
G G
S
v relative =
v1 + v 2
=
( ) ( )
3 c+ 3 c
5 5 Seeta clocks 4 years in
v v
1+ 122
c 1+
( ) ( )
3 c 3 c
5 5 her clock and shoots off
2
c
toward Geeta at a speed
( 5 )c
6
( 5) ( )
15
v relative = = 6 c × 25
34
= 15
17
c of c , and in as
1+ 9 17
25 much time (4-Seeta-yrs),
Again, Geeta would age by 10
catches up with Geeta.
years, and Seeta by 8 years.
PCD-09 No paradox! 40
In some (published)
comments on the twin-
paradox, resolution has
been sought by invoking
Seeta’s acceleration when
she would begin Geeta’s
chase after 4 years
(previous slide).
One would then argue that we are no longer considering
inertial frames of references, since acceleration of Seeta’s
frame would come into the analysis. However, such an
‘explanation’ is not called for.
PCD-09 41
One can do away completely with Seeta’s
acceleration by considering in our thought-
experiment a third observer Jayalalitha
( weren’t there a set of triplets, rather than
mere twins? ) who does not undergo any
acceleration, but passes by Seeta and sets her
clock with her as she (J) breezes past her (S).
PCD-09 42
Jayalalitha
would then
catch up with
Geeta and
compare her
clock with
Geeta’s as
she crosses
her, and then send that information back to Seeta.
PCD-09 45
Unit 8: Physical examples of fields. Potential
energy function. Gradient, Directional Derivative,
Learning goals:
unstable
stable
neutral
stable
neutral
a c e d b
stable
3
PCD-09
z Primary, elementary entities in Physics:
r
ψ = ψ ( r , t ) = ψ ( x, y , z , t ) = ψ ( r , θ , ϕ , t ) = ψ ( ρ , ϕ , z , t )
ur ur r ur ur ur
A = A(r , t ) = A( x, y, z , t ) = A(r , θ , ϕ , t ) = A( ρ , ϕ , z , t )
6
PCD-09
DIRECTIONAL DERIVATIVE
is a SCALAR QUANTITY
which has a DIRECTIONAL ATTRIBUTE.
7
PCD-09
r r r r
dψ δψ ψ ( r + δ r ) −ψ ( r ) ψ (r)=ψ ( x, y, z )
= lim = lim
ds δ s →0 δ s δ s →0 δs r
ψ (r)=ψ ( ρ , ϕ , z )
r
∂ψ ∂ψ ∂ψ ψ (r)=ψ (r ,θ , ϕ )
δψ = δx+ δy+ δz
∂x ∂y ∂z Cartesian Coordinate
System
∂ψ ∂ψ ∂ψ
= δρ + δϕ + δz Cylindrical Polar
∂ρ ∂ϕ ∂z
Coordinate System
∂ψ ∂ψ ∂ψ
= δr + δθ + δϕ Spherical Polar
∂r ∂θ ∂ϕ
Coordinate System
8
PCD-09
r
The directional derivative ψ (r)= ψ ( x , y , z )
r r r r
dψ ψ ( r + δ r ) −ψ ( r ) ψ (r)= ψ ( ρ , ϕ , z )
= lim
ds δ s →0 δs r
ψ (r)= ψ ( r , θ , ϕ )
dψ ∂ψ dx ∂ψ dy ∂ψ dz Cartesian Coordinate
= + +
ds ∂x ds ∂y ds ∂z ds System
dψ ∂ψ d ρ ∂ψ dϕ ∂ψ dz Cylindrical Polar
= + +
ds ∂ρ ds ∂ϕ ds ∂z ds Coordinate System
dψ ∂ψ dr ∂ψ dθ ∂ψ dϕ Spherical Polar
= + +
ds ∂r ds ∂θ ds ∂ϕ ds Coordinate System
These expressions can be written very nicely by using the
uur
9
PCD-09 expressions for the displacement dr in various coordinate systems.
uur
dr = ? Cartesian Unit vectors cylindrical polar coordinate
are constant vectors; ∂eˆρ ∂eˆρ
= 0, = eˆϕ ,
Unit vectors of the ∂ρ ∂ϕ
cylindrical polar and the ∂eˆϕ ∂eˆϕ
spherical polar coordinate = 0, = −eˆρ
∂ρ ∂ϕ
systems are not!
∂eˆr r spherical polar coordinate
=0 ∂eˆθ r ∂eˆϕ r
∂r =0 =0
∂r ∂r
∂eˆr ∂eˆϕ r
= eˆθ ∂eˆθ =0
∂θ = −eˆr ∂θ
∂θ
∂eˆr ∂eˆϕ
= sin θ eˆϕ ∂eˆθ ∂ϕ
= − cos θ eˆθ − sin θ eˆr
∂ϕ = cos θ eˆϕ
PCD-09
∂ϕ 10
Position and Displacement vectors in various coordinate systems
r r
r = x eˆ x + y eˆ y + z eˆ z d r =eˆ x dx + eˆ y dy + eˆ z dz
r
r
r = ρ eˆ ρ + z eˆ z d r =(d ρ )eˆ ρ + ρ ( d eˆ ρ ) + ( dz )eˆ z
r
d r =eˆ ρ d ρ + eˆ ϕ ρ d ϕ + eˆ z dz
r r
r =reˆ r dr = (dr )eˆr + r (deˆr )
In order to avoid making
careless mistakes, always
r ⎡ ∂ eˆr ∂ eˆr ⎤
dr = ( dr ) eˆr + r ⎢ dθ + dϕ ⎥
try to write unit vectors ⎣ ∂θ ∂ϕ ⎦
first, differential r
elements last! dr = dreˆr + rd θ eˆθ + r sin θ d ϕ eˆϕ
Example: r
dr = eˆr dr + eˆθ rd θ + eˆϕ r sin θ d ϕ
11
PCD-09
r r r r
dψ δψ ψ ( r + δ r ) −ψ ( r ) ψ (r)=ψ ( x, y, z )
= lim = lim
ds δ s →0 δ s δ s →0 δs r
ψ (r)=ψ ( ρ , ϕ , z )
∂ψ ∂ψ ∂ψ r
δψ = δx+ δy+ δz ψ (r)=ψ (r ,θ , ϕ )
∂x ∂y ∂z
δψ ∂ψ δ x ∂ψ δ y ∂ψ δ z Cartesian Coordinate
= + +
δs ∂x δ s ∂y δ s ∂z δ s System
dψ δψ r
= lim d r =eˆ x dx + eˆ y dy + eˆ z dz
ds δ s →0 δ s
∂ψ δ x ∂ψ δ y ∂ψ δ z
= lim + +
δ s →0 ∂x δ s ∂y δ s ∂z δ s
uur
dψ dr ⎡ ∂ψ ∂ψ ∂ψ ⎤
= • ⎢eˆ x + eˆ y + eˆ z ⎥
PCD-09
ds ds ⎣ ∂x ∂y ∂z ⎦ 12
Gradient in the Cartesian Coordinate System
r r r
dψ ψ ( r + δ r ) −ψ ( r )
= lim
ds δ s →0 δs
dψ r
= uˆ • ∇ψ
dψ ∂ψ dx ∂ψ dy ∂ψ dz ds
= + + uur uur
ds ∂x ds ∂y ds ∂z ds δ r dr
uˆ = lim =
r δ s →0 δ s ds
d r =e x dx + e y dy + e z dz
ˆ ˆ ˆ uur
ur ∂ ∂ ∂ δ s = δ r , tiny
∇ =eˆ x + eˆ y + eˆ z
∂x ∂y ∂z
increament
uur
ur ⎡ ∂
∇ψ = ⎢eˆ x + eˆ y
∂ ∂⎤
+ eˆ z ⎥ψ ds = dr , differential
⎣ ∂x ∂y ∂z ⎦
PCD-09 increament
13
r
The GRADIENT ∇ r of a
scalar point function ψ ( r )
yields a vector point function
such that the component of
the resultant vector along dψ r
any direction (given by a unit
= uˆ • ∇ψ
ds
vector û ) gives the uur uur
δ r dr
DIRECTIONAL DERIVATIVE uˆ = lim =
δ s →0 δ s
dψ of the scalar ds
ds function in the
direction of that unit vector.
14
PCD-09
δψ = Cylindrical Polar Coordinate System
ur
∂ψ ∂ψ ∂ψ How should we express ∇ψ
= δρ + δϕ +δ z
∂ρ ∂ϕ ∂z such that :
uur ur
uur ur δψ = δ r • ∇ψ
= δ r • ∇ψ
ur where :
( eˆ ρδρ + eˆϕ ρδϕ + eˆ zδ z ) • ∇ψ r
δ r =eˆ ρ δρ + eˆ ϕ ρδϕ + eˆ zδ z
Following form of the Note how the
gradient operator will work ! 1
ur ∂ 1 ∂ ∂ ρ cancels
∇ =eˆ ρ + eˆ ϕ + eˆ z ρ
∂ρ ρ ∂ϕ ∂z
δψ =
⎛ ∂ 1 ∂ ∂ ⎞
( eˆ ρδρ + eˆϕ ρδϕ + eˆ zδ z ) • ⎜ eˆ ρ ∂ρ + eˆϕ ρ ∂ϕ + eˆ z ∂z 15⎟ψ
PCD-09 ⎝ ⎠
δψ = Spherical Polar Coordinate System
∂ψ ∂ψ ∂ψ ur
=δr + δθ + δϕ How should we express ∇ψ
∂r ∂θ ∂ϕ
uur ur
= δ r • ∇ψ such that :
ur uur ur
⎡⎣eˆ r (δ r ) + eˆ θ ( rδθ ) + eˆ ϕ ( r sin θ δϕ ) ⎤⎦ • ∇ψ δψ = δ r • ∇ψ
where :
Following form of the r
δ r =eˆ r (δ r ) + eˆ θ ( rδθ ) + eˆ ϕ ( r sin θ δϕ )
gradient operator will work !
ur ∂ 1 ∂ 1 ∂
∇ =eˆ r + eˆ θ + eˆ ϕ Note the cancellation of
∂r r ∂θ r sin θ ∂ϕ the factors that are
circled
δψ =
⎛ ∂ 1 ∂ 1 ∂ ⎞
⎡⎣eˆ r (δ r ) + eˆ θ ( rδθ ) + eˆ ϕ ( r sin θ δϕ ) ⎤⎦ • ⎜ eˆ r + eθ
ˆ + eϕ
ˆ ⎟ψ
⎝ ∂r r ∂θ r sin θ ∂ϕ ⎠
16
PCD-09
dψ r
Consolidated expressions for = uˆ • ∇ψ
the GRADIENT ds
uur uur
Cartesian Coordinate System δ r dr
ur u = lim =
∂ ∂ ∂ ˆ
δ s →0 δ s
∇ =eˆ x + eˆ y + eˆ z uur
ds
∂x ∂y ∂z δ s = δ r , tiny
Cylindrical Polar Coordinate System increament
uur
ur ∂ 1 ∂ ∂ ds = dr , differential
∇ =eˆ ρ + eˆ ϕ + eˆ z
∂ρ ρ ∂ϕ ∂z increament
ur ∂ ∂ ∂
∇ =eˆ x + eˆ y + eˆ z
∂x ∂y ∂z
r
∇T = eˆx + eˆy
18
PCD-09
T = 100 + x y
ur ∂ ∂ ∂
∇ =eˆ x + eˆ y + eˆ z
∂x ∂y ∂z
Brighter regions
have
larger values of T
r
∇T = eˆx y + eˆy x
19
PCD-09
T = 100 + x y
r
∇T = xeˆy + yeˆx
20
PCD-09
fields -- gradient
Example:
Brighter regions
have
larger values of T
21
PCD-09
fields -- gradient
Equipotential surface
T = x 2 y 2 = 100
T = x 2 y 2 = 400
Example
Φ(x,y,z) = x2 +y2+z2
23
PCD-09
Example 1 Example 2 : uniform gradient
r
a) T = 20 + x + y ∇T = eˆx + eˆy
r
b) T= 100+xy ∇T = xeˆy + yeˆx
c) T = x2y2 r
∇T = 2 xy ex + 2 x yeˆy
2
ˆ 2
25
PCD-09
T = 15 + ρ cos φ
r
∇T ?
26
PCD-09
Geostationary Orbits ⊂ Geosynchronous Orbits
Geostationary Orbit:
The satellite orbits the earth at exactly the same speed as
the earth turns and at zero latitude (equator). It appears
at the same zenith point from the equatorial point
underneath it. Geosynchronous Orbit:
The satellite’s orbit is synchronized with the earth's
rotation. The plane of the satellite’s orbit may be tilted
with respect to the equatorial plane. The satellite will
then appear to move around along a latitude even if it
will remain at the same longitude.
Many advantages! We need detailed maps of ‘potentials’
and ‘fields’ to launch such objects in the sky. 27
PCD-09
Potentials & their Gradients
Force / Field
http://map.gsfc.nasa.gov/media/990533/index.html
Colonizing space?
At Lagrange points?
Movie available, showing
trajectory 32
PCD-09
L1, L2, L3: Unstable ; L4, L5: Stable
Derivation of the L1, L2, L3 points
Reference:
http://www-istp.gsfc.nasa.gov/stargaze/Slagrang.htm
Derivation of the L4 and L5 points,
based on
"When Trojans and Greeks Collide"
by I. Vorobyov,
"Quantum," p. 16-19, Sept-Oct. 1999.
Uses rotating frames of reference.
Reference:
http://www-istp.gsfc.nasa.gov/stargaze/Slagrng3.htm
33
PCD-09
Interesting Reading on ‘Lagrangian Points’
http://www.aeiveos.com/~bradbury/Authors/Engineering/ONeill-GK/TCoS.html
Colonies in Space
by T. A. Heppenheimer
Reference:
http://www.nss.org/settlement/ColoniesInSpace/colonies_chap02.html
34
PCD-09
October 15, 2009
http://www.youtube.com/watch?v=TSNYu86CqSA
Planck Surveyor – European Space Agency
Launched on 14 May 2009
http://www.rssd.esa.int/index.php?project=PLANCK&page=dev_news
20th October 2009, 05:20
35
PCD-09
http://www.aerith.net/comet/catalog/1998J1/pictures.html
Roughly eighty-five percent of the SOHO discoveries, and also this one,
are fragments from a once great comet that split apart in a death plunge
around the Sun, probably many centuries ago. The fragments are known
as the Kreutz group and now pass within 1.5 million kilometres of the
Sun's surface when they return from deep space.
Downloaded on
20th Oct. 2009
08:00am
Pictures by SOHO@L1
36
PCD-09
http://soho.nascom.nasa.gov/hotshots/2008_06_23/1500thcomet.tif
Downloaded on 20th
Oct. 2009, 08:00am
Pictures by
STEREO
http://stereo.gsfc.nasa.gov/
Other solar observatories:
STEREO (Solar TErrestrial
RElations Observatory) -
NASA's Solar Terrestrial
Probes program - Employs
two nearly identical space-
based observatories - one
ahead of Earth in its orbit, the
other trailing behind
- to provide the first-ever stereoscopic measurements to
study the Sun and the nature of its CME (coronal mass
39
PCD-09 ejections)
Understanding gradient
Electric field and potential inside a conductor
r kp
r Thus E (r , θ , ϕ ) = 3 (eˆr 2 cos θ + eˆθ sin θ )
since p = eˆr p cos θ − eˆθ p sin θ r
r r r
3( p ⋅ eˆr )eˆr − p Ref: Introduction to
E=k 3 Electrodynamics, Griffiths
PCD-09
r
Understanding gradient
Electric field and potential due to an infinitely long wire
r 2k λ
E = eˆρ where λ is the uniform charge density per unit length
ρ
ur ⎡ ∂ 1 ∂ ∂⎤
∇U= ⎢eˆ ρ + eˆ ϕ + eˆ z ⎥ U
⎣ ∂ρ ρ ∂ϕ ∂z ⎦
r r
E = −∇U ; comparing coefficients of eˆρ and eˆϕ
∂U 2k λ ∂U
=− ; =0
∂ρ ρ ∂ϕ
⎧1⎫
hence U = 2k λ ln ⎨ ⎬ + some constant (indep. of ρ )
⎩ρ ⎭
⎧c⎫
i.e. U = 2k λ ln ⎨ ⎬
PCD-09
⎩ρ ⎭ 42
Understanding gradient
Obtaining the field for a given
potential
⎧ ⎛ 2a ρ sin ϕ ⎞ ⎫⎪
⎪π
consider an electrostatic potential U = U 0 ⎨ + tan −1 ⎜ 2 ⎟⎬
⎪⎩ 2 ⎜ ρ − a ⎟⎪2
ur ⎡ ∂ 1 ∂ ∂⎤ ⎝ ⎠⎭
∇U= ⎢eˆ ρ + eˆ ϕ + eˆ z ⎥ U
⎣ ∂ρ ρ ∂ϕ ∂z ⎦
⎧π −1 ⎛ 2 a sin ϕ ⎞ ⎫
U = U 0 ⎨ + tan ⎜ ⎟ ⎬ at ρ >> a
⎩2 ⎝ ρ ⎠⎭
U ( x) = 2 Ax − 3Bx
2 4
PCD-09
44
Slightly magnified view between -2 and +2 on the Y-axis.
Note that
the force
goes to
zero at
both the
points of
‘unstable’
equilibrium,
and also at
the point of
‘stable’
equilibrium
45
PCD-09
Given: Force experienced by a particle is
ur
F ( ρ , ϕ , z ) = −eˆρ ρ cos 2ϕ + eˆϕ ρ sin 2ϕ + eˆz z
ur r ∂U 1 ∂U ∂U
F = −∇ U = − eˆρ − eˆϕ − eˆz Obtain the potential
∂ρ ρ ∂ϕ ∂z for the given field.
∂U ρ2
i.e. = ρ cos 2ϕ ⇒ U= cos 2ϕ + f (ϕ ,z )
∂ρ 2
∂U ∂f
This implies = − ρ sin 2ϕ +
2
.
∂ϕ ∂ϕ
ur 1 ∂U
The azimuthal component of − ∇U is − = ρ sin 2ϕ ,
ρ ∂ϕ
∂f
hence : = 0 which gives f = f ( z ) + c
∂ϕ
∂U ∂f ∂U −z2
Now = ; and = −z ⇒ f =
∂z ∂z ∂z 2
ρ2 z2
∴ U (ρ ,ϕ , z) = cos 2ϕ − + k
PCD-09 2 2 46
Next Unit is ‘Unit 9’: Gauss’ Law;
Equation of Continuity. Hydrodynamic and
Electrodynamic illustrations.
You will find that the gradient operator has very many
important applications in physics.
http://en.wikipedia.org/wiki/Image:Carl_Friedrich_Gauss.jpg
Learning goals:
The divergence theorem is a mathematical expression of
a conservation principle. In the absence of any source or
sink, the density of matter in a well-defined volume
element can change if and only of matter flows in, or out,
of that region across the surface that bounds that
volume region. 1
PCD-09
The result is equally consequential
with regard to fields just as well for
matter.
ur ur r ur r uur
∫∫∫
volume
dτ ⎡⎣∇ • A(r ) ⎤⎦ = ∫∫
surface
A(r ) • da
region enclosing
that
region 2
PCD-09
Field Lines
r r 1 q
E (r ) = eˆ − − − −(1)
4πε 0 r 2 r
ur r uur
Flux crossing a surface = ∫∫
surface
A(r ) • da
ur r
A(r ) : A vector
r
nˆ (r ) : unit vector point function
r
normal to the surface at the point r
6
PCD-09
Flux crossing a surface
ur r uur Z Y
= ∫∫ A(r ) • da
surface
⎡ ∂Az ⎤ ⎡ ∂Az ⎤
=⎢ δ zδ xδ y =⎢ ⎥ δV
⎥
⎣ ∂z ⎦ ( x0, y0 , z0 ) ⎣ ∂z ⎦ ( x0, y0 , z0 )
Flux crossing all
the six surface elements Adding the flux through all faces, total flux
that enclose the cell r r r ⎡ ∂Ax ∂Ay ∂Az ⎤
∫ ∫ ⎢⎣ ∂x + ∂y + ∂z ⎥⎦ dV
ur r uur
= ∫∫ A(r ) • da A(r ).dS =
surface
Cube 7
PCD-09
Flux crossing a surface
ur r uur Adding the flux through Z Y
= ∫∫ A(r ) • da all faces, total flux
surface
ur ur r ur r uur
Take the limit of
∫∫∫
volume
dτ ⎡⎣∇ • A(r ) ⎤⎦ = ∫∫
surface
A ( r ) • d a
the ratio of
region enclosing
that total flux over δs
region
to δV
ur r
∫∫ A(r) • ndS
ˆ
flux per unit volume,
r enclosing r r
divA = lim = ∇• A
surface
δτ →0 δτ at that point
remember: flux is defined through a SURFACE,
whereas divergence is defined at a POINT
Flux is a scalar quantity. It is not a scalar field; it is not a
local quantity – not a ‘point function’.
Divergence is a scalar field; it is a scalar point function,
defined at each point of space 9
PCD-09
Gauss’s Divergence Theorem
If a volume V is bounded by a surface S, then, for vector A,
Physical Meaning:
Integration of the faucets
(source of vector field) over
a volume is equal to the
flux flowing out through the
surface enclosing the
volume. 11
PCD-09
ur ur r ur r uur r r r r
∫∫∫
volume
dτ ⎡⎣∇ • A(r ) ⎤⎦ = ∫∫
surface
A(r ) • da
∫ divA(r )dV =
V
∫ A(r ).nds
S
ˆ
region enclosing
that
region
r r q
E (r ) = eˆr
4πε 0 r 2
12
PCD-09
ur ur r ur r uur
∫∫∫
volume
dτ ⎡⎣∇ • A(r ) ⎤⎦ = ∫∫
surface
A(r ) • da
the electric flux
region enclosing
that
through a spherical
region
surface is given as
uur
dS r = r sin θ dθ dϕ eˆr
2
π 2π
r r ⎛ q ⎞
∫ E (r ) • nds
ˆ = ∫θ φ∫ ⎜ 2 r ⎟
eˆ • r 2
( )
sin θ dθ dφ eˆr
s =0 =0 ⎝ 4πε 0 r ⎠
Solving the integral we get
r r q
∫
Gauss’s law in electrostatics
E ( r ). ˆ
nds = ( in the integral form)
S ε0
The total flux emanating from the closed surface is independent
of radius r, in fact it is also independent of the shape of the
surface chosen!
Exercise: Find the flux of electric field due to a
point charge through a cube of length L. 13
PCD-09
if we have charge distribution instead of a point charge,
r r q 1
∫∫ E (r ).nds
S
ˆ = =
ε0 ε0 ∫∫∫
volume
ρ dV
of the charge
distribution
r r
∇• A = ⎡ ∂ ⎤
⎢eˆ ρ ∂ρ ⎥ • ⎡⎣eˆ ρ Aρ ( ρ , ϕ , z ) + eˆ ϕ Aϕ ( ρ , ϕ , z ) + eˆ z Az ( ρ , ϕ , z ) ⎤⎦ +
⎣ ⎦
⎡ 1 ∂ ⎤
⎢ eˆ ϕ ρ ∂ϕ ⎥ • ⎡⎣eˆ ρ Aρ ( ρ , ϕ , z ) + eˆ ϕ Aϕ ( ρ , ϕ , z ) + eˆ z Az ( ρ , ϕ , z ) ⎤⎦ +
⎣ ⎦
⎡ˆ ∂ ⎤ ˆ
⎢⎣e z ∂z ⎥⎦ • ⎡⎣e ρ Aρ ( ρ , ϕ , z ) + eˆ ϕ Aϕ ( ρ , ϕ , z ) + eˆ z Az ( ρ , ϕ , z ) ⎤⎦
15
PCD-09
r r
∇• A = ⎡ ∂ ⎤
⎢eˆ ρ ∂ρ ⎥ • ⎡⎣eˆ ρ Aρ ( ρ , ϕ , z ) + eˆ ϕ Aϕ ( ρ , ϕ , z ) + eˆ z Az ( ρ , ϕ , z ) ⎤⎦ +
⎣ ⎦
⎡ 1 ∂ ⎤
⎢ eˆ ϕ ρ ∂ϕ ⎥ • ⎡⎣eˆ ρ Aρ ( ρ , ϕ , z ) + eˆ ϕ Aϕ ( ρ , ϕ , z ) + eˆ z Az ( ρ , ϕ , z ) ⎤⎦ +
⎣ ⎦
⎡ˆ ∂ ⎤ ˆ
⎢⎣e z ∂z ⎥⎦ • ⎡⎣e ρ Aρ ( ρ , ϕ , z ) + eˆ ϕ Aϕ ( ρ , ϕ , z ) + eˆ z Az ( ρ , ϕ , z ) ⎤⎦
r r ⎧∂ ⎫
∇ • A = e ρ • ⎨ ∂ρ ⎬ ⎡⎣eˆ ρ Aρ ( ρ , ϕ , z ) + eˆ ϕ Aϕ ( ρ , ϕ , z ) + eˆ z Az ( ρ , ϕ , z ) ⎤⎦ +
ˆ
⎩ ⎭
⎧1 ∂ ⎫
ˆeϕ • ⎨ ⎬ ⎡⎣eˆ ρ Aρ ( ρ , ϕ , z ) + eˆ ϕ Aϕ ( ρ , ϕ , z ) + eˆ z Az ( ρ , ϕ , z ) ⎤⎦ +
⎩ ρ ∂ϕ ⎭
⎧ ∂⎫
eˆ z • ⎨ ⎬ ⎡⎣eˆ ρ Aρ ( ρ , ϕ , z ) + eˆ ϕ Aϕ ( ρ , ϕ , z ) + eˆ z Az ( ρ , ϕ , z ) ⎤⎦16
PCD-09 ⎩ ∂z ⎭
r r ⎧∂ ⎫
∇ • A = eˆ ρ • ⎨ ⎬ ⎡⎣ eˆ ρ Aρ ( ρ , ϕ , z ) + eˆ ϕ Aϕ ( ρ , ϕ , z ) + eˆ z Az ( ρ , ϕ , z ) ⎤⎦ +
∂ρ⎩ ⎭
⎧1 ∂ ⎫
eˆ ϕ • ⎨ ⎬ ⎡⎣eˆ ρ Aρ ( ρ , ϕ , z ) + eˆ ϕ Aϕ ( ρ , ϕ , z ) + eˆ z Az ( ρ , ϕ , z ) ⎤⎦ +
⎩ ρ ∂ϕ ⎭
⎧ ∂⎫
eˆ z • ⎨ ⎬ ⎡⎣eˆ ρ Aρ ( ρ , ϕ , z ) + eˆ ϕ Aϕ ( ρ , ϕ , z ) + eˆ z Az ( ρ , ϕ , z ) ⎤⎦
⎩ ∂z ⎭
∂eˆρ ∂eˆρ
= 0, = eˆϕ ,
∂ρ ∂ϕ
∂eˆϕ ∂eˆϕ
= 0, = −eˆρ
∂ρ ∂ϕ
r r ∂ 1
∇• A = Aρ ( ρ , ϕ , z ) + Aρ ( ρ , ϕ , z ) +
∂ρ ρ
1 ∂ ∂
+ Aϕ ( ρ , ϕ , z ) + Az ( ρ , ϕ17 , z )
PCD-09 ρ ∂ϕ ∂z
Expression for ‘divergence’ in spherical polar coordinate system
∂eˆθ r ∂eˆϕ r
∂eˆr r =0 =0
=0 ∂r
∂r ∂r
∂eˆr ∂eˆθ ∂eˆϕ r
= eˆθ = −eˆr =0
∂θ ∂θ ∂θ
∂eˆr ∂eˆθ ∂eˆϕ
= sin θ eˆϕ
∂ϕ = cos θ eˆϕ = − cos θ eˆθ − sin θ eˆr
∂ϕ ∂ϕ
r r
∇• A =
⎧ ∂ 1 ∂ 1 ∂ ⎫
⎨eˆ r + eˆθ + eˆϕ ⎬ • ⎡⎣eˆ r Ar (r ,θ , ϕ ) + eˆθ Aθ (r ,θ , ϕ ) + eˆϕ Aϕ (r,θ , ϕ )⎤⎦
⎩ ∂r r ∂θ r sin θ ∂ϕ ⎭
=
1 ∂ 2 1 ∂ 1 ∂
⎡
⎣ r Ar (r,θ , ϕ )⎦ +
⎤ [ Aθ (r,θ ,ϕ )sin θ ] + Aϕ (r,θ ,ϕ )
r ∂r
2
r sin θ ∂θ r sin θ ∂ϕ
18
PCD-09
Examples for solenoidal and nonsolenoidal fields
r ⎡ ∂Ax ∂Ay ∂A ⎤ r r
div A = ⎢ + + z ⎥ = ∇. A
⎣ ∂x ∂y ∂z ⎦
r
r r
∇.A = 0
A = (x -y)eˆ x + ( x + y ) eˆ y
r r
Influx balances the outflux ∇ .A = 2
Solenoidal Example: B
The divergence of a vector field is the extent to which the flux lines
behave like being near a source or a sink at a given point.
19
PCD-09
Flux crossing a surface
ur r uur Is there any net accumulation of
= ∫∫ A(r ) • da
surface
the flux in a volume element?
Source
Sources and
Sinks may be
present in the
region !
Sink
r r
∫∫ J (r ).nds
S
ˆ = I , the current through the surface element.
23
PCD-09
Divergence theorem: Conservation principle
r
ρ
Conservation of mass or charge r , t ) represents mass/charge density
( r
r
J ( r , t ) : mass/charge current density
What shall we get if we integrate the flux
emanating from all the six enclosing surfaces?
ur ur r ur r uur
r r ∫∫∫ dτ ⎡⎣∇ • J (r ) ⎤⎦ = ∫∫ J (r ) • da
∫∫ J (r ).nds
S
ˆ = I. volume
region
surface
enclosing
that
region
r r ∂qtotal ∂ ⎧ ∂ρ ⎫
i.e. ∫∫ J (r ).nds
ˆ =− = − ∫∫∫ ρ dV = ∫∫∫ ⎨- ⎬ dV
S
∂t ∂t V V ⎩
∂t ⎭
Negative sign:
Outward flux is at the expense of the charge inside!
Integral and
⎧ur ur r ∂ρ ⎫
∫∫∫ dV ⎨∇• J (r ) + ⎬ = 0
⎩ ∂t ⎭
Differential forms of
volume
region the equation of
continuity:
⎧ ur ur r ∂ρ ⎫
⎨∇ • J (r ) + ⎬ = 0 conservation
⎩ ∂t ⎭ principle 27
PCD-09
Divergence theorem: Conservation principle Equation of Continuity
ur ur r ur r uur
∫∫∫ {
volume
} ∫∫
dτ ∇• J (r ) =
surface
J (r) • da In the absence of any
region enclosing source or sink, the
that
region density of matter in a
well-defined volume
ur ur r ∂ρ
∇ • J (r ) = − element can change if
∂t and only of matter flows
ur ur r ∂ρ
∇ • J (r ) + =0 in, or out, of that region
∂t across the surface that
bounds that volume
region.
⎧ur ur r ∂ρ ⎫
∫∫∫
volume
dV ⎨∇• J (r ) + ⎬ = 0
⎩ ∂t ⎭
region 28
PCD-09
FLUID MECHANICS
We consider an incompressible fluid.
A fluid ‘flows’
Deformation of solids, fluid flow: “rheology”
Non-Newtonian fluids --- example paints, foams, molten
plastics…..
30
PCD-09
The pressure applied to an enclosed fluid is transmitted undiminished to
every portion of the fluid and to the walls of the container.
The pressure is the same in every direction. The shape of the container
does not matter.
S does not depend on the direction of uˆ N .
Pascal ' s law.
Stress at the point P is S .
δA
P uˆN =
δA
Sindhu Indus
Vitastha Jhelum
Asikini Chenab
Airavati Ravi
Satadru Sutlej
Rivers
in
Jammu
and
FLUID DYNAMICS Kashmir36
PCD-09
r rr
How do we track density ρ (r) & velocity v(r) at a point in a river?
Alakananda
Bhagirathi
At 2700 feet, Deoprayag
("Divine confluence") is the
meeting place for the Bhagirathi
River (left) and the Alakananda
River, the two major Himalayan
tributaries of the Ganges.
How do we have
TRIWENI
SANGAM at
Prayag?
http://www.google.co.in/imgres?imgurl=http://personal.carthage.edu/jlochtefeld/picturepages/pilgrimage/deoprayag05.jpg&im
grefurl=http://personal.carthage.edu/jlochtefeld/picturepages/pilgdeoprayag.html&h=481&w=700&sz=87&tbnid=ZQt0BRpsi0E
J::&tbnh=96&tbnw=140&prev=/images%3Fq%3Dbhagirathi%2Briver%2Bpicture&hl=en&usg=__899jmh0E8OAKvcBWqaIJPd 37
ON_k8=&sa=X&oi=image_result&resnum=1&ct=image&cd=1
PCD-09
r rr
How do we track density ρ (r) & velocity v(r) at a point in a river?
Bhagirathi Alakananda
Jamuna
Saraswathi
Is merely the position
How did the waters Ganga vector of a particular
point in space, or is it
of SARASWATHI the position vector of a
ever get to the Prayag moving molecule of
TRIWENI SANGAM? Triweni water that has come
Sangam from either Bhagirathi
38 or
PCD-09 from Alakananda?
r rr
How do we track density ρ (r) & velocity v(r) at a point in a river?
ONLY A ‘LAGRANGIAN’
TRACKING OF A WATER
MOLECULE WOULD REVEAL
HOW WATERS OF
SARASWATHI MADE IT TO
PRAYAG! The waters of
Yamuna would mix with the
waters of Saraswathi and bring
them to Prayag!
http://www.google.co.in/imgres?imgurl=http://personal.carthage.edu/jlochtefeld/picturepages/pilgrimage/deoprayag05.jpg&im
grefurl=http://personal.carthage.edu/jlochtefeld/picturepages/pilgdeoprayag.html&h=481&w=700&sz=87&tbnid=ZQt0BRpsi0E
J::&tbnh=96&tbnw=140&prev=/images%3Fq%3Dbhagirathi%2Briver%2Bpicture&hl=en&usg=__899jmh0E8OAKvcBWqaIJPd 39
ON_k8=&sa=X&oi=image_result&resnum=1&ct=image&cd=1
PCD-09
Equation of motion for fluids
two basic approaches
Lagrangian Approach:
Follow the motion of some
particle of the fluid; this must be
done for all particles of the fluid
Joseph-Louis Lagrange
1736 - 1813
Eulerian Approach:
⎡ ∂J x ⎤
⎢ ∂x ⎥ δxδyδz
Amount of mass of fluid ⎣ ⎦P
crossing face ABCD in unit
time
⎧ ⎡ ∂J x ⎤ ⎛ δx ⎞ ⎫ ⎡ ∂J x ⎤
= ⎨ J x (r ) + ⎢ ⎜ ⎟ ⎬δ y δ z =⎢ ⎥ δV
⎥ ⎣ ∂ x ⎦ P 42
⎩ ⎣ ∂x ⎦ P ⎝ 2 ⎠ ⎭
PCD-09
Net OUTWARD flux
through the two faces
orthogonal to x-axis =
G C
⎡ ∂J x ⎤ H D
⎢ ∂x ⎥ δxδyδz P
⎣ ⎦P F
B
⎡ ∂J x ⎤ E
r
=⎢ ⎥ δV A
⎣ ∂x ⎦ P êx
This quantity
Net OUTWARD flux through the whole
must be
parallelepiped, per unit volume =
equal to
⎧⎪ ⎡ ∂ J x ⎤
⎨⎢ ⎥
⎡ ∂J y ⎤
+⎢
⎡ ∂ J z ⎤ ⎫⎪
⎥ +⎢ ⎥ ⎬= ∇•J
⎪⎩ ⎣ ∂ x ⎦ P ⎣ ∂ y ⎦ P ⎣ ∂ z ⎦ P ⎪⎭
[ ]
P −
∂ρ
∂t
The choice of the term ‘DIVERGENCE’ ∂ρ
∇ • J = − 43
PCD-09
is thus well justified. ∂t
Equation of Continuity
Conservation of Matter
∂ρ (r , t )
∇ • J (r , t ) = −
∂t
J (r , t ) • uˆ = mass flux
in the direction
of uˆ
∂ρ (r , t )
∇ • J (r , t ) + =0
∂t
44
PCD-09
ê y Face ‘1’
of the Face ‘2’
parallelepiped of the
ê x region parallelepiped
region
ê z
Ideal fluid: stress at any point is
essentially one of COMPRESSION.
45
PCD-09
⎧ ⎡ ∂p ⎤ δx ⎫
F (1) = ⎨ p ( r ) − ⎢ ⎥ ⎬(δyδz )(+ eˆx ) → 1 2 ←
⎩ ⎣ ∂x ⎦ P 2 ⎭ ê x
⎧ ⎡ ∂p ⎤ δx ⎫
⎬(δyδz )(− eˆx )
Ideal fluid: stress
F ( 2) = ⎨ p ( r ) + ⎢ ⎥ at any point is
⎩ ⎣ ∂x ⎦ P 2 ⎭ essentially one of
COMPRESSION.
⎡ ∂p ⎤ ⎡ ∂p ⎤
F (1) + F ( 2) = − ⎢ ⎥ δx (δyδz )(+ eˆx ) = − ⎢ ⎥ δV eˆx
⎣ ∂x ⎦ P ⎣ ∂x ⎦ P
6
∑ F (i) = −∇ p δV
i =1
HYDROSTATIC force. Of course, there
may be some additional external force
acting, such as gravity.
Net force per unit volum e = − ∇ p = negative
gradient of
HYDROSTATIC force
pressure 46
PCD-09
Net HYDROSTATIC force acting on
the parallelepiped per unit volume = − ∇p → ←
External force (such as gravity) Total {hydrostatic +
acting on the parallelepiped external (gravity)} force
per unit volume acting on the parallelepiped
F external per unit volume
= lim
δV → 0 δV δm d v
lim
F external δm δV → 0 δ V dt
= lim
δV → 0 δm δV dv
ρ (r ) = −∇ p + g ρ ( r )
F external dt
= lim ρ (r )
δV → 0 δm Mass x Acceleration
“Cause-Effect”
= g ρ (r ) Newton’s law:
Equation of Motion 47
PCD-09
Mass x Acceleration / “Cause-Effect” δm d v
lim
Newton’s law: Equation of Motion δV → 0 δ V dt
dv
ρ (r ) = −∇ p + g ρ ( r )
dt
The position vector r in this equation must represent the
position of a moving/flowing material particle of the fluid, just
as it is used the Newton’s law. It is the ‘LAGRANGIAN’
position vector of a moving/flowing fluid ‘particle/molecule’, not
the EULERIAN position vector of a fixed point in space.
r Lagrangian = r (t ) This is a function of time
⎡ ∂⎤ d v(r , t ) − ∇ p
⎢⎣ v • ∇ + ∂t ⎥⎦ v ( r , t ) = = − ∇φ + F viscous
dt ρ (r )
Hydrodynamic Viscous, frictional,
term dissipative term.
This terms makes
“dry water wet”
Next: Unit 10 - Feynman 50
PCD-09
Unit 10: Fluid Flow, Bernoulli’s Principle
Steady flow.
Bernoulli’s equation/principle, some illustrations.
1
PCD-09
Unit 10: Fluid Flow, Bernoulli’s Principle.
Equation of motion for fluid flow. Definition of
curl, vorticity, Irrotational flow and circulation.
Steady flow. Bernoulli’s principle, some
illustrations. Introduction to applications of Gauss’
law and Stokes’ theorem in Electrodynamics.
Learning goals:
Learn that both the divergence and the curl of a
vector field are involved (along with the boundary
conditions) in determining its properties. Learn
how a rigorous treatment of the velocity field is
necessary to explain quantitatively the observed
phenomena in fluid dynamics. Get ready for a
theory of electrodynamics. 2
PCD-09
Recall the discussion on directional derivative
dψ r
= uˆ • ∇ψ
ds Gradient: direction in
uur uur which the function varies
δ r dr fastest / most rapidly.
uˆ = lim = ur r
δ s →0 δ s ds F = −∇ψ
uur
δ s = δ r , tiny Force: Negative gradient of
the potential
increament
uur ‘negative’ sign is the result
ds = dr , differential of our choice of natural
motion as one occurring
increament from a point of ‘higher’
potential to one at a ‘lower’
potential. 3
PCD-09
ur r ur r
Is it obvious that the ‘force’ defined by these F = −∇ψ
F = ma two equations is essentially the same?
Consistency in
Is it obvious that the ‘force’ defined by
these relations
these two equations is essentially the exists only for
same? ‘conservative’
forces.
ur uur
The compatibility of the two expressions ∫ F • dr = 0
emerges if, and only if, b ur uur
the potential ψ is defined in such way ∫ F • dr is
a
that the work done by the force given by INDEPENDENT
r
−∇ψ in displacing the object on which this
of the path
force acts, is independent of the path
a to b
along which the displacement occurred. PATH INTEGRAL4
PCD-09 “CIRCULATION”
ur uur
∫ F • dr = 0 It is only when the line integral of the work done is
path independent that the force is conservative
bur uur
∫ F • dr is and accounts for the acceleration it generates
a when it acts on a particle of mass m through the
INDEPENDENT causality relation of ur r
of the path Newtonian mechanics:
F = ma
a to b
ur
This alternative expression employsurwhaturis known as CURL of a
VECTOR FIELD F , denoted as ∇× F .
5
PCD-09
Definition of Curl of a vector:
ur ur
ur r ∇ × F is a vector point function
of the vector field F ( r ) such that for an orthonormal basis set
r
of unit vectors
{uˆ (r ), i = 1, 2,3} ,
i
Mean (average) ‘circulation’
per unit area taken at the
ur r uur
r ur ur r
uˆi (r ) • ∇× F (r ) = lim
∫ F (r ) • dr ,
point when the elemental
area becomes infinitesimally
Δs
ΔS →0
small.
where the path integral is taken over a closed path C, taken
over a tiny closed loop C which bounds an elemental vector
uuur r
surface area ΔS = ΔS uˆi (r ).
The direction of the unit vector is to be taken such that
a right-hand screw would propogate forward when turned
along the sense in which the path integral is determined.
6
PCD-09
ur r uur In general, the unit
r ur ur r
uˆi (r ) • ∇× F (r ) = lim
∫ F (r ) • dr vectors may depend
ΔS →0 Δs
r on the particular point
{ }
uˆi (r ), i = 1, 2,3 under discussion, and
Cartesian unit vectors, of course, hence written as
do not change from point to point. r r
functions uˆi (r ) of r.
They are constant vectors.
9
PCD-09
curl of a vector field at a point represents the net circulation of
the field around that point.
ur r r
If ∇×F=0 in a region then there would be no
curliness/rotation, and the field is called irrotational.
10
PCD-09
r r r
r ∫ A ( r ).dl
r r
( curl A ).nˆ = lim C
= ( ∇ × A ).nˆ
δ S→0 δS
Remember!
The criterion that a force field is conservative is that its path
integral over a closed loop (i.e. “circulation”) is zero. This is
ur r r
equivalent to the condition that ∇×F=0
ur r r
If ∇×F=0 in a region, then there would be no curliness
(rotation), and the field is called irrotational.
r r
V ( x, y ) = − yeˆx + xeˆ y A(x,y)=(x-y)eˆ x + ( x + y )eˆy
r r r r
∇ × V = 2eˆz ∇ × A = 2eˆ z
x
v = xeˆ y
r
∇ × v = êz
16
PCD-09
What is the DIVERGENCE and the CURL of the following
vector field?
eˆ x eˆ y eˆ z
∂ ∂ ∂
∇ × ∇φ = The final result will be
∂x ∂y ∂z independent of the
∂φ ∂φ ∂φ coordinate system.
∂x ∂y ∂z
⎛d ⎞ r ⎛d ⎞ r ur r
⎜ ⎟ r =⎜ ⎟ r + ω x r
⎝ dt ⎠ I ⎝ dt ⎠ R
⎛d ⎞ r r ⎛ d ⎞ r ur r
When ⎜ ⎟ r =0, ⎜ ⎟ r =ω x r
⎝ dt ⎠ R ⎝ dt ⎠ I
⎛d ⎞ r r ur r
⎜ ⎟ r = vI = ω x r
⎝ dt ⎠ I
24
PCD-09
⎛d ⎞ r r ur r r r r
⎜ ⎟ r = vI = ω x r
⎝ dt ⎠ I
v =ω×r
r r r ur r eˆx eˆy eˆz
r
∇ × v = ∇ × (ω × r ) = ∇ × ωx ω y ωz
x y z
r
= ∇ × ⎡⎣(ω y z − ωz y )iˆ + (ωz x − ωx z )eˆy + (ωx y − ω y x)eˆz ⎤⎦
eˆx eˆy eˆz The ‘curl’ of
the linear
∂ ∂ ∂
= velocity gives
∂x ∂y ∂z a measure of
(ω y z − ωz y ) (ω z x − ω x z ) (ω x y − ω y x) (twice) the
angular
r r r velocity; thus
∇ × v = 2(ωx eˆx + ω y eˆy + ωz eˆz ) = 2ω justifying the
term ‘curl’. 25
PCD-09
Remember:
ur r uur
r ur ur r
uˆi (r ) • ∇× F (r ) = lim
∫ F (r ) • dr
ΔS →0 Δs
∫ A ( r ) • dl = ∫∫ ( ∇ × A ) • dS (1824-1907)
Note!
It is STOKES’
THEOREM
not STOKE’S THEOREM
0
K temperature
This theorem is named after George Gabriel Stokes
(1819–1903), although the first known statement of
George Gabriel the theorem is by William Thomson (Lord Kelvin) and
Stokes
(1819–1903) appears in a letter of his to Stokes in July 1850.
27
PCD-09 Reference: http://www.123exp-math.com/t/01704066342/
Proof of Stokes’ theorem follows from the very definition of the curl:
r r r
r ∫ A ( r ) • dl
r r
D efinition : ( curl A ) • nˆ = lim C
= ( ∇ × A ) • nˆ
δ S→0 δS
For a path δC, which binds an area δS,
r r r
∫δ
C
A • dl = δ S (curl A) • nˆ
r r
∇× A =
⎡ ∂ 1 ∂ ∂⎤
⎢eˆ ρ ∂ρ + eˆ ϕ ρ ∂ϕ + eˆ z ∂z ⎥ × ⎡⎣ eˆ ρ Aρ ( ρ , ϕ , z ) + eˆ ϕ Aϕ ( ρ , ϕ , z ) + eˆ z Az ( ρ , ϕ , z ) ⎤⎦
⎣ ⎦
r r
∇ × A = ⎡⎢eˆ ρ ∂ ⎤⎥ × ⎡eˆ ρ Aρ ( ρ , ϕ , z ) + eˆ ϕ Aϕ ( ρ , ϕ , z ) + eˆ z Az ( ρ , ϕ , z ) ⎤ +
∂ρ ⎣ ⎦
⎣ ⎦
⎡ 1 ∂ ⎤
⎢ eˆ ϕ ρ ∂ϕ ⎥ × ⎡⎣ eˆ ρ Aρ ( ρ , ϕ , z ) + eˆ ϕ Aϕ ( ρ , ϕ , z ) + eˆ z Az ( ρ , ϕ , z ) ⎤⎦ +
⎣ ⎦
⎡ˆ ∂ ⎤ ˆ
⎢⎣e z ∂z ⎥⎦ × ⎡⎣e ρ Aρ ( ρ , ϕ , z ) + eˆ ϕ Aϕ ( ρ , ϕ , z ) + eˆ z Az ( ρ , ϕ , z ) ⎤⎦
33
PCD-09
r r
∇× A =
⎡ ∂
ˆ
e + ˆ
1 ∂
⎢ ρ ∂ρ ϕ ρ ∂ϕ
e + ˆ
e z
∂⎤
∂ ⎥ × ⎡⎣eˆ ρ Aρ ( ρ , ϕ , z ) + eˆ ϕ Aϕ ( ρ , ϕ , z ) + eˆ z Az ( ρ , ϕ , z ) ⎤⎦
{eˆ , eˆ , eˆ }
ρ ϕ z
⎣ z ⎦
⎡ ∂ ⎤
r r ⎢ ρ ∂ρ ⎥ × ⎡⎣eˆ ρ Aρ ( ρ , ϕ , z ) + eˆ ϕ Aϕ ( ρ , ϕ , z ) + eˆ z Az ( ρ , ϕ , z ) ⎤⎦ +
ˆ
e
∇× A = ⎣ ⎦
⎡ 1 ∂ ⎤
⎢ eˆ ϕ ρ ∂ϕ ⎥ × ⎡⎣ eˆ ρ Aρ ( ρ , ϕ , z ) + eˆ ϕ Aϕ ( ρ , ϕ , z ) + eˆ z Az ( ρ , ϕ , z ) ⎤⎦ +
⎣ ⎦
⎡ˆ ∂ ⎤ ˆ
⎢e z ∂z ⎥ × ⎡⎣e ρ Aρ ( ρ , ϕ , z ) + eˆ ϕ Aϕ ( ρ , ϕ , z ) + eˆ z Az ( ρ , ϕ , z ) ⎤⎦
⎣ ⎦
r r ∂
∇ × A = ⎡⎣ eˆ ρ ⎤⎦ × ⎡⎣ eˆ ρ Aρ ( ρ , ϕ , z ) + eˆ ϕ Aϕ ( ρ , ϕ , z ) + eˆ z Az ( ρ , ϕ , z ) ⎤⎦ +
∂ρ
⎛1 ∂ ⎞
⎡⎣ eϕ ⎤⎦ × ⎜
ˆ ⎟ ⎡⎣ eˆ ρ Aρ ( ρ , ϕ , z ) + eˆ ϕ Aϕ ( ρ , ϕ , z ) + eˆ z Az ( ρ , ϕ , z ) ⎤⎦ +
⎝ ρ ∂ϕ ⎠
∂
[eˆ z ] × ⎡⎣eˆ ρ Aρ ( ρ , ϕ , z ) + eˆ ϕ Aϕ ( ρ , ϕ , z ) + eˆ z Az ( ρ , ϕ , z ) ⎤⎦ 34
PCD-09 ∂z
Expression for ‘curl’ in cylindrical polar coordinate system
{eˆ , eˆ , eˆ }
ρ ϕ z
r r ∂
⎡ ⎤
∇ × A = ⎣e ρ ⎦ ×
ˆ ⎡ eˆ ρ Aρ ( ρ , ϕ , z ) + eˆ ϕ Aϕ ( ρ , ϕ , z ) + eˆ z Az ( ρ , ϕ , z ) ⎦⎤ +
⎣
∂ρ
⎛1 ∂ ⎞
⎣⎡ eˆ ϕ ⎦⎤ × ⎜ ρ ∂ϕ ⎟ ⎣⎡ eˆ ρ Aρ ( ρ , ϕ , z ) + eˆ ϕ Aϕ ( ρ , ϕ , z ) + eˆ z Az ( ρ , ϕ , z ) ⎦⎤ +
⎝ ⎠
∂
[e z ] × ⎡⎣eˆ ρ Aρ ( ρ , ϕ , z ) + eˆ ϕ Aϕ ( ρ , ϕ , z ) + eˆ z Az ( ρ , ϕ , z ) ⎤⎦
ˆ
∂z
r r
∇ × A=
⎛ ∂ 1 ∂ ∂ ⎞
⎟ × ( eˆ r Ar (r , θ , ϕ ) + eˆ θ Aθ (r ,θ , ϕ ) + eˆ ϕ Aϕ (r ,θ , ϕ ) )
1
ˆ
= ⎜ er + eθ
ˆ + eϕ
ˆ
⎝ ∂r r ∂θ r sin θ ∂ϕ ⎠
r r ⎧ ∂ ∂Aθ ⎫
⎨ ( sin θ Aϕ ) −
1
∇ × A=eˆ r ⎬
r sin θ ⎩ ∂θ ∂ϕ ⎭
1 ⎧ 1 ∂Ar ∂ ⎫
+eˆ θ ⎨ − ( rAϕ ) ⎬
r ⎩ sin θ ∂ϕ ∂r ⎭
1⎧ ∂ ∂Ar ⎫
+ eˆ ϕ ⎨ ( rAθ ) − ⎬
r ⎩ ∂r ∂θ ⎭ 36
PCD-09
Do you like the ear-ring? The Möbius strip used to be
common in belt drives (like a
car fan belt). With an ordinary
belt only the inside of the belt
was in contact with the
wheels, so it would wear out
before the outside did. Since
a Möbius strip has only one
side, the wear and tear on
It is not
the at
beltall
wasobvious….
spread out more
evenly and they would last
longer. However, modern
It is MOBIUS !!! belts are made from several
layers of different materials,
with a definite inside and
outside, and do not have a
twist. 37
PCD-09 http://mathssquad.questacon.edu.au/mobius_strip.html
Maurits
Cornelius
M.C. Escher
Escher
1898 - 1972
http://www.mcescher.com/
“The laws of mathematics are not merely human inventions or creations. They
simply 'are'; they exist quite independently of the human intellect. The most that
any(one) ... can do is to find that they are there and to take cognizance of them. ”
38
PCD-09
An important identity: divergence of a curl is zero
r r r r r r
By
Gauss’ divergence theorem
∫ ∇ ⋅ (∇ × A)dτ = ∫ (∇ × A) ⋅ dS
V S
S2
Volume
S1 c2
Surface enclosing a volume c1
r r r r
= ∫ A ⋅ dl + ∫ A ⋅ dl = 0 r r r
∇ ⋅ (∇ × A) = 0
C1 C2 39
PCD-09
The curl of a vector is an important quantity.
40
PCD-09
2 sons of Posed the
Nicolaus Bernoulli brachistrocrone
Bernoulli’s brothers problem
Bernoulli Bernoulli’s
Family Principle
Math/Phys Tree
http://www.york.ac.uk/depts/maths/histstat/people/bernoulli_tree.htm
http://library.thinkquest.org/22584/temh3007.htm
Daniel Bernoulli
1700 - 1782
http://www-groups.dcs.st-and.ac.uk/~history/PictDisplay/Bernoulli_Daniel.html 42
PCD-09
Result of Page ⎡ ∂⎤ d v(r , t ) − ∇ p
50, Unit 9
⎢⎣ v • ∇ + ∂t ⎥⎦ v ( r , t ) = = − ∇φ
dt ρ (r )
( )
Use now the
following ∇ A• B =
( ) ( )A + A × (∇ × B )+ B × (∇ × A )
vector
identity: A•∇ B + B •∇
∇ (v • v ) =
(v • ∇ )v + (v • ∇ )v + v × (∇ × v )+ v × (∇ × v )
i.e. ∇ (v • v ) = (v • ∇ )v + v × (∇ × v )
1
2
1
2
( )
∇ v• v − v× ∇× v + (
∂v − ∇ p
∂t
= )
ρ (r )
− ∇φ
43
PCD-09
1
2
( )
∇ v• v − v× ∇× v + (
∂v − ∇ p
∂t
=
ρ (r )
)− ∇φ
∇ p(r ) ⎧⎪ p (r ) ⎫⎪
Now, ≈ ∇⎨ ⎬
ρ (r ) ⎪⎩ ρ ⎪⎭
1
2
( )
∇ v• v − v× ∇× v +(∂v
∂t
= −∇ ⎨)
⎧⎪ p ( r ) ⎫⎪
⎪⎩ ρ ⎪⎭
⎬ − ∇φ
i.e.,
1
2
( )
∇ v• v − v× ∇× v +(∂v
∂t
= −∇ ⎨)
⎧⎪ p ( r )
⎪⎩ ρ
⎫⎪
+φ⎬
⎪⎭
44
PCD-09
1
2
( ) (
∇ v• v − v× ∇× v +
∂v
∂t
) ⎧⎪ p ( r )
= −∇ ⎨
⎪⎩ ρ
⎫⎪
+φ⎬
⎪⎭
⎛d ⎞ r ⎛d ⎞ r ur r
Recall that : ⎜ ⎟ r = ⎜ ⎟ r + ω x r
⎝ dt ⎠ I ⎝ dt ⎠ R
r r ur r r
v I = v R + ω x r R , where v I is just the
velocity that is employed in the equation
of motion for the fluid.
ur r ur r ur ur r
∴ ∇ × vI = ∇ × vR + ∇ × ω x r R { }
ur ur r
{ }
To determine ∇ × ω x r R we now use another vector
Identity, for the curl of cross-product of two vectors:
( ) ( ) ( )
∇ × A × B = B • ∇ A − A • ∇ B + A(∇ • B ) − B (∇ • A)
45
PCD-09
ur r ur r ur ur r
{
∇ × vI = ∇ × vR + ∇ × ω x r R }
( ) ( ) ( )
∇ × A × B = B • ∇ A − A • ∇ B + A (∇ • B ) − B (∇ • A)
( ) ( ) ( )
∇ × ω × rR = rR • ∇ ω − ω • ∇ rR + ω (∇ • rR ) − rR (∇ • ω )
∇ × (ω × r ) = − (ω • ∇ )r + ω (∇ • r ) = 2ω
R R R
∇ × v I = ∇ × v R + 2ω
1 2 ∂v ⎧⎪ p ( r ) ⎫⎪
∇ v − v× χ + = −∇ ⎨ +φ⎬ For ‘STEADY STATE’
2 ∂t ⎪⎩ ρ ⎪⎭
∂v ⎧⎪ p ( r ) ⎫⎪ 1 ∂v
=0
2
− v × χ = −∇ ⎨ +φ⎬− ∇ v
∂t ⎪⎩ ρ ⎪⎭ 2 ∂t
⎧ 2
⎫
⎧ 2
⎫ ⎪ p(r ) v ⎪
∂v ⎪ p(r ) v ⎪ − v × χ = −∇ ⎨ +φ + ⎬
− v × χ = −∇ ⎨ +φ + ⎬ ⎪ ρ 2 ⎪
∂t ⎪ ρ 2 ⎪ ⎩ ⎭
⎩ ⎭
⎧ v
2
⎫
⎪ p(r ) ⎪
Hence, 0 = v • ∇ ⎨ +φ + ⎬
⎪ ρ 2 ⎪ 47
PCD-09 ⎩ ⎭
⎧ v
2
⎫
⎪ p(r ) ⎪
0 = v • ∇⎨ +φ + ⎬
⎪ ρ 2 ⎪
⎩ ⎭
⎧ v
2
⎫
⎪ p (r ) ⎪
⇒ ∇⎨ +φ + ⎬ must be ORTHOGONAL to v,
⎪ ρ 2 ⎪
⎩ ⎭
⎧ v
2
⎫
⎪ p(r ) ⎪
i.e., ∇ ⎨ +φ + ⎬ must be ORTHOGONAL to STREAMLINES
⎪ ρ 2 ⎪
⎩ ⎭
⇒ ∇Ψ must be ORTHOGONAL to streamlines,
2
p(r ) v
where Ψ = +φ + .
ρ 2
2
p (r ) v
⇒Ψ= +φ + = constant for a given streamline
ρ
48
PCD-09 2 Daniel Bernoulli’s Theorem
2
p (r ) v
⇒Ψ= +φ + = constant for a given streamline
ρ 2
We derived the above result for a ‘STEADY STATE’ and made use of the relation
⎧ 2
⎫
⎪ p(r ) v ⎪ ∂v
− v × χ = −∇ ⎨
ρ
+φ + ⎬ =0
⎪
⎩
2 ⎪
⎭ ∂t
2
∇× v = χ = 0
p(r ) v
⇒Ψ= +φ +
ρ 2 Daniel Bernoulli’s Theorem
49
PCD-09 is constant for the entire velocity field in the liquid.
Mass Current Density Vector J ( r , t ) = ρ (r , t )v(r , t )
Dimensions : ML − 2T −1
G C
H D
For Steady State Flow,
P
ρvA = constant , F
B
E
A : cross - sectional area r A
since,
ur ur r ur r uur ∂ρ
∫∫∫ { } ∫∫
dτ ∇ • J (r ) = J (r ) • dS = 0, for STEADY STATE as
∂t
=0
volume surface
region enclosing
that
region
Work done on the fluid by the pressure Net work done on the fluid in the
that the fluid exerts on Face 1 is: parallelepiped by the pressure
δ W 1 = F1δ s = p 1 A1δ s = p 1 A1 v 1δ t that the fluid exerts on Faces 1
& 2 is:
Work done by the fluid on Face 2 is:
δW1 − δW2 = p1 A1v1δt − p2 A2 v 2δt
δ W 2 = F2 δ s = p 2 A 2 δ s = p 2 A 2 v 2 δ t 50
PCD-09
Net work done on the fluid in the
parallelepiped by the pressure δW1 −δW2 = p1 A1v1δt − p2 A2v2δt
that the fluid exerts at Faces 1 & 2 :
[ p1 A1 v 1 − p 2 A2 v 2 ]δ t ⎡1 ⎤ ⎡1 ⎤
= ⎢ v 2 + φ + U int ernal ⎥ − ⎢ v 2 + φ + U int ernal ⎥
ρ (δ sA ) ⎣2 ⎦2 ⎣2 ⎦1
51
PCD-09
[ p1 A1 v 1 − p 2 A2 v 2 ]δ t ⎡1 2 ⎤ ⎡1 2 ⎤
= ⎢ v + φ + U int ernal ⎥ − ⎢ v + φ + U int ernal ⎥
ρ (δ sA ) ⎣2 ⎦2 ⎣2 ⎦1
⎡ p1 A1 v 1 p 2 A2 v 2 ⎤ ⎡1 2 ⎤ ⎡1 2 ⎤
⎢ − ⎥ δ t = ⎢⎣ 2 v + φ + U int ernal ⎥ − ⎢ v + φ + U int ernal ⎥
⎣ ρ ( v 1δ t ) A1 ρ ( v 2 δ t ) A 2 ⎦ ⎦2 ⎣ 2 ⎦1
⎡ p1 p 2 ⎤ ⎡ 1 2 ⎤ ⎡1 2 ⎤
⎢ρ − =
⎥ ⎢2 v + φ + U int ernal ⎥ − ⎢2 v + φ + U int ernal ⎥
⎣ ρ ⎦ ⎣ ⎦2 ⎣ ⎦1
⎡1 p⎤ ⎡1 p⎤
0 = ⎢ v 2 + φ + U int ernal + ⎥ − ⎢ v 2 + φ + U int ernal + ⎥
⎣2 ρ ⎦2 ⎣2 ρ ⎦1
1 2 p From page 48 :
i .e. v + φ + U internal + = constant
2 ρ 2
Daniel Bernoulli’s Theorem p(r ) v
Ψ= +φ +
ρ 2
52
PCD-09 is constant for the entire velocity field in the liquid.
2
p (r ) v
+φ + = constant
ρ 2
The swing of a ball is
governed by Bernoulli's
theorem.
A swing bowler rubs only
one side of the ball. The ball
is then more rough on one Ishant Sharma
side than on the other. Inswing / Outswing
bowler
A white ball has a thin lacquer that is applied to its surface to avoid discoloring
the ball. During play, the shiny surface of the white ball remains shinier than
that of a red ball, which has a rougher surface to begin with.
The difference between the rough and shiny surface of a white ball is
much more, and thus it swings more than the red ball.
53
PCD-09
ρ James Clerk Maxwell
∇•E = 1831-1879
ε0
Divergence
∇•B = 0 and Curl of
∇× E = −
∂B (E, B )
∂t
∂J 1 ∂E
∇ × B = μ0 + 2
∂t c ∂t
1
c=
μ 0ε 0
Helmholtz Theorem
PH102: PHYSICS II 54
PCD-09
Revision Topic
ur ur ur ur ⎛ d ⎞ r ur ur r
F R = F I − F ω& − 2mω × ⎜ ⎟ r - mω × ω × r
⎝ dt ⎠ R
( )
‘Leap second’ term
PCD-09
‘Coriolis force’ ‘Centrifugal force’
r
Earth’s ω êv Choose eˆz = eˆv
rotation axis eˆr eˆx = eˆy × eˆz
θ −δ δ
eˆv = eˆvertical π
eˆy = eˆNorth = eˆEast
2 r
eˆh = eˆhorizontal Fcf = Fcf eˆρ
Frame of
r reference on the
Fg θ θ − δ eˆρ rotating earth to
θ −δ study centrifugal
êh and/or Coriolis
effects. The
1. Note earth’s axis of rotation. ‘vertical’ is not
r
2. Recognize that Fcf = Fcf eˆρ Note
r
that:
r
along the radial
3. Recognize that eˆ ≠ eˆ . ( Fg + Fcf ) • eˆh = 0 line, nor along
v r the axis of
4. Choose eˆz = eˆv . ∠ ( ˆ
eρ , ˆ
eh ) = θ − δ
earth’s rotation!
5. Choose eˆ y = eˆNorth , pointing toward North.
ω = (ω • eˆy ) eˆy + (ω • eˆz ) eˆz
r r r
6. Now, eˆx = eˆy × eˆz = eˆEast
PCD-09
ur ur ur ur ⎛ d ⎞ r ur ur r
F R = F I − F ω& − 2mω × ⎜ ⎟ r - mω × ω × r
⎝ dt ⎠ R
( )
ω = (ω • eˆ ) eˆ + (ω • eˆ ) eˆ
ur ur ⎛ d ⎞ r r r r
F Coriolis = −2mω × ⎜ ⎟ r y y z z
⎝ ⎠R
dt
= −2m ⎣(ω • eˆy ) eˆy + (ω • eˆz ) eˆz ⎦⎤ × ⎡⎣ v x eˆx +v y eˆy + v z eˆz ⎤⎦
r r
⎡
r
FR = m ( aRx ex + aRy ey + aRz ez )
) ) )
INTEGRATION
ω
approximations are often employed, such
as ignoring higher powers of
aRx = −2ω ( v Rz cos λ − v Ry sin λ )
v Ry = −2ω ( x sin λ ) + y& R 0
v Rz = − gt + 2ω ( x cos λ ) + z&R 0
integrating :
x&R = ω gt 2 cos λ − 2ωt ( z&R 0 cos λ − y& R 0 sin λ ) + x&R 0
1
xR = ω gt 3 cos λ − ωt 2 ( z&R 0 cos λ − y& R 0 sin λ ) + x&R 0t + xR 0
3
Ω = ω sin λ
r&& r r ur ⎛ d ⎞ r
Foucault Pendulum mrR = mg + S − 2mω × ⎜ ⎟ r
⎝ dt ⎠ R
xR = g cos α + 2Ωy&
&&
Coupled differential equations
yR = g cos β − 2Ωx&
&&
Ω = ω sin λ