Measurement, which is a circuit on a circuit

The resistance that depends on the reading

Three of a shunt's errors are free of the current being measured, which is the whole content of the burden-voltage optimum. The fourth is not: the shunt dissipates, warms, and its resistance rises — so the divisor the reading uses is a function of the reading. Solved as a fixed point it agrees with R₀/(1 − αθI²R₀) to 2×10⁻¹⁶, and along the optimum, where the dissipation is I·u* rather than I²R, the error is the FIRST power of the current: 7.75 ppm at an ampere, 775 at a hundred, fitted exponent 1.0007.

Assumes: The ammeter that is a resistor · The loss that depends on what it causes

The burden-voltage optimum is a statement about three errors and its most quoted property is what it leaves out. A shunt’s insertion error is the burden over the supply; the amplifier’s offset error is the offset over the burden; the shunt’s own Johnson noise is four orders below both. None of the three contains the current being measured. The ammeter that is a resistor makes a good deal of that, and it is right to: ten milliamps and a hundred amps want the same 7.75 millivolts of burden, and only the resistance moves.

There is a fourth error and it has the current in it, twice.

A shunt dissipates I2RI^2R. Its temperature rises by that times a thermal resistance to whatever it is mounted on. Its resistance rises with its temperature, by whatever its coefficient is. And the current is computed by dividing the measured voltage by the resistance — so the divisor is a function of the dividend, and the reading is the solution of an equation rather than the result of a division.

The shunt's resistance as a function of what it is measuringcomputed by solving, not by drawing, as a fixed point: the shunt dissipates I²R, its temperature rises by 20 K per watt, and at 50 ppm/K its resistance rises with its temperature — so the resistance the reading is divided by depends on the reading. Iterated to convergence it agrees with the closed form R₀/(1 − αθI²R₀) to 2.2e-16. Along the burden-voltage optimum, where R = u*⁄I, the dissipation is I·u* rather than I²R, so the temperature rise is 155 mK per ampere and the error is the FIRST power of the current — fitted exponent 1.0007 over five decades. That is the only one of the shunt's errors with the current in it, and it puts a term in I² into the reading, which is a curvature no single-current calibration removes. The upper curve is a shunt of fixed resistance, where the error is quadratic. The fixed point stops existing at 129 kA and never at a current a shunt will see.10n100n10µ100µ1m10m100m10m100m1101001kcurrent being measured (amperes)error from self-heating, as a fraction of the readingslope 1.001: the burden optimumslope 2: a shunt of fixed resistancea hundred parts per millioncoefficient50 ppm/Kthermal resistance20 K/Wburden held at7.746 mVrise per ampere154.9 mK/Aat 1 A7.746 ppmat 100 A775.2 ppmexponent1.0007fixed point fails at129.1 kAsolved, then checked — the reading is in the divisorlinear along the optimum, quadratic off it
Fig. 1 A shunt held at the best burden voltage, solved as a fixed point: the resistance sets the dissipation, the dissipation sets the temperature, the temperature sets the resistance. The lower curve is the resulting error along the optimum; the upper one is a shunt of fixed resistance, where the dissipation is quadratic. The slider is the temperature coefficient.

A fixed point, and the closed form that agrees with it

The three statements are

P=I2R,ΔT=θP,R=R0(1+αΔT)P = I^2R, \qquad \Delta T = \theta P, \qquad R = R_0\left(1 + \alpha \Delta T\right)

and they are circular. Substituting the first two into the third gives R=R0(1+αθI2R)R = R_0(1 + \alpha\theta I^2 R), an equation in RR whose solution is

R=R01αθI2R0R = \frac{R_0}{1 - \alpha\theta I^2 R_0}

The figure does not use that. It iterates the three statements to convergence — resistance to dissipation to temperature to resistance, four hundred times or until the change is below a part in 101510^{15} — and then compares the result against the closed form. They agree to 2×10162\times10^{-16} across five decades of current, which is the arithmetic’s floor.

That pairing is worth the cost. A fixed point solved by iteration and a fixed point solved in closed form are two computations of one number that share no steps, and the closed form is the kind of expression whose sign is easy to get wrong — the denominator is 1αθI2R01 - \alpha\theta I^2 R_0 and a first version of almost anybody’s derivation writes 1+1 + there, which gives an error of the right size and the wrong direction.

