magnetic field

Calculating the magnetostatic field from the moment

November 14, 2016 math and physics play , , , , , ,

[Click here for a PDF of this post with nicer formatting]

The vector potential, to first order, for a magnetostatic localized current distribution was found to be

\begin{equation}\label{eqn:magneticFieldFromMoment:20}
\BA(\Bx) = \frac{\mu_0}{4 \pi} \frac{\Bm \cross \Bx}{\Abs{\Bx}^3}.
\end{equation}

Initially, I tried to calculate the magnetic field from this, but ran into trouble. Here’s a new try.

\begin{equation}\label{eqn:magneticFieldFromMoment:40}
\begin{aligned}
\BB
&=
\frac{\mu_0}{4 \pi}
\spacegrad \cross \lr{ \Bm \cross \frac{\Bx}{r^3} } \\
&=
-\frac{\mu_0}{4 \pi}
\spacegrad \cdot \lr{ \Bm \wedge \frac{\Bx}{r^3} } \\
&=
-\frac{\mu_0}{4 \pi}
\lr{
(\Bm \cdot \spacegrad) \frac{\Bx}{r^3}
-\Bm \spacegrad \cdot \frac{\Bx}{r^3}
} \\
&=
\frac{\mu_0}{4 \pi}
\lr{
-\frac{(\Bm \cdot \spacegrad) \Bx}{r^3}
– \lr{ \Bm \cdot \lr{\spacegrad \inv{r^3} }} \Bx
+\Bm (\spacegrad \cdot \Bx) \inv{r^3}
+\Bm \lr{\spacegrad \inv{r^3} } \cdot \Bx
}.
\end{aligned}
\end{equation}

Here I’ve used \( \Ba \cross \lr{ \Bb \cross \Bc } = -\Ba \cdot \lr{ \Bb \wedge \Bc } \), and then expanded that with \( \Ba \cdot \lr{ \Bb \wedge \Bc } = (\Ba \cdot \Bb) \Bc – (\Ba \cdot \Bc) \Bb \). Since one of these vectors is the gradient, care must be taken to have it operate on the appropriate terms in such an expansion.

Since we have \( \spacegrad \cdot \Bx = 3 \), \( (\Bm \cdot \spacegrad) \Bx = \Bm \), and \( \spacegrad 1/r^n = -n \Bx/r^{n+2} \), this reduces to

\begin{equation}\label{eqn:magneticFieldFromMoment:60}
\begin{aligned}
\BB
&=
\frac{\mu_0}{4 \pi}
\lr{
– \frac{\Bm}{r^3}
+ 3 \frac{(\Bm \cdot \Bx) \Bx}{r^5} %
+ 3 \Bm \inv{r^3}
-3 \Bm \frac{\Bx}{r^5} \cdot \Bx
} \\
&=
\frac{\mu_0}{4 \pi}
\frac{3 (\Bm \cdot \ncap) \ncap -\Bm}{r^3},
\end{aligned}
\end{equation}

which is the desired result.

Magnetostatic force and torque

October 18, 2016 math and physics play , , , , , , ,

[Click here for a PDF of this post with nicer formatting]

In Jackson [1], the following equations for the vector potential, magnetostatic force and torque are derived

\begin{equation}\label{eqn:magnetostaticsJacksonNotesForceAndTorque:20}
\Bm = \inv{2} \int \Bx’ \cross \BJ(\Bx’) d^3 x’
\end{equation}
\begin{equation}\label{eqn:magnetostaticsJacksonNotesForceAndTorque:40}
\BF = \spacegrad( \Bm \cdot \BB ),
\end{equation}
\begin{equation}\label{eqn:magnetostaticsJacksonNotesForceAndTorque:60}
\BN = \Bm \cross \BB,
\end{equation}

where \( \BB \) is an applied external magnetic field and \( \Bm \) is the magnetic dipole for the current in question. These results (and a similar one derived earlier for the vector potential \( \BA \)) all follow from
an analysis of localized current densities \( \BJ \), evaluated far enough away from the current sources.

For the force and torque, the starting point for the force is one that had me puzzled a bit. Namely

\begin{equation}\label{eqn:magnetostaticsJacksonNotesForceAndTorque:80}
\BF = \int \BJ(\Bx) \cross \BB(\Bx) d^3 x
\end{equation}

This is clearly the continuum generalization of the point particle Lorentz force equation, which for \( \BE = 0 \) is:

