Metal Detector Coil Calculator and PI Pulse Timing

Size a pulse-induction coil, estimate resonance and damping, compare recovery with target decay, and explore calibrated air range versus coil diameter.

#air-core coil #eddy currents #pulse induction #RLC damping #Wheeler #wire packing

This metal detector coil calculator helps engineers, students and coil builders compare an air-core search winding with the timing needed by a pulse-induction receiver. It combines Wheeler inductance, copper resistance, lumped self-resonance, parallel damping and a target-decay model. Use it to plan a winding and an oscilloscope test before committing to a coil former or switching circuit.

Start with the mean winding diameter, turn count and the two dimensions of the winding bundle. The radial build is the difference between outer and inner radii; it is not the outer diameter. Enter bare copper diameter and the combined length of both leads. The initial resistance estimate applies at 20 °C. After winding, select measured inductance and resistance so the timing calculation uses the actual coil and cable.

For pulse induction metal detector design, enter supply voltage, pulse width, repetition rate and the additional resistance of the switch and supply. The calculator reports end-of-pulse current and stored magnetic energy. Total shunt capacitance includes the winding, cable, transistor and receiver. Together with inductance it sets the lumped self-resonant frequency and the critical parallel damping resistor. A larger resistor can allow ringing; a much smaller resistor can produce a slow overdamped tail.

The recovery result is a conservative planning allowance. It combines a constant-clamp ramp estimate with a linear-network settling envelope, then compares that allowance with your measured receiver recovery time. It does not simulate MOSFET avalanche, a voltage-dependent clamp or amplifier overload. Confirm the real waveform and component stresses on an oscilloscope. Magnetic energy per pulse is not a damping-resistor power rating.

The target section uses a solid, non-magnetic copper or aluminium sphere to show the fundamental eddy-current decay time. Press play to see the circulating-current illustration fade and follow its position on the decay curve. Marker speed is illustrative. The displayed sample and gate fractions describe the late exponential mode, not the total early-time response or a probability of detection. Coins, foil, rings and steel can behave very differently.

Depth requires an empirical anchor. The default reference follows a Clone PI-W builder's reported coin air test with a 210 mm coil; replace it with your own repeatable test. The diameter curve holds target, turns, peak current, sampling and receiver threshold fixed. It cannot predict how a larger coil changes noise, ground response or drive current. Air range is not guaranteed burial depth.

Before field use, check landowner consent and local detecting permits. Monument-protection rules can restrict searching as well as excavation. In Germany these rules vary by state; the linked North Rhine-Westphalia guidance is an example checked on 27 October 2026.

How it is calculated

Coil and transmit pulse

Wheeler's rectangular-section air-core approximation is

LW[μH]=0.8r2N26r+9b+10c,L_W[\mu H]=\frac{0.8r^2N^2}{6r+9b+10c},

with mean radius rr, axial height bb and radial build cc in inches. The coefficient is explicitly converted for millimetre inputs. See Wheeler (1928), DOI 10.1109/JRPROC.1928.221309; worked examples are attached in CREI 740 pp.33–34 and Kesler, attached p.26. A thin detector bundle is outside Wheeler's original near-equal-denominator-term 1% accuracy claim; measure the completed coil.

Copper length is ℓ=πDN+ℓleads\ell=\pi DN+\ell_{leads} and R20=ρℓ/AR_{20}=\rho\ell/A, with annealed-copper resistivity 1/58  Ω mm2/m1/58\;\Omega\,mm^2/m at 20 °C (IEC 60028 convention used by the platform material data). Measured L and R override these estimates for pulse timing. With the damping resistor permanently connected, the supply seen by the coil is Vth=VbRd/(Rd+Rs)V_{th}=V_bR_d/(R_d+R_s) and Rth=RsRd/(Rs+Rd)R_{th}=R_sR_d/(R_s+R_d). Thus Ip=Vth[1−e−ton(Rc+Rth)/L]/(Rc+Rth)I_p=V_{th}[1-e^{-t_{on}(R_c+R_{th})/L}]/(R_c+R_{th}) and W=LIp2/2W=LI_p^2/2.

Discrete wire packing

