



An introduction to the art of Solenoid inductance calculation By David W Knight
Solenoid
inductance calculation (v 0.20, unfinished, 4th
Feb. 2016) 
On
the Geometrical Mean Distance of two figures in a plane.
J C Maxwell, Trans Roy. Soc. Edinburgh, Vol 26, 1872. [Maxwell
1872 GMD] Calculation of the selfinductance of 1layer coils. E B Rosa, BBS Vol 2, No.2, p161187 [BS Sci. 31]. This is the original paper giving Rosa's method for solenoid inductance, including derivation of the external GMD for currentsheet segments. The Self and Mutual Inductances of Linear Conductors. E B Rosa. BBS Vol 4, No 2, 1908 [BS Sci. 80]. The Inductance Coefficients of Solenoids. H Nagaoka. J. Coll. Sci. Tokyo, Vol 27, 6, 1909 (1.7 MiB). Note that all of Nagaoka's papers can be obtained from the University of Tokyo repository . Formulas and Tables for the Calculation of Mutual and Self Induction. E B Rosa & F W Grover, revised 1916 3rd edition, 1948. [BS Sci. 169]. (15 MiB)  corrected version of BBS Vol. 8, No. 1 (1911 & 1912) (same page numbering). Additions to the formulas for the calculation of mutual and self inductance. F W Grover. BBS, Vol 14, No 4, July 12, 1919 [BS Sci. 320] (1 MiB). Formula for the inductance of a helix made with wire of any section. C. Snow. [BS Sci. 537], 1926 (5.74 MiB). Methods for the Derivation and Expansion of Formulas for the Mutual Inductance of Coaxial Circles and the Inductance of SingleLayer Solenoids. F W Grover, NBS J. Research, Vol 1, 1928. [BS RP16]. p487511 Comparison of Formulas for the Calculation of Inductance of Coils and Spirals Wound with Wire of Large CrossSection. F W Grover, JBS Vol 3. 1929. [BS RP90]. A simplified Precision Formula for the Inductance of a Helix with Corrections for the Leadin Wires. C Snow, JBS Vol 9. 1932, p419426. (0.5 MiB). Formulas for Computing Capacitance and Inductance. C Snow. NBS Circular 544, 1954. [BS Circ. 544]. (1.1 MiB). 
The following easytouse routine for calculating solenoid inductance by the RosaNagaoka method has been added to the macro library in the spreadsheet Lcalcs.ods . The three functions called are also in the macro library and are discussed in the solenoid incuctance calculation article. Note that the calculation is for external partial inductance only. Internal inductance and lead inductance should be added to the result 
Function
LxRosa(byval h as double, D as double, N as double, dw as double, p as
double) 'External inductance in henrys of a roundwire solenoid, using the NagaokaRosa method. 'D W Knight. v 1.00, 20160111 'Calls functions W82W(), KMGO() and Rosaks90() 'h = coil length / m 'D = coil diameter / m 'N = number of turns 'dw = wire diameter / m 'p = winding pitch / m Dim kNagaoka as double, Lsheet as double, Rosacorr as double kNagaoka = W82W(D/h) Lsheet = 1e7*pi()*pi()*D*D*N*N*kNagaoka/h Rosacorr = KMGO(N) + Rosaks90( p/dw , p/D , 0) LxRosa =Lsheet 2E7*pi()*D*N*Rosacorr end function 
Material not yet
transferred. The text below is in the process of being updated and extended. The information will at some point be transferred to the main solenoid inductance article above. 
B1. Solenoid inductance
calculation vs. measurement: Using the techniques outlined in the previous sections, it is possible to calculate the lowfrequency external inductance of solenoids to an accuracy of a few parts in 1000 using a micrometer or a good set of engineer's callipers and a personal computer. If a scientific calculator is used, the computation can be simplified by using one of the Wheeler unrestricted formulae with only a small loss of accuracy. In either case, the outcome is usually better than can be obtained using generalpurpose measuring equipment. All solenoid inductance calculations are dependent on the effective currentsheet diameter (D) or radius (r=D/2). In section 3 it was noted that at low frequencies D=D0, where D0≈Da and Da is the coil diameter measured from wirecentre to wirecentre, whereas at very high frequencies D=D∞, where D∞ is less than D0 but always greater than the diameter at the inside of the conducting cylinder (i.e., D0d). Consequently, the measured radiofrequency inductance of a coil (substantially below the SRF at least) will always lie between the low frequency inductance L0, which is the value calculated using D=D0; and the minimum inductance L∞, which is the value which obtains when D=D∞ (presuming that we can determine D∞). The inductance of a solenoid without a magnetic core, presuming that the wire diameter is << the solenoid diameter, is given within a few parts per thousand by:
where D=2r, k_{L} is Nagaoka's coefficient , k_{s(e)} Rosa's roundwire self inductance correction, for external inductance only, and k_{m} is Rosa's roundwire mutual inductance correction.. The internal inductance L_{i} , generally < 1% of the total inductance, can be calculated to within 0.016% using the LiPACAML approximation, see Zint.pdf. (also function Lintern() described on the Zint index page) From the previous discussion, the inductance of a solenoid, neglecting the proximity effect, is:
where k_{L(0)} is Nagaoka's coefficient calculated using /D_{0}. The worstcase minimum inductance, due to reduction in the effective diameter is:
where k_{L(∞)} is Nagaoka's coefficient calculated using /D_{∞}. 
Example B1a R G Medhurst, on page 42 of his 1947 paper "H.F Resistance and SelfCapacitance of SingleLayer Solenoids", gives accurate dimensions of a testcoil (#31) used in a study of Q and selfcapacitance. The coil is wound on a grooved former using 40 turns of 20SWG bare wire (d=0.9144 mm), with a diameter of 51.9 mm measured at the outside of the winding. The length of the coil is given as 70.1 mm, but the winding pitch is stated to be: p=(70.1d)/(N1) from which we may deduce that the measurement was made to the outside of the wire on the side away from the terminations, (as shown in the diagram on the right). Hence the stated length is not the same as the equivalent current sheet length, which is defined as: 

= Np Thus it appears that Medhurst uses a definition of coil length that differs from that intended by Rosa and Grover (as discussed in section 2) but is an understandable misinterpretation of the explanation given [Sci169, p119]. Ambiguities of this type can introduce errors of several % in inductance calculations, and there is no use in claiming accuracy in parts per 1000 unless care is taken to avoid them. In this case the equivalent currentsheet length is: = 40p = 40(70.10.9144)/39 = 70.9596 mm. Using this value, with Nagaoka's coefficient kL calculated using Lundin's formula , we have:
The roundwire correction ks(e) is: ks(e) = (3/2)  ln(2p/d)
and k_{m} =0.3142. Medhurst measured the inductance of this coil by averaging readings taken over the range 780 to 860 kHz. Taking the average frequency as 820 kHz, the internal inductance is 101.8 nH. The results of the calculation are shown below, where the 'greater than' symbol (>) is used for the L_{∞} value because the coil inside diameter (D_{0 } d) was used instead of the slightly larger but undefined D_{∞}. L_{0} = (μ_{0} π D_{0}² N² k_{L(0)} /4) [μ_{0} D_{0} N (k_{s(e)} + k_{m} )/2] + L_{i} = 43.7331  0.5873 + 0.1018 μH
L_{∞} = (μ_{0} π D_{∞}² N² k_{L(∞)} /4)[μ_{0} D_{∞} N(k_{s(e)} +k_{m} )/2] + L_{i} > 42.3694  0.5768 + 0.1018 μH
Hence we expect the measured inductance to lie in the range: L = 43.248, +0, 1.354 μH 
Medhurst gives nine inductance measurements for this coil, taken at 10 kHz intervals over the range 780 to 860 kHz and corrected for self capacitance. These measurements are scattered about the mean value, whereas if the inductance were changing significantly with frequency in this region we would expect them to diminish progressively as the frequency increased. Hence we can infer that deviations from the mean are principally due to experimental error, and we can use this information to estimate the standard deviation of the mean value. The measurements are tabulated below: 

/ kHz 
L_{k} / μH 
L_{k}  L_{mean} 
(L_{k}  L_{mean})² 













































Lmean = 42.3744 


The mean (average) of the
measurements is given by:
Where 'Σ' (capital Sigma) means "the sum of" (over the range specified below and above), Lk is the 'kth' measurement of L, and in this case the number of measurements n=9. The square of the standard deviation (known as the variance) can be estimated from:
where we divide the sum of the squares of the deviations by n1 rather than n because the use of the sample mean rather than the population mean (i.e.the unknown 'true' value of L) removes a degreeoffreedom from the data [see Data Analysis]. Thus, for Medhurst's measurements: σ = √(0.02482222 / 8) = 0.0557 Hence:
Comparing this with our calculations (43.248 > L > 41.894 [μH] ), we can observe that at 820 ±40 kHz, the inductance of the coil is lower than that predicted from physical dimensions, but within the range allowed by variation of effective diameter. It is tempting to think that the highfrequency drop in inductance has been quantified by this comparison (and possibly it has), but such is not necessarily the case. The standard deviation calculated from the scatter in the data is a measure of precision not accuracy, and unknown systematic errors, or simple uncertainties in the coil dimensions, prevent us from drawing firm conclusions. 
Example B1b Once the principle has been established, accurate inductance calculation does not have to be a laborious process.We might note from the earlier discussion however; that the available literature on the subject contains mistakes, with the consequence that existing computer programs and models cannot be trusted without verification of the numbers that they produce. It is a relatively straightforward matter however, to set up a spreadsheet to do the calculation, and once the entered formulae have been verified, new rows can be added and old ones deleted at will. A suitable calculation template is provided by the OpenDocument Spreadsheet file: L_calcs.ods (sheet 2), which can be downloaded and amended as required. The spreadsheet as provided is filled with example inductance calculations for seven coils described in the academic literature. Medhurst's coil #31 (M31) from the previous example is there , along with another coil (M32) described in the same paper (these are the only test coils for which Medhurst gave full details). There are also: a coil described in reference [10] (GKMR1), and four coils described in reference [11] (MKG14). Interest in the last five coils lies in the fact that the inductance measurements are accompanied by calculations using a method different to the one described here. The results are summarised in the table below: 








































































































L_{0} /μH 


[10] 
[11] 
[11] 
[11] 
[11] 
The main point to note is that for
coils MKG14, lowfrequency (L0)
inductances calculated using the modified National Bureau of Standards
(NBS) method (i.e., the method used here) are in exact agreement with
the method used in ref. [11],
and that for a measurement frequency of 1 kHz, the internal inductance
is the same as the DC value. This provides a vindication, not of the
NBS method which is not in doubt, but of the formula verification
procedures applied here. The calculation of ref. [10] however, is not
in agreement with the NBS method, despite the statement that it was
carried out using the same method as in ref. [11]. Although very
close to the measured value, the calculation of ref. [10] does not
correspond to L_{0} or L_{0}L_{i},
and appears for some reason to be in error. What we find here is that
the measured inductance is sandwiched between our L_{0}
and L_{∞}>
values as in the case of coils M31 and M32, and that the very low p/d
ratio of 1.02 (and probable use of copper tubing rather than wire) has
given rise to an early onset of the highfrequency regime. The results
for coils MKG14 show nothing more onerous than a measurement standard
deviation of around 1%, and that the HF regime is not evident for these
coils at a frequency of 1 kHz. If there is anything to be inferred from these comparisons between measurement and calculation, it is that the skin effect and the proximity (effective diameter) effect are linked. All of the coils measured at 1 kHz have an internal inductance component that is the same as the DC value, and show no sign of the diameter effect even though some have low p/d ratios. All of the coils measured at frequencies above 1 kHz have an internal inductance component substantially less than the DC value, and the inductances lie between the L_{0} and L_{∞} values. Logically, this is to be expected, because the proximity effect modifies the current distribution within the wire and therefore changes the internal inductance. The proximity effect however, is not simply a perturbation of the internal inductance, because the total variation it causes is greater than the internal inductance. >>>>>> Calculation methods have been refined since 2008. More accurate formulae for effective current sheet diameter have been derived. New spreadsheet: Lcalcs.ods (sheet 1, still a work in progress). >>>>>> 
B2. Additional sources
of deviation: Points to note when comparing inductance measurements with calculations: The presence of conducting material in the vicinity of the coil disposed in such a way as to form a shortedturn coupled to the coil, and especially of any metal screeningcan or box, will cause the effective inductance to be less than that predicted for an isolated coil. The presence of any nonconductive ferromagnetic material in the vicinity of the coil will cause the inductance to be greater than that predicted by formulae which assume that μ=μ_{0}. Solenoid inductance is paradoxical: the reason being that inductance is only completely defined when the terminals of the inductor are coincident in space, whereas the terminals of a solenoid are by definition separated by a distance . All measured solenoid inductances therefore contain a contribution from the connecting wires. This added inductance, usually a few 10s of nH, can be estimated using the methods outlined in Comps & materials 2, (section 6) although the rough guide "add about 20 nH per inch" (8 nH/cm) will probably suffice in many instances. This implies about 12.5 cm of connecting wire before a change of 1 will occur in the first decimal place of an inductance expressed in microHenrys (and is about right for 1 mm diameter wire). The lumped inductance model breaks down when a coil is operated close to, at, or above its selfresonance frequency. 
[10]
"Stray Capacitances of SingleLayer Solenoid AirCore
Inductors", G. Grandi, M K Kazimierczuk, A Massarini, U
Reggiani. IEEE Transactions on Industry Applications, Vol 35, No. 5,
Sept/Oct 1999, p11621168. Available from: Classic Tesla at time of writing. www.classictesla.com/download/ia99.pdf . [11] "Lumped Parameter Models for Single and MultipleLayer Inductors", A Massarini, M K Kazimierczuk, G Grandi. Proc PESC '96, June 1996, p295301. Available from: Classic Tesla at time of writing. www.classictesla.com/download/pesc96.pdf . 