\begin{equation}\label{eqn:magnetostaticsJacksonNotesForceAndTorque:100}
\BF = q \Bv \cross \BB
\end{equation}

For the point particle, this is the force on the particle when it is in the external field \( BB \). i.e. this is the force at the position of the particle. My question is what does it mean to sum all the forces on the charge distribution over all space.
How can a force be applied over all, as opposed to a force applied at a single point, or against a surface?

In the special case of a localized current density, this makes some sense. Considering the other half of the force equation \( \BF = \ddt{}\int \rho_m \Bv dV \), where \( \rho_m \) here is mass density of the charged particles making up the continuous current distribution. The other half of this \( \BF = m\Ba \) equation is also an average phenomena, so we have an average of sorts on both the field contribution to the force equation and the mass contribution to the force equation. There is probably a centre-of-mass and centre-of-current density interpretation that would make a bit more sense of this continuum force description.

It’s kind of funny how you can work through all the detailed mathematical steps in a book like Jackson, but then go right back to the beginning and say “Hey, what does that even mean”?

Force

Moving on from the pondering of the meaning of the equation being manipulated, let’s do the easy part, the derivation of the results that Jackson comes up with.

Writing out \ref{eqn:magnetostaticsJacksonNotesForceAndTorque:80} in coordinates

\begin{equation}\label{eqn:magnetostaticsJacksonNotesForceAndTorque:320}
\BF = \epsilon_{ijk} \Be_i \int J_j B_k d^3 x.
\end{equation}

To first order, a slowly varying (external) magnetic field can be expanded around a point of interest

\begin{equation}\label{eqn:magnetostaticsJacksonNotesForceAndTorque:120}
\BB(\Bx) = \BB(\Bx_0) + \lr{ \Bx – \Bx_0 } \cdot \spacegrad \BB,
\end{equation}

where the directional derivative is evaluated at the point \( \Bx_0 \) after the gradient operation. Setting the origin at this point \( \Bx_0 \) gives

\begin{equation}\label{eqn:magnetostaticsJacksonNotesForceAndTorque:340}
\begin{aligned}
\BF
&= \epsilon_{ijk} \Be_i
\lr{
\int J_j(\Bx’) B_k(0) d^3 x’
+
\int J_j(\Bx’) (\Bx’ \cdot \spacegrad) B_k(0) d^3 x’
} \\
&=
\epsilon_{ijk} \Be_i
\Bk_0 \int J_j(\Bx’) d^3 x’
+
\epsilon_{ijk} \Be_i
\int J_j(\Bx’) (\Bx’ \cdot \spacegrad) B_k(0) d^3 x’.
\end{aligned}
\end{equation}

We found

earlier
that the first integral can be written as a divergence

\begin{equation}\label{eqn:magnetostaticsJacksonNotesForceAndTorque:140}
\int J_j(\Bx’) d^3 x’
=
\int \spacegrad’ \cdot \lr{ \BJ(\Bx’) x_j’ } dV’,
\end{equation}

which is zero when the integration surface is outside of the current localization region. We also found

that

\begin{equation}\label{eqn:magnetostaticsJacksonNotesForceAndTorque:160}
\int (\Bx \cdot \Bx’) \BJ
= -\inv{2} \Bx \cross \int \Bx’ \cross \BJ = \Bm \cross \Bx.
\end{equation}

so
\begin{equation}\label{eqn:magnetostaticsJacksonNotesForceAndTorque:180}
\begin{aligned}
\int (\spacegrad B_k(0) \cdot \Bx’) J_j
&= -\inv{2} \lr{ \spacegrad B_k(0) \cross \int \Bx’ \cross \BJ}_j \\
&= \lr{ \Bm \cross (\spacegrad B_k(0)) }_j.
\end{aligned}
\end{equation}

This gives

\begin{equation}\label{eqn:magnetostaticsJacksonNotesForceAndTorque:200}
\begin{aligned}
\BF
&= \epsilon_{ijk} \Be_i \lr{ \Bm \cross (\spacegrad B_k(0)) }_j \\
&= \epsilon_{ijk} \Be_i \lr{ \Bm \cross \spacegrad }_j B_k(0) \\
&= (\Bm \cross \spacegrad) \cross \BB(0) \\
&= -\BB(0) \cross (\Bm \cross \lspacegrad) \\
&= (\BB(0) \cdot \Bm) \lspacegrad – (\BB \cdot \lspacegrad) \Bm \\
&= \spacegrad (\BB(0) \cdot \Bm) – \Bm (\spacegrad \cdot \BB(0)).
\end{aligned}
\end{equation}