The iteration also makes the failure visible, which the closed form does not. When αθI2R0\alpha\theta I^2R_0 reaches one the denominator is zero and the closed form returns infinity; the iteration simply fails to converge, walking upward for ever. Both are describing thermal runaway, and only one of them looks like it.

The shunt's resistance as a function of what it is measuring. computed by solving, not by drawing, as a fixed point: the shunt dissipates I²R, its temperature rises by 20 K per watt, and at 400 ppm/K its resistance rises with its temperature — so the resistance the reading is divided by depends on the reading. Iterated to convergence it agrees with the closed form R₀/(1 − αθI²R₀) to 2.2e-16. Along the burden-voltage optimum, where R = u⁄I, the dissipation is I·u rather than I²R, so the temperature rise is 155 mK per ampere and the error is the FIRST power of the current — fitted exponent 1.0056 over five decades. That is the only one of the shunt's errors with the current in it, and it puts a term in I² into the reading, which is a curvature no single-current calibration removes. The upper curve is a shunt of fixed resistance, where the error is quadratic. The fixed point stops existing at 16.1 kA and never at a current a shunt will see.
Fig. 2 Four hundred parts per million a kelvin, which is a cheap thick-film chip resistor rather than a purpose-made shunt alloy. The error at a hundred amps is 6200 ppm — 0.62 per cent — which is five times the entire error the burden optimum was designed to minimise.

Along the optimum the heat is linear in the current

Here is the part that is specific to a shunt and is not obvious from the expression.

The dissipation is I2RI^2R, which is quadratic. But a shunt at the burden optimum has R=u/IR = u^*/I, so

P=I2uI=IuP = I^2\cdot\frac{u^*}{I} = I\,u^*

The dissipation is linear in the current and the resistance has cancelled out of it. Every shunt sized at the optimum dissipates 7.75 millivolts times its current, whatever that current is: 77.5 milliwatts at ten amps, 775 at a hundred, 7.75 watts at a thousand.

The temperature rise follows: ΔT=θIu\Delta T = \theta I u^*, which is 154.9 millikelvin per ampere at a thermal resistance of twenty kelvin per watt, and contains no resistance either. The fractional error is αθIu\alpha\theta I u^* and is therefore the first power of the current rather than the second.

current dissipation rise error at 50 ppm/K
0.1 A 0.775 mW 15.5 mK 0.775 ppm
1 A 7.75 mW 155 mK 7.75 ppm
10 A 77.5 mW 1.55 K 77.5 ppm
100 A 775 mW 15.5 K 775 ppm
1000 A 7.75 W 155 K 7750 ppm

Fitted over five decades the exponent comes out at 1.0007 rather than exactly one, and the departure is the fixed point feeding back on itself: a hotter shunt is a larger shunt is more dissipation. It grows with the coefficient — 1.0000 at ten parts per million a kelvin and 1.0146 at a thousand — which is a measurement of the feedback’s strength and is checked as such.

A shunt of fixed resistance, by contrast, has an exponent of exactly two, and the figure draws both. The difference between the two curves is the whole of what the burden optimum does to the thermal problem: it converts a quadratic error into a linear one, which is a much better position to be in at high current and a slightly worse one at low.

The shunt's resistance as a function of what it is measuring. computed by solving, not by drawing, as a fixed point: the shunt dissipates I²R, its temperature rises by 20 K per watt, and at 10 ppm/K its resistance rises with its temperature — so the resistance the reading is divided by depends on the reading. Iterated to convergence it agrees with the closed form R₀/(1 − αθI²R₀) to 2.2e-16. Along the burden-voltage optimum, where R = u⁄I, the dissipation is I·u rather than I²R, so the temperature rise is 155 mK per ampere and the error is the FIRST power of the current — fitted exponent 1.0001 over five decades. That is the only one of the shunt's errors with the current in it, and it puts a term in I² into the reading, which is a curvature no single-current calibration removes. The upper curve is a shunt of fixed resistance, where the error is quadratic. The fixed point stops existing at 645 kA and never at a current a shunt will see.
Fig. 3 Ten parts per million a kelvin, which is a good manganin or Zeranin shunt. The error is 1.55 ppm at an ampere and 155 at a hundred, the exponent is 1.0000 to four decimals, and the runaway current is 645 kiloamps. This is the part a precision measurement uses and the reason it exists.

