CM M01 3

From QED

< PlasmaWiki | Prelims(Link to this page as [[PlasmaWiki/Prelims/CM M01 3]])
Jump to: navigation, search

A particle of mass m and charge q moves freely in a gravitational field \mathbf{B}=B\hat{k}. At time t=0 , the particle is released from the origin O with no initial velocity. It traces a curve in the x , y plane.

a. Find the parametric equations x=x(t), y=y(t) describing the curve. Sketch the curve on an x, y diagram.

The above motion is idealized, because two effects have been ignored: air drag and radiation damping

b. Now assume that the particle also feels a drag force due to the surrounding atmosphere, F = − βv. Derive the motion of the particle. What is its final velocity?

c. Instead of air drag, suppose we include the damping effect caused by the electromagnetic radiation emitted during its motion. Describe, qualitatively, how this modifies the motion found in part a. What is the final velocity of the particle?

We have the equations of motion:

m\ddot{x}=q\dot{y}B
m\ddot{y}=-q\dot{x}B+mg

Dotting these, dividing by m, and substituting:

\dot{\ddot{x}}=\frac{q}{m}\ddot{y}B=\frac{q}{m}B\left(-\frac{q}{m}\dot{x}B+g\right)
\dot{\ddot{y}}=-\frac{q}{m}\ddot{x}B=-\frac{q}{m}B\frac{q}{m}\dot{y}B

From these, we find the solutions have the form:

\dot{x}=A\cos\left(\frac{qB}{m}t\right)+C
\dot{y}=D\sin\left(\frac{qB}{m}t\right)+E

Solving for the constants, using \dot{x}(t=0)=\dot{y}(t=0)=0, we get:

\dot{x}=\frac{gm}{qB}\cos\left(\frac{qB}{m}t\right)-\frac{gm}{qB}
\dot{y}=\frac{gm}{qB}\sin\left(\frac{qB}{m}t\right)

Integrating:

x=\frac{gm^{2}}{q^{2}B^{2}}\sin\left(\frac{qB}{m}t\right)-\frac{gm}{qB}t+A
y=-\frac{gm^{2}}{q^{2}B^{2}}\cos\left(\frac{qB}{m}t\right)+C

Solving for the constants such that x(t=0)=y(t=0)=0 :

<math>x=\frac{gm^{2}}{q^{2}B^{2}}\sin\left(\frac{qB}{m}t\right)-\frac{gm}{qB}t</math>
y=-\frac{gm^{2}}{q^{2}B^{2}}\cos\left(\frac{qB}{m}t\right)+\frac{gm^{2}}{q^{2}B^{2}}

Which looks like a drifting loop:

The total force by the drag will be:

|F|=-\beta\sqrt{\dot{x}^{2}+\dot{y}^{2}}

which means that in the x and y directions:

F_{x}=-\beta\dot{x};\qquad F_{y}=-\beta\dot{y}

So we get the equations of motion:

<math>m\ddot{x}=q\dot{y}B-\beta\dot{x}</math>
m\ddot{y}=-q\dot{x}B+mg-\beta\dot{y}

We can find the final velocity by finding \dot{x} such that these equations are zero:

0=q\dot{y}B-\beta\dot{x}
0=-q\dot{x}B+mg-\beta\dot{y}

Plugging the first equation into the second:

0=-q\frac{\dot{y}qB}{\beta}B+mg-\beta\dot{y}
\dot{y}\left(\frac{q^{2}B^{2}}{\beta}+\beta\right)=mg
\dot{y}=\frac{mg\beta}{q^{2}B^{2}+\beta^{2}}

Giving:

\dot{x}=\frac{mgqB}{q^{2}B^{2}+\beta^{2}} v_{term}=\sqrt{\dot{x}^{2}+\dot{y}^{2}}=\frac{mg}{\sqrt{q^{2}B^{2}+\beta^{2}}}

The radiation will only slow the motion perpendicular to the motion of the guiding center, since we can consider that in the frame of the guiding center it is moving in a circle and will come to rest. Once this radial motion has all been emitted, it will simply move with the guiding center velocity:

v_{gc}=\frac{gm}{qB}

This page was recovered in October 2009 from the Plasmagicians page on Prelim_M01_Mech3 dated 21:20, 20 December 2005.

Personal tools