Analytic Solution

Purpose

Here, Maxwell’s equations are solved for a harmonic electrical current dipole source. This is accomplished by using the method of Schelkunoff potentials, as shown in Ward and Hohmann ([WH88]). Analytic expressions for the electric field, the magnetic field and the corresponding vector potential are provided. Numerical modeling tools for visualizing the fields are provided after the asymptotics section.

For an electrical current source (\(\mathbf{J_e^s}\)), Maxwell’s equations in the frequency domain are given by:

(197)\[\nabla \times \mathbf{E_e} + i\omega \mu \mathbf{H_e} = 0\]
(198)\[\nabla \times \mathbf{H_e} - (\sigma + i\omega \varepsilon) \mathbf{E_e} = \mathbf{J}_e^s\]

where subscripts \(_e\) remind us that we are considering an electric source. For an electrical current source (\(\mathbf{J_e^s}\)), the electric and magnetic fields can be expressed in terms of the Schelkunoff vector potential (\(\mathbf{A}\)), where:

(199)\[\mathbf{H_e} \equiv \nabla \times \mathbf{A}\]

and

(200)\[\mathbf{E}_e = -i\omega\mu\mathbf{A} + \frac{1}{(\sigma + i\omega\varepsilon)} \nabla (\nabla \cdot \mathbf{A})\]

Eq. (199) can be obtained simply by taking the divergence of Eq. (197). Eq. (200) is obtained by manipulating Eqs. (197), (198) and (199), and choosing an appropriate Gauge. We can see from Eqs. (199) and (200) that \(\mathbf{A}\) contains all the information corresponding to the electric and magnetic fields. Therefore, Maxwell’s equations will be manipulated to solve for \(\mathbf{A}\); which can then be used to obtain \(\mathbf{E_e}\) and \(\mathbf{H_e}\).

By manipulating Eqs. (197), (198) and (199) and choosing an appropriate Gauge, we find that \(\mathbf{A}\) can be expressed using the Helmholtz equation:

(201)\[\nabla^2 \mathbf{A} + k^2 \mathbf{A} = - \mathbf{J}_e^s, \ \ \ \ \text{where} \ \ k^2 = \omega^2\mu\epsilon -i\omega\mu\sigma\]

The Helmholtz equation with boundary conditions can be solved to generate \(\mathbf{A}\). For infinite media, the boundary condition is such that \(\mathbf{A} \rightarrow 0\) as \(r \rightarrow \infty\). From the Helmholtz equation, we can see that \(\mathbf{A}\) will only have a component along the direction of \(\mathbf{J_e^s}\). The scalar Green’s function for the Helmholtz equation is:

(202)\[G(r) = \frac{e^{-ikr}}{4\pi r}.\]

and hence the vector potential for an arbitrary electric current source is:

(203)\[\mathbf{A}(\mathbf{r}) = \int_{V^\prime} \frac{e^{-ik|\mathbf{r}-\mathbf{r}'|}}{4\pi |\mathbf{r}-\mathbf{r}'|} \mathbf{J_e^s}(\mathbf{r}') dV^\prime\]

where \(\mathbf{r}\) is the observation location, \(\mathbf{r^\prime}\) refers to locations within the source region and \(V^\prime\) is the volume of the source region. For an electric current dipole oriented in the \(\mathbf{\hat{x}}\) direction, the source term is given by:

(204)\[\mathbf{J_e^s} = \mathbf{\hat{x}} I ds \delta(x) \delta(y) \delta(z)\]

and the solution to Eq. (203) is:

(205)\[\mathbf{A} = \frac{I ds}{4\pi r} e^{-ikr} \mathbf{\hat{x}}\]

Recall the \(\mathbf{A}\) can be used to obtain the electric and magnetic field according to Eqs. (199) and (200). Thus the electric field for an electrical current dipole in the \(\mathbf{\hat x}\) direction is:

\[\mathbf{E_e} = \frac{I ds}{4 \pi (\sigma + i \omega \varepsilon)} \left[ \left( k^2 + \frac{\partial^2}{\partial x^2} \right) \mathbf{\hat{x}} + \frac{\partial^2}{\partial x \partial y} \mathbf{\hat{y}} + \frac{\partial^2}{\partial x \partial z} \mathbf{\hat{z}} \right] \frac{e^{-ikr}}{r}\]

which is equal to:

(206)\[\begin{split}\begin{split} \mathbf{E_e} = \frac{I ds}{4 \pi (\sigma + i \omega \varepsilon) r^3} e^{-ikr} \Bigg [ \Bigg ( \frac{x^2}{r^2} \mathbf{\hat{x}} + & \frac{xy}{r^2} \mathbf{\hat{y}} + \frac{xz}{r^2} \mathbf{\hat{z}} \Bigg ) ... \\ &\big ( -k^2 r^2 + 3ikr +3 \big ) + \big ( k^2 r^2 - ikr -1 \big ) \mathbf{\hat{x}} \Bigg ] . \end{split}\end{split}\]

The magnetic field is:

\[\mathbf{H_e} = \frac{I ds}{4 \pi} \left[ \frac{\partial}{\partial z} \mathbf{\hat{y}} - \frac{\partial}{\partial y} \mathbf{\hat{z}} \right] \frac{e^{-ikr}}{r}\]

which is equal to:

(207)\[\mathbf{H_e} = \frac{I ds}{4 \pi r^2} \left( ikr + 1 \right) e^{-ikr} \left( -\frac{z}{r} \mathbf{\hat{y}} + \frac{y}{r} \mathbf{\hat{z}} \right) .\]

On the following page, we show how Eqs. (206) and (207) can be simplified for various cases.