The Nonzero External Magnetic Field of Long Solenoids ()

Shiwei Ye^{}

Beijing Computational Science Research Center, Beijing, China.

**DOI: **10.4236/jemaa.2022.142002
PDF
HTML XML
183
Downloads
856
Views
Citations

Beijing Computational Science Research Center, Beijing, China.

A widespread assertion has existed for a long time, believing the external field of an infinitely long solenoid should be zero, but it is proofed to be wrong in this work. The components of magnetic flux density of current-carrying, closely wound cylindrical solenoids are calculated. At a distant field point, the external field definitely has a nonzero component, being equal to that of a straight wire of equal length. Since this equivalence is length-independent, it still holds true for ideal solenoids having infinite length. Hence the incorrect and still spreading inference about long solenoids should be rectified. Furthermore, theoretical and experimental discussions involving solenoids should be reviewed again carefully.

Share and Cite:

Ye, S. (2022) The Nonzero External Magnetic Field of Long Solenoids. *Journal of Electromagnetic Analysis and Applications*, **14**, 19-30. doi: 10.4236/jemaa.2022.142002.

1. Introduction

Scientific research related to magnetic field not only has a long history [1] [2], but also covers a wide range of fields [3]. Solenoids and toroids are widely known as basic circuit structures that could generate magnetic field. It is believed that the external magnetic field generated by an infinitely long, closely wound, current-carrying ideal cylindrical solenoid is zero [4]. Moreover, this inference is used as a basic knowledge in further studies, such as the early research works related to the Aharonov-Bohm effect [1] [2]. However, as will be elucidated in this work, the inference about solenoids is wrong. Besides, a similar incorrect inference about the external magnetic field of closely wound toroids also exists [5]. Therefore, the studies involving solenoids and toroids might be affected by these erroneous inferences.

This paper is devoted to correcting this long-standing and widespread fallacy about long solenoids which is still spreading in textbooks [4] [5], by giving the correct results of the external magnetic field distribution of the ideal solenoids. In this work, the components of the magnetic flux density generated by solenoids are numerically calculated using the Biot-Savart Law. By comparing with the results of numerical calculations, one could notice that there is a universal equivalence between the magnetic field of a solenoid and that of an equally long straight wire located along the solenoid axis, no matter whether the length is infinite. With the length increasing, the calculated field approaches that of an infinitely long solenoid. In this extreme case, the field distribution at distant points turns out to be exactly the same as that of an infinite straight wire.

2. Calculations and Results

2.1. Theoretical Explanations

A simple transforming process could explain intuitively why the external field of current-carrying solenoids might not be zero. For simplicity, the values are calculated and discussed without any other magnetic media nearby. As shown in Figure 1(a), a highly flexible thin wire is closely wound around a cylinder with cross-section radius *R*, forming a helical line with pitch *H*_{0}. The diameter of the wire, *D*_{0}, is smaller than *H*_{0}. A closely wound solenoid requires
${D}_{0}<{H}_{0}\ll R$. When carrying a steady current with intensity *I*, the helical wire forms a solenoid

Figure 1. The structure of a solenoid and the transforming process. (a) The parameters of the solenoid. A flexible thin wire is wound around a cylinder with radius *R*, forming a solenoid with total length *L* and pitch *H*_{0}. (b) The transforming process from a solenoid to a straight wire. The diameter of the wire, *D*_{0}, is smaller than *H*_{0}. The steady current intensity in each wire is *I*. *P* is a distant field point at (0, *d*, *z*_{0}).