The second term is killed by the magnetic Gauss’s law, leaving to first order

\begin{equation}\label{eqn:magnetostaticsJacksonNotesForceAndTorque:220}
\BF = \spacegrad \lr{\Bm \cdot \BB}.
\end{equation}

Torque

For the torque we have a similar quandary at the starting point. About what point is a continuum torque integral of the following form

\begin{equation}\label{eqn:magnetostaticsJacksonNotesForceAndTorque:240}
\BN = \int \Bx’ \cross (\BJ(\Bx’) \cross \BB(\Bx’)) d^3 x’?
\end{equation}

Ignoring that detail again, assuming the answer has something to do with the centre of mass and parallel axis theorem, we can proceed with a constant approximation of the magnetic field

\begin{equation}\label{eqn:magnetostaticsJacksonNotesForceAndTorque:260}
\begin{aligned}
\BN
&= \int \Bx’ \cross (\BJ(\Bx’) \cross \BB(0)) d^3 x’ \\
&=
-\int (\Bx’ \cdot \BJ(\Bx’)) \BB(0) d^3 x’
+\int (\Bx’ \cdot \BB(0)) \BJ(\Bx’) d^3 x’ \\
&=
-\BB(0) \int (\Bx’ \cdot \BJ(\Bx’)) d^3 x’
+\int (\Bx’ \cdot \BB(0)) \BJ(\Bx’) d^3 x’.
\end{aligned}
\end{equation}

Jackson’s trick for killing the first integral is to transform it into a divergence by evaluating

\begin{equation}\label{eqn:magnetostaticsJacksonNotesForceAndTorque:280}
\begin{aligned}
\spacegrad \cdot \lr{ \BJ \Abs{\Bx}^2 }
&=
(\spacegrad \cdot \BJ) \Abs{\Bx}^2
+
\BJ \cdot \spacegrad \Abs{\Bx}^2 \\
&=
\BJ \cdot \Be_i \partial_i x_m x_m \\
&=
2 \BJ \cdot \Be_i \delta_{im} x_m \\
&=
2 \BJ \cdot \Bx,
\end{aligned}
\end{equation}

so

\begin{equation}\label{eqn:magnetostaticsJacksonNotesForceAndTorque:300}
\begin{aligned}
\BN
&=
-\inv{2} \BB(0) \int \spacegrad’ \cdot \lr{ \BJ(\Bx’) \Abs{\Bx’}^2 } d^3 x’
+\int (\Bx’ \cdot \BB(0)) \BJ(\Bx’) d^3 x’ \\
&=
-\inv{2} \BB(0) \oint \Bn \cdot \lr{ \BJ(\Bx’) \Abs{\Bx’}^2 } d^3 x’
+\int (\Bx’ \cdot \BB(0)) \BJ(\Bx’) d^3 x’.
\end{aligned}
\end{equation}

Again, the localized current density assumption kills the surface integral. The second integral can be evaluated with \ref{eqn:magnetostaticsJacksonNotesForceAndTorque:160}, so to first order we have

\begin{equation}\label{eqn:magnetostaticsJacksonNotesForceAndTorque:360}
\BN
=
\Bm \cross \BB.
\end{equation}

References

[1] JD Jackson. Classical Electrodynamics. John Wiley and Sons, 2nd edition, 1975.

Magnetic moment for a localized magnetostatic current

October 13, 2016 math and physics play , , , , , , , , ,

[Click here for a PDF of this post with nicer formatting]

Motivation.

I was once again reading my Jackson [2]. This time I found that his presentation of magnetic moment didn’t really make sense to me. Here’s my own pass through it, filling in a number of details. As I did last time, I’ll also translate into SI units as I go.

Vector potential.

The Biot-Savart expression for the magnetic field can be factored into a curl expression using the usual tricks