The curvature no calibration removes

The reading is V=IR0(1+αθIu)V = I\,R_0(1 + \alpha\theta I u^*), which expands to

V=IR0+αθuR0I2V = I R_0 + \alpha\theta u^* R_0\,I^2

A term in I2I^2 in a measurement that is supposed to be linear. That is a different kind of error from a gain error or an offset, and the difference matters: a gain error is removed by a single-point calibration, an offset by a zero, and a quadratic term by neither.

Calibrating the gain at full scale makes the reading exact there and leaves it low everywhere below, by c(IfsI)c(I_{\text{fs}} - I) as a fraction of the reading. Expressed the way an instrument’s integral non-linearity is specified — as a fraction of full scale — that is cI(IfsI)/Ifsc\,I(I_{\text{fs}}-I)/I_{\text{fs}}, which is zero at both ends and peaks at mid-scale at exactly a quarter of the full-scale term. At fifty parts per million a kelvin and a hundred amps full scale the full-scale term is 775 ppm, so the worst departure is 194 ppm of full scale — which for a four-and-a-half-digit measurement is two counts and for a six-digit one is two hundred.

The quarter is worth keeping because it is exact and because it is the signature of a quadratic. A term in I2I^2 calibrated out at one point always leaves a parabola whose peak is a quarter of what was removed, and an instrument whose non-linearity is a quarter of its full-scale gain error has a quadratic in it whatever else is said about the mechanism.

So the honest specification of a shunt-based instrument has an integral non-linearity in it, that non-linearity is thermal in origin, and it depends on the mounting rather than on the part. The same shunt on a large copper plane has a third of the thermal resistance and a third of the non-linearity, and nothing about the electrical design has changed.

This is exactly the shape of the sensor inside its own answer, where a sensor’s own excitation changes the quantity it is measuring, and of the bridge that is linear near one point, where the linearity is a property of an operating region rather than of a circuit. A measurement that dissipates has a non-linearity, and the only questions are how large and what shape.

The shunt's resistance as a function of what it is measuring. computed by solving, not by drawing, as a fixed point: the shunt dissipates I²R, its temperature rises by 20 K per watt, and at 100 ppm/K its resistance rises with its temperature — so the resistance the reading is divided by depends on the reading. Iterated to convergence it agrees with the closed form R₀/(1 − αθI²R₀) to 2.2e-16. Along the burden-voltage optimum, where R = u⁄I, the dissipation is I·u rather than I²R, so the temperature rise is 155 mK per ampere and the error is the FIRST power of the current — fitted exponent 1.0014 over five decades. That is the only one of the shunt's errors with the current in it, and it puts a term in I² into the reading, which is a curvature no single-current calibration removes. The upper curve is a shunt of fixed resistance, where the error is quadratic. The fixed point stops existing at 64.5 kA and never at a current a shunt will see.
Fig. 4 A hundred parts per million a kelvin, which is an ordinary metal-film or a poorly chosen alloy. The error is 15.5 ppm at an ampere and 1550 at a hundred, and the runaway current has fallen to 64.5 kiloamps — still unreachable, which is the finding of the next section.

The runaway that exists and cannot be reached

A fixed point of this kind has a condition for existing, and it has turned up here before: αθI2R0<1\alpha\theta I^2R_0 < 1. Past that there is no solution and the resistance runs away.

Along the burden optimum, substituting R0=u/IR_0 = u^*/I gives the condition αθuI<1\alpha\theta u^* I < 1, so the critical current is

Icrit=1αθu=150×106207.746×103=129 kAI_{\text{crit}} = \frac{1}{\alpha\theta u^*} = \frac{1}{50\times10^{-6}\cdot 20\cdot 7.746\times10^{-3}} = 129\ \text{kA}

A hundred and twenty-nine kiloamperes. At a thousand parts per million a kelvin it is 6.45 kA and at ten it is 645 kA. None of those is reachable by a shunt, which would have vaporised at a thousandth of any of them — 129 kA at a 7.75 mV burden is a kilowatt of dissipation in a part sized for milliwatts.

That is a genuine and slightly unusual result here, where the usual finding is that a boundary is closer than expected. Here the boundary is real, it is the same mathematics that produces thermal runaway in the ratio that does not walk to one and in a bipolar transistor’s second breakdown, and for a shunt it is six orders of magnitude away from anything that can happen.