with total length *L*. Then the winding turns per unit length of this solenoid is, *n* = 1/*H*_{0}. The axis of the cylinder is set as the *Z* axis. The distance from a field point, *P*, to *Z*-axis is *d*. The perpendicular line from *P* to *Z*-axis is parallel to the *Y* axis.

As shown in Figure 1(b), for an infinitely long solenoid, let the radius *R *decrease gradually and slowly. The solenoid keeps shrinking and getting thinner and thinner. Finally, when its radius *R* becomes less than the radius of the wire (*R* < *D*_{0}/2), the solenoid can be recognized as a straight wire. During this transforming process, the external magnetic field of the solenoid at *P* changes continuously and finally becomes the field of an infinite straight wire, which is nonzero apparently. This indicates that the characteristic of the external field of a solenoid is closely related to that of a straight wire. Especially, it inspires us to speculate that the external field of an infinite solenoid might also be nonzero.

In general, the detailed characteristics of the wire-shape localized within the scale of *R*, should not play an important role in determining the magnetic field at a far away point*P* (
$d\gg R$ ). The overall magnetic field distribution of a very thin solenoid should be very close to that of a simple straight wire. On the other hand, when a straight wire changes its structure (such as bending or twisting) around a certain spot localized within the scale of *R*, the corresponding change of its field could not be reflected significantly at those distant points.

Theoretically, each elementary current *I*d*l* of the wire could be decomposed into three components, d*l** _{z}* (paralleling to

$\text{d}l=\text{d}{l}_{z}+\text{d}{l}_{\tau}$ (1)

The infinitesimal, d*B*, of the field generated by *I*d*l* at a point *P* is given by the Biot-Savart Law,

$\text{d}B=\left(\frac{{\mu}_{0}}{4\pi}\right)I\frac{\text{d}l\times r}{{r}^{3}}$ (2a)

Then the magnetic flux density *B* generated by the current-carrying wire, would be the integration of d*B* along the wire,

$B={\displaystyle \int \text{d}B}=\left(\frac{{\mu}_{0}}{4\pi}\right)I{\displaystyle \int \frac{\text{d}{l}_{z}\times r}{{r}^{3}}}+\left(\frac{{\mu}_{0}}{4\pi}\right)I{\displaystyle \int \frac{\text{d}{l}_{\tau}\times r}{{r}^{3}}}$ (2b)

At a distant point *P*, the first term is approximately equal to the field of a straight wire with the same length of the solenoid, while the second term is approximately equal to that of a series of circular currents with radius *R*, arranged at equal intervals along the *Z* axis. If the total length *L *is infinite, then the second term is zero, while the first term is not. Therefore, the result of *B* at *P* is approximately equal to that generated by an infinitely long straight wire.

As for toroids, a similar assertion exists, believing that their external magnetic field is zero [5]. Based on a similar analysis as the transforming process, it is easy to know the truth that at distant points, the field is similar to that of a magnetic dipole. At nearer points, the field is close to that of a circular current or a donut-shaped current distribution. In this work, we focus on solenoids mainly, while toroids will be discussed elsewhere.

In the following section, the magnetic field distribution of cylindrical solenoids will be discussed with detailed calculations and analyses.

2.2. Parameters and Descriptions of the Calculations

Let the expression of the helical wire be

$\begin{array}{l}x=R\mathrm{cos}\left(\alpha \right)\\ y=R\mathrm{sin}\left(\alpha \right)\\ z=\left(\frac{{H}_{0}}{2\pi}\right)\alpha \end{array}$ (3)

where *α *is an independent variable for determining the position of each elementary current, *I*d*l*. An *I*d*l* could be expressed by the differentials of the coordinates,

$I\text{d}l=I\left(\text{d}x,\text{d}y,\text{d}z\right)=\left(-R\mathrm{sin}\left(\alpha \right),R\mathrm{cos}\left(\alpha \right),\frac{{H}_{0}}{2\pi}\right)I\text{d}\alpha $ (4)

The position of the field point *P* is at (*x _{P}*,

$\begin{array}{l}{x}_{P}=0\\ {y}_{P}=d\\ {z}_{P}={z}_{0}\end{array}$

where *z*_{0} is another variable for determining the position of *P*.

The following calculations are based on the Biot-Savart Law. In Equation (2a), *I* is the current intensity in the wire, *r* is the position vector pointing from elementary current *I*d*l* to *P*, and *r* is the module of *r*.

$r={r}_{P}-{r}_{I\text{d}l}=\left(-R\mathrm{cos}\left(\alpha \right),d-R\mathrm{sin}\left(\alpha \right),{z}_{0}-\frac{{H}_{0}}{2\pi}\alpha \right)$

$r=\left|r\right|=\sqrt{{\left[-R\mathrm{cos}\left(\alpha \right)\right]}^{2}+{\left[d-R\mathrm{sin}\left(\alpha \right)\right]}^{2}+{\left({z}_{0}-\frac{{H}_{0}}{2\pi}\alpha \right)}^{2}}$

The result of the outer product could be expanded as,

$\begin{array}{c}\text{d}l\times r=u\left(\alpha \right)\text{d}\alpha \\ ={\left[\begin{array}{c}i\\ j\\ k\end{array}\right]}^{\text{T}}\left[\begin{array}{c}-\frac{{H}_{0}d}{2\pi}-\frac{\alpha {H}_{0}R\mathrm{cos}\left(\alpha \right)}{2\pi}+R{z}_{0}\mathrm{cos}\left(\alpha \right)+\frac{{H}_{0}R\mathrm{sin}\left(\alpha \right)}{2\pi}\\ -\frac{{H}_{0}R\mathrm{cos}\left(\alpha \right)}{2\pi}-\frac{\alpha {H}_{0}R\mathrm{sin}\left(\alpha \right)}{2\pi}+R{z}_{0}\mathrm{sin}\left(\alpha \right)\\ {R}^{2}-Rd\mathrm{sin}\left(\alpha \right)\end{array}\right]\text{d}\alpha \end{array}$ (5)

where *u*(*α*) is a vector function of variable *α*. The magnetic flux density at *P* is given by the integration of d*B*,

$B={\displaystyle \int \text{d}B}=\left(\frac{{\mu}_{0}}{4\pi}\right)I{\displaystyle \int \frac{\text{d}l\times r}{{r}^{3}}}={\displaystyle {\int}_{{\alpha}_{1}}^{{\alpha}_{2}}f\left(\alpha \right)\text{d}\alpha}$, (6)

where the vector function *f*(*α*) could be expressed as,

$f\left(\alpha \right)=\left(\frac{{\mu}_{0}}{4\pi}\right)I\frac{u\left(\alpha \right)}{r{\left(\alpha \right)}^{3}}$.

In a practical calculation, the integration could be replaced by the following summation with a properly selected Δ*α*,

$B={\displaystyle \underset{\alpha ={\alpha}_{1}}{\overset{{\alpha}_{2}}{\sum}}f\left(\alpha \right)\Delta \alpha}$. (7)

2.3. Results and Comparisons

In the following calculations, the steady current in the wire is known as *I*. The unit of length is selected as the pitch of the helical wire *H*_{0}, while the unit of *B* is
$\frac{{\mu}_{0}I}{4\pi {H}_{0}}$. Then *B* is related to several variables.

$B=B\left(L,d,R,{z}_{0}\right)$. (8)

The meanings of the variables are listed earlier in Section 2.1. In this section, the relation between *B* and each variable is discussed respectively.

The *B _{X}* component of the external field of a solenoid is always compared to that of a straight wire with the same length and current intensity. Since

${B}_{L}=\left|{B}_{L}\right|=\left|{B}_{Lx}\right|=\left(\frac{{\mu}_{0}}{4\pi}\right)\frac{I}{d}\left[\frac{L/2+{z}_{0}}{\sqrt{{\left(L/2+{z}_{0}\right)}^{2}+{d}^{2}}}+\frac{L/2-{z}_{0}}{\sqrt{{\left(L/2-{z}_{0}\right)}^{2}+{d}^{2}}}\right]$, (9)

where *L* is also the length of the wire, while other variables are the same as those mentioned earlier. The theoretical values of *B _{Lx}* are used as the standards to compare with the calculated values of

2.3.1. The Relationship between *B* and *L*, the Total Length of the Solenoid

Figure 2 shows the calculated *B _{x}* of a solenoid and the

It is verified by the calculation that, at a fixed distant field point *P* (
$d=100R\gg R$ ), the *B _{X}* of the solenoid is always equal to the

Figure 2. The components of the external magnetic flux density of a solenoid at a distant point, as functions of the solenoid’s length *L*. (The parameters are as follows, *R* = 100*H*_{0}, *d* = 10,000*H*_{0} = 100*R*, and *z*_{0} = 0. The red line represents the theoretical values of the corresponding straight wire, *B _{Lx}*, while the black marks show the calculated values of

what the value of *L* is. The relative difference between them is less than 1 × 10^{−3}. Moreover, when *L* increases, the absolute value of *B _{X}*

${B}_{L\to \infty}=\frac{{\mu}_{0}I}{2\pi d}=\left(\frac{{\mu}_{0}I}{4\pi {H}_{0}}\right)\frac{2{H}_{0}}{d}=\left(\frac{{\mu}_{0}I}{4\pi {H}_{0}}\right)\times 0.0002$.

For those long enough solenoids (with *L* > 150*R* = 15,000*H*_{0}), the absolute value of *B _{z}* keeps decreasing when

*B _{y}* approximately remains zero at a point with

2.3.2. The Relationship between *B* and *d*, the Distance from Point *P* to the Solenoid Axis

Figure 3 shows the calculated components of *B*, as functions of distance *d*. The fixed parameters are as follows, *R* = 100*H*_{0}, *L* = 20,000*H*_{0} = 200*R*, and *z*_{0} = 0.

The value of *B _{X}* of the solenoid still fits very well to the corresponding

Figure 3. The relation between the external *B* of a solenoid and the distance *d.* (The inset of Figure 3 shows the comparison of *B _{Z}* with corresponding

the corresponding *B _{Lx}* quite well with any

The magnitudes of *B _{X}* and

Again, the magnetic induction lines *projected* on the *YOZ*-plane of a solenoid indicates that *B _{Y}* = 0 at a point with

2.3.3. The Components of *B* as Functions of *z*_{0}

Figure 4 shows the relation between *B* and *z*_{0}, the *z* coordinate of the field point. The fixed parameters are as follows, *R* = 100*H*_{0}, *L* = 20,000*H*_{0} = 200*R*, and *d* = 10,000*H*_{0} = 100*R*.

It is clear that *z*_{0} can affect all three components of *B* significantly. *B _{X}* and

Again, at all the distant field points, the values of *B _{X}* of the solenoid consistently match those of

2.3.4. The Components of *B* with Different Radius *R*

The fixed parameters are as follows, *z*_{0} = 0, *L* = 20,000*H*_{0}, and *d* = 10,000*H*_{0}. As shown in Figure 5, the tightness of winding, represented by (*R/H*_{0}), hardly affects

Figure 4. The components of *B* of a solenoid, as functions of *z*_{0}. (The parameters are as follows, *R* = 100*H*_{0}, *L* = 20,000*H*_{0} = 200*R*, and *d* = 10,000*H*_{0} = 100*R*. The red line represents the theoretical values of the corresponding straight wire, *B _{Lx}*. The relative difference between

Figure 5. The relation between *B* and the solenoid’s radius *R*. (The parameters are as follows, *z*_{0} = 0, *L* = 20,000*H*_{0}, and *d* = 10,000*H*_{0}.)

the value of *B _{X}* at a distant point, though

On the contrary, the radius *R* of the solenoid affects *B _{Z}* more significantly, because the field distribution on the

2.3.5. Summary of the Results and Comparisons

Since adequate discussions have been made for solenoids with finite length, it is possible to infer the results for a very long solenoid.

It has been verified that, no matter how long the solenoid is, the *B _{X}* of its external field always equals that of a straight wire having an equal length and carrying an equal current. This equivalence itself is independent of the value of length

$\left|{B}_{x}\right|=\frac{{\mu}_{0}I}{2\pi d}$ (10)

On the other hand, according to the symmetry of the infinite solenoid, *B _{Y}* equals zero. Since the length

$B=\left|{B}_{\text{ext}}\right|=\left|{B}_{x}\right|=\frac{{\mu}_{0}I}{2\pi d}$. Naturally, the external *B*_{ext} satisfies the Ampere’s law,

$\oint {B}_{\text{ext}}\cdot \text{d}l}={\mu}_{0}I$, (11)

where *I *is exactly the flux of the current density through an arbitrary cross-section of the solenoid, *i.e.*, the current enclosed by the Amperian loop parallel to the *XOY*-plane.

3. The Substitution of Solenoids

We have proofed that, no matter how long it is, a solenoid always has an external nonzero *B _{X}* component. So the problem is finding a substitution of solenoid. The answer is simple. A pair of coaxial but chiral-opposite solenoids can build a field-free area enclosing a nonzero amount of magnetic flux.

As for the Aharonov-Bohm effect [1], Figure 6(a) shows the schematic experiment observing the inference of two beams of electron passing through a field-free area (with *B* = 0) enclosing a nonzero flux amount. Many experiments have been established to measure the results of the Aharonov-Bohm effect, *i.e.*, the change of inference fringes [2] [6]. In fact, to distinguish the effects of *A* and *B*, one can just make these two variables change conversely, in the area where the beams pass through.

Figure 6. A substitution of a solenoid. (a) Schematic experiment to demonstrate the Aharonov-Bohm effect. (b) A pair of two coaxial solenoids with opposite chirality.

3.1. Coaxial, Chiral-Opposite Solenoids Carrying Equal Steady Currents

As shown in Figure 6(b), suppose there are two coaxial solenoids, with opposite chirality and different cross-section radius, carrying steady current with intensity,

${I}_{1}={I}_{2}=I$. (12)

If the *Z*-components of their current density are in opposite directions, then their internal fields are mutually reinforcing, while their external fields are canceling each other out. Thus, the structure of an external field-free area where the beams passing through with nonzero magnetic flux enclosed between the two beams is ideally realized.

3.2. Coaxial, Chiral-Opposite Solenoids Carrying Slowly Varying Currents

The device mentioned in Section 3.1 can be used as a substitution of a solenoid in certain cases, such as the experiments related to the Aharonov-Bohm effect. Though the leakage of field of the solenoid has been discussed in other works [2] [7], only the *B _{Z}* component is considered, while the

In this section, we don’t need *B* to be zero on the routes of the beams. On the contrary, we let both *B* and *A* exist together in the same area, changing conversely. It is only necessary to distinguish which variable can affect the phase difference, the magnetic flux in the area surrounded by the routes related to vector *A*, or the fields along the routes related to vector *B*. Therefore, the existence of the field-free areas is not a necessity, if the changing trend of *B* is opposite to that of *A*.

In the first step, let the current intensity in solenoid-1 be zero, *I*_{1} = 0, and that in solenoid-2, *I*_{2}, increase very slowly. Then the magnetic flux increase gradually and may cause a changing trend *T*_{1} of the interference fringes.

In the second step, fix *I*_{2} and let *I*_{1} increase very slowly. The magnetic flux generated by both currents increases gradually, while the field outside solenoid-2 decreases gradually. The changing trend of the interference fringes in this step is noted as *T*_{2}.

By comparing *T*_{1} and *T*_{2}, it could be distinguished whether the change of interference fringes is caused by the magnetic flux or by the magnetic field itself. In these steps, the changing trends of the flux are the same, while those of the field intensity are opposite. If *T*_{1} and *T*_{2} always show the same tendency, then the change of flux dominates this process. Otherwise, the change of field intensity may have played an important role in this process.

4. Summary

The assertion that the external field of a very long solenoid is zero has been proved to be wrong. Although it directly violates the Ampere’s law and the Biot-Savart Law, this fallacy is still spreading widely.

The occurrence and spread of this fallacy might be due to the following reasons. A possible theoretical reason might be the improper substitution of a series of coaxial circular currents for the original solenoid, which changes the symmetric properties of the system and totally ignores the contribution of the d*l _{z}*-term to

${B}_{\text{ext}}=\frac{{\mu}_{0}I}{2\pi d}<\frac{{\mu}_{0}I}{2\pi R}\ll \frac{{\mu}_{0}I}{2\pi {H}_{0}}<\frac{{\mu}_{0}I}{{H}_{0}}={B}_{\mathrm{int}}$

The difference of intensity would become more significant when ferromagnetic medium is inserted inside the solenoid. Besides, only the magnetic induction lines *projected* on the *YOZ*-plane of solenoids are concerned and shown in experiments, while those *projected* on the *XOY*-plane are ignored.

In this work, the external magnetic field flux density of solenoids is calculated numerically. The results of *B _{X}* are compared with the

Acknowledgements

The author would like to thank Professor Bang-jiao Ye for the valuable guidance and insightful comments, and Professor Si-zhao Jin for the long-term, tireless help and support.

Conflicts of Interest

The author declares no conflicts of interest regarding the publication of this paper.

[1] |
Aharonov, Y. and Bohm, D. (1959) Significance of Electromagnetic Potentials in the Quantum Theory, Physical Review, 115, 485-491.
https://doi.org/10.1103/PhysRev.115.485 |

[2] |
Chambers, R.G. (1960) Shift of an Electron Interference Pattern by Enclosed Magnetic Flux. Physical Review Letters, 5, 3-5. https://doi.org/10.1103/PhysRevLett.5.3 |

[3] |
Drochon, A., Robin, V., Fokapu, O. and Rodriguez, D. (2016) Stationary Flow of Blood in a Rigid Vessel in the Presence of an External Magnetic Field: Considerations about the Forces and Wall Shear Stresses. Applied Mathematics, 7, 130-136.
https://doi.org/10.4236/am.2016.72012 |

[4] | Feynman, R.P., Leighton, R.B. and Sands, M.L. (2010) The Feynman Lectures on Physics, Vol. 2. Chapter 14-4. Addison-Wesley, Reading. |

[5] | Griffiths, D.J. (2012) Introduction to Electrodynamics. Addison-Wesley, New York, Chapter 5. |

[6] |
Tonomura, A., et al. (1986) Evidence for Aharonov-Bohm Effect with Magnetic Field Completely Shielded from Electron Wave. Physical Review Letters, 56, 792-795.
https://doi.org/10.1103/PhysRevLett.56.792 |

[7] |
Roy, S.M. (1980) Condition for Nonexistence of Aharonov-Bohm Effect. Physical Review Letters, 44, 111-114. https://doi.org/10.1103/PhysRevLett.44.111 |

Journals Menu

Contact us

+1 323-425-8868 | |

customer@scirp.org | |

+86 18163351462(WhatsApp) | |

1655362766 | |

Paper Publishing WeChat |

Copyright © 2023 by authors and Scientific Research Publishing Inc.

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.