5.7: Gauss’ Law - Differential Form
- Page ID
- 24254
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)The integral form of Gauss’ Law is a calculation of enclosed charge \(Q_{encl}\) using the surrounding density of electric flux:
\[\oint_{\mathcal S} {\bf D}\cdot d{\bf s} = Q_{encl} \label{m0045_eGLIF} \]
where \({\bf D}\) is electric flux density and \({\mathcal S}\) is the enclosing surface. It is also sometimes necessary to do the inverse calculation (i.e., determine electric field associated with a charge distribution). This is sometimes possible using Equation \ref{m0045_eGLIF} if the symmetry of the problem permits; see examples in Section 5.5 and 5.6. If the problem does not exhibit the necessary symmetry, then it seems that one must fall back to the family of techniques presented in Section 5.4 requiring direct integration over the charge, which is derived from Coulomb’s Law.
However, even the Coulomb’s Law / direct integration approach has a limitation that is very important to recognize: It does not account for the presence of structures that may influence the electric field. For example, the electric field due to a charge in free space is different from the electric field due to the same charge located near a perfectly-conducting surface. In fact, these approaches do not account for the possibility of any spatial variation in material composition, which rules out their use in many engineering applications.
To address this broader scope of problems, we require an alternative form of Gauss’ Law that applies at individual points in space. That is, we require Gauss’ Law expressed in the form of a differential equation, as opposed to an integral equation. This facilitates the use of Gauss’ Law even in problems that do not exhibit sufficient symmetry and that involve material boundaries and spatial variations in material constitutive parameters. Given this differential equation and the boundary conditions imposed by structure and materials, we may then solve for the electric field in these more complicated scenarios. In this section, we derive the desired differential form of Gauss’ Law. Elsewhere (in particular, in Section 5.15) we use this equation as a tool to find electric fields in problems involving material boundaries.
There are in fact two methods to develop the desired differential equation. One method is via the definition of divergence, whereas the other is via the divergence theorem. Both methods are presented below because each provides a different bit of insight. Let’s explore the first method:
Derivation via the Definition of Divergence
Let the geometrical volume enclosed by \({\mathcal S}\) be \({\mathcal V}\), which has volume \(V\) (units of m\(^3\)). Dividing both sides of Equation \ref{m0045_eGLIF} by \(V\) and taking the limit as \(V \to 0\):
\[\lim_{V \to 0} \frac{ \oint_{\mathcal S}{ {\bf D} \cdot d{\bf s} } }{V} = \lim_{V \to 0} \frac{ Q_{encl} }{V} \nonumber \]
The quantity on the right hand side is the volume charge density \(\rho_v\) (units of C/m\(^3\)) at the point at which we converge after letting the volume go to zero. The left hand side is, by definition, the divergence of \({\bf D}\), indicated in mathematical notation as “\(\nabla \cdot {\bf D}\)” (Section 4.6). Thus, we have Gauss’ Law in differential form:
\[\boxed{ \nabla \cdot {\bf D} = \rho_v } \label{m0045_eGLDF} \]
To interpret this equation, recall that divergence is simply the flux (in this case, electric flux) per unit volume.
Gauss’ Law in differential form (Equation \ref{m0045_eGLDF}) says that the electric flux per unit volume originating from a point in space is equal to the volume charge density at that point.
Derivation via the Divergence Theorem
Equation \ref{m0045_eGLDF} may also be obtained from Equation \ref{m0045_eGLIF} using the Divergence Theorem, which in the present case may be written:
\[\int_{\mathcal V} \left( \nabla \cdot {\bf D} \right) dv = \oint_{\mathcal S} {\bf D}\cdot d{\bf s} \nonumber \]
From Equation \ref{m0045_eGLIF}, we see that the right hand side of the equation may be replaced with the enclosed charge:
\[\int_{\mathcal V} \left( \nabla \cdot {\bf D} \right) dv = Q_{encl} \nonumber \]
Furthermore, the enclosed charge can be expressed as an integration of the volume charge density \(\rho_v\) over \({\mathcal V}\): \[\int_{\mathcal V} \left( \nabla \cdot {\bf D} \right) dv = \int_{\mathcal V}\rho_v dv \nonumber \] The above relationship must hold regardless of the specific location or shape of \({\mathcal V}\). The only way this is possible is if the integrands are equal. Thus, \(\nabla \cdot {\bf D}=\rho_v\), and we have obtained Equation \ref{m0045_eGLDF}.
The electric field intensity in free space is \[{\bf E}({\bf r}) = \hat{\bf x}A x^2 + \hat{\bf y}B z + \hat{\bf z}C x^2 z \nonumber \] where \(A=3 \: ~\mbox{V/m}^3\), \(B=2 \: ~\mbox{V/m}^2\), and \(C=1 \: ~\mbox{V/m}^4\). What is the charge density at \({\bf r} = \hat{\bf x}2 - \hat{\bf y}2\) m?
Solution
First, we use \({\bf D} = \epsilon {\bf E}\) to get \({\bf D}\). Since the problem is in free space, \(\epsilon = \epsilon_0\). Thus we have that the volume charge density is
\[\begin{aligned} \rho _ { v } & = \nabla \cdot \mathbf { D } \\ & = \nabla \cdot \left( \epsilon _ { 0 } \mathbf { E } \right) = \epsilon _ { 0 } \nabla \cdot \mathbf { E } \\ & = \epsilon _ { 0 } \left[ \frac { \partial } { \partial x } \left( A x ^ { 2 } \right) + \frac { \partial } { \partial y } ( B z ) + \frac { \partial } { \partial z } \left( C x ^ { 2 } z \right) \right] \\ & = \epsilon _ { 0 } \left[ 2 A x + 0 + C x ^ { 2 } \right] \end{aligned} \nonumber \]
Now calculating the charge density at the specified location \({\bf r}\):
\[\begin{array} { l } { \epsilon _ { 0 } \left[ 2 \left( 3 \: \mathrm { V } / \mathrm { m } ^ { 3 } \right) ( 2 \: \mathrm { m } ) + 0 + \left( 1 \: \mathrm { V } / \mathrm { m } ^ { 4 } \right) ( 2 \: \mathrm { m } ) ^ { 2 } \right] } \\ { = \epsilon _ { 0 } ( 16 \: \mathrm { V } / \mathrm { m } ) } \\ { = 142 \: \mathrm { pC } / \mathrm { m } ^ { 3 } } \end{array} \nonumber \]
To obtain the electric field from the charge distribution in the presence of boundary conditions imposed by materials and structure, we must enforce the relevant boundary conditions. These boundary conditions are presented in Sections 5.17 and 5.18. Frequently, a simpler approach requiring only the boundary conditions on the electric potential (\(V(\bf{r})\)) is possible; this is presented in Section 5.15.
Furthermore, the reader should note the following. Gauss’ Law does not always necessarily fully constrain possible solutions for the electric field. For that, we might also need Kirchoff’s Voltage Law; see Section 5.11.
Before moving on, it is worth noting that Equation \ref{m0045_eGLDF} can be solved in the special case in which there are no boundary conditions to satisfy; i.e., for charge only, in a uniform and unbounded medium. In fact, no additional electromagnetic theory is required to do this. Here’s the solution:
\[{\bf D}({\bf r}) = \frac{1}{4\pi} \int_{\mathcal V} { \frac{{\bf r}-{\bf r}'}{\left|{\bf r}-{\bf r}'\right|^3}~\rho_v({\bf r}')~dv} \nonumber \]
which we recognize as one of the results obtained in Sections 5.4 (after dividing both sides by \(\epsilon\) to get \({\bf E}\)). It is reasonable to conclude that Gauss’ Law (in either integral or differential form) is fundamental, whereas Coulomb’s Law is merely a consequence of Gauss’ Law.