The copper-area inequality is only a necessary screening bound; area alone cannot establish winding feasibility. A separate discrete check constructs aligned and staggered layouts of equal bare round wires in the rectangular bundle, trying both axial and radial row stacking.

Aligned rows hold ⌊b/d⌋⌊c/d⌋\lfloor b/d\rfloor\lfloor c/d\rfloor wires. Half-diameter staggered rows have pitch 3d/2\sqrt{3}d/2; alternating row counts are ⌊w/d⌋\lfloor w/d\rfloor and ⌊w/d−1/2⌋\lfloor w/d-1/2\rfloor. The first row consumes one diameter of height, so the row count is 1+⌊(h/d−1)/(3/2)⌋1+\lfloor(h/d-1)/(\sqrt{3}/2)\rfloor. We also test equal-count rows with offset s=min⁡[d/2,w−d⌊w/d⌋]s=\min[d/2,w-d\lfloor w/d\rfloor] and pitch d2−s2\sqrt{d^2-s^2}. These are geometric constructions derived here from circle contact, not an empirical fill factor.

The largest tested layout supplies the discrete capacity. When either bundle dimension is exactly one wire diameter, all wire centres lie on a line and capacity reduces to the floor of the other dimension divided by diameter. Thus a 0.63 mm by 5 mm bundle holds at most seven bare 0.63 mm wires, and nine turns need at least 5.67 mm along the row.

This is a conservative regular-row construction check, not a search for every possible irregular circle packing. All counts use bare copper diameter; passing does not account for enamel, winding tolerances or manufacturability. Allow additional space for those effects.

Resonance, damping and recovery allowance

The turn-off network is a series Rc,LR_c,L branch in parallel with RdR_d and total shunt C. Kirchhoff's laws give

v′′+2αv′+ωn2v=0,2α=Rc/L+1/(RdC),ωn2=(1+Rc/Rd)/(LC).v''+2\alpha v'+\omega_n^2v=0,\quad 2\alpha=R_c/L+1/(R_dC),\quad \omega_n^2=(1+R_c/R_d)/(LC).

The lossless lumped resonance is 1/(2πLC)1/(2\pi\sqrt{LC}). Setting α=ωn\alpha=\omega_n gives

Rcrit=L/C2+Rc/L/C.R_{crit}=\frac{\sqrt{L/C}}{2+R_c/\sqrt{L/C}}.

These equations and the following planning allowance are derived here, rather than copied from a detector design.

The assumed clamp releases at Ie=min⁡(Ip,Vc/Rd)I_e=\min(I_p,V_c/R_d). Its constant-voltage ramp estimate is tf=(L/Rc)ln⁡[(Vc+RcIp)/(Vc+RcIe)]t_f=(L/R_c)\ln[(V_c+R_cI_p)/(V_c+R_cI_e)]. The tail starts from the conservative voltage V0=min⁡(Vc,IpRd)V_0=\min(V_c,I_pR_d). Let λ=α\lambda=\alpha for an underdamped network, otherwise λ=ωn2/(α+α2−ωn2)\lambda=\omega_n^2/(\alpha+\sqrt{\alpha^2-\omega_n^2}). A normalized zero-slope voltage tail is bounded by (1+αt)e−λt(1+\alpha t)e^{-\lambda t}. For ϵ=Vready/V0\epsilon=V_{ready}/V_0, use the deliberately generous allowance

tb=max⁡ ⁣[0,2λln⁡1+α/λϵ].t_b=\max\!\left[0,\frac{2}{\lambda}\ln\frac{1+\alpha/\lambda}{\epsilon}\right].

The inequality (1+cx)e−x≤(1+c)e−x/2(1+cx)e^{-x}\le(1+c)e^{-x/2} for x≥0x\ge0, c≥1c\ge1, supplies this bound. Add 1/λ1/\lambda for peak arrival. Suggested sample start is max⁡(tf+1/λ+tb,trx)\max(t_f+1/\lambda+t_b,t_{rx}). The plotted coil curve starts its decay after the ramp and peak-arrival allowance and caps the settling envelope at unity. It is an allowance curve, not a simulated clamp waveform.

Target decay and calibrated air range

