Current, Current Density and Charge Conservation
~8 min read · Tier toggle changes depth (T1 plain · T2 standard · T3 rigour)
You should be comfortable with vector calculus operators, particularly the divergence and the divergence theorem, and understand the concept of charge density \(\rho\) from electrostatics.
Why this matters
Current density provides the microscopic link between moving charges and the macroscopic currents that power circuits.
While electrostatics deals with stationary charges, real-world electrical systems involve charges in motion. To describe this motion rigorously, we move beyond the simple scalar current \(I\) to the vector field \(\mathbf{J}\), which tells us the direction and magnitude of charge flow at every point in space. This transition is crucial because it allows us to apply the powerful tools of vector calculus to dynamic systems. Furthermore, the requirement that charge is neither created nor destroyed imposes a strict mathematical constraint on how \(\mathbf{J}\) and \(\rho\) can vary in time and space, leading to the continuity equation, a cornerstone of field theory.
The Flow of Charge
Imagine a river. The total amount of water passing a cross-section per second is the flow rate. In electricity, the total charge passing a cross-section per second is the current \(I\). However, just as the river might be wider in some places and narrower in others, the "density" of the flow varies. Current density \(\mathbf{J}\) is the electrical analogue of mass flux in fluid dynamics. It is a vector field defined such that its magnitude at a point is the amount of charge crossing a unit area perpendicular to the flow per unit time. If the flow is uniform and perpendicular to a surface, the total current is simply the current density multiplied by the area. If the flow is angled or non-uniform, we must integrate the component of \(\mathbf{J}\) normal to the surface over the entire area. This local description is essential because it allows us to analyze how current behaves inside conductors, not just at their boundaries.
Defining Current Density
Consider a small surface element \(d\mathbf{a}\) with an outward normal vector. The current \(dI\) passing through this element in time \(dt\) is the charge \(dq\) that crosses it. If the charge density is \(\rho\) and the velocity of the charges is \(\mathbf{v}\), the volume of charge crossing the surface in time \(dt\) is \(dV = (\mathbf{v} \cdot d\mathbf{a}) dt\). Thus, \(dq = \rho (\mathbf{v} \cdot d\mathbf{a}) dt\). Dividing by \(dt\), we find the differential current \(dI = \rho (\mathbf{v} \cdot d\mathbf{a})\). By definition, the current density \(\mathbf{J}\) is the current per unit area perpendicular to the flow, leading to the fundamental relation:
For a general surface \(S\), the total current \(I\) is the flux of \(\mathbf{J}\) through that surface:
This integral form is the bridge between the microscopic particle motion and the macroscopic circuit current. It highlights that current is not a property of a wire alone, but a property of the field of moving charges within the volume.
The Continuity Equation
Charge is a conserved quantity. It cannot appear or disappear; it can only move from one place to another. Consider a fixed volume \(V\) bounded by a closed surface \(S\). The rate of change of the total charge \(Q\) inside \(V\) must equal the net current flowing into the volume. The net current flowing *out* is the flux integral \(\oint_S \mathbf{J} \cdot d\mathbf{a}\). Therefore, the rate of change of charge inside is the negative of the outward flux:
Using the divergence theorem, we can convert the surface integral into a volume integral of the divergence of \(\mathbf{J}\):
Since this must hold for *any* arbitrary volume \(V\), the integrands must be equal at every point in space. This yields the local form of the continuity equation, which expresses local charge conservation:
This equation is a fundamental constraint on any physical solution in electromagnetism. It ensures that the fields we calculate respect the conservation of charge.
Ohm's Law in Microscopic Form
In many materials, such as metals, the current density is proportional to the electric field driving it. This is the microscopic form of Ohm's law. We define the conductivity \(\sigma\) (the reciprocal of resistivity \(\rho_{\text{res}}\)) such that:
Combining this with the continuity equation for a steady state (where \(\partial \rho / \partial t = 0\)), we get \(\nabla \cdot \mathbf{J} = 0\). Substituting Ohm's law, we find \(\nabla \cdot (\sigma \mathbf{E}) = 0\). If the material is homogeneous (\(\sigma\) is constant), this simplifies to \(\nabla \cdot \mathbf{E} = 0\). In a region with no free charge (\(\rho = 0\)), Gauss's law also states \(\nabla \cdot \mathbf{E} = 0\). Thus, in a homogeneous ohmic conductor with no net charge accumulation, the electric field is divergence-free. This consistency check reinforces the validity of linking the transport properties (conductivity) with the field equations.
The derivations
Continuity Equation and Charge Conservation
For any distribution of electric charge described by a charge density \(\rho(\mathbf{r},t)\) and current density \(\mathbf{J}(\mathbf{r},t)\), local conservation of charge is expressed by the continuity equation \(\displaystyle \nabla\cdot\mathbf{J} + \frac{\partial \rho}{\partial t} = 0\): at every point and instant, the divergence of the current density is exactly balanced by the local rate of decrease of charge density. Equivalently, the net current out of any closed surface equals the rate at which the enclosed charge falls.
Ohm's Law, Resistance and Kirchhoff's Rules for DC Circuits
For a uniform conductor of length \(L\), cross-sectional area \(A\), and resistivity \(\rho\), the potential difference \(V\) across its ends is proportional to the current \(I\) flowing through it, with proportionality constant \(R = \rho L / A\). In a closed DC network of such ohmic resistors and ideal voltage sources, the current at any junction is conserved (Kirchhoff's junction rule) and the algebraic sum of potential differences around any closed loop is zero (Kirchhoff's loop rule), allowing the complete solution of branch currents and node potentials.
A common misconception is that \(\mathbf{J} = \rho \mathbf{v}\) implies that current density is always parallel to the electric field. This is only true for isotropic, ohmic materials. In anisotropic crystals or plasmas, the relationship between \(\mathbf{J}\) and \(\mathbf{E}\) can be tensorial, meaning the current may flow at an angle to the applied field. Always verify the constitutive relation for the specific medium before assuming collinearity. C
Worked example
Problem. A cylindrical copper wire of radius \(a = 1.0 \times 10^{-3} \ \mathrm{m}\) carries a steady current \(I = 5.0 \ \mathrm{A}\). The current density is uniform across the cross-section. Calculate the magnitude of the current density \(J\) and the average drift velocity \(v_d\) of the conduction electrons, given the electron number density \(n = 8.5 \times 10^{28} \ \mathrm{m^{-3}}\) and the elementary charge \(e = 1.6 \times 10^{-19} \ \mathrm{C}\).
- Calculate the cross-sectional area \(A\) of the wire: \(A = \pi a^2 = \pi (1.0 \times 10^{-3})^2 = \pi \times 10^{-6} \ \mathrm{m^2}\).
- Determine the current density magnitude \(J\) using \(J = I/A\): \[ J = \frac{5.0}{\pi \times 10^{-6}} \approx 1.59 \times 10^6 \ \mathrm{A/m^2} \]
- Relate current density to drift velocity using \(J = n e v_d\). Note that \(J\) is the magnitude of the current flow, while \(v_d\) is the speed of the electrons. \[ v_d = \frac{J}{n e} \]
- Substitute the values: \[ v_d = \frac{1.59 \times 10^6}{(8.5 \times 10^{28})(1.6 \times 10^{-19})} \]
- Calculate the denominator: \(8.5 \times 1.6 \times 10^{9} = 13.6 \times 10^9 = 1.36 \times 10^{10}\). \[ v_d = \frac{1.59 \times 10^6}{1.36 \times 10^{10}} \approx 1.17 \times 10^{-4} \ \mathrm{m/s} \]
Units check. \(J\) is \(\mathrm{A/m^2}\). \(v_d\) is \(\mathrm{m/s}\). The result is physically reasonable; drift velocities are typically very slow compared to the speed of light, which is why the electric field propagates through the wire almost instantaneously while the electrons themselves move slowly.
Three quick tests
- If the charge density \(\rho\) in a region is increasing with time, what must be true about the divergence of the current density \(\mathbf{J}\) in that region?
Answer
Since \(\partial \rho / \partial t > 0\), the continuity equation \(\nabla \cdot \mathbf{J} + \partial \rho / \partial t = 0\) implies that \(\nabla \cdot \mathbf{J}\) must be negative. This means there is a net current flowing *into* the region. - Can the current density \(\mathbf{J}\) be non-zero in a region where the electric field \(\mathbf{E}\) is zero?
Answer
Yes. In a superconductor, current can flow without an electric field. Also, in a neutral plasma, currents can exist due to pressure gradients or magnetic forces, not just \(\mathbf{E}\). Ohm's law \(\mathbf{J} = \sigma \mathbf{E}\) is not a universal law of nature, but a constitutive relation for specific materials. - Why is the continuity equation considered a local law rather than a global one?
Answer
It is derived by considering an arbitrary infinitesimal volume element. It asserts that charge conservation holds at every single point in space and time, not just for the total charge in the universe or a specific circuit. This local form is necessary to be compatible with the finite speed of propagation of electromagnetic influences.
Problem set for PU-102 →
Each derivation in this unit ends with five graded problems with full solutions.
Lecture notes drafted with a local model and checked line by line by the editor; every number in the worked example was recomputed.
Derivations cited: D-049 · D-427 · Threads: force, energy, fields