5. Steady Electric Currents, Conduction and Grounding

Introduction and context

Electrostatics describes electric fields produced by stationary charges. In practical electrical systems, however, charges frequently move through conductors, semiconductors, electrolytes, insulation, soil, and biological materials. When the macroscopic quantities associated with this motion do not vary with time, the system is described as a steady-current field. The analysis of such fields connects electric-field theory with resistance, leakage, earthing, insulation coordination, power dissipation, and electrothermal design.

Steady-current problems are governed by a particularly useful combination of Maxwell’s equations, the constitutive relation between electric field and current density, and boundary conditions at material interfaces. Although the familiar circuit equation V = IR is often sufficient for a simple conductor, field methods are needed when geometry is nonuniform, materials are anisotropic, current spreads through a volume, or current enters the earth. Grounding systems, cable insulation, contact resistance, electrochemical electrodes, and semiconductor devices are examples in which a circuit-only description may conceal important physical effects.

This chapter develops the field theory of conduction, derives resistance expressions from current density, examines continuity and boundary conditions, and applies these ideas to grounding, leakage, insulation, and electrothermal problems.

Learning objectives

After studying this chapter, the reader should be able to:

1. distinguish steady electric-current fields from electrostatic fields and time-varying electromagnetic fields;

2. derive and apply the continuity equation for charge conservation;

3. use Ohm’s law in point, integral, and circuit forms;

4. relate conductivity, resistivity, mobility, and carrier concentration;

5. calculate resistance for uniform, cylindrical, spherical, and nonuniform geometries;

6. apply boundary conditions for current density and electric field at material interfaces;

7. analyze leakage current, insulation resistance, grounding resistance, and touch or step voltage;

8. assess Joule heating and coupled electrothermal behavior;

9. identify modelling assumptions, practical hazards, and design limitations in grounding and conduction systems.

The physical basis of conduction

Charge conservation and the continuity equation

Electric current is the rate at which charge crosses a surface. If a volume contains charge and charge cannot be created or destroyed, the rate of decrease of charge within the volume must equal the net outward current. In integral form,

∮S J · dS = − d/dt ∫V ρv dV,

where J is current density in A/m² and ρv is volume charge density in C/m³. Applying the divergence theorem gives the differential continuity equation,

∇ · J = − ∂ρv/∂t. (5.1)

For a steady-current field, the charge distribution is time-invariant, so ∂ρv/∂t = 0 and

∇ · J = 0. (5.2)

This does not mean that current density is uniform. It means that current has no net accumulation at any point. Current may spread, converge, or change direction, but every closed control volume has zero net current inflow imbalance. Equation (5.2) is the field form of Kirchhoff’s current law.

The continuity equation is essential when analysing cable junctions, current crowding, grounding electrodes, and conduction through composite media. A proposed current-density distribution that violates ∇ · J = 0 cannot represent a physically possible steady state unless charge is accumulating somewhere.

Ohm’s law in point form

For an isotropic, linear conducting material,

J = σE, (5.3)

where σ is conductivity in S/m and E is electric field intensity in V/m. The reciprocal quantity ρ = 1/σ is resistivity, measured in Ω·m. Consequently,

E = ρJ. (5.4)

The word “Ohm’s law” can refer either to this local constitutive relation or to the circuit equation V = IR. The latter follows from the former only after the field and geometry have been considered. In anisotropic materials, conductivity is a tensor and the more general relation is

J = [σ]E.

The vectors J and E then need not be parallel. This is important in crystalline materials, laminated structures, carbon composites, and engineered conductive media.

For a homogeneous conductor in steady conditions, E is conservative in the usual low-frequency approximation, so

∇ × E = 0,

and a scalar potential V may be defined by E = −∇V. Combining this expression with Ohm’s law and continuity gives

∇ · (σ∇V) = 0. (5.5)

For constant conductivity, equation (5.5) becomes Laplace’s equation, ∇²V = 0. Thus many steady-conduction problems have the same mathematical form as electrostatic potential problems, although their physical interpretation and boundary conditions differ.

Microscopic interpretation of conductivity

In a simple electronic conductor, an applied electric field produces a small average drift velocity superimposed on the random thermal motion of carriers. A basic carrier model gives