For a solid non-magnetic sphere, the Smythe/Wait magnetic-diffusion idealization has fundamental time constant τ=μ0σa2/π2\tau=\mu_0\sigma a^2/\pi^2; Geonics TN32 Eq.2 states this result. Material resistivity is from the platform's 20 °C metal data. The late-mode fraction is e−t/τe^{-t/\tau}; the gate result is its exact rectangular-gate average. No absolute target amplitude is inferred.

For an on-axis small compact target, a circular coil's coupling is proportional to Nr2/(r2+z2)3/2Nr^2/(r^2+z^2)^{3/2}. Reciprocity makes the transmit-receive signal proportional to its square at fixed current. Equating it to an empirical reference gives

z2=(D/Dref)4/3[(Dref/2)2+zref2]−(D/2)2.z^2=(D/D_{ref})^{4/3}\left[(D_{ref}/2)^2+z_{ref}^2\right]-(D/2)^2.

A negative right side returns zero range with a message. The default anchor is DesAlex's 18 June 2009 Clone PI-W report: 21 cm coil and about 25 cm coin air range. That coin measurement is separate from the sphere timing surrogate.

Assumptions

  • Circular, air-core, rectangular winding bundle; all turns carry current in the same sense. Lead inductance and distributed winding details are not calculated.
  • Discrete packing checks use bare-wire aligned or staggered rows in both bundle orientations. The area check alone does not establish feasibility; enamel and winding tolerances need additional space. Irregular optimized packing is outside this construction model.
  • Copper resistance estimate at 20 °C. Enamel, joints, heating, skin and proximity effects are excluded; use measured resistance at operating temperature when these matter. The material constants follow the platform's IEC 60028 copper and IEC 60889 conductor-aluminium data; check against the standard in force for procurement.
  • Transmit current starts from zero. The RL rise ignores the short capacitive switching-edge transient; the damping resistor's loading through source resistance is included.
  • The reverse-clamp voltage is an effective design assumption across the coil. Turn-off is approximated by a clamp ramp and a separately budgeted linear RLC tail. Supply-referenced clamp topology, clamp capacitance, voltage-dependent MOSFET capacitance and amplifier recovery are not simulated.
  • The sphere calculation concerns only the fundamental late exponential mode of a homogeneous, non-magnetic solid conductor. It is not a model of a coin, steel target or hollow shell.
  • The depth comparison is calibrated to one empirical air test. The calibration test's target, orientation, turns, peak current, sampling and receiver threshold stay fixed. Changing electrical inputs does not automatically rescale this independent diameter comparison.
  • Drawings are schematic. Moving markers indicate current circulation; their speed is illustrative and their fading follows the sphere mode. Chart axes explicitly retain their stated mm and µs units; physical input and result cards support the viewer's metric/US conversion.

Where it applies

Use for preliminary sizing and timing comparisons of circular mono-coils, followed by measurement. Wheeler's original approximately 1% claim applies only to suitably proportioned rectangular windings, not every thin detector bundle. The self-resonance result requires a supplied lumped capacitance; it does not predict that capacitance from wire spacing. Discrete packing acceptance only establishes a tested bare-wire row layout; verify insulated wire diameter, winding clearance and the finished bundle.

The recovery budget is a conservative engineering allowance within its assumed linear network, not a guaranteed safe receiver delay or a semiconductor stress calculation. Verify the actual peak voltage, avalanche or clamp energy, resistor pulse rating and receiver overload recovery. The energy-rate result excludes transmit copper losses, damping-resistor on-time loss and the rest of the electronics. No thermal rating is selected.

Sphere late-mode retention is meaningful for comparison after higher modes have decayed. At early times, finite current turn-off, multiple diffusion modes and pulse history affect the total signal. A mode fraction is not a metal-discrimination result or a detection threshold.

Air-range extrapolation assumes a small target on the coil axis in a nearly uniform applied field across that target. Keep coil diameter within 0.5–2 times the measured reference and keep the winding thin. Ground mineralization, target geometry, receiver noise, electromagnetic interference and sweep speed can dominate actual depth. The Clone PI-W anchor is a builder report with unspecified test uncertainty; it is not independent validation of the extrapolation law.

Detecting and recovery of finds can require permits. The NRW portal describes permission under §15 DSchG NRW and separate landowner consent, checked 27 October 2026. Other German states and countries differ. Check the rules currently in force before searching or digging.