\begin{equation}\label{eqn:magneticMomentJackson:20}
\begin{aligned}
\BB
&= \frac{\mu_0}{4\pi} \int \frac{\BJ(\Bx’) \cross (\Bx – \Bx’)}{\Abs{\Bx – \Bx’}^3} d^3 x’ \\
&= -\frac{\mu_0}{4\pi} \int \BJ(\Bx’) \cross \spacegrad \inv{\Abs{\Bx – \Bx’}} d^3 x’ \\
&= \frac{\mu_0}{4\pi} \spacegrad \cross \int \frac{\BJ(\Bx’)}{\Abs{\Bx – \Bx’}} d^3 x’,
\end{aligned}
\end{equation}

so the vector potential, through its curl, defines the magnetic field \( \BB = \spacegrad \cross \BA \) is given by

\begin{equation}\label{eqn:magneticMomentJackson:40}
\BA(\Bx) = \frac{\mu_0}{4 \pi} \int \frac{J(\Bx’)}{\Abs{\Bx – \Bx’}} d^3 x’.
\end{equation}

If the current source is localized (zero outside of some finite region), then there will always be a region for which \( \Abs{\Bx} \gg \Abs{\Bx’} \), so the denominator yields to Taylor expansion

\begin{equation}\label{eqn:magneticMomentJackson:60}
\begin{aligned}
\inv{\Abs{\Bx – \Bx’}}
&=
\inv{\Abs{\Bx}} \lr{1 + \frac{\Abs{\Bx’}^2}{\Abs{\Bx}^2} – 2 \frac{\Bx \cdot \Bx’}{\Abs{\Bx}^2} }^{-1/2} \\
&\approx
\inv{\Abs{\Bx}} \lr{ 1 + \frac{\Bx \cdot \Bx’}{\Abs{\Bx}^2} } \\
&=
\inv{\Abs{\Bx}} + \frac{\Bx \cdot \Bx’}{\Abs{\Bx}^3}.
\end{aligned}
\end{equation}

so the vector potential, far enough away from the current source is
\begin{equation}\label{eqn:magneticMomentJackson:80}
\BA(\Bx)
=
\frac{\mu_0}{4 \pi} \int \frac{J(\Bx’)}{\Abs{\Bx}} d^3 x’
+\frac{\mu_0}{4 \pi} \int \frac{(\Bx \cdot \Bx’)J(\Bx’)}{\Abs{\Bx}^3} d^3 x’.
\end{equation}

Jackson uses a sneaky trick to show that the first integral is killed for a localized source. That trick appears to be based on evaluating the following divergence

\begin{equation}\label{eqn:magneticMomentJackson:100}
\begin{aligned}
\spacegrad \cdot (\BJ(\Bx) x_i)
&=
(\spacegrad \cdot \BJ) x_i
+
(\spacegrad x_i) \cdot \BJ \\
&=
(\Be_k \partial_k x_i) \cdot\BJ \\
&=
\delta_{ki} J_k \\
&=
J_i.
\end{aligned}
\end{equation}

Note that this made use of the fact that \( \spacegrad \cdot \BJ = 0 \) for magnetostatics. This provides a way to rewrite the current density as a divergence

\begin{equation}\label{eqn:magneticMomentJackson:120}
\begin{aligned}
\int \frac{J(\Bx’)}{\Abs{\Bx}} d^3 x’
&=
\Be_i \int \frac{\spacegrad’ \cdot (x_i’ \BJ(\Bx’))}{\Abs{\Bx}} d^3 x’ \\
&=
\frac{\Be_i}{\Abs{\Bx}} \int \spacegrad’ \cdot (x_i’ \BJ(\Bx’)) d^3 x’ \\
&=
\frac{1}{\Abs{\Bx}} \oint \Bx’ (d\Ba \cdot \BJ(\Bx’)).
\end{aligned}
\end{equation}

When \( \BJ \) is localized, this is zero provided we pick the integration surface for the volume outside of that localization region.

It is now desired to rewrite \( \int \Bx \cdot \Bx’ \BJ \) as a triple cross product since the dot product of such a triple cross product has exactly this term in it

\begin{equation}\label{eqn:magneticMomentJackson:140}
\begin{aligned}
– \Bx \cross \int \Bx’ \cross \BJ
&=
\int (\Bx \cdot \Bx’) \BJ

\int (\Bx \cdot \BJ) \Bx’ \\
&=
\int (\Bx \cdot \Bx’) \BJ

\Be_k x_i \int J_i x_k’,
\end{aligned}
\end{equation}

so
\begin{equation}\label{eqn:magneticMomentJackson:160}
\int (\Bx \cdot \Bx’) \BJ
=
– \Bx \cross \int \Bx’ \cross \BJ
+
\Be_k x_i \int J_i x_k’.
\end{equation}

