Calculating Electric Flux Φ w/o typical Gaussian surface

AI Thread Summary
The discussion revolves around calculating electric flux through a uniform cube with a charge at its center, contrasting the simpler calculation for a sphere using Gauss's Law. The participant expresses difficulty in constructing the appropriate integral and understanding the electric field's behavior at different points on the cube's surface. Suggestions include using Cartesian coordinates to express the electric field and calculating the flux through one face of the cube before multiplying by six for the total flux. The importance of mastering vector calculus and the divergence theorem is emphasized as essential for understanding these concepts. Overall, the thread highlights the challenges of applying theoretical principles to complex geometries in electrostatics.
chawk
Messages
3
Reaction score
0
Hello all, new here.

In the past few weeks, I have been trying to gain a basic understanding of classical electricity and magnetism through the fantastic lectures at http://ocw.mit.edu, specifically the physics course http://ocw.mit.edu/courses/physics/8-02-electricity-and-magnetism-spring-2002/ taught by Dr. Walter Lewin.

In his explanation of electric flux (lecture 3), he derives Φ for an arbitrary sphere of radius R containing charge Q, and comes up with Φ = Q/ε0, as expected per Gauss's Law. He then explains that any closed surface containing charge Q will have that same flux.

Exploiting symmetry with a sphere, the calculation is fairly easy, so I wanted to challenge myself and see if I could calculate the flux for a uniform cube with a charge Q at its center. I am finding the problem unexpectedly difficult and I'm not even sure what kind of integral to construct.

Is this too ambitious to calculate given my little experience in this area? My calculus skills are rusty, but I am eager to learn and enjoy a challenge.

My only reasoning so far is that the total flux through the cube would be 6 times the flux through a single plane, so I am approaching the problem by trying to calculate the flux through a square on the xy-plane with 1 corner at the origin and charge Q under it. While the normal at every point is the same, simply constructing the correct expression for the changing E-field vector and the correct integral is kicking my butt :P

Any tips or suggestions to point me in the right direction are appreciated!
 
Last edited by a moderator:
Physics news on Phys.org
The fact that the electrostatic flux through any closed surface is equal to the charge enclosed divided by epsilon-naught can be easily proven using the divergence theorem and Coulomb's law (or Gauss' law in differential form).

Start with Coulomb's law:

\textbf{E}(\textbf{x})=\frac{1}{4\pi\epsilon_0}\int \rho(\textbf{x}')\frac{\textbf{x}-\textbf{x}'}{|\textbf{x}-\textbf{x}'|^3}d^3x'

Take the divergence of both sides of the equation, use the Leibniz integral rule, and the fact that

\mathbf{\nabla}\cdot\frac{\textbf{x}-\textbf{x}'}{|\textbf{x}-\textbf{x}'|^3}=4\pi\delta(\textbf{x}-\textbf{x}')

and you find that

\mathbf{\nabla}\cdot\textbf{E}(\textbbf{x})=\frac{1}{\epsilon_0}\int \rho(\textbf{x}')\delta(\textbf{x}-\textbf{x}')d^3x'=\frac{\rho(\textbf{x})}{\epsilon_0}

Which is Gauss' law in differential form (one of Maxwell's equations). From there, just apply the divergence theorem, along with the definition of charge density and you find

\begin{aligned}\Phi &= \oint_{\mathcal{S}}\textbf{E}\cdot d\textbf{a} \\ &= \int_{\mathcal{V}}(\mathbf{\nabla}\cdot\textbf{E})d^3x \\ &= \frac{1}{\epsilon_0}\int_{\mathcal{V}}\rho(\textbf{x})d^3x \\ &= \frac{Q_{\text{enc}}}{\epsilon_0}\end{aligned}

If this kind of vector calulus is too advanced for your current level, and you simply want to convince yourself with your cube example, just begin by expressing the electric field of the point charge Q at the origin in Cartesian coordinates by using the fact that the radial unit vector is

\mathbf{\hat{r}}=\frac{x\mathbf{\hat{x}}+y\mathbf{\hat{y}}+z\mathbf{\hat{z}}}{\sqrt{x^2+y^2+z^2}}

From there, just compute \textbf{E}\cdot d\textbf{a} for one face of the cube and integrate over that face.
 
You want to calculate the flux through the sides of a cube of the field of a charge placed at its center.
You were right while reasoning that the total flux through the cube would be 6 times the flux through a single side.
So, you have got to calculate the flux through one side.

Take, say, the top surface.
Divide it into rectangular strips.
Let x be constant for each strip. Let each strip be of thickness, dx.
(the length of the strip = a is along y and the width = dx is along x).

Choose one strip. Let it be at position 'x'.
You'll have to further divide this strip (differential elements with sides dy & dx).
Choose one differential element. Let it be at the point (x,y).
Find out the value of E at this point. To make calculations simple choose the center of the square to be the origin. To calculate the flux, you'll need to find the normal component of E at this point. You should get {(\frac{{a}^{2}}{4} + {x}^{2} + {y}^{2}})}^{\frac{3}{2}} in the denominator.
Integrate this w.r.t y to find the flux through the strip.
Integrate this again w.r.t x to find the flux through the side.
 
Thank you both for the replies!

To gabbagabbahey, yeah that vector calculus is a bit beyond my skill, though I did still look into it and it gives me an idea of what math I can look forward to learning :)

I realize now this thread really is more mathematical than physics, so I'll definitely be re-working my approach.

Thank you graphene for that push in the right direction. Seems like I'm almost there!
 
If you take a look online for 'divergence of a vector' you'll find that the case of flux from a cube is a fundamental concept called 'divergence', which is pretty self-explanatory.

It's not a difficult bit of math (because the cube considered is infinitesimal you can forget about all the complication of taking strips and sub-elements - just consider a simple cube with an average flux through the centre of each face. All the complications vanish in the limit)

If you can master that concept you'll be a long way down the road.
 
This is from Griffiths' Electrodynamics, 3rd edition, page 352. I am trying to calculate the divergence of the Maxwell stress tensor. The tensor is given as ##T_{ij} =\epsilon_0 (E_iE_j-\frac 1 2 \delta_{ij} E^2)+\frac 1 {\mu_0}(B_iB_j-\frac 1 2 \delta_{ij} B^2)##. To make things easier, I just want to focus on the part with the electrical field, i.e. I want to find the divergence of ##E_{ij}=E_iE_j-\frac 1 2 \delta_{ij}E^2##. In matrix form, this tensor should look like this...
Thread 'Applying the Gauss (1835) formula for force between 2 parallel DC currents'
Please can anyone either:- (1) point me to a derivation of the perpendicular force (Fy) between two very long parallel wires carrying steady currents utilising the formula of Gauss for the force F along the line r between 2 charges? Or alternatively (2) point out where I have gone wrong in my method? I am having problems with calculating the direction and magnitude of the force as expected from modern (Biot-Savart-Maxwell-Lorentz) formula. Here is my method and results so far:- This...
Back
Top