Proof of Mean Value Property

  • Thread starter yungman
  • Start date
  • #1
4,261
44

Main Question or Discussion Point

Attached is a copy of p181 of Strauss Partial Differential Equation. This is part of proof of Mean Value Property using Green's 1st identity:
[tex]\int_Dv\nabla u+\nabla v\cdot\nabla u \;dV=\oint_A v\frac {\partial u}{\partial r}dA[/tex]

Let ##v=1## and let ##u(r,\theta,\phi)## be a harmonic function ie ##\nabla^2 u=0## and ##u## has continuous 1st and 2nd partial derivatives.
[tex]\Rightarrow\;\oint_A \frac {\partial u}{\partial r}dA=0[/tex]
This is (5) on the top of the copy of the book.

If you follow the steps, the next step shows integration of
[tex]\int_0^{2\pi}\int_0^{\pi}\frac {\partial u(a,\theta,\phi)}{\partial r}\;a^2\;\sin\theta\; d\theta \;d\phi\;=0,\hbox{ Where the integration is on surface at }\;r=a[/tex]

My question is in the very next step when the book claimed since it is equal to zero, then ##u(a,\theta,\phi)## is not a function of ##r##, therefore you can move the ##\frac{\partial}{\partial r}## outside the integral. That's all reasonable. But then the next line you see the equation using ##u(r,\theta,\phi)##, which is function of ##r##!! How can you do that?

Please explain to me why, Thanks

Alan
 

Attachments

Answers and Replies

  • #2
612
23
[...]since it is equal to zero, then ##u(a,\theta,\phi)## is not a function of ##r##, therefore you can move the ##\frac{\partial}{\partial r}## outside the integral. [...]
(Bolding by me.)

Riddle me this: is that why you can pull out the ##\frac{\partial}{\partial r}##?
 
  • #3
4,261
44
(Bolding by me.)

Riddle me this: is that why you can pull out the ##\frac{\partial}{\partial r}##?
No, the book pull it out as shown in the copy. You see the second equation using ##u_r(a,\theta,\phi)##? That's the ##\frac{\partial}{\partial r}## inside the integral.

Thanks
 
  • #4
612
23
No, the book pull it out as shown in the copy.

Thanks
I'm not familiar with the book, but I didn't see any indication of "therefore" for pulling it out. I'm betting the referenced section ("Section A.3") is something about differentiation under the integral sign.

At least, that's what I'm seeing. I haven't had my chocolate-covered espresso beans tonight, so trusting my judgement might be a little precarious. :tongue:
 
  • #5
4,261
44
I'm not familiar with the book, but I didn't see any indication of "therefore" for pulling it out. I'm betting the referenced section ("Section A.3") is something about differentiation under the integral sign.

At least, that's what I'm seeing. I haven't had my chocolate-covered espresso beans tonight, so trusting my judgement might be a little precarious. :tongue:
Did you see in the copy, right above the 3rd equation, it said "Then we pull ##\partial /\partial r## outside the integral ( see Section A.3),........

Tell me about it, the doctor tell me to cut out diet Coke.....Caffeine!!! this is my second day!!!! Who said caffeine is not addictive?!!
 
  • #6
612
23
Did you see in the copy, right above the 3rd equation, it said "Then we pull ##\partial /\partial r## outside the integral ( see Section A.3),........

Tell me about it, the doctor tell me to cut out diet Coke.....Caffeine!!! this is my second day!!!! Who said caffeine is not addictive?!!
Yes. Did you refer to Section A.3 in the book? It probably gives the justification for the pulling out ##\frac{\partial}{\partial r}## (differentiation under the integral sign).

Paul Erdős claimed a mathematician was a machine for turning coffee into theorems. :biggrin:
 
  • #7
4,261
44
Yes. Did you refer to Section A.3 in the book? It probably gives the justification for the pulling out ##\frac{\partial}{\partial r}## (differentiation under the integral sign).

Paul Erdős claimed a mathematician was a machine for turning coffee into theorems. :biggrin:
I think I found it reason. There might be a typo:

[tex]\int_0^{2\pi}\int_0^{\pi}\frac {\partial u(a,\theta,\phi)}{\partial r}\;a^2\;\sin\theta\; d\theta \;d\phi\;=0[/tex]
should be
[tex]\int_0^{2\pi}\int_0^{\pi}\frac {\partial u(r,\theta,\phi)}{\partial r}\;a^2\;\sin\theta\; d\theta \;d\phi\;=0[/tex]

