⬇ Download Markdown  ·  ⬇ Download LaTeX  ·  ⬇ original  ·  view original

Self inductance of a wire loop as a curve integral

arXiv:1204.1486v2 [physics.class-ph] 15 Sep 2013

R. Dengler

Rohde& Schwarz GmbH & Co KG, 1GP1 Mühldorfstr. 15, 81671 Munich, P.O.B. 801469.

September 17, 2013

Abstract

It is shown that the self inductance of a wire loop may be written as a curve
integral akin to the Neumann formula for the mutual inductance of two wire loops.
The only difference is that contributions where the two integration variables get too
close to each other must be excluded from the curve integral and evaluated in detail.
The contributions of these excluded segments depend on the distribution of the
current in the cross section of the wire. They add to a simple constant proportional
to the wire length. The error of the new expression is of first order in the wire radius
if there are sharp corners and of second order in the wire radius for smooth wire
loops.

1 Introduction

Electrical inductance plays a crucial role in power plants, transformers and electronic devices. The coefficients of self and mutual inductance required to quantitatively describe inductance belong to the field of magnetostatics. Calculating inductance coefficients with analytic techniques is, however, impossible except in simple cases. The mathematical reason for the difficulty is that the Laplace equation allows analytic solutions only for some symmetric constellations. There thus are only a few closed-form expressions for these coefficients. In practice one often is forced to use approximations, finite element methods or other numerical techniques. The situation simplifies when the current flows in thin wires. This situation is analogous to an electrostatic system of point charges, where electric field and electrostatic energy directly follow from the given charge distribution, while in a generic system charge and current distributions also are unknown at the outset. The purpose of this article is to derive a new expression for the self inductance of a wire loop, giving self inductance as a curve integral similar to the Neumann formula for mutual inductance. The starting point is the expression W = µ0 8π ˆ j (x) · j (x′ ) |x − x′ | d3 xd3 x′ (1) for the magnetic field energy of a system with current density j (x), where µ0 is the magnetic constant.[1] This expression essentially was already given by Neumann in 1845.[2] 1

It resembles the expression for gravitational or electrostatic potential energy, the only new ingredient is the scalar product between the current elements. For a current density j (x) = Im j m (x) corresponding to N separate current loops with currents Im and normalized current densities jm it follows W = µ0 8π N m,n=1 Im I n ˆ jm (x) · jn (x′ ) |x − x′ | d3 xd3 x′ != 1 2 N m,n=1 Lm ,nI m I n . (2) If the currents flow in thin wires, then the integrals become curve integrals, and one immediately reads off the Neumann expression for mutual inductance of two (filamentary) current loops[2] L1 ,2 = µ0 4π ˛ dx1 · dx2 |x1 − x2 | . (3) It is plausible that there exists a similar expression for the self inductance of a wire loop, but we were not able to find any hint in the literature. Formally one might read off from equation (2) an expression similar to equation (3), where the two closed curves coincide. But this cannot be correct, because |x − x′ | now vanishes and the integral isn’t defined. Instead we will prove L = µ0 4π ˛ dx · dx′ |x − x′ | |s−s′ |>a/2 + µ0 4π lY + ... (4) where a denotes the wire radius and l the length of the wire. The variable s measures the length along the wire axis. The constant Y depends on the distribution of the current in the cross section of the wire: Y = 0 if the current flows in the wire surface, Y = 1/2 when the current is homogeneous across the wire. The ellipses represents terms like O (µ0 a ) and O (µ0 a 2 /l), which are negligible for l ≫ a. In hindsight it is completely natural to use a cutoff of order a in the curve integral. In fact, the exact value of this cutoff is arbitrary, because the contribution proportional to lY also depends on this cutoff. The simplest way to determine Y would be to compare the expression with the self inductance of a long rectangle.

2 Simple derivation