References

  • Wheeler, H. A. (1928), Simple Inductance Formulas for Radio Coils, Proceedings of the IRE 16(10), 1398–1400. https://doi.org/10.1109/JRPROC.1928.221309
  • CREI 740, Receivers, Radio Frequency Amplification, printed pp.28–29 and p.26; attached file pp.33–34 and 31. https://www.worldradiohistory.com/Archive-Courses/CREI/CREI-TV-740.pdf
  • Kesler, R., Multilayer air-cored coils, Electronics World, September 1997, pp.752–753 (attached p.26). https://worldradiohistory.com/UK/Wireless-World/90s/Electronic-World-1997-09.pdf
  • Wait, J. R. (1951), A conducting sphere in a time varying magnetic field, Geophysics 16(4), 666–672. https://doi.org/10.1190/1.1437716
  • Geonics Technical Note TN32 (2002), sphere time constant, Eq.2, and practical time-domain EMI limitations. https://geonics.com/pdfs/technicalnotes/tn32.pdf
  • Lammert Bies, Pulse induction metal detector with DSP: search-coil current and flyback worked examples. https://www.lammertbies.nl/electronics/pi-metal-detector
  • DesAlex, Clone PI-W thread, coil guidance 5 May 2009 and 21 cm coil / 25 cm coin air-test report 18 June 2009. https://forum.cxem.net/topic/47662-clone-pi-w-своими-руками/
  • North Rhine-Westphalia official Bauportal, permission for metal detecting under §15 DSchG NRW; checked 27 October 2026. https://bauportal.nrw/antraege/antrag-auf-erlaubnis-zur-sondengaengerei

Inputs

InputSymbolUnitDefault
Mean winding diameter
Diameter through the centre of the winding bundle; inner diameter = mean diameter − radial build.
D_coilmm250
Number of turnsN_turn25
Axial winding height
Thickness perpendicular to the coil plane. Wheeler rectangular-section bundle model.
b_axmm5
Radial winding build
Outer radius minus inner radius, measured across the bundle.
c_radmm5
Bare copper wire diameter
Excludes enamel. Resistance is calculated at 20 °C; packing checks use bare diameter, so allow extra space for insulation and winding tolerances.
d_wiremm0.63
Total copper lead length
Sum of both leads, assumed to use the same conductor diameter as the winding.
l_leadm1
Use measured inductance and resistance
Recommended after building. Both resistance values include the leads; the switch and supply resistance below is additional.
measuredfalse
Measured coil inductance
Measure at a frequency well below self-resonance.
L_measuH400
Measured coil and lead resistance
DC resistance at the operating temperature, excluding the switch and battery.
R_measohm1.5
Total shunt capacitance
Coil, cable, switch and input capacitance together. Estimate from measured self-resonance when possible.
C_parpF250
Parallel damping resistor
Across the series coil-and-lead resistance and inductance, in parallel with the total capacitance.
R_dampohm680
Switch and supply resistance
Additional series resistance during the transmit pulse.
R_driveohm0.15
Applied transmit voltage
Voltage available to the series RL circuit during the pulse.
V_batV12
Transmit pulse width
Pulse current starts from zero; incomplete reset is flagged.
t_onus75
Pulse repetition rate
Used for the timing budget and magnetic energy rate.
f_repHz500
Effective reverse clamp voltage
Reverse voltage across the coil during fast turn-off, excluding the supply. This is a timing assumption, not a MOSFET rating.
V_clampV400
Coil tail voltage for settling budget
Conservative envelope threshold after the clamp stage. Receiver overload recovery can take longer.
V_readymV1
Measured receiver recovery time
Time from turn-off until the receiver is usable; confirm on an oscilloscope.
t_rxus20
Sample start after turn-off
Compare with the conservative recovery budget and the target time constant.
t_delayus25
Sample gate width
The target result averages its late exponential mode over this gate.
t_gateus10
Non-magnetic sphere material at 20 °C
Solid non-magnetic sphere only. Coins, foil, rings and steel require a measured target decay time or a different model.
target_metalCopper, annealed (IACS)
Sphere diameter for timing model
This sphere is a timing surrogate. The separate depth calibration can represent a coin or another small compact target.
d_targetmm5
Air-test reference mean coil diameter
Default follows DesAlex's 21 cm Clone PI-W coil report; replace with your own measurement.
D_refmm210
Measured reference air range
Default: a reported coin air test, not this sphere. Keep the same target, turns, peak current, sampling and receiver threshold for the diameter comparison.
z_refmm250
Time after turn-off ▶
Current markers show circulation schematically; their speed is illustrative. Fading follows the sphere's late mode.
t_animus25

