Rajeev’s Almanack » Math/Physics

Archive for the 'Math/Physics' Category

A Random Pursuit

Saturday, December 1st, 2007

A missile is launched when an airplane is directly overhead, at a height h. The missile moves at a constant speed v always heading directly towards the aircraft, which is moving along a straight line at constant velocity u. What is the shape of the missile’s trajectory? How far will the aircraft fly before it is hit by the missile? What if the plane takes evasive action by randomly changing its direction, but heading in the original direction on the average? What is the best strategy for the missile to maximize the probability of a hit?

Pursuit
(more…)

Complex Time in Quantum Tunneling

Sunday, October 21st, 2007

Perhaps the most spectacular early prediction of quantum mechanics was tunneling: that particles can do things that are forbidden in Newton’s mechanics, although with a small probability. (more…)

The Geometry of Thermodynamics

Monday, October 8th, 2007

Discriminant

Thermodynamics is the study of heat. Originally developed to understand steam engines and such, it led to a revolution in physics. It showed that time has a preferred direction. Also, that physics is not fully deterministic: the best we can do for large systems is to predict averages of physical quantities and probabilities of events. But with the even greater revolutions of quantum mechanics and relativity that happened soon after , thermodynamics lost some of its original wonder. Nowadays it is thought of a staid old field, barely taught in physics departments anymore ( except as a preparation for a Stat Mech course). This is a pity, because thermodynamics is perhaps the most remarkable of all physical theories. We have none other than Albert Einstein vouching for this1: (more…)

Reduction or Emergence

Wednesday, September 5th, 2007

Earnest Rutherford used to say that all science is either physics or stamp collecting. This could have been a dig at the biologists of his time, who were still collecting samples and classifying species. He probably would have thought more highly of modern molecular biology, which is a lot like his physics in outlook: everything is determined by the DNA. It is said that Rutherford’s worst insult for a student who had done something stupid was–Chemist. The chemists had the last laugh though: Rutherford was awarded the Nobel Prize not in Physics but in Chemistry for having achieved the transmutation of elements.

Should we understand the world bottom up or top down? Which is the proper scientific view? (more…)

Matrix Diagonalization by Sampling

Saturday, September 1st, 2007

It is hard to think of a problem more ubiquitous than the diagonalization of a matrix.I will discuss today a statistical approximation method for finding the eigenvalues of a symmetric matrix. (more…)

Macroscopic Fluid Mechanics 3: References

Friday, June 1st, 2007

Part 1 Part 2 SIAM Talk
A recent review,closest to our point of view :


B. Khesin Topological fluid dynamics. Notices of the AMS, 52:1
(2005), 9-19

(more…)

Macroscopic Fluid Mechanics 2

Tuesday, May 29th, 2007