Consider equation (2) with N = 1 for a thin wire with circular cross section, radius a and length l. Let s denote the length along the axis of the wire. The planes perpendicular to the wire axis then define a projection from the bulk of the wire onto the axis, x → s (x) . Selecting a length scale b satisfying a ≪ b ≪ l allows to write L = (µ0 / 4π) L + L with L = ˆ j (x) j (x′ ) |x − x′ | d3 xd3 x′ |s−s′ |>b , (5) L = ˆ j (x) j (x′ ) |x − x′ | d3 xd3 x′ |s−s′ |<b . The second part contains contributions from all point pairs {x, x′ } with a distance along the axis smaller than b, the first the complement of this set (s is a cyclic quantity). For 2

Figure 1: A section of a wire with radius a, with a segment of length 2b, and a plane perpendicular to the wire axis at the center of the segment. given x the planes at s (x) ± b delimitate the points x′ contributing to the first or second integral, see figure (1). L now approximately becomes a curve integral and L essentially consists of cylinders of length 2b. The strategy then is to replace L with the curve integral and to explicitly evaluate the contribution of the cylinders in L. The cylinders are long in comparison to the radius because of a ≪ b and straight (at least most of them) be cause of b ≪ l. Actually the only requirement for the lengths is a ≪ l, the length b = √ al then satisfies a ≪ b ≪ l. The approximation thus is exact in the limit a ≪ l except in special cases. Inserting L0 for a straight segment from equation (A.4) in the appendix thus leads to L = µ0 4π ˛ dx · dx′ |x − x′ | |s−s′ |>b + µ0 l 2π ln 2b a + Y 2 + ... (6) This expression cannot depend on the (more or less) arbitrary length scale b. The curve integral thus is1 L (b) = const − µ0 l 2π ln (2b/b0 ) . But b is the only “short” length scale in the curve integral, and L (b) thus also is valid for b = a/2. The expression (6) therefore doesn’t change if one formally sets b = a/2. Equation (6) now agrees with equation (4), but some questions remain. First of all, how accurate is formula (4)? The curve integral is a purely geometric quantity with dimension “length” and order of magnitude l. Plausible expressions for the order of magnitude of the relative error are a/l, (a/l)2 and (a/R)2 , with R a typical curvature radius of the wire loop. Errors of this order normally are negligible (and also occur in the Neumann formula for mutual inductance). But the derivation of formula (4) is not as straightforward, and so what are the actual limits or exceptions?

3 Examples and comparison with exact self induc-

tance

To get an impression of the accuracy we have compared self-inductances calculated with the curve integral with the result of a numeric evaluation of the nominally 6-dimensional integral in equation (2). This integral becomes 4-dimensional if the currents flow in the wire surface (skin effect, Y = 0), and the results below correspond to this skin effect case. The order of magnitudes of the error terms identified here are corroborated below in a more detailed derivation of formula (4). 1 This also can easily be verified with an explicit calculation. 3

Straight segment

The first example is a straight segment with length c and complete skin effect. This of course isn’t a closed circuit, but it might be an edge of a rectangle. The orthogonal edges of the rectangle don’t interact with the segment because of the scalar product j · j′ . What is missing for a rectangle are the interaction terms with the opposite edges (and the small contributions from the corners). In this case the volume integral (2) may even be evaluated analytically with the result L = µ0 4π 2c ln 2c a − 1 + 8a/π − a2 /c + ... , while the curve integral (4) leads to Lc (c) = µ0 4π 2c ln 2c a − 1 + a . (7) The difference is of order O (µ0 a ), much smaller than µ0 c for c ≫ a.

Circular loop

The next example is a ring with radius R. The curve integral (4) gives Lc = µ0 R (ln (8R/a) − 2 + Y /2) + µ0 O a2 /R . This expression also may be found in the literature, derived with the help of elliptic functions and some approximations in a much more complicated way. The table displays the ratio of the exact inductance and Lc for some values of R/a, R/a 1 2 3 5 10 20 L/Lc 5.971 1.238 1.092 1.031 1.00756 1.00188 The expression Lc is more accurate than one might expect. It gives a reasonable approxi- mation already for R = 3a, and the error roughly decays like O (a2 /R2 ).

Rectangle

This case is more complicated because in principle also the shape of the corners comes into play (curvature radius?). But the simplest thing to do is to evaluate the curve integral (4) for a rectangle with edges of length c and d. Orthogonal edges decouple because of the scalar product j · j′ , and the first contribution are the terms (7) for the four edges by themselves. The second contribution are the parts of the curve integral (4) with x on one edge and x′ on the opposite one. The condition |s − s′ | > a/2 is irrelevant for these cross terms, and one easily obtains for parallel edges of length c and distance d Lc (c, d) = µ0 4π 4 √ c2 + d2 − 4d − 4c asinh (c/d) , and the sum together with the Y -term of equation (4) is Lc = µ0 π { c ln 2c a + d ln 2d a − (c + d) (2 − Y /2) (8) +2 √ c2 + d2 − c asinh (c/d) − d asinh (d/c) + a }. 4

This expression also may be found in the literature, with sometimes a factor 2 at the a-term.[3] The table displays the ratio of the numerically evaluated self-inductance L of a square with border length c and corners with curvature radius a and the curve integral Lc for different border length c, c/a 5 10 20 40 80 160 L/Lc 1.168 1.056 1.0205 1.00776 1.00345 1.00144 (L − Lc ) /µ0 a 0.501 0.562 0.5871 0.57826 0.63559 0.63185 L/Le xact c 1.053 1.021 1.0080 1.0029 1.00147 1.00061 The curvature radius a is minimal in that the centre of curvature lies on the inner border of the wire. It is remarkable that the absolute error nearly remains constant. The last row of the table displays the ratio of the exact self inductance and the exact curve integral (4) (with round corners), also evaluated numerically. This expression is a better approximation for small c/a, where the square with round corners degenerates to a ring.

Equilateral triangle

The curve integral (6) for an equilateral triangle with edge length c consists of three times the expression (4) for the edges by themselves and three times the interaction energy Lc (c, c, 120) of adjacent edges (with s on one edge and s′ on the other, see appendix C. There is no such interaction for rectangles because of the scalar product), Lc = µ0 2π 3c ln c a − 1 − ln 3 2 . The table displays the ratio of the exact self-inductance L of an equilateral triangle with border length c and corners with curvature radius a and the curve integral Lc for different border length c, c/a 5 10 20 40 80 160 L/Lc 2.983 1.2534 1.0743 1.0264 1.0102 1.0043 (L − Lc ) /µ0 a 0.966 1.0856 1.1288 1.1501 1.1599 1.2017 L/Le xact c 1.376 1.0931 1.0302 1.0110 1.0043 1.0019 The absolute error (L − Lc ) /µ0 is nearly constant also here. The last row again is the ratio of the exact self inductance and the curve integral (4) with round corners, evaluated numerically.

Parallel wires

For a loop consisting of infinitely long parallel wires the condition a ≪ l is perfectly met and the curve integral gives the exact self inductance even for minimal distance d = 2a, Lc = µ0 l π ln d a + Y /2 . This expression is the limiting case of the expression (8) for a long rectangle. The point is, that the contribution of the corners becomes negligible for a long rectangle, and that the replacement of the (circular symmetric) current by a line current doesn’t change the 5

Figure 2: A wire loop and distances relevant for calculating the energy of a disc at s = 0 relative to other points in the wire. magnetic field according to Ampere’s law. Of course, the assumption of a circular sym- metric current distribution gets wrong in the skin effect case if the wires are close to each other because of additional screening currents. To summarize, formula (6) is rather accurate even for circuits with a linear extension as small as 20 times the wire radius, even if the circuit contains sharp corners.

4 Error estimation