To get of this second term, the next sneaky trick is to consider the following divergence

\begin{equation}\label{eqn:magneticMomentJackson:180}
\begin{aligned}
\oint d\Ba’ \cdot (\BJ(\Bx’) x_i’ x_j’)
&=
\int dV’ \spacegrad’ \cdot (\BJ(\Bx’) x_i’ x_j’) \\
&=
\int dV’ (\spacegrad’ \cdot \BJ)
+
\int dV’ \BJ \cdot \spacegrad’ (x_i’ x_j’) \\
&=
\int dV’ J_k \cdot \lr{ x_i’ \partial_k x_j’ + x_j’ \partial_k x_i’ } \\
&=
\int dV’ \lr{J_k x_i’ \delta_{kj} + J_k x_j’ \delta_{ki}} \\
&=
\int dV’ \lr{J_j x_i’ + J_i x_j’}.
\end{aligned}
\end{equation}

The surface integral is once again zero, which means that we have an antisymmetric relationship in integrals of the form

\begin{equation}\label{eqn:magneticMomentJackson:200}
\int J_j x_i’ = -\int J_i x_j’.
\end{equation}

Now we can use the tensor algebra trick of writing \( y = (y + y)/2 \),

\begin{equation}\label{eqn:magneticMomentJackson:220}
\begin{aligned}
\int (\Bx \cdot \Bx’) \BJ
&=
– \Bx \cross \int \Bx’ \cross \BJ
+
\Be_k x_i \int J_i x_k’ \\
&=
– \Bx \cross \int \Bx’ \cross \BJ
+
\inv{2} \Be_k x_i \int \lr{ J_i x_k’ + J_i x_k’ } \\
&=
– \Bx \cross \int \Bx’ \cross \BJ
+
\inv{2} \Be_k x_i \int \lr{ J_i x_k’ – J_k x_i’ } \\
&=
– \Bx \cross \int \Bx’ \cross \BJ
+
\inv{2} \Be_k x_i \int (\BJ \cross \Bx’)_j \epsilon_{ikj} \\
&=
– \Bx \cross \int \Bx’ \cross \BJ

\inv{2} \epsilon_{kij} \Be_k x_i \int (\BJ \cross \Bx’)_j \\
&=
– \Bx \cross \int \Bx’ \cross \BJ

\inv{2} \Bx \cross \int \BJ \cross \Bx’ \\
&=
– \Bx \cross \int \Bx’ \cross \BJ
+
\inv{2} \Bx \cross \int \Bx’ \cross \BJ \\
&=
-\inv{2} \Bx \cross \int \Bx’ \cross \BJ,
\end{aligned}
\end{equation}

so

\begin{equation}\label{eqn:magneticMomentJackson:240}
\BA(\Bx) \approx \frac{\mu_0}{4 \pi \Abs{\Bx}^3} \lr{ -\frac{\Bx}{2} } \int \Bx’ \cross \BJ(\Bx’) d^3 x’.
\end{equation}

Letting

\begin{equation}\label{eqn:magneticMomentJackson:260}
\boxed{
\Bm = \inv{2} \int \Bx’ \cross \BJ(\Bx’) d^3 x’,
}
\end{equation}

the far field approximation of the vector potential is
\begin{equation}\label{eqn:magneticMomentJackson:280}
\boxed{
\BA(\Bx) = \frac{\mu_0}{4 \pi} \frac{\Bm \cross \Bx}{\Abs{\Bx}^3}.
}
\end{equation}

Note that when the current is restricted to an infintisimally thin loop, the magnetic moment reduces to

\begin{equation}\label{eqn:magneticMomentJackson:300}
\Bm(\Bx) = \frac{I}{2} \int \Bx \cross d\Bl’.
\end{equation}

Refering to [1] (pr. 1.60), this can be seen to be \( I \) times the “vector-area” integral.

References

[1] David Jeffrey Griffiths and Reed College. Introduction to electrodynamics. Prentice hall Upper Saddle River, NJ, 3rd edition, 1999.

[2] JD Jackson. Classical Electrodynamics. John Wiley and Sons, 2nd edition, 1975.

Electric and magnetic fields at an interface

October 9, 2016 math and physics play , , ,

[Click here for a PDF of this post with nicer formatting]