Part 1
To understand the origin of this non-commutativity, let us again consider the example of a hurricane. It is an extended object, whose radius is of the order of 100 km. It wouldn’t make sense to have two such objects within a 100 km of each other: the two hurricanes will interact strongly with each other and combine into a single one. (This phenomenon of a `reverse cascade’ can been seen clearly in some simulations.) Thus there is a limit to the resolution of the co-ordinates of a hurricane, given by the area of a hurricane. This is reminiscent of the uncertainty principle of quantum mechanics, except that there it is the co-ordinates of phase space that is fuzzy: the analogue of the area is Plank’s constant.
(more…)

Macroscopic Fluid Mechanics 1

Tuesday, May 29th, 2007

Introduction to Talk at SIAM DS07 Conference Snowbird Utah May 28th-June 1


The equations of motion of a fluid are obtained by averaging over the equations of motion of the large number of molecules that occupy even a small volume: we are not for the most part interested in the details of the motion of individual molecules. The equations of a fluid so obtained ( Euler or Navier-Stokes) are quite different from those of particle mechanics, being partial differential equations. Nevertheless the fundamental symmetries (translation and rotation invariance) of particle mechanics are preserved in this reformulation.
The conservation laws (energy, momentum,angular momentum) are preserved for ideal flow (Euler). In the next approximation, the effect of the transport of these conserved quantities to molecular scales are incorporated (viscosity). Even higher order corrections from molecular scales can be incorporated ( Chapman and Enskog) but are rarely needed.
(more…)

Fuzzy Fluid Mechanics

Sunday, May 13th, 2007

Terrence Tao has made some deep observations on why the regularity of three dimensional Navier-Stokes is such a hard problem. Although he has gone on to many other equally interesting topics, I remain fascinated by his main point there: that Navier-Stokes is supercritical. The nonlinearities become stronger at small distance scales, making it impossible to know (using present techniques) whether solutions remain smooth for all time. Thus, it is crucial to understand the scale dependence of non-linearities in fluid mechanics.
(more…)

Dynamics of Coin Tossing II

Sunday, April 29th, 2007

It appears that the instability that we described in the last article is indeed new; at least not known to experts such as Persi Diaconis and Richard Mongomery. At their encouragement I describe in more detail the effect of small pertrubations on a coin. We can look at this as a problem in Riemannian geometry or one in Classical Mechanics. Each brings its own insights and informs the other. Let us begin with the geometric point of view.

The Riemannian Submersion S^3\to S^2

In the orthogonal basis used above,

<e_1,e_1>=<e_2,e_2>=1, <e_3,e_3>=a^2

where we choose units such that

Q_1=1, Q_3=a^2. Since e_3 is a Killing vector whose orbits are circle, we can quotient by it to get a manifold S^2=SU(2)/SO(2). We work in the covering space SU(2)=S^3 of the orthogonal group for convenience.

There is an induced metric on S^2; it sectional curvature is given by O’Neill’s formula for Riemannian submersions ( Petersen p. 58)

S(e_1,e_2)=R_{1212}+3\left|{1\over 2}[e_1,e_2]\right|^2=1-{3\over 4}a^2+{3\over 4}a^2=1.

Thus the induced metric on the sphere is the standard one. In fact S^3 in this point of view is the circle bundle of S^2. But it is not the circle of unit tangent vectors; their length is a.

Geodesics on S^3 do not go project to geodesics in the base. The bundle U(1)\to SU(2)\to S^2 carries a natural connection; the horizontal vectors are orthogonal to the vertical vectors. Note that this connection is independent of a: the horizontal component is the same no matter what the radius of the fiber is. Thus the induced curves on the sphere are not geodesics, but are the same as what one would get for an isotropic rigid body.

The fact that a may not be one only affects the lift of these curves to the total space; that is where the negative sectional curvature can come in. When a>1 we are `stretching’ the vectors in the vertical direction. If this stretching is too large, it leads to an instability.

More coming…

Classical Mechanics

Landau Lifshitz sections 33 and 35 is a standard discussion of the rigid body.

The Lagrangian is, in Euler co-ordinates,

L={1\over 2}Q_1[\dot\phi^2+\sin^2\theta+\dot\theta^2]+{1\over 2}Q_3[\dot\phi\cos\theta+\dot\psi]^2.

This leads to the hamiltonian

H={1\over 2Q_1}\left[p_\theta^2 +{p_\phi^2-2p_\phi p_\psi\cos\theta +p_\psi^2\over \sin^2\theta}\right]+{p_\psi^2\over 2Q_1}[1-\eps],\ \eps=1-{I_1\over I_3}.

p_\phi and p_\psi are conserved. So is the total angular momentum:

L^2=p_\theta^2 +{p_\phi^2-2p_\phi p_\psi\cos\theta +p_\psi^2\over \sin^2\theta}.

The equation for \theta is not affected by the fact that Q_1\neq Q_3. This is what we saw geometrically also.

More coming …

Tags:

Dynamics of Coin Tossing

Monday, April 23rd, 2007

The standard example of a random event we give in probability courses is the tossing of a coin. On the other hand, in mechanics courses we teach that rigid body mechanics is one of the few examples of an integrable system. Not only are its equations deterministic, there is an explicit solution in terms of elliptic functions. In other words, a rigid body is one of the most predictable of all mechanical systems.
(more…)