According to equation (5) the self inductance may be written as L = (µ0 / 4π) L + L , where L contains the short segments and L the complement. The self inductance L of course doesn’t depend on the arbitrary segment length b. Let us now introduce some notation. We use a coordinate system {s, r, φ} in the wire where the length s along the wire axis is cyclic with period l, and the coordinates {r, φ} describe planes perpendicular to the wire axis. The intersections of the planes and the wire are assumed to be circular, 0 ≦ r ≦ a and 0 ≦ φ ≦ 2π. The volume element reads dV = (1 + r cos φ/R (s)) rdrdφds, where R (s) denotes the curvature radius of the wire, and φ = 0 at the outer border of the wire (this is possible at least locally). The volume element also may be written as dV = dsdA with dA = (1 + r cos φ/R (s)) dA and area element dA = rdrdφ. The coordinates become cylindrical coordinates for straight wire segments, i.e. for R = ∞. The current density j is normalized, that is ´ dA |j| = ´ dA |j| = 1. We will also need the radial moments an = ⟨rn ⟩ = ˆ dA rn |j| of the current distribution. In the skin effect case of course an = an . One quantity of interest then is L (s) = ˆ ds′ θ (|s′ − s| − b) ˆ dA dA ′ j (s, r, φ) · j (x′ ) |x (s, r, φ) −x (s′ , r′ , φ′ ) | , (9) the energy of the current in the plane at s = 0 with respect to the current at |s′ − s| > b, see figure (2). The symbol θ denotes the step function. To obtain L from L (s) requires 6

to integrate over s. The curve integral

to integrate over s. The curve integral Lγ (s) = ˛ ds′ θ (|s′ − s| − b) cos (j (s) , j (s′ )) |x (s, 0, 0) −x (s′ , 0, 0) | , (10) is an approximation for L (s). Similarly we write the short segment around the plane at s = 0 as L (s) = ˆ ds′ θ (b − |s′ − s|) ˆ dA dA ′ j (s, r, φ) · j (s′ , r′ , φ′ ) |x (s, r, φ) −x (s′ , r′ , φ′ ) | , (11) These definitions allow to write 4π µ0 L (s) = L (s) + L (s) = (12) = Lγ + Lγ + L − Lγ + ˆP0 + L − L0 s , where we have added and subtracted the curve integral Lγ (s), the segment integral L0 (s) for a straight segment from equation (A.4) and the approximation Lγ (s) = 2 (ln (2b/a) + Y /2) for L0 (s). The first bracket in equation (12) now is formula (6), and the two other brackets thus represent the error. The second bracket contains the difference of the volume and the curve integral (for |s′ − s| > b) plus the power series ˆP0 = L0 − Lγ from equation (A.4), the third bracket is the difference of the actual segment integral (11) and the segment integral for a straight segment. It is evident that the error becomes small in suitable limits, and we now want to determine the order of magnitude of the error.

Smooth current loops

We first consider smooth current loops, that is current loops with a minimal curvature radius R comparable with the system size. We also assume that the loop returns imme- diately and doesn’t touch itself anywhere in between. Such complications are considered below. The coordinates x and x′ in the integral (9) above may be expanded like x = x (s, 0, 0)+ x1 (s, r, φ), and thus x − x′ = xs ,s′ + x1 − x′1 , (x − x′ )2 = x2s ,s′ + 2µxs ,s′ + ν 2 with xs ,s′ = x (s, 0, 0) − x (s′ , 0, 0) the distance of the projections onto the axis and |x1 | and |x′1 | of order O (a). The abbreviations are µ = xs s′ · x1 + xs ′ s · x′1 , ν 2 = (x1 − x′1 ) 2 . The procedure now is to use the multipole expansion 1 |x − x′ | = 1 xs ,s′ − µ x2 s,s′ + 1 2x3 s,s′ 3µ2 − ν 2 − 1 2x4 s,s′ 5µ3 − 3µν 2 (13) + 1 8x5 s,s′ 35µ4 − 30µ2 ν 2 + 3ν 4 + O a5 x− 6 s,s′ . 7