As pointed out in [1] the fields at an interface that is not a perfect conductor on either side are related by

\begin{equation}\label{eqn:fieldsAtInterface:20}
\begin{aligned}
\ncap \cdot \lr{ \BD_2 – \BD_1 } &= \rho_{es} \\
\ncap \cross \lr{ \BE_2 – \BE_1 } &= -\BM_s \\
\ncap \cdot \lr{ \BB_2 – \BB_1 } &= \rho_{ms} \\
\ncap \cross \lr{ \BH_2 – \BH_1 } &= \BJ_s.
\end{aligned}
\end{equation}

Given the fields in medium 1, assuming that boths sets of media are linear, we can use these relationships to determine the fields in the other medium.

\begin{equation}\label{eqn:fieldsAtInterface:40}
\begin{aligned}
\ncap \cdot \BE_2 &= \inv{\epsilon_2} \lr{ \epsilon_1 \ncap \cdot \BE_1 + \rho_{es} } \\
\ncap \wedge \BE_2 &= \ncap \wedge \BE_1 -I \BM_s \\
\ncap \cdot \BB_2 &= \ncap \cdot \BB_1 + \rho_{ms} \\
\ncap \wedge \BB_2 &= \mu_2 \lr{ \inv{\mu_1} \ncap \wedge \BB_1 + I \BJ_s}.
\end{aligned}
\end{equation}

Now the fields in interface 2 can be obtained by adding the normal and tangential projections. For the electric field

\begin{equation}\label{eqn:fieldsAtInterface:60}
\begin{aligned}
\BE_2
&=
\ncap (\ncap \cdot \BE_2 )
+ \ncap \cdot (\ncap \wedge \BE_2) \\
&=
\inv{\epsilon_2} \ncap \lr{ \epsilon_1 \ncap \cdot \BE_1 + \rho_{es} }
+
\ncap \cdot (\ncap \wedge \BE_1 -I \BM_s).
\end{aligned}
\end{equation}

Note that this manipulation can also be done without Geometric Algebra by writing \( \BE_2 = \ncap (\ncap \cdot \BE_2 ) – \ncap \cross (\ncap \cross \BE_2) \)).
Expanding \( \ncap \cdot (\ncap \wedge \BE_1) = \BE_1 – \ncap (\ncap \cdot \BE_1) \), and \( \ncap \cdot (I \BM_s) = -\ncap \cross \BM_s \), that is

\begin{equation}\label{eqn:fieldsAtInterface:80}
\boxed{
\BE_2
=
\BE_1
+ \ncap (\ncap \cdot \BE_1) \lr{ \frac{\epsilon_1}{\epsilon_2} – 1 }
+ \frac{\rho_{es}}{\epsilon_2}
+ \ncap \cross \BM_s.
}
\end{equation}

For the magnetic field

\begin{equation}\label{eqn:fieldsAtInterface:100}
\begin{aligned}
\BB_2
&=
\ncap (\ncap \cdot \BB_2 )
+
\ncap \cdot (\ncap \wedge \BB_2) \\
&=
\ncap \lr{ \ncap \cdot \BB_1 + \rho_{ms} }
+
\mu_2 \ncap \cdot \lr{ \lr{ \inv{\mu_1} \ncap \wedge \BB_1 + I \BJ_s} },
\end{aligned}
\end{equation}

which is

\begin{equation}\label{eqn:fieldsAtInterface:120}
\boxed{
\BB_2
=
\frac{\mu_2}{\mu_1} \BB_1
+
\ncap (\ncap \cdot \BB_1) \lr{ 1 – \frac{\mu_2}{\mu_1} }
+ \ncap \rho_{ms}
– \ncap \cross \BJ_s.
}
\end{equation}

These are kind of pretty results, having none of the explicit angle dependence that we see in the Fresnel relationships. In this analysis, it is assumed there is only a transmitted component of the ray in question, and no reflected component. Can we do a purely vectoral treatment of the Fresnel equations along these same lines?

References

[1] Constantine A Balanis. Advanced engineering electromagnetics. Wiley New York, 1989.

Green’s function inversion of the magnetostatic equation

September 27, 2016 math and physics play , , , , , , , , ,

[Click here for a PDF of this post with nicer formatting]

A previous example of inverting a gradient equation was the electrostatics equation. We can do the same for the magnetostatics equation, which has the following Geometric Algebra form in linear media