The reason is the linearity of the previous section. A device whose dissipation is quadratic in its drive has a runaway condition that scales as 1/αθR01/\sqrt{\alpha\theta R_0} and can be met; the burden optimum makes the dissipation linear, so the condition scales as 1/αθu1/\alpha\theta u^* and is fixed by a burden voltage of a few millivolts. Holding the burden constant has not merely improved the error — it has moved the instability out of reach entirely.

The shunt's resistance as a function of what it is measuring. computed by solving, not by drawing, as a fixed point: the shunt dissipates I²R, its temperature rises by 20 K per watt, and at 1000 ppm/K its resistance rises with its temperature — so the resistance the reading is divided by depends on the reading. Iterated to convergence it agrees with the closed form R₀/(1 − αθI²R₀) to 2.2e-16. Along the burden-voltage optimum, where R = u⁄I, the dissipation is I·u rather than I²R, so the temperature rise is 155 mK per ampere and the error is the FIRST power of the current — fitted exponent 1.0146 over five decades. That is the only one of the shunt's errors with the current in it, and it puts a term in I² into the reading, which is a curvature no single-current calibration removes. The upper curve is a shunt of fixed resistance, where the error is quadratic. The fixed point stops existing at 6.45 kA and never at a current a shunt will see.
Fig. 5 A thousand parts per million a kelvin — a copper shunt, which is what an improvised measurement uses because a piece of wire is at hand. The error is 155 ppm at an ampere and 1.57 per cent at a hundred, the fitted exponent has risen to 1.0146 as the feedback starts to bite, and the runaway is at 6.45 kiloamperes.

Which of the three thermal resistances matters

The single θ\theta above hides three paths and only one of them is a design variable.

The element’s own path to its terminations is set by the part and is typically a few kelvin per watt. The terminations’ path into the board is set by the pad geometry and the copper attached to it, and is where most of the twenty comes from. The board’s path to the ambient is the largest of the three and the slowest, and it is shared with everything else on the board.

Only the middle one is the designer’s, and it is worth a factor of three or four: a shunt on two small pads sees perhaps thirty kelvin per watt into the board, and the same shunt on pads with an acre of copper and a dozen thermal vias sees eight. That is a factor of nearly four on every error in the table, bought with board area and nothing else — no extra part, no calibration, no bandwidth cost, which is unusual enough in this field to be worth stating on its own.

It also means the error is not a property of the shunt. Two identical parts on two boards have different non-linearities, and a specification written from the data sheet alone is not a specification of the instrument. The quantity to measure is the rise per ampere, which is θu\theta u^* and has no resistance in it, so it can be measured once on the assembly with a thermocouple and a current source and used for every range the instrument has.

The refusal, and what it is protecting

With the temperature coefficient set to zero the loop is not a fixed point at all. One pass through it returns the cold resistance and so does every subsequent pass, at every current, and the curvature in the reading is absent rather than small.

That is worth claiming because the alternative — a very small coefficient — looks the same in a figure and is not the same statement. A zero-coefficient shunt is exactly linear, so a single-point calibration is exactly right, and an instrument built on one has no thermal non-linearity to specify. A ten-parts-per-million shunt has one of 155 ppm at a hundred amps and needs it in the specification.

The distinction is the difference between a term that is negligible and a term that is not there, and it is kept deliberately. It is the same distinction as a reactance is not warm — where a lossless capacitor’s noise is zero rather than small — and it matters for the same reason: a term that is absent stays absent when the conditions change, and a term that is small does not.

What one thermal resistance does not describe

That twenty kelvin per watt is the right thermal resistance. It is a reasonable figure for a surface-mount shunt on an ordinary board and it is the one thing in this analysis that belongs to the assembly rather than to the part. On a heavy copper plane with thermal vias it might be five; hanging in free air it might be a hundred. Every temperature rise and every error above scales in direct proportion, and the runaway current inversely.

That the thermal path is a single resistance. It is not — it is a network with several time constants, which is what two loops, and one heatsink measures — and the steady state computed here is the limit of that network after long enough. A current pulse shorter than the shunt’s own thermal time constant produces a smaller rise and a smaller error, and the figure’s numbers are the worst case.

That the temperature coefficient is a constant. Shunt alloys are chosen because their coefficient is small and because it is roughly parabolic with a maximum near room temperature, so the effective coefficient over a rise is smaller than the specified one at one end of the range and larger at the other. The linear model here is the first term of that and is honest about being one.

That the current is the only heat source. A shunt on a board near a power device is warmed by it, and that rise produces exactly the same resistance error with none of the current dependence — an offset in the gain rather than a curvature. The two are distinguishable only by measuring at two currents.

A fixed point iterated and in closed form, and the exponent fitted

The fixed point is solved twice. By iteration to a part in 101510^{15}, and by the closed form R0/(1αθI2R0)R_0/(1 - \alpha\theta I^2R_0), agreeing to 2×10162\times10^{-16} over five decades of current.

The exponent is fitted rather than claimed, over five decades, and required to be one within two per cent — and separately checked never to be below one, because the fixed point feeds its own output back and any departure must be upward.

The temperature rise per ampere is checked against θu\theta u^* to a part in a thousand, which is the statement that the resistance has cancelled and is the whole reason the error is linear.

The runaway current is checked to be above a kiloampere at every coefficient on the slider, so that “the fixed point never fails for a shunt” is a measurement rather than an expectation.

And the loop is refused at zero coefficient, where it returns the cold resistance at every current and the curvature is absent.

The one error with the current in it

Four errors, and the first three do not know what current is flowing.

That is a strange property for a current measurement to have and it is the reason the burden-voltage result is worth carrying: it means one number sizes every shunt, from a milliamp to a kiloamp, and the design does not have to be redone for each range. The three errors it balances are properties of the instrument rather than of the measurement.

The fourth is a property of the measurement, and its arrival changes what kind of specification the instrument has. Three current-free errors give a fixed accuracy — 0.129 per cent, everywhere. A fourth error proportional to the current gives an accuracy that degrades across the range, and one proportional to its square would give a curvature that no calibration removes. The optimum turns the second of those into the first, which is the best available outcome short of removing the term.

What it cannot do is remove it, because the term is a consequence of measuring a current by dissipating some of it. Every instrument that works that way has it. A current transformer does not — it dissipates nothing in the primary, which is the whole of its appeal and which the ammeter that is not in the circuit prices in a different currency entirely.

So the three essays after the first one have found three quantities the optimum handed back: a bandwidth that falls with the current, a burden that scales with the rail, and a non-linearity that grows with the current. The optimum is still the right place to start and it is still resistance-free and current-free. What it is not is the whole design, and the pattern across the three is the same: a constraint that eliminates a variable from one quantity puts it into the others.

Still open: the transient rise, the alloy whose coefficient is not linear, and the four-wire connection that does not help

The pulse shorter than the thermal time constant. Everything here is a steady state. A shunt has a thermal time constant of some tens of milliseconds to its board and perhaps a millisecond to its own element, so a current pulse of a hundred microseconds produces a rise set by the element’s heat capacity alone and much smaller than the steady figure. Solving the thermal network rather than the single resistance would give the error as a function of pulse length, which is what a fault-current measurement actually needs.

The coefficient that is parabolic. A shunt alloy’s resistance is specified as a coefficient at twenty degrees and a curvature, and the effective coefficient over a fifteen-kelvin rise from an ambient of twenty-five is not the specified one. Substituting the real curve into the fixed point would say how much of the 775 ppm at a hundred amps is a real prediction and how much is an artefact of linearising a parabola about the wrong point.

And the four-wire connection, which repairs a different defect. Two terminals measure the leads as well removes the leads’ resistance from a measurement, and a Kelvin-connected shunt has the sense terminals inside the current terminals for exactly that reason. It does nothing whatever about self-heating: the element between the sense points is the part that is warming, and no connection scheme puts the measurement outside it. Whether a shunt could be built with a sense element that carries no current — a second, unheated element in the same thermal environment, differenced against the first — is a construction question this analysis makes precise.

Part 5 on current sensing

One argument about Current sensing, and one of 5 essays on it so far, each part numbered by how much of the idea it assumes. What sits either side of it:

The objects named here

The third axis, after the field and the idea: the things themselves, and every essay that touches each one.

Current sensingFixed pointModel rangeNonlinearitySelf heatingTemperature coefficientThermal resistance