J = nqμE,

where n is carrier concentration, q is carrier charge magnitude, and μ is mobility. Hence,

σ = nqμ. (5.6)

In metals, conductivity is high because the concentration of mobile electrons is large. In semiconductors, conductivity depends strongly on temperature, doping, and the mobilities of electrons and holes. In electrolytes, ions carry current and chemical reactions may occur at electrodes. In gases, conduction may require ionization and can become nonlinear, particularly near breakdown.

This microscopic picture explains why conductivity is not a universal constant. Temperature, mechanical stress, frequency, moisture, impurities, electric-field magnitude, and material ageing can all modify it. For many metals over a moderate temperature range,

ρ(T) ≈ ρ0[1 + α(T − T0)],

where α is the temperature coefficient of resistivity. Metals generally have positive α, whereas some semiconducting materials have a negative temperature coefficient over practical ranges.

From current density to resistance

The current through a surface is

I = ∫S J · dS. (5.7)

For a uniform conductor of length L and constant cross-sectional area A, with J parallel to the area normal,

I = JA = σEA.

Since the voltage drop is V = EL, the resistance is

R = V/I = L/(σA) = ρL/A. (5.8)

This result is a consequence of field theory, not an independent law. It shows why long, thin conductors have greater resistance than short, thick conductors, and why the same material may have very different resistance in different geometrical forms.

Nonuniform conductors

If area varies along a one-dimensional path, the differential resistance is

dR = ρ(x) dx/A(x),

so that

R = ∫0L ρ(x)/A(x) dx. (5.9)

For radial conduction, the area changes with position and must be included in the integration. For a cylindrical shell carrying current radially between radii a and b and having length l,

R = ρ/(2πl) ln(b/a). (5.10)

For spherical radial conduction between radii a and b,

R = ρ/(4π) (1/a − 1/b). (5.11)

The spherical expression is central to the idealized analysis of a point or hemispherical grounding electrode. If current flows from a hemispherical electrode of radius a into a uniform half-space, the approximate spreading resistance is

R_g ≈ ρ/(2πa). (5.12)

Actual grounding resistance depends on electrode shape, depth, soil stratification, moisture, temperature, frequency, and neighbouring electrodes. Equation (5.12) is therefore a useful estimate, not a universal design value.

Resistance, conductance, and networks

Conductance is G = 1/R and, for a uniform conductor, G = σA/L. Conductors in series carry the same current and their resistances add. Conductors in parallel have the same voltage and their conductances add. These circuit rules are valid when each region can be represented by a lumped resistance. Where current distribution is two- or three-dimensional, a field solution or numerical method is more appropriate.

Boundary conditions in conduction

Consider an interface between two conducting media with conductivities σ1 and σ2. Let n be the unit normal from medium 1 to medium 2. Charge conservation applied to a thin pillbox gives the normal-current condition

n · (J2 − J1) = −∂ρs/∂t.

Under steady conditions, if no surface charge is accumulating,

J1n = J2n. (5.13)

The tangential electric field is continuous in the electroquasistatic approximation:

E1t = E2t. (5.14)

For isotropic media, J = σE. Thus the normal electric fields may be discontinuous when conductivities differ:

E2n/E1n = σ1/σ2,

provided the normal current density is continuous. A highly conductive region can therefore carry the same normal current with a smaller electric-field gradient than a poorly conductive region.

At an insulating boundary, Jn = 0. At a perfect conductor, the internal electric field is ideally zero and the potential is constant. At an electrode or contact, the condition may instead be specified as a potential, a total current, or a contact law such as

Jn = gc(Velectrode − Vsurface),

where gc is an interfacial conductance per unit area. Contact resistance is often decisive in busbars, battery terminals, grounding connections, and power-electronic assemblies.

Leakage, insulation, and dielectric conduction

No practical insulator has infinite resistance. Leakage current may pass through the bulk material, along a contaminated surface, across joints, or through imperfect interfaces. For a uniform insulation layer of thickness d and area A,

Ileak = σA V/d,

and the insulation resistance is Rins = d/(σA). In a cable, the equivalent radial resistance must be calculated using the cylindrical expression in equation (5.10).

Leakage is influenced by humidity, contamination, temperature, ageing, voids, mechanical damage, and applied voltage. Surface leakage may dominate bulk leakage in outdoor equipment because moisture and deposits create a conductive film. Insulation testing therefore requires controlled test conditions and careful interpretation. A high resistance at one moment does not prove long-term reliability if the material has strong temperature or moisture dependence.

Polarization and absorption can also affect insulation measurements. When a voltage is first applied, displacement current and polarization current may accompany conduction current. A simple resistance measurement taken immediately may therefore differ from a timed insulation-resistance or polarization-index test. The test voltage must remain within the equipment’s rating, and stored charge must be safely discharged afterward.

Grounding and earth-current fields

Grounding provides a controlled reference and a low-impedance path for fault, lightning, static, and leakage currents. Its purposes are not identical in every installation. Protective grounding seeks to limit dangerous touch voltage and operate protective devices; functional grounding establishes a reference for signal or power-system operation; lightning grounding manages impulse currents; and static-control grounding drains charge gradually.

Grounding resistance and potential distribution

When current I enters the earth through an electrode, the surrounding soil develops a potential distribution. For a point source in a homogeneous infinite medium,

V(r) = ρI/(4πr).

For a hemispherical source at the earth surface, the coefficient becomes approximately ρI/(2πr). Potential decreases with distance, but it is not uniform near the electrode. The voltage between two points separated by a distance is called a step voltage when a person’s feet occupy those points. Touch voltage is the potential difference between a grounded structure and a point on the earth surface touched by a person.

Ground potential rise is the electrode potential relative to remote earth during current injection:

GPR = Iground Rg.

A low resistance reduces GPR, but safety cannot be judged by resistance alone. Soil resistivity, fault duration, surface layers, current division, bonding, geometry, and permissible body current are also important. A ground grid can reduce surface voltage gradients by distributing current over a larger area and by bonding accessible metallic objects.

Soil resistivity and measurement

Soil is heterogeneous and its resistivity varies with moisture, temperature, salt content, compaction, and stratification. Four-terminal methods, such as the Wenner arrangement, are commonly used to estimate apparent soil resistivity while minimizing electrode-contact errors. Interpretation of measurements may require multilayer-soil models rather than a single resistivity value.

Grounding design normally includes electrode selection, conductor sizing, corrosion assessment, fault-current paths, equipotential bonding, inspection access, and verification by measurement. Parallel rods do not always provide the resistance predicted by ideal parallel-resistor formulas because their soil current fields overlap. Separation and geometry matter.

Electrothermal coupling

Whenever current flows through a resistive medium, electrical energy is converted into heat. The volumetric power density is

p = J · E = σE² = J²/σ (W/m³). (5.15)

For a lumped resistor, P = VI = I²R = V²/R. The temperature field is governed by the heat equation, which may be written in a simple isotropic form as

ρm cp ∂T/∂t = ∇ · (k∇T) + J · E,

where ρm is mass density, cp is specific heat, and k is thermal conductivity. Because resistivity commonly changes with temperature, the electrical and thermal problems are coupled. In a metal conductor, increasing temperature generally increases resistance, which can reduce current under voltage-source operation but increase losses under some current-controlled conditions.

Electrothermal analysis is needed for busbars, cable ampacity, fuses, resistive heaters, transformer windings, battery interconnects, semiconductor packages, and high-current contacts. Design must consider steady temperature, transient overload, thermal runaway, cooling, insulation temperature class, and local hot spots caused by current crowding.

Methods for solving steady-conduction problems

A sound solution begins by defining the domain, material properties, sources, and boundary conditions. Symmetry should be used whenever justified. The potential formulation, ∇ · (σ∇V) = 0, is particularly convenient for regions of known conductivity. After solving for V, calculate E = −∇V, J = σE, total current, resistance, and power density. The result should be checked against conservation of current and dimensional consistency.

Analytical methods are effective for uniform bars, coaxial structures, spherical or cylindrical symmetry, and simple layered media. Separation of variables and conformal methods can handle certain two-dimensional geometries. Numerical methods, especially the finite-element method, are preferred for irregular electrodes, anisotropic materials, complex grounding grids, current crowding, and coupled thermal analysis. Finite-difference and boundary-element approaches are also useful in suitable domains.