This expansion converges for xs ,s′ ≧ |x′1 − x1 | , the coefficients are the coefficients of the Legendre polynomials. Inserting the leading monopole term 1/xs ,s′ of equation (13) into equation (9) reproduces the curve integral Lγ (s). The higher multipole terms describe the difference between the volume integral L (s) and the curve integral. The s′ -integral in the multipole terms converges and the difference thus is mainly a local quantity. The no minal order of magnitude of the multipole terms Ln is an /bn , n ≧ 1. With a length b = √ aR the order of magnitude becomes (a/R)n /2 , and the expansion up to the hexadecupole (n = 4) is needed to get the (a/R)2 approximation. The volume element is (1 + r cos φ/R (s)) rdrdφds = 1 + x1 · R (s) /R2 (s) rdrdφds, where R (s) is the local curvature radius vector. Integrals over φ and φ′ then may be evaluated with the help of ⟨x1 ⟩ = 0, ⟨(x1 ) m (x1 ) n⟩ = r2 Pm ,n/ 2, where Pm ,n is the projection operator projecting onto the plane perpendicular to the wire axis. This implies P · R = R. Inserting now the multipole expansion (13) into the volume integral (9) generates an expansion of the difference of the volume and the curve integral in the region |s′ | > b. Dipole The voltage drops by the same amount along inner and outer border along curved parts of the loop. Electric field and current density thus are larger at the inner border, and the current density thus depends on r and ϕ. This ϕ-dependence of the current density compensates the ϕ-dependence 1 + r cos (ϕ) /R from the volume element. The simple result then is: there is no dipole contribution. The factors cos ϕ and cos ϕ′ from µ give 0 after integration over the angles. There are in fact short transition regions between straight and curved parts of the loop where the current distribution changes from uniform to non-uniform, but this only leads to corrections of higher order. Quadrupole After integration over the angles φ and the radial coordinates r there remains L2 (s) = a2 ˛ ds′ θ (|s′ − s| − b) 3 4 2 − (n · xs ,s′ )2 − (n′ · xs ,s′ )2 − 1 cos α x3 s,s′ , where n denotes a unit vector in the direction of the wire axis and the expression P = 1−nn for the projection operator was used. In L2 (s) one may recognize 1 − (n · xs ,s′ )2 = sin2 ψ and 1 − (n′ · xs ,s′ )2 = sin2 ψ′ , where ψ denotes the angle between the distance vector xs ,s′ and the wire axis. With ψ′ = ψ = 2α = (s − s)′ /2R for a smooth current loop there remains an er ror (a2 /R2 ) ln (R/b). The −a2 / b2 term (from the −1) is peculiar. With the choice b = √ aR it would be of order a/R, but it gets cancelled against the leading term of the power series ˆP0 of equation (12). This of course is the reason for combining ˆP0 with the multipole expansion in equation (12): for a straight segment the formula is exact, and the third bracket in equation (12) vanishes. The multipole expansion together with ˆP0 thus also vanishes for R = ∞. 8