The reason is this is a surface integral at ##r=a## where ##dA=a^2\;\sin\theta\; d\theta \;d\phi\;##. The integral does not involve ##r##, so even if ##u(r,\theta,\phi)## is function of ##r##, you can pull ##\frac{\partial}{\partial r}## out of the integral. Does that make sense.

Maybe it's the caffeine......or the lack of it, I should have seen this.

Thanks
 
  • #8
4,261
44
But still there is something that does not make sense. because even if you pull ##\frac {\partial}{\partial r}## out,
[tex]\int_0^{2\pi}\int_0^{\pi} u(r,\theta,\phi)\;a^2\;\sin\theta\; d\theta \;d\phi[/tex]
Still a function of ##r##. It just mean:
[tex]\frac{\partial}{\partial r}\left[\int_0^{2\pi}\int_0^{\pi} u(r,\theta,\phi)\;a^2\;\sin\theta\; d\theta \;d\phi\right]=0[/tex]
 
  • #9
34
2
Your confusion is caused just by a matter of notation.

Strauss wrote [itex]u(r,\theta,\phi)[/itex] because function [itex]u[/itex] depends on [itex]r[/itex].
What is meant is that the expression

[tex]\frac{1}{4\pi}\int_0^{2\pi}\!\!\!\int_0^\pi u(r,\theta,\phi)\,d\theta\,d\phi[/tex]

does not depend on [itex]r[/itex]. Notice it doesn't depend on [itex]\theta[/itex] and [itex]\phi[/itex], either, since we are integrating with respect to those variables. So what he actually proves is that the above integral is a constant. That is, whatever value of [itex]r[/itex]
we take in the above integration, we will obtain the same value. Now let [itex]r=0[/itex] (or more rigorously [itex]r \to 0[/itex]), and the rest of the proof follows easily.
 
  • #10
34
2
Just another detail: when [itex]\frac{\partial f}{\partial r}=0[/itex], this means that a function doesn't depend on variable [itex]r[/itex]. So, if it is a multivariate function, it will depend only on the other variables.
 
  • #11
4,261
44
Just another detail: when [itex]\frac{\partial f}{\partial r}=0[/itex], this means that a function doesn't depend on variable [itex]r[/itex]. So, if it is a multivariate function, it will depend only on the other variables.
I understand, but ##u(r,\theta,\phi)## is specified to be a function of ##r##. There lies the contradiction. The integration cannot be zero if ##u(r,\theta,\phi)## is a function of ##r##.

Put it in another way, if ##u## is not a function of ##r##, the usefulness of this Mean Value is quite useless as not too many function is not a function of ##r## in EM.
 
  • #12
34
2
It is not the integral that is zero, but the derivative of the integral divided by [itex]4\pi[/itex]. Thus, the integral is a constant (not necessarily zero). Also, [itex]u[/itex] IS a function of [itex]r[/itex]. The mean value IS NOT a function of [itex]r[/itex]
 
  • #13
4,261
44
It is not the integral that is zero, but the derivative of the integral divided by [itex]4\pi[/itex]. Thus, the integral is a constant (not necessarily zero). Also, [itex]u[/itex] IS a function of [itex]r[/itex]. The mean value IS NOT a function of [itex]r[/itex]
Thanks for the reply, I still don't get if the ##u## inside the integral is a function of ##r##, how can the integral of ##u## not a function of ##r## if the integral does the touch the ##r## part of ##u##.
 
  • #14
34
2
That's the beauty of the mean value property! It doesn't matter if you're calculating the average of a harmonic function in huge or in a tiny ball -- the average is always the same and equals the harmonic function evaluated at the center of the ball. This has a lot of consequences, besides having interesting physical interpretations, too. The mean value property is perhaps the reason why everybody loves the qualitative theory of solutions of Laplace's equations.
 
  • #15
612
23
I understand, but ##u(r,\theta,\phi)## is specified to be a function of ##r##. There lies the contradiction. The integration cannot be zero if ##u(r,\theta,\phi)## is a function of ##r##.
Consider a function $$c:\mathbb{R}^5\to\mathbb{R},~(x,y,z,u,v)\mapsto c.$$ Call it a constant function of 5 real variables.

What is ##\frac{\partial c}{\partial x}##?
 
  • #16
4,261
44
That's the beauty of the mean value property! It doesn't matter if you're calculating the average of a harmonic function in huge or in a tiny ball -- the average is always the same and equals the harmonic function evaluated at the center of the ball. This has a lot of consequences, besides having interesting physical interpretations, too. The mean value property is perhaps the reason why everybody loves the qualitative theory of solutions of Laplace's equations.
Thanks
Can you give an example of ##u## is a function of ##r##, but after the integration, it becomes independent of ##r## so the ##\partial /\partial r## =0?
 
  • #17
4,261
44
Consider a function $$c:\mathbb{R}^5\to\mathbb{R},~(x,y,z,u,v)\mapsto c.$$ Call it a constant function of 5 real variables.

What is ##\frac{\partial c}{\partial x}##?
Thanks
But you now limits the function ##u## to be a constant function. That will not be useful for electromagnetics. All the potentials and fields depend on distance ##r##. None are constant.

Green's function is part of the advanced EM material, I just don't think it meant to be a constant function.
 
  • #18
34
2
Can you give an example of u is a function of r, but after the integration, it becomes independent of r so the ∂/∂r =0?
Take [itex]u(r,\theta,\phi)=r^2\sin^2 \theta\cos 2\phi[/itex]. The function [itex]u[/itex] is not constant, depends on [itex]r[/itex] and is unbounded. Now calculate the surface integral of [itex]u[/itex] on a sphere of radius [itex]r[/itex]. You'll see that it equals zero, independently of [itex]r[/itex]. The same is true of the average of [itex]u[/itex] on that sphere.(Remark: in accordance to the notation of Strauss, take [itex]\theta \in [0,\pi][/itex] (polar angle or colatitude) and [itex]\phi \in [0,2\pi][/itex] (azimuthal angle or longitude)).