A numerical model is only as reliable as its assumptions. Mesh refinement is important near sharp electrode edges and narrow contacts, where current density may become strongly concentrated. Material properties should be temperature- and frequency-appropriate. The model should be validated with hand calculations, limiting cases, or measurements.

Illustrative calculations

Consider a copper conductor 20 m long with cross-sectional area 25 mm². Taking ρ = 1.72 × 10⁻⁸ Ω·m,

R = ρL/A = (1.72 × 10⁻⁸)(20)/(25 × 10⁻⁶) ≈ 0.0138 Ω.

At 100 A, the voltage drop is approximately 1.38 V and the heat generation is I²R ≈ 138 W. The calculation is valid only if the conductor temperature is close to the temperature associated with the selected resistivity and if contact resistances are negligible.

As a grounding estimate, suppose soil resistivity is 100 Ω·m and a hemispherical electrode has radius 0.10 m. Equation (5.12) gives

Rg ≈ 100/(2π × 0.10) ≈ 159 Ω.

This large value demonstrates why a small isolated electrode may be inadequate in relatively resistive soil. Increasing effective electrode area, using multiple well-spaced electrodes, installing a grid, or improving the soil-electrode interface may be necessary. Chemical treatment must be assessed for environmental impact and corrosion.

Illustrative case study: grounding of a medium-voltage substation

Background

The following is an illustrative engineering case. A medium-voltage substation is installed on soil whose measured apparent resistivity ranges from 80 to 180 Ω·m. The site contains a steel fence, transformer tank, cable trenches, and control equipment. A fault may inject 8 kA into the grounding system for a short duration.

Situation and problem

A preliminary design uses four isolated ground rods connected only to the transformer neutral. Calculations indicate a ground resistance of approximately 12 Ω. During a fault, the estimated ground potential rise could therefore approach 96 kV. Although the fault duration is short, the fence and remote metallic services may remain at different potentials, creating dangerous touch and transfer voltages.

Analysis

The issue is not simply that the resistance exceeds an arbitrary target. The electrode arrangement is small, the current paths are poorly bonded, and the potential gradient near the transformer and fence is steep. The design also fails to control voltage differences between simultaneously accessible conductive parts. Soil layering may further invalidate a single-resistivity calculation.

Approach and solution

The design team develops a buried grounding grid extending across the substation, bonds the transformer tank, switchgear, fence, cable sheaths, and structural steel to the grid, and provides controlled connections to neutrals and surge-protection equipment. A crushed-rock surface layer is included to increase the resistance between a person and the earth. The grid spacing is reduced near equipment and perimeter regions where gradients are most significant. Soil measurements are interpreted with a layered model, and fault-current division between overhead shield wires, cable screens, and the grid is included.

A finite-element or equivalent grounding analysis is used to map surface potentials and calculate touch and step voltages for credible fault locations. Connections are specified for mechanical strength, corrosion resistance, inspection, and continuity. The completed installation is tested, and the results are compared with the design model.

Outcome

The improved grid may not eliminate ground potential rise, but it substantially equalizes accessible metalwork and reduces local surface gradients. Protective clearing time, bonding quality, surface resistivity, and current division become part of the safety solution rather than relying solely on a low measured resistance.

Learning points

Grounding is a distributed field problem. Resistance is an important parameter, but touch and step voltage, fault duration, geometry, bonding, and soil structure determine practical safety. Measurements and models should be used together, and the installation must be maintained because corrosion, construction changes, and damaged bonds can invalidate the original design.

Challenges, risks, and common mistakes

Several errors recur in steady-conduction analysis. Treating every conductor as a lumped resistor can conceal current crowding and local heating. Using resistivity at an incorrect temperature produces inaccurate voltage-drop and loss estimates. Confusing soil resistance with electrode resistance, or assuming rods behave as ideal parallel resistors, can lead to unsafe grounding conclusions. Another mistake is to calculate ground resistance without evaluating touch and step voltage.