Oktupole The oktupole contributes at most a term of order a3 /b3 . But the cos φ and sin ψ from the odd power of µ make the actual contribution smaller than the expected ∼ a2 /R2 . Hexadecupole The hexadecupole term is of order a4 /b4 ∼ a2 /R2 and its leading part must be kept. The µ factors contain a factor sin ψ ∼ b/R and may be dropped. The ν 4 leads to L4 (s) = 3 8 ˛ ds′ θ (|s′ − s| − b) 2a4 + a22 3 + cos2 α cos α x5 s,s′ . The factors cos α may be replaced with 1 because of α = (s′ − s) /R is small. The integral converges and contributes an error like a2 /R2 . But this term gets cancelled by the second term of the power series ˆP0 from equation (A.4). For smooth current loops it finally follows 4π µ0 L (s) = Lγ + Lγ + O a2 /R ln R b + L − L0 , (14) where the first bracket on the r.h.s is formula (6), evaluated with b = √ aR. The remaining segment integral L for a segment with curvature radius R is evaluated in appendix B. The result (B.1) contributes another logarithmic error (a2 /R2 ) ln (R/b), and a larger error of order b2 /R2 ∼ a/R. We now drop all errors of order O ((a2 /R) ln (R/b)) and O (a2 /R) from the multipole expansion and the curved segment . As in the simple derivation of the formula above one may notice that there is no small lenght scale in the formula ¯Lγ + Lγ and its b-dependence is under control for all b < √ aR and given by (B.2). It contains a “large” b2 /R2 ∼ a/R term. The essential point now is that this b-dependence of the formula exactly absorbs the remaining large error from L (the cancellation comes about, because b is the lower limit in the curve integral and the upper limit in the volume integral). Nothing therefore changes if we now set b = a/2, except that there only remains formula (6) evaluated for b = a/2 and small terms of order O ((a2 /R) ln (R/a)) and O (a2 /R). For a circular loop the leading error terms can even be evaluated in closed form. The results for the self-inductances in the skin effect and volume current case are LS = µ0 R 1+ 3 4 a2 R2 ln 8R a − 2 − 3 2 a2 R2 , LV = µ0 R 1+ 3 8 a2 R2 ln 8R a − 1.75 − 2 3 a2 R2 . The corrections perfectly agree with a numeric evaluation of the 6- or 4-dimensional in- tegrals, but strongly disagree with expressions in the literature.[4] The reason appears to be that these calculations assume a constant current distribution in the surface or cross section of the wire, wrong just in the curved parts of the loop.

Current loops with sharp corners

The errors originate from the curved parts of the current loop, and for current loops with
sharp corners one may expect larger errors. A simple way to construct such loops is to
9

insert straight segments into a circular loop. It was shown above that the absolute error for a circular loop with radius R is of order O (µ0 a 2 /R). This estimation is valid even for minimal curvature radius R = a (the inductance has dimension length × µ0 and a is the only available length for such loops). Since the loop with inserted straight segments is better approximated by the curve integral the absolute error can only be smaller. The formal reason is that a filamentary straight segment generates the same magnetic field as the actual axially symmetric current distribution. The inductance of a loop of extension l ≫ a generally is of order O (µ0 l · ln (l/a)). The ratio with the absolute error leads to a generic estimation of the relative error of formula (4) ∆L/L = n O a2 lRn , where the index n enumerates the corners and the logarithmic factor is neglected. We have allowed here corners with different curvature radius Rn . The special case with curvature radii of order O (l) leads back to the estimation of the error for smooth current loops above. Sharp corners with curvature radius R = O (a) contribute a relative error of order O (a/l). Corners with a small angle of course should get a smaller weight in the sum. This estimation could be made more rigorous, but the procedure only is circumstantial and of no interest here.

Current loops not closing immediately

The error estimations above fail if the current loop comes close to itself before it actually closes, that is for tight spirals, coils or something like that (the estimation xs ,s′ ∼ |s′ − s| gets invalid). In this case the replacement of the actual current with a filamentary cur- rent generates additional errors at these positions. But this isn’t specific for equation (4), exactly the same (small) errors are contained in the Neumann formula for mutual inductance. The error may easily be estimated, because the condition |s′ − s| > b is irrelevant if s′ is on one winding and s on another. A simple possibility is to consider as a worst case scenario two current loops of radius R a distance d > 2a on top of each other. The error comes from the multipole expansion (13) valid everywhere now, and it is a simple matter to estimate the integrals. The order of magnitude of the relative error is O (a2 /Rd), small except for small curvature radius and small distance (both of order O (a)).

5 Conclusions

The curve integral (4) for the self inductance of a wire loop is only a little bit more complicated than the Neumann formula for the mutual inductance of two wire loops. The exact expression for self inductance is a 6-dimensional integral with a logarithmic divergence and several length scales. Nevertheless clear statements follow for the accuracy of formula (4), for loops consisting of straight segments as well as for smooth loops. The error originates from the curved parts of the loop, and is of order µ0 a or µ0 a 2 /l, negligible for most practical purposes. The leading error presumably even might be given as a sum over the corners for loops consisting of straight segments or as additional curve integrals 10

for smooth current loops (see the quadrupole contribution above). The techniques used for error estimation also may be used for the Neumann formula. The self inductance curve integral can be evaluated analytically in many cases, for instance for current loops consisting of coplanar straight segments (not a new result; see also appendix C and D). But equation (4) is valid for abitrary curves, and the numerical evaluation of two-dimensionals integrals with a computer program is a breeze with appro- priate numerical libraries. The information for the self inductance is contained in the curve spanned by the current loop, and any self inductance calculation at least requires a double integral along the curve. Simpler methods or estimations based on “partial inductance” or only the magnetic flux miss this point and only may work in special cases. The fact that two coinciding points cause problems in self inductance calculations is well known and has been circumvented in several ways, for instance by distributing the current onto two filamentary loops. But no systematic approximations can be obtained in this way. Current distributions which are not circular symmetric also lead to formula (4), with a cutoff and a constant Y depending on the current distribution. An example are circuits consisting of coplaner flat strips of width w. The self inductance of such circuits is given in the accuracy described above by formula (4) with a = w and Y = 3.

A Contribution from straight segments of length 2b

The contribution L to the self inductance in equation (5) is due to the interaction of the current in the plane s with the current in all planes s′ with |s′ − s| < b. This value depends on the current distribution in the wire and on the wire geometry within the segment [s − b,s + b], but may be evaluated if the segment is straight or slightly curved. This s-dependent value of course still is to be integrated over all s. To get an approximation for L in the straight wire case use cylinder coordinates with a length s along the axis and area element dA = rdrdφ (see figure (1)). This leads to

˛

L0 = dsL (s) , (A.1) L0 (s) = ˆ j (r) j (r′ ) |x (s, r, φ) −x′ | ds′ dA′ dA |s(x′ )−s|<b . In the latter integral x extends over the plane through the centre of a cylindrical segment, x′ extends over the complete segment. The integral L0 (s) of course is independent of s. The integral over s′ (from −b to b) may be performed using |x − x′ |2 = N 2 + s′ 2 , N 2 = r2 + r′ 2 − 2rr′ cos ( φ − φ′ ),

ˆ

L0 (0) ∼ = 2 ˆ dAdA′ j (r) j (r′ ) asinh (b/N )

ˆ

= 2 ˆ dAdA′ j (r) j (r′ ) {ln (2b/a) − ln (N/a) + A1 (N/b) + ...} In the second line the expansion asinh (x) = ln (2x) + A1 (1/x) (A.2) A1 (x) = ∞ n=1 1 · 3... (2n − 1) 2 · 4...2n (−1)n +1 2n x2 n = x2 /4 − 3x4 /32 + ... 11

was used. The expansion converges because of N = O (a) ≪ b. It doesn’t matter which φ′ occurs in the φ-integral and thus we now set φ′ = 0. Because of ´ dA |j| = 1 the leading term simply becomes 2 ln (2b/a). The second term follows from ln (N/a) = 1 2 ln (ρ2 + ρ′ 2 − 2ρρ′ cos φ) with ρ = r/a and ρ′ = r′ /a and 1 2π ˆ 2π 0 ln ρ2 + ρ′ 2 − 2ρρ′ cos φ dφ = 2 ln (ρ> ) , (A.3) where ln (ρ> ) = θ (ρ − ρ′ ) ln ρ + θ (ρ′ − ρ) ln ρ′ . This term thus vanishes in the skin effect case where the current differs from 0 only for ρ = ρ′ = 1. The current density in the constant current case is j (r) = 1/ (πa2 ) and the second term becomes − (2π)2 π2 2 ˆ 1 0 dρρ ˆ 1 0 dρ′ ρ′ ln (ρ> ) = 1/2. A rapidly convergent expansion approximation for b ≫ a thus is L0 (0) = Lγ (0) + P0 (0) , (A.4) Lγ (0) = 2 ln (2b/a) + Y, P0 (0) = 2A1 N b = a2 b2 − 3 8b4 a4 + 2a22 + O a6 b6 with Y = 1/2 for a constant current distribution and Y = 0 in the skin effect case.