Try to understand why this happens in terms of the mean value property. (Hint: show that [itex]u[/itex] satisfies Laplace's equation.)

If you have trouble with the details, do not hesitate to let me know.
 
  • #19
4,261
44
Take [itex]u(r,\theta,\phi)=r^2\sin^2 \theta\cos 2\phi[/itex]. The function [itex]u[/itex] is not constant, depends on [itex]r[/itex] and is unbounded. Now calculate the surface integral of [itex]u[/itex] on a sphere of radius [itex]r[/itex]. You'll see that it equals zero, independently of [itex]r[/itex]. The same is true of the average of [itex]u[/itex] on that sphere.(Remark: in accordance to the notation of Strauss, take [itex]\theta \in [0,\pi][/itex] (polar angle or colatitude) and [itex]\phi \in [0,2\pi][/itex] (azimuthal angle or longitude)).

Try to understand why this happens in terms of the mean value property. (Hint: show that [itex]u[/itex] satisfies Laplace's equation.)

If you have trouble with the details, do not hesitate to let me know.
I worked through ##u=r^2\sin^2\theta\cos 2\phi##, ##\nabla^2 u \;\neq \;0##.

I am using
[tex]\nabla ^2 u=\frac {1}{r^2}\frac{\partial}{\partial r}\left(r^2\frac{\partial u}{\partial r}\right)+\frac{1}{r^2\sin\theta}\frac{\partial}{\partial\theta}\left(\sin\theta\frac{\partial u}{\partial \theta}\right)+\frac{1}{r^2\sin^2\theta}\frac{\partial^2u}{\partial \phi^2}[/tex]
I have
[tex]\nabla^2 u=-4\cos 2\phi\;\neq\;0[/tex]
So this is not a Harmonic function. This Mean Value property require ##u## being a Harmonic function.
 
Last edited:
  • #20
4,261
44
I have even more confusion: ##\oint_A \frac {\partial u}{\partial r}dA=0## is based on ##\nabla^2 u=0##.

Let ##u=\frac{1}{r}##
[tex]\nabla ^2 u=\frac {1}{r^2}\frac{\partial}{\partial r}\left(r^2\frac{\partial u}{\partial r}\right)+\frac{1}{r^2\sin\theta}\frac{\partial}{\partial\theta}\left(\sin\theta\frac{\partial u}{\partial \theta}\right)+\frac{1}{r^2\sin^2\theta}\frac{\partial^2u}{\partial \phi^2}[/tex]
[tex]\Rightarrow\;\nabla^2 u=\frac {1}{r^2}\frac{\partial}{\partial r}\left(r^2\frac{\partial \left(\frac{1}{r}\right)}{\partial r}\right)+0+0=\frac {1}{r^2}\frac{\partial}{\partial r}(-1)=0[/tex]

to check 1st and 2nd partial derivatives are continuous:
[tex]\frac {\partial u}{\partial r}+\frac {\partial u}{\partial \theta}+\frac {\partial u}{\partial \phi}=\frac {\partial u}{\partial r}=-\frac{1}{r^2}[/tex]
[tex]\frac {\partial^2 u}{\partial r^2}=2\frac{1}{r^3}[/tex]
Both are continuous for ##r\neq 0##. Therefore ##u=1/r## is a Harmonic Function.

Green's first identity:
[tex]\int_Dv\nabla u+\nabla v\cdot\nabla u \;dV=\oint_A v\frac {\partial u}{\partial r}dA[/tex]
Let ##v=1##
[tex]v=1\Rightarrow\;\int_Dv\nabla^2 u+\nabla v\cdot\nabla u \;dV=0[/tex]
But
[tex]\int_A\frac{\partial u}{\partial r}dA=-\int_A\frac{1}{r^2} a^2\sin\theta\;d\theta\;d\phi\;\neq\;0[/tex]

I can't even get the Green's first identity.
[tex]\int_Dv\nabla u+\nabla v\cdot\nabla u \;dV=\oint_A v\frac {\partial u}{\partial r}dA[/tex]

Please help.

Thanks
 
Last edited:
  • #21
34
2
Well, first I'm glad you enjoyed my little example. I'm going to show (I don't know if I'm allowed) how to prove that [itex]u=r^2 \sin^2 \theta \cos 2 \phi[/itex] is harmonic. There are two ways to do it.

First method (easier): Remember that [itex]\cos 2 \phi=\cos^2 \phi - \sin^2 \phi[/itex] so that
[itex]u=(r\sin \theta \cos \phi)^2-(r\sin \theta \sin \phi)^2[/itex]. But, using the formula for spherical coordinates (in the notation of Strauss), we get [itex]u=x^2-y^2[/itex], which is easily seen to have zero laplacian.

Second method (harder): Uses the laplacian formula in spherical coordinates, as you tried to do.

[tex]\frac{\partial u}{\partial r}=2r \sin^2 \theta \cos 2 \phi[/tex]
[tex]r^2\frac{\partial u}{\partial r}=2r^3 \sin^2 \theta \cos 2 \phi[/tex]
[tex]\frac{\partial}{\partial r}\left(r^2\frac{\partial u}{\partial r}\right)=6r^2 \sin^2 \theta \cos 2 \phi[/tex]
[tex]\frac{1}{r^2}\left(\frac{\partial}{\partial r}r^2\frac{\partial u}{\partial r}\right)=6 \sin^2 \theta \cos 2 \phi[/tex]

[tex]\frac{\partial u}{\partial \theta}=2r^2 \sin \theta \cos \theta \cos 2 \phi[/tex]
[tex]\sin \theta \frac{\partial u}{\partial \theta}=2r^2 \sin^2 \theta \cos \theta \cos 2 \phi[/tex]
[tex]\frac{\partial}{\partial \theta}\left(\sin \theta \frac{\partial u}{\partial \theta}\right)=2r^2 (2\sin \theta \cos^2 \theta -\sin^3 \theta)\cos 2 \phi[/tex]
[tex]\frac{1}{r^2 \sin \theta}\frac{\partial}{\partial \theta}\left(\sin \theta \frac{\partial u}{\partial \theta}\right)=2 (2\cos^2 \theta -\sin^2 \theta)\cos 2 \phi=4\cos^2 \theta\cos 2\phi-2\sin^2 \theta\cos 2 \phi[/tex]

[tex]\frac{\partial^2 u}{\partial \phi^2}=-4r^2\sin^2 \theta \cos 2\phi[/tex]
[tex]\frac{1}{r^2 \sin^2 \theta}\frac{\partial^2 u}{\partial \phi^2}=-4 \cos 2\phi[/tex]

Putting it all together, we have

[tex]\Delta u=4\sin^2 \theta \cos 2 \phi+4\cos^2 \theta\cos 2\phi-4 \cos 2\phi=[/tex]
[tex]=4(\sin^2 \theta+ \cos^2 \theta) \cos 2 \phi-4 \cos 2\phi=4 \cos 2\phi-4 \cos 2\phi=0.[/tex]

I'd like you to check my calculations and tell me what you think. If you disagree in any point, please warn me.

But the point here was that you could calculate the surface integral of [itex]u[/itex] (and also its average) directly, since it's in the form [itex]u=F(r)G(\theta)H(\phi)[/itex] and the region of integration is a rectangle (in the [itex](\theta,\phi)[/itex] plane). Try to do it and let me know if you got that it is zero, independently of [itex]r[/itex]. So this is the example you asked for.

Now, about the function [itex]w=1/r[/itex], your reasoning is not working because [itex]D[/itex] is the ball centered at the origin with radius [itex]a[/itex] and and [itex]w[/itex] is not defined at the origin. In fact, you'd better think of [itex]w[/itex] as satisfying not Laplace's, but Poisson's equation with forcing term equal to [itex]-4 \pi \delta[/itex], where [itex]\delta[/itex] stands for Dirac's delta. Now try to apply Green's identities with that in mind.
 
  • #22
4,261
44
Thanks, I got all the step except forgetting ##4(\cos^2\theta+\sin^2\theta)=4## to cancel the -4.
 
  • #23
4,261
44
Now, about the function [itex]w=1/r[/itex], your reasoning is not working because [itex]D[/itex] is the ball centered at the origin with radius [itex]a[/itex] and and [itex]w[/itex] is not defined at the origin. In fact, you'd better think of [itex]w[/itex] as satisfying not Laplace's, but Poisson's equation with forcing term equal to [itex]-4 \pi \delta[/itex], where [itex]\delta[/itex] stands for Dirac's delta. Now try to apply Green's identities with that in mind.
In this part, Another book use ##v(x,y)=\frac{1}{2}\ln[(x-x_0)^2+(y-y_0)^2]=\ln(r)## and proof ##v## is a Harmonic function for all points where ##(x,y)\neq(x_0,y_0)##. The situation is the same as the example I gave.

I can scan a copy if you want. It is in the book Partial Differential Equation 2nd edition by Nakhle' Asmar p623.
 
Last edited:
  • #24
34
2
Thanks, but I already have a copy of that book. In fact, I even like it better than Strauss. Anyway I think now you understand there was nothing wrong in the proof by Strauss, except that he used a partial derivative when he should've used an ordinary one while proving that the integral does not depend on [itex]r[/itex]. That integral could not depend on [itex]\theta[/itex] and [itex]\phi[/itex], anyway, since those were variables of integration.

By the way, do you want another (nicer) example of an integral that doesn't depend on [itex]r[/itex]?
 
  • #25
4,261
44
Thanks, but I already have a copy of that book. In fact, I even like it better than Strauss. Anyway I think now you understand there was nothing wrong in the proof by Strauss, except that he used a partial derivative when he should've used an ordinary one while proving that the integral does not depend on [itex]r[/itex]. That integral could not depend on [itex]\theta[/itex] and [itex]\phi[/itex], anyway, since those were variables of integration.

By the way, do you want another (nicer) example of an integral that doesn't depend on [itex]r[/itex]?
Thanks, I like to see another example. I actually like Strauss a little better. I compare the two books, you are right about ##u## is not harmonic at the center and Strauss has to cut that out. Asmar never said that. The explanation in Strauss are very short, but it is very to the point. Asmar tends to go around and around. Strauss is not the easiest book to read, but if you follow and write it out step by step, it is really to the point.

Also I was not happy with Asmar in Chapter 3 on D'Alembert formulas.

Thanks for all your help.

Alan
 

Related Threads on Proof of Mean Value Property

  • Last Post
Replies
4
Views
4K
  • Last Post
Replies
4
Views
896
Replies
2
Views
3K
Replies
2
Views
2K
Replies
2
Views
12K
  • Last Post
Replies
0
Views
935
  • Last Post
Replies
2
Views
788
Replies
1
Views
1K
Replies
1
Views
2K
Replies
3
Views
2K
Top