Boundary conditions are often oversimplified. A current-injection boundary must have a physically consistent return path, and an insulating boundary must have zero normal current rather than necessarily zero electric field. Singular fields at sharp corners may be mathematical idealizations; real electrodes have finite radii, contact impedances, and material limits.

Grounding and insulation work also involves serious hazards: induced voltages, transferred potentials, stored charge, arc flash, and unexpected parallel current paths. Testing should be performed by competent personnel using isolation, discharge, barriers, appropriate instruments, and documented procedures. Protective conductors must not be disconnected casually, and grounding arrangements must comply with applicable electrical codes and utility requirements.

Professional and emerging practice

Modern engineering increasingly combines field simulation, automated design checks, distributed sensing, and condition monitoring. Thermal cameras, leakage-current monitors, online partial-discharge systems, and earth-current sensors can identify degradation before failure. Digital models can compare grounding-grid configurations, estimate seasonal soil effects, and couple fault-current calculations with thermal and mechanical limits.

Data-driven methods and machine learning may assist anomaly detection in insulation and grounding systems, but they do not replace physical models, calibrated measurements, or engineering judgement. Training data may not represent rare faults, and a prediction system must not be treated as a protective device unless its reliability has been demonstrated and accepted by the relevant safety framework. Emerging applications include conductive concrete, corrosion-resistant grounding materials, high-voltage direct-current earth-return studies, battery-grounding diagnostics, and electrochemical energy systems. These developments make accurate treatment of interfaces, transients, corrosion, and environmental effects increasingly important.

Best-practice checklist

A professional analysis should:

• define whether the problem is steady, quasi-static, or transient;
• identify all conducting regions, interfaces, sources, and return paths;
• use conductivity or resistivity values appropriate to temperature, frequency, moisture, and material condition;
• enforce charge continuity and verify total current balance;
• distinguish bulk, surface, contact, and electrode resistance;
• evaluate current density and power density, not only total resistance;
• assess touch, step, and transferred voltages in grounding designs;
• include thermal limits and temperature dependence where losses are significant;
• validate numerical models through convergence studies, limiting cases, and measurements;
• document assumptions, test conditions, maintenance requirements, and applicable standards.

Chapter Summary

Steady electric-current fields describe charge motion for which there is no continuing accumulation of charge. Their central conservation law is ∇ · J = 0, while the constitutive relation for an isotropic conductor is J = σE. Combining these with E = −∇V produces the potential equation ∇ · (σ∇V) = 0. Resistance expressions for bars, cylindrical layers, spherical regions, and grounding electrodes follow by integrating current density over the appropriate geometry.

At material interfaces, normal current is continuous when no surface charge accumulates, while tangential electric field is continuous under the electroquasistatic approximation. These conditions explain field redistribution in composite conductors and insulation. Leakage current depends on bulk and surface conduction and is strongly affected by environmental conditions. Grounding is a distributed conduction problem in which ground potential rise, touch voltage, step voltage, soil resistivity, bonding, and fault duration are as important as ground resistance. Joule heating couples electrical and thermal fields through p = J · E. Analytical methods, numerical simulation, measurement, and professional safety practice must be integrated for reliable design.

Key Terms

Steady-current field: A conduction field in which macroscopic charge density and current distribution do not vary with time.

Current density: Electric current per unit area normal to the direction of flow, represented by J.

Continuity equation: The differential statement of charge conservation, ∇ · J = −∂ρv/∂t.

Conductivity: A material property σ relating current density to electric field in J = σE.

Resistivity: The reciprocal of conductivity, measured in Ω·m.

Mobility: Average carrier drift velocity per unit applied electric field.

Drift velocity: The average directed velocity acquired by charge carriers in an electric field.

Ohm’s law: The constitutive relation J = σE for a linear isotropic conductor, or V = IR in lumped form.

Spreading resistance: Resistance caused by current diverging or converging through a conducting medium.

Contact resistance: Resistance associated with an interface between contacting conductors or an electrode and its medium.

Leakage current: Unintended or weak conduction through insulation, surfaces, interfaces, or equipment structures.

Insulation resistance: Effective resistance opposing current through an insulation system.

Ground potential rise: The potential of a grounding system relative to remote earth during current injection.