Results

With the default inputs above, the calc gives:

ResultValue
Coil inductance used364 uH
End-of-pulse current2.17 A
Lumped self-resonance527.6 kHz
Critical parallel damping resistor603 Ω
Damping at the selected resistorUnderdamped
Suggested earliest sample20 us
Sphere fundamental time constant46.2 us
Gate-average late-mode fraction (0–1)0.523
Magnetic energy per pulse0.857 mJ
Reference air range, fixed turns and current278 mm

Checked against sources

Each case below was calculated by this calc and agrees with the value from its source (within the tolerance where one is shown).

CREI worked rectangular coil: 204 turns, 444 µH

Source: CREI 740, attached p.34 (printed p.29): r = 0.372 in, l = 0.35 in, h = 0.493 in, 204 turns; printed result 444 µH. Substitution of the page's own rounded dimensions gives 446.78 µH, a 0.63% arithmetic/rounding discrepancy. The correct Wheeler expression is retained; 3 µH tolerance compares with the printed number while explicitly recording this discrepancy. This is a published calculation, not a coil measurement.

Inputs: Mean winding diameter = 18.898 mm; Number of turns = 204; Axial winding height = 8.89 mm; Radial winding build = 12.522 mm; Bare copper wire diameter = 0.635 mm; Total copper lead length = 0 m; the others as above.

ResultSourceThis calc
Wheeler rectangular-bundle inductance444 ± 3 uH446.8 uH

Kesler worked multilayer coil: 450 turns, 2 mH

Source: Robert Kesler, Electronics World September 1997, printed pp.752–753, attached page 26: initial former inner diameter 12 mm, width 16 mm, radial build 6 mm, mean diameter 18 mm, N = 450 for 2 mH. Check uses this initial geometry before the article changes wire size and radial build. The printed coefficient 7.87 nH/mm is rounded; N is rounded to a whole turn. It does not use the article's inconsistent first design-example turn count.

Inputs: Mean winding diameter = 18 mm; Number of turns = 450; Axial winding height = 16 mm; Radial winding build = 6 mm; Bare copper wire diameter = 0.35 mm; Total copper lead length = 0 m; the others as above.

ResultSourceThis calc
Wheeler rectangular-bundle inductance2,000 ± 12 uH2,002 uH

CREI tuned circuit: 206 µH, 40 pF, 1.75 MHz

Source: CREI 740, attached page 31 (printed p.26), Eq.24 gives approximately 206 µH with 40 pF at 1,750 kc/s. This independently checks the lumped LC resonance conversion; it is not a prediction of distributed self-capacitance.

Inputs: Use measured inductance and resistance = yes; Measured coil inductance = 206 uH; Total shunt capacitance = 40 pF; the others as above.

ResultSourceThis calc
Lumped self-resonance1,750 ± 5 kHz1,753 kHz

Bies PI example: 9 V, 400 µH, 2 Ω, 250 µs → 3.2 A

Source: Lammert Bies, Pulse induction metal detector with DSP, Search coil design: 400 µH assumed, 2 Ω circuit resistance, 9 V minimum supply and 250 µs pulse give about 3.2 A. https://www.lammertbies.nl/electronics/pi-metal-detector . Source is an author design calculation, not an experimental measurement.

Inputs: Use measured inductance and resistance = yes; Measured coil and lead resistance = 2 Ω; Switch and supply resistance = 0 Ω; Applied transmit voltage = 9 V; Transmit pulse width = 250 us; the others as above.

ResultSourceThis calc
End-of-pulse current3.2 ± 0.02 A3.211 A

Related Resources

Responses

Your email address will not be published. Required fields are marked *

L ading...