\begin{equation}\label{eqn:biotSavartGreens:20}
\spacegrad I \BB = – \mu \BJ.
\end{equation}

The Green’s inversion of this is
\begin{equation}\label{eqn:biotSavartGreens:40}
\begin{aligned}
I \BB(\Bx)
&= \int_V dV’ G(\Bx, \Bx’) \spacegrad’ I \BB(\Bx’) \\
&= \int_V dV’ G(\Bx, \Bx’) (-\mu \BJ(\Bx’)) \\
&= \inv{4\pi} \int_V dV’ \frac{\Bx – \Bx’}{ \Abs{\Bx – \Bx’}^3 } (-\mu \BJ(\Bx’)).
\end{aligned}
\end{equation}

We expect the LHS to be a bivector, so the scalar component of this should be zero. That can be demonstrated with some of the usual trickery
\begin{equation}\label{eqn:biotSavartGreens:60}
\begin{aligned}
-\frac{\mu}{4\pi} \int_V dV’ \frac{\Bx – \Bx’}{ \Abs{\Bx – \Bx’}^3 } \cdot \BJ(\Bx’)
&= \frac{\mu}{4\pi} \int_V dV’ \lr{ \spacegrad \inv{ \Abs{\Bx – \Bx’} }} \cdot \BJ(\Bx’) \\
&= -\frac{\mu}{4\pi} \int_V dV’ \lr{ \spacegrad’ \inv{ \Abs{\Bx – \Bx’} }} \cdot \BJ(\Bx’) \\
&= -\frac{\mu}{4\pi} \int_V dV’ \lr{
\spacegrad’ \cdot \frac{\BJ(\Bx’)}{ \Abs{\Bx – \Bx’} }

\frac{\spacegrad’ \cdot \BJ(\Bx’)}{ \Abs{\Bx – \Bx’} }
}.
\end{aligned}
\end{equation}

The current \( \BJ \) is not unconstrained. This can be seen by premultiplying \ref{eqn:biotSavartGreens:20} by the gradient

\begin{equation}\label{eqn:biotSavartGreens:80}
\spacegrad^2 I \BB = -\mu \spacegrad \BJ.
\end{equation}

On the LHS we have a bivector so must have \( \spacegrad \BJ = \spacegrad \wedge \BJ \), or \( \spacegrad \cdot \BJ = 0 \). This kills the \( \spacegrad’ \cdot \BJ(\Bx’) \) integrand numerator in \ref{eqn:biotSavartGreens:60}, leaving

\begin{equation}\label{eqn:biotSavartGreens:100}
\begin{aligned}
-\frac{\mu}{4\pi} \int_V dV’ \frac{\Bx – \Bx’}{ \Abs{\Bx – \Bx’}^3 } \cdot \BJ(\Bx’)
&= -\frac{\mu}{4\pi} \int_V dV’ \spacegrad’ \cdot \frac{\BJ(\Bx’)}{ \Abs{\Bx – \Bx’} } \\
&= -\frac{\mu}{4\pi} \int_{\partial V} dA’ \ncap \cdot \frac{\BJ(\Bx’)}{ \Abs{\Bx – \Bx’} }.
\end{aligned}
\end{equation}

This shows that the scalar part of the equation is zero, provided the normal component of \( \BJ/\Abs{\Bx – \Bx’} \) vanishes on the boundary of the infinite sphere. This leaves the Biot-Savart law as a bivector equation

\begin{equation}\label{eqn:biotSavartGreens:120}
I \BB(\Bx)
= \frac{\mu}{4\pi} \int_V dV’ \BJ(\Bx’) \wedge \frac{\Bx – \Bx’}{ \Abs{\Bx – \Bx’}^3 }.
\end{equation}

Observe that the traditional vector form of the Biot-Savart law can be obtained by premultiplying both sides with \( -I \), leaving

\begin{equation}\label{eqn:biotSavartGreens:140}
\BB(\Bx)
= \frac{\mu}{4\pi} \int_V dV’ \BJ(\Bx’) \cross \frac{\Bx – \Bx’}{ \Abs{\Bx – \Bx’}^3 }.
\end{equation}

This checks against a trusted source such as [1] (eq. 5.39).

References

[1] David Jeffrey Griffiths and Reed College. Introduction to electrodynamics. Prentice hall Upper Saddle River, NJ, 3rd edition, 1999.