General Central Force 1. Definition and Notations: Here vectors will be displayed in bold blue, unit vectors in bold plum, and scalars (and components) in regular black. Since force is a vector, F should be a vector, F. If the force depends on position only, we can express force as a function of position as F(r), which in three dimensions means F(x,y,z) or F(r,,) . A central force is one where the force is directed along the radial direction, r, and not in the other two directions, or . Its magnitude also depends only on the radial distance, r, and not on or . This can be expressed as F = F(r) r . Be aware, however, that the radial direction, r, does depend on and : F = F(r) r(,). r(x,y,z) = xx + yy + zz = r sin() cos() x + r sin() sin() y + r cos() z = rr where r = r(,) = sin()cos() x + sin() sin() y + cos() z . 2. Is a central force (always, sometimes, never) a conservative force? To answer this, recall Stokes Theorem: (F) n dS = closed loop F dl (where S is an area) lf so that if F = 0 everywhere, then closed loop F dl = 0, and so li F dl is independent of path, and so a potential energy can be defined: rf PE = - ri F dl . Let’s now calculate F for a central force (using spherical coordinates): F = [r /r + (1/r) / + (1/{r sin()}) / ] F(r) r = (F(r)/r) (r r) + F(r) (r r/r) + (1/r) (F(r)/) ( r) + (F(r)/r) ( (r/)) + (1/{r sin()}) (F(r)/) ( r) + (1/{r sin()}) F(r) ( r/) . There are six terms, so let’s look at these one by one. The first term is zero since (r r) = 0 (since any vector crossed itself is zero). The second term is zero since r/r = 0 (since r depends only on and ). The third term is zero since F(r)/ = 0 (since F(r) does not depend on ). The fourth term is zero since r/ = , and ( ) = 0. The fifth term is zero since F(r)/ = 0 (since F(r) does not depend on ). The sixth term is zero since r/ = sin() , and ( ) = 0. Therefore, F = 0 for F = F(r)r, so a central force is always a conservative force! We can then define a potential energy function: r PE(r) = - r-standard F(r) dr . 3. Torque and Angular Momentum for a Central Force Starting with the definition of torque as ≡ r F, for a central force, F = F(r)r, we have (recalling that r = rr): = r F(r)r = r F(r) r r = 0. Central forces can create no torque. Now we recall the relation between torque and angular momentum (L): = dL/dt . Since torque under a central force is zero, the angular momentum must remain constant. Let’s now get an expression for angular momentum in polar coordinates. We can use 2-D polar (actually 3-D cylindrical) coordinates since a central force has no component of either motion or force perpendicular to the (r v) plane, i.e., the z plane. By definition, L ≡ r p = r mv . In cylindrical form, v = dr/dt = d(rr)/dt = (dr/dt)r + r(dr/dt) = r'r + r' (where r' = dr/dt, and ' = d/dt) . We will use this ' notation to only denote time derivatives, not any other derivatives. Therefore, L = rr m(r'r + r') = mrr'(r r) + mr2'(r ) = 0 + mr2' z = constant. Note that if we know the initial conditions, ro and 'o, we can calculate the value of L. If we just consider the magnitude, since the direction will not change, L = mr2' = constant. Eq. 1 This equation provides a relation between the radius, r and the angular speed, ' (which we sometimes have called ). We can write this as: ' = L/mr2 . Eq. 1a 4. Newton’s Second Law for a Central Force F = ma is a vector equation, which means it is one equation for each dimension. For a central force, F = F(r)r, in polar coordinates we have Fr = F(r) , and F = 0 . The acceleration expressed 2 2 2 2 in polar form is: ar = r''– r' (where r'' = d r/dt is the usual linear acceleration in the radial direction, and r' 2 2 2 = r is the circular acceleration); and a = r'' + 2r'' (where r'' is r(d /dt ), the acceleration in the circular direction, and 2r'' is that strange Coriolis acceleration). Therefore, in component form, Newton’s Law becomes: r component: F(r) = m(r'' – r'2) , or F(r) + mr'2 = mr'' Eqs. 2 2 2 [where the term mr’ = m r can be considered the “centrifugal” force] component: F = 0, so 0 = m(r'' + 2r'') . Notice that both equations have both r and so they are simultaneous differential equations. However, if we use the angular momentum equation, ' = L/mr2 in the radial equation, we get: F(r) + mr(L/mr2)2 = mr'', or F(r) + L2/mr3 = mr'' Eq. 2a which is an expression that is equivalent to a one-dimensional problem since there is no mention of . If we can solve this for r(t), we can then use the component equation to solve for '(t) and then q(t). We could also use ' = L/mr2 with r(t) and both L and m constant to solve for '(t) and then q(t). 5. Conservation of Energy Since a central force is a conservative force, we can define a potential energy. If there are no other forces, then energy is conserved and we can write this in an equation: KE + PE = E = constant 2 2 2 2 2 v = vrr + v , and v = vr + v = (r') + (r') so we get: ½ m v2 + PE(r) = ½ m r'2 + ½ m r2 '2 + PE(r) = E . If we use the angular momentum expression, ' = L/mr2, this becomes: ½ m r'2 + ½ m r2 (L/mr2)2 + PE(r) = E = ½ m r'2 + ½L2/mr2 + PE(r) . If we treat the term ½L2/mr2 as an additional potential energy which is equivalent to treating the centripetal acceleration as a centrifugal force (since it only depends on r), we can define an effective potential energy: 2 2 PEeff(r) = PE(r) + L /2mr . Eq. 3 We now have an equation that is essentially one-dimensional: 2 ½ m r' + PEeff(r) = E = constant . Eq. 4 As always, this is a first order differential equation whose solution is r(t). Then we can use ' = L/mr2 to solve for q'(t) and then q(t). 6. Straightforward solutions a) Use the Conservation of Energy (1-D) equation with the equivalent potential energy, and use the technique we used earlier in the course to first get r(t) and then to get (t): 2 2 2 ½ mr' + PEeff(r) = E where PEeff(r) = PE(r) + L /2mr ½ dr/dt = r' = [2(E-PEeff(r))/m ] t r ½ 0 dt = ro {1 / [2(E-PEeff(r))/m ] } dr Assuming we can actually do the integration over r, we get an expression for t(r). Then we must find the inverse to get r(t). This is sometimes very difficult. If we can get r(t), we can then use the angular momentum expression: L = mr2' to get an expression for ': '(t) = L/mr(t)2 . We can then use ' = d/dt to get: t t 2 0 d = 0 ' dt = 0 [L/mr(t) ] dt . If we can do the integration, we then have an expression for (t) and the problem is solved. b) Use Newton’s Second Law to try to get the trajectory r() [instead of r(t) and (t)]. While this will lose some information (time information), it has the advantage of giving us the orbits directly. We start with the equivalent 1-D relation (see part 4 above): F(r) + L2/mr3 = mr'' . Eq. 2a Since we want r(), we still have to deal with dr/dt, but we note that (recalling that ' = L/mr2) : r' = dr/dt = (dr/d)*(d/dt) = (dr/d)*' = (dr/d)*(L/mr2) . r'' = d(r')/dt = (dr'/d)*(d/dt) = d[(dr/d)*(L/mr2)]/d * ' = (d2r/d2)*(L/mr2)*' + (dr/d)*(L/m)*(d[1/r2]/d)*' note: d[1/r2]/d = (d[1/r2]/dr)*(dr/d) = (-2/r3)*(dr/d) r'' = (d2r/d2)*(L/mr2)*' + (dr/d)*(L/m)*(-2/r3)*(dr/d)*' (again using ' = L/mr2 ) r'' = (d2r/d2)*(L2/m2r4) - (2L2/m2r5)*(dr/d)2 . Using this, we can now write Newton’s 2nd Law (r component equation, Eq. 2a) as: F(r) + L2/mr3 = (d2r/d2)*(L2/mr4) - (2L2/mr5)*(dr/d)2 . Eq. 5 This doesn’t look all that easy to solve for r(q). 7. Tricky way to solve for r() Sometimes things look nicer in one form than in another. We will try a trick that involves substitution: let’s define u = 1/r (or r = 1/u). L = mr2' = m'/u2, or ' = Lu2/m r' = dr/dt = (dr/du)*(du/dt) = (dr/du)*(du/d)*(d/dt) = (d(1/u)/du)*(du/d)*' = (-1/u2)*(du/d)*' = (-1/u2)*(du/d)*(Lu2/m) = (-L/m)*(du/d). r'' = dr’/dt = (dr'/d)*(d/dt) = d[(-L/m)*(du/d)]/d * ' = (-L/m)*(d2u/d2)*(Lu2/m) = (-L2u2/m2)*d2u/d2 . Newton’s Second Law (r component equation) F(r) + L2/mr3 = mr’'' Eq. 2a thus becomes: F(1/u) + L2u3/m = m*(-L2u2/m2)*d2u/d2 , or (-m/L2u2)*F(1/u) + (-m/L2u2)*(L2u3/m) = d2u/d2 , or (-m/L2u2)*F(1/u) + -u = d2u/d2 .