Electric Flux from Distributed Charge

From Charge Density to Electric Flux

Electric fields can be produced by charge spread continuously through a region, rather than concentrated at a single point. A volume charge density describes how much charge occupies each small piece of space. The divergence theorem turns that volume description into a statement about the total electric flux through the region’s boundary.

This connection is useful when the charge density is known and the desired quantity is the net outward flux through a closed surface. In particular, the total flux depends on the charge enclosed, even when the boundary is not spherical. A spherical boundary offers an additional advantage: if the charge distribution and field are spherically symmetric, the flux can also determine the field’s magnitude on the boundary.

The electric flux through an oriented surface \(S\) is the surface integral of the normal component of the electric field. For a closed surface, we use the outward-pointing unit normal \(\mathbf n\), so positive net flux means more field exits the surface than enters it.

Here \(\mathbf E\) is measured in newtons per coulomb, \(\mathrm{N/C}\), and flux has units \(\mathrm{N\,m^2/C}\). The surface must be closed for the enclosed-charge form of Gauss’s law used below.

The Divergence Theorem and Gauss’s Law

The divergence theorem applies to a vector field \(\mathbf E\) with continuous first partial derivatives on a region containing a bounded solid \(V\) whose boundary \(S=\partial V\) is piecewise smooth. With the outward orientation, it states that the flux through the boundary equals the integral of the divergence throughout the solid.

Gauss’s law for the electric field states that its divergence is the charge density divided by the vacuum permittivity \(\varepsilon_0\), where \(\varepsilon_0\) is measured in \(\mathrm{C^2/(N\,m^2)}\). For a smooth charge density and a field satisfying this local law throughout \(V\), the divergence theorem gives the integral law.

The implication follows by substituting \(\nabla\cdot\mathbf E=\rho/\varepsilon_0\) into the volume integral in the divergence theorem. Since \(\varepsilon_0\) is constant, it can be taken outside the integral. This is the central calculation: integrate the density over the enclosed volume, then divide by \(\varepsilon_0\).

This calculation gives total flux, not necessarily the electric field at each point of the boundary. To infer a field magnitude from the flux, additional information such as symmetry is needed. Also, the smoothness conditions matter: a singular point source inside the region cannot simply be treated as a smooth density there. The theorem must be applied only when its hypotheses hold, or with a suitable treatment of the singularity.

Worked Example: A Smooth Radial Density in a Sphere

Worked Example: A Rectangular Volume with Varying Density

Worked Example: A Linear Density on a Centered Sphere

Why Symmetry Helps—and What It Does Not Do

Symmetry can simplify the charge integral. In the centered-sphere example, reflection symmetry made the integral of \(x\) vanish. In a spherical density, spherical coordinates reduce the angular integrals to a simple factor. These are ways to compute \(Q_{\mathrm{enc}}\); they do not replace the divergence theorem.

When the entire electric field is spherically symmetric, it is also constant in magnitude on a sphere and points in the radial direction. Then \(\mathbf E\cdot\mathbf n=E(R)\) everywhere on a sphere of radius \(R\), so the flux is \(4\pi R^2E(R)\). Combined with Gauss’s law, this gives \(E(R)=Q_{\mathrm{enc}}/(4\pi\varepsilon_0R^2)\). This field conclusion requires spherical symmetry of the field, not merely a spherical boundary.

A Common Pitfall: Mixing Up Density and Charge

A frequent error is to use a density value as though it were already a charge. For example, if one writes \(Q_{\mathrm{enc}}=\rho_0\) for a uniform density, the units on the two sides are \(\mathrm{C/m^3}\) and \(\mathrm C\), respectively. They do not match. The density must be integrated over the volume:

for uniform \(\rho_0\). The resulting charge has units \((\mathrm{C/m^3})(\mathrm{m^3})=\mathrm C\). Another error is to reverse the normal: the divergence theorem uses the outward normal. Reversing it changes the sign of the surface flux, while the enclosed charge integral remains unchanged.

It is also essential to keep dimensions consistent inside a field formula. Terms added as components of \(\mathbf E\) must have the same units. In the last example, the coefficients were explicitly given units so each bracketed term had units \(\mathrm{C/m^2}\); division by \(\varepsilon_0\) then produced \(\mathrm{N/C}\). The same dimensional check should be made before differentiating or integrating a proposed physical field.