B Contribution from a curved segment

The goal is to evaluate the integral L (0) from equation (11) for a segment of length 2b
and constant curvature radius R,

LR (0) = ˆ ds′ dA dA ′ |x (0, r, φ) −x′ | θ (b − |s′ |) j (r) j (r′ ) cos s′ R . The distance up to order O (R− 2 ) follows from (x − x′ )2 ∼ = s′ 2 + N 2 + q, N 2 = r2 + r′ 2 − 2rr′ cos (φ − φ′ ) , q = s′ 2 (r′ cos φ′ + r cos φ) R− 1 + rr′ R− 2 cos φ cos φ′ − s′ 4 / 12R2 + s′ 6 / 360R4 . Expanding in q gives 1 |x − x′ | = 1 (s′ 2 + N 2 )1 /2 − q 2 (s′ 2 + N 2 )3 /2 + 3q2 8 (s′ 2 + N 2 )5 /2 + ... The N in the denominators is negligible for b ≫ a and the elementary integrals lead to LR (0) − L0 (0) = − 11 24 b2 R2 + O a2 R2 ln b a . (B.1) 12

Figure 3: Two non-adjacent coplanar straight segments. The coordinates m, m′ and n, n′ measure the distance of the end points from the intersection of the segment extensions. The term of order O (a2 / R2 ) contributes to the error of the self inductance formula with the expected order of magnitude. It is essential however, that the O (b2 /R2 ) term, which is of order a/R because of b = √ aR, cancels against a contribution from the curve integral Lγ (0) from equation (10). To verify this cancellation start with d db Lγ (0) = − cos (j (0) , j (b)) |x (0) −x (b) | − − cos (j (0) , j (−b)) |x (0) −x (−b) | . Inserting |x − x′ | = 2R sin b/2R and cos (j (0) , j (b)) = cos b/r for a section with curvature radius R gives d db Lγ (0) = −2 b 1 − 11 24 b R 2 + ... , with integral Lγ (0) = const − 2 ln 2b a + 11 24 b2 R2 + ... (B.2)

Curve integral for adjacent straight segments

For completeness we display here the curve integral contribution to the self inductance from adjacent straight segments of length c and d with an angle α between the currents, Lγ (c, d, α) = µ0 2π cos (α) {c asinh d + c cos α c sin α + d asinh c + d cos α d sin α − (c + d) asinh cos α sin α − 2b 1 (1 − cos α) asinh 1 − cos α sin α }. For each corner such a term is to be added to the contribution (7) of the (straight) segments by themselves. The b-term is of order O (a) for b = a/2 and normally may be neglected.

D Curve integral for non-adjacent coplanar straight

segments

What then is missing for calculating the self inductance of a loop consisting of arbitrary
coplanar straight segments is the mutual contribution from non-adjacent straight segments,
13

see figure (3). L (m, m′ , n, n′ , α) = µ0 2π {A (m′ , n′ , n, α) + A (n′ , m′ , m, α) + A (m, n, n′ , α) + A (n, m, m′ , α)} , A (w, u, v, α) = w asinh u + w cos α w sin α − w asinh v + w cos α w sin α cos α. This leads to an unwieldy expression already for a hexagon, but the calculation of the self inductance of such loops is a matter of algebra and geometry (not a new result).

References

[1] J.D. Jackson, ”Classical Electrodynamics”, p. 261. Wiley 1975.
[2] F. E. Neumann, “Allgemeine Gesetze der inducirten elektrischen Str¨ome”. Abhand-
lungen der K¨oniglichen Akademie der Wissenschaften zu Berlin, aus dem Jahre 1845:
1-87, 1847.
[3] E.B. Rosa, The self and mutual inductances of linear conductors. Bulletin of the Bureau
of Standards 4 (2) (1907) 301-344.
[4] E.B. Rosa and F.W. Grover, Formulas and tables for the calculation of mutual and self-
inductance. Scientific papers of the bureau of standards (1916). Government printing
office, Washington 1948.
14