Touch voltage: Voltage between accessible grounded metal and the earth point on which a person stands.

Step voltage: Voltage between two points on the earth surface separated by a human step.

Grounding grid: An interconnected buried conductor network used to distribute fault current and control surface potentials.

Joule heating: Heat generation caused by electrical conduction, with density p = J · E.

Current crowding: Nonuniform concentration of current near edges, contacts, constrictions, or geometric transitions.

Electrothermal coupling: Mutual dependence of electrical conduction and temperature through material properties and heat generation.

Apparent soil resistivity: Measured effective resistivity representing the response of a possibly nonuniform soil structure.

Review Questions

Short Answer Questions

1. State the continuity equation and explain its steady-state form.

2. Distinguish conductivity from resistivity.

3. Why can electric field be discontinuous at an interface between two conductors?

4. Define touch voltage and step voltage.

5. State the expression for Joule power density.

Long Answer Questions

1. Derive R = ρL/A from the point form of Ohm’s law.

2. Explain the microscopic origin of conductivity in metals and semiconductors.

3. Derive the resistance of a cylindrical insulation layer carrying radial current.

4. Discuss the role of soil resistivity, bonding, and grid geometry in grounding design.

5. Explain the assumptions and limitations of the hemispherical-electrode grounding-resistance formula.

Analytical / Application Questions

1. A conductor has length 15 m, area 10 mm², and resistivity 1.8 × 10⁻⁸ Ω·m. Calculate its resistance and heat loss at 80 A.

2. A grounding electrode injects 2 kA into soil with an estimated resistance of 4 Ω. Determine the ground potential rise and explain why this value alone does not establish safety.

3. Two conducting layers carry the same normal current density, but the second has one-quarter the conductivity of the first. Compare their normal electric fields and discuss the implication for power density.

22. Multiple Choice Questions with Answers

Question 1

For a steady conduction field with no charge accumulation, which equation must hold?

A. ∇ × J = 0

B. ∇ · J = 0

C. ∇ · E = 0 in every material

D. ∂J/∂t = ρv

Correct Answer: B. ∇ · J = 0

Question 2

A uniform conductor has its length doubled and cross-sectional area halved. Its resistance becomes

A. one-fourth of the original

B. one-half of the original

C. twice the original

D. four times the original

Correct Answer: D. four times the original

Question 3

At a steady interface with no surface charge accumulation, the normal component of current density is

A. zero in both materials

B. continuous across the interface

C. inversely proportional to area only

D. always discontinuous when conductivities differ

Correct Answer: B. continuous across the interface

Question 4

For an isotropic material, which relation gives local conduction?

A. J = σE

B. J = E/σ²

C. E = σJ

D. J = ρvE

Correct Answer: A. J = σE

Question 5

The approximate resistance of a hemispherical electrode of radius a in uniform soil of resistivity ρ is

A. ρa/(2π)

B. 2πρa

C. ρ/(2πa)

D. ρ/(4πa²)

Correct Answer: C. ρ/(2πa)

Question 6

If conductivity is reduced while normal current density remains unchanged, the local Joule power density J²/σ will

A. decrease

B. remain unchanged

C. increase

D. become zero

Correct Answer: C. increase

Question 7

Which quantity is most directly associated with the potential difference between a person’s two feet during a fault?

A. Touch voltage

B. Step voltage

C. Contact conductance

D. Insulation resistance

Correct Answer: B. Step voltage

Question 8

For radial conduction through a cylindrical insulation layer, the resistance varies with

A. the difference of radii only

B. the logarithm of the ratio of outer to inner radius

C. the square of the outer radius

D. the inverse of insulation length

Correct Answer: B. the logarithm of the ratio of outer to inner radius

Question 9

Why is a grounding grid often preferable to a single small ground rod for a substation?

A. It eliminates all fault current

B. It increases soil resistivity uniformly

C. It distributes current and reduces dangerous surface gradients

D. It makes bonding unnecessary

Correct Answer: C. It distributes current and reduces dangerous surface gradients

Question 10

Which expression represents electrical power density in a conducting material?

A. p = E/J

B. p = J · E

C. p = ∇ · E

D. p = ρvJ

Correct Answer: B. p = J · E