- #1
Christopher314
- 6
- 0
Electrostatic conundrum: How "brute force" and symmetry arg. give different answers!
Hello everyone, this is my first post -- and an interesting one (or so I think).
**Note to Mods: this is NOT me asking for a solution to a HW type problem -- I know and have the solution already -- rather, this is an investigation on how two valid methods to answering a physics problem can have two complete different answers.**
Background: I do private tutoring, and a college student client of mine wanted help on this problem. I looked up the problem in my own copy of Griffiths and solved it and got the right answer. The thing is: I couldn't find anything wrong with how she did the problem. The problem is from Griffiths E&M book, page 106 problem 2.41: "Find E at a height z above the center of a square sheet (side a) carrying a uniform surface charge s. Answer: (s/2*e0)[(4/Pi)Arctangent[Sqrt[1+(a^2/(2z^2))]]]."
Here's the methods:
"Elegant": We can use a previous result of the field due to a finite line of constant line charge, except, now, we give it a constant surface charge s and an infinitesimal thickness. The point is to construct the square from thin square "frames." To do this we divide the square into four equal regions along the diagonals using the functions y = x and y = -x. So we have, essentially, line charges of varying lengths. And we use that by definition a square has all sides that are equal, so that if we vary x by dx, we inadvertently vary y by dy, i.e. since x = y, then dx = dy. Integrate and voila. This was my way and led to the correct answer.
**Edit: Using my client's notes here's her way exactly**
"Brute Force": My client, on the other hand, decided to go the "brute force" way and use the integral equation for the field. Since our field point is at z in the z-direction, r = (0, 0, z). On the x-y plane, any element of area dA is located at r' = (x, y, 0). So our |R| = |r - r'| = Sqrt[x^2 + y^2 + z^2]. Our element of area is dA = dxdy, a tiny square. By symmetry, the field points in the z-direction, so R = z cos Theta, where cos Theta = |r|/|R|. So then E = (zs/4 Pi e0)*Double Integral[dxdy/(x^2 + y^2 + z^2)^(3/2), where the limits of integration are -a/2<x<+a/2, and -a/2<y<+a/2 since we're integrating over a rectangle. This looks fine, right? Well if you compute this using Mathematica (or integral tables) you get a TOTALLY different answer!
The Plot Thickens: What's MORE interesting is that if, instead of integrating over the whole square and you just consider the contribution from 1/4 of the square cut along the diagonals (as in my derivation), where the limits of integration change to -a/2<x<+a/2, -x<y<+x and you multiply back in the factor of 4, then you DO get the correct answer.
The conundrum is: How can the solution change simply by changing the manner in which the surface integral is being calculated? Since this is a finite region, both integrals should be equal...but aren't.
Let's see what the forum's thoughts are. And if you'd like to compute the integrals yourself but don't have Mathematica, use this free Mathematica tool: http://integrals.wolfram.com/index.jsp"
Hello everyone, this is my first post -- and an interesting one (or so I think).
**Note to Mods: this is NOT me asking for a solution to a HW type problem -- I know and have the solution already -- rather, this is an investigation on how two valid methods to answering a physics problem can have two complete different answers.**
Background: I do private tutoring, and a college student client of mine wanted help on this problem. I looked up the problem in my own copy of Griffiths and solved it and got the right answer. The thing is: I couldn't find anything wrong with how she did the problem. The problem is from Griffiths E&M book, page 106 problem 2.41: "Find E at a height z above the center of a square sheet (side a) carrying a uniform surface charge s. Answer: (s/2*e0)[(4/Pi)Arctangent[Sqrt[1+(a^2/(2z^2))]]]."
Here's the methods:
"Elegant": We can use a previous result of the field due to a finite line of constant line charge, except, now, we give it a constant surface charge s and an infinitesimal thickness. The point is to construct the square from thin square "frames." To do this we divide the square into four equal regions along the diagonals using the functions y = x and y = -x. So we have, essentially, line charges of varying lengths. And we use that by definition a square has all sides that are equal, so that if we vary x by dx, we inadvertently vary y by dy, i.e. since x = y, then dx = dy. Integrate and voila. This was my way and led to the correct answer.
**Edit: Using my client's notes here's her way exactly**
"Brute Force": My client, on the other hand, decided to go the "brute force" way and use the integral equation for the field. Since our field point is at z in the z-direction, r = (0, 0, z). On the x-y plane, any element of area dA is located at r' = (x, y, 0). So our |R| = |r - r'| = Sqrt[x^2 + y^2 + z^2]. Our element of area is dA = dxdy, a tiny square. By symmetry, the field points in the z-direction, so R = z cos Theta, where cos Theta = |r|/|R|. So then E = (zs/4 Pi e0)*Double Integral[dxdy/(x^2 + y^2 + z^2)^(3/2), where the limits of integration are -a/2<x<+a/2, and -a/2<y<+a/2 since we're integrating over a rectangle. This looks fine, right? Well if you compute this using Mathematica (or integral tables) you get a TOTALLY different answer!
The Plot Thickens: What's MORE interesting is that if, instead of integrating over the whole square and you just consider the contribution from 1/4 of the square cut along the diagonals (as in my derivation), where the limits of integration change to -a/2<x<+a/2, -x<y<+x and you multiply back in the factor of 4, then you DO get the correct answer.
The conundrum is: How can the solution change simply by changing the manner in which the surface integral is being calculated? Since this is a finite region, both integrals should be equal...but aren't.
Let's see what the forum's thoughts are. And if you'd like to compute the integrals yourself but don't have Mathematica, use this free Mathematica tool: http://integrals.wolfram.com/index.jsp"
Last edited by a moderator: