* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
Download Small Current-Loops [ [ ].
Friction-plate electromagnetic couplings wikipedia , lookup
Maxwell's equations wikipedia , lookup
Magnetometer wikipedia , lookup
Giant magnetoresistance wikipedia , lookup
Neutron magnetic moment wikipedia , lookup
Earth's magnetic field wikipedia , lookup
Magnetotactic bacteria wikipedia , lookup
Superconducting magnet wikipedia , lookup
Magnetic monopole wikipedia , lookup
Electric dipole moment wikipedia , lookup
Multiferroics wikipedia , lookup
Relativistic quantum mechanics wikipedia , lookup
Magnetoreception wikipedia , lookup
Magnetotellurics wikipedia , lookup
Electromagnetism wikipedia , lookup
Electromagnet wikipedia , lookup
History of geomagnetism wikipedia , lookup
Magnetochemistry wikipedia , lookup
Magnetohydrodynamics wikipedia , lookup
Electromotive force wikipedia , lookup
Electromagnetic field wikipedia , lookup
Mathematical descriptions of the electromagnetic field wikipedia , lookup
Ferromagnetism wikipedia , lookup
PHY2206 (Electromagnetic Fields) Small Current Loops Small Current-Loops Vector potential due to a current loop The magnetic vector potential at a distant point r due to a current I flowing round a small loop C, is found by substituting J(r′ )d 3r ′ → Idr ′ in the definition of the potential so A( r ) = µ 0 I ⌠ dr ′ 4 π ⌡C r − r ′ (1) which can be approximated by a series in much the same way as was described on the previous Multipole Expansions handout. 2 2( r ⋅ r ′ ) r′ 1 1 = 1 + − r − r′ r r r 2 −1/ 2 (2) and since r >> r ′ a simple binomial series will converge rapidly 1 1 (r ⋅ r ′ ) = 1 + 2 + L r − r′ r r (3) where only terms to first order in r ′ /r have been retained, an approximation which requires that the origin of the primed coordinate system is close to the centre of the loop. Substituting into equation 1 A( r ) = µ0 I 1 ⌠ 1⌠ dr ′ + 3 (r ⋅ r ′ ) dr ′ + L 4 π r ⌡C r ⌡C (4) and the first term of this obviously vanishes. The second term can be rewritten with the help of a vector identity (VAF-2) which states that (r′ × dr′) × r = − r′(r ⋅ dr′) + dr′(r′ ⋅ r) (5) The small change in r′(r′ ⋅ r) due to a small change dr ′ in r ′ is d[r ′(r′ ⋅ r)] = r ′(r ⋅ dr ′ ) + dr ′(r ′ ⋅ r) (6) which is an exact differential. After adding equation 6 to equation 5 we find dr′(r′ ⋅ r) = 1 2 (r′ × dr′) × r + 12 d[r′(r′ ⋅ r)]. (7) The second term vanishes on integration because it is an exact differential and its integral between two points is therefore independent of the path, a closed loop in this case, leaving 1 PHY2206 (Electromagnetic Fields) Small Current Loops A( r ) = µ0I ⌠ 1 r 2 (r ′ × dr ′ ) × 3 + L 4 π ⌡C r (8) so, if the magnetic dipole moment m of the circuit is defined and the higher order terms are dropped A( r ) = µ0 m × r 4π r3 m= where I⌠ r′ × d r′ . 2 ⌡C (9) Field due to a current loop Having found the vector potential its curl can be used to find the field B(r) due to the current loop described above. We start by using a vector identity (VAF-15) B(r) = ∇ × A(r) = µ r m×r r ∇ × 3 = 0 −( m ⋅ ∇ ) 3 + m ∇ ⋅ 3 r r r 4π 4π µ0 (10) where the two terms in the identity involving ∇ ⋅ m have been dropped because m doesn’t depend on the coordinates. The first term is best dealt with by writing it out in component form, which in index notation and summing over all indices ∂ r −3 1 ∂ xk r ∂ xk x̂k x̂ = −m = −m + 3 j j k xk r3 ∂x j r3 ∂ x j r ∂ x j (11) r 3(m ⋅ r)r m −5 + r −3δ jk = − 3. 3 = −m j x̂k −3xk x j r r r5 r (12) −( m ⋅ ∇ ) and therefore −( m ⋅ ∇) [ ] The second term of equation 10 is zero, which can be shown by using (VAF-9) r 1 1 −3r 3 m ∇ ⋅ 3 = m r ⋅ ∇ 3 + 3 ∇ ⋅ r = m r ⋅ 5 + 3 = 0 r r r r r (13) so, combining all these results B(r) = µ 0 3(m ⋅ r)r 4 π r 5 − m r 3 which is known as the magnetic dipole field because of its similarity to the electric dipole field. 2 (14) PHY2206 (Electromagnetic Fields) Small Current Loops Scalar potential due to a current loop In regions of space where the current density is zero the curl of the magnetic field must also be zero so it can be described by a magnetic scalar potential φm B = − µ 0 ∇φ m . (15) Since the divergence of B is also zero φm satisfies Laplace’s equation which means that many results derived for electrostatics can be reused for magnetostatics. Things are not entirely straightforward as φm is often not single valued and getting boundary conditions right can be tricky. A simple example is the scalar potential outside a wire carrying current I, in cylindrical coordinates, φm = − Iθ . 2π (16) By comparing equation 14 with its electrostatic equivalent, the scalar potential of the magnetic dipole moment m is φ m (r ) = m⋅r . 4π r3 (17) This quantity is useful when calculating the field due to large current-loops, which can be represented as an array of many small loops, and when dealing with problems involving magnetic materials. Magnetic forces on a small current-loop To calculate the forces due to an inhomogeneous magnetic field B(r) on a loop carrying a current I consider a small rectangular loop of side δx and δy lying in the x–y plane. The net force components are ∂ Bz ∂B δx = m z ∂x ∂x ∂B ∂B Fy = Iδ xBz ( y + δ y ) − Iδ xBz ( y ) = Iδ x z δ x = m z ∂y ∂y Fx = Iδ yBz ( x + δ x ) − Iδ yBz ( x ) = Iδ y (18) ∂ B ∂ By ∂B ∂B = −m(∇ ⋅ B) + m z = m z Fz = −m x + ∂y ∂z ∂z ∂x where mẑ is the magnetic dipole moment of the loop. Adding these components together and generalising for a magnetic moment pointing in an arbitrary direction the net force is F = m∇(B ⋅ m̂) = ∇(B ⋅ m) − (B ⋅ m̂)∇m. and in the case that m doesn’t depend on the position of the loop, this simplifies and becomes 3 (19) PHY2206 (Electromagnetic Fields) Small Current Loops F = ∇(m ⋅ B) = (m ⋅ ∇)B + m × (∇ × B) (20) Note: Since an electric field E has no curl the equivalent expression for the force on an electrostatic dipole moment p is F = ∇(p ⋅ E) = (p ⋅ ∇)E. (21) The torque on the loop when it lies in the x–y plane has components Γ x = + δ xIδ yBx = +mz Bx Γ y = − δ yIδ xBy = −mz By Γz = 0 (22) with similar expressions obtained (from symmetric permutations of the coordinates) for the cases when it lies in the x–z and y–z planes. By considering each of these thre cases as the projection of an arbitrarily orientated loop the results can be summed to obtain the general expression Γ x = my Bz − mz By Γ y = mz Bx − mx Bz Γ z = mx By − my Bx i.e. Γ = m × B. (23) The results in this section apply to loops of any shape because these can be approximated to arbitrary accuracy by superpositions of smaller square loops. Potential energy of a current loop Since the definition of the potential energy V is that it satisfies F = −∇V and equation 20 has exactly this form the potential energy of the current loop is simply V m = −m ⋅ B (24) This potential must be used with care as it is not the total energy of the current loop because it was derived subject to the assumption that m is constant and this often not be the case, for example if we were to move the loop to a position where B was different the Lentz’s law current set up would change m, and energy would be needed to counteract this and keep m constant. © Copyright CDH Williams Exeter 1996, CW960430/2 4