Subsections

### 2.2 The Poisson equation in infinite space

(Book, section 8.2)

#### 2.2.1 Overview

This section works out the Green’s function idea for the Poisson equation

in infinite space. In two dimensions, this can be understood to be heat conduction in an infinite plate, with the temperature and the heat added to the plate per unit surface area. In three dimensions, the heat could be generated internally, due to chemical or nuclear reactions, say, or be due to absorbed radiation passing through the body. Or can be understood to be charge density in electro-statics, with being the potential. In the fluid dynamics of viscous unidirectional flows, would be the pressure gradient. In various numerical schemes for viscous incompressible flows, could be pressure and velocity terms, or in two dimensions could be streamfunction and vorticity.

In this section, the Green’s function in infinite two-dimensional space will be derived. For discussion purposes, the problem will be assumed to be heat conduction in a plate, but the mathematical solution does not depend on what the physical meaning is. In the homework you will derive the Green’s function for the Poisson equation in infinite three-dimensional space; the analysis is similar but the result will be quite different.

First of all, a Green’s function for the above problem is by definition a solution when function is a delta function. A delta function is an infinitely narrow spike that integrates to one. We will write for the point at which the temperature is desired. Further is the position of the delta function.

In two dimensions and . Also, the two-dimensional delta function can be written in terms of one-dimensional ones as

In any number of dimensions, the Green’s function solution of the Poisson equation with a delta function as right hand side will be written as

(Since it is a constant coefficient problem, the Green’s function will simplify to , and even further to because the problem is rotationally symmetric.)

In terms of the Green’s function, the solution to the Poisson equation with an arbitrary right hand side can be written as

 (2.3)

Here stands for a “volume” integral over all components of vector . In two dimensions that becomes
 (2.4)

The “volume” is here really an area.

It will be found that the Green’s function for the two-dimensional infinite-domain Poisson problem is:

 (2.5)

The physical meaning of is the distance between the point where the heat is added and the point at which the temperature is desired.

Similarly, you will find in the homework that the Green’s function for the three-dimensional infinite-domain Poisson problem is:

 (2.6)

Therefore the two-dimensional temperature distribution corresponding to a general distribution of added heat is:

 (2.7)

while the three-dimensional temperature distribution is
 (2.8)

Given an added-heat distribution , you can find the temperature distribution by doing the corresponding integral. Especially if you only want the temperature at a few points, this can be quite effective. (If you want the temperature at essentially all points, a “multigrid” numerical method that directly solves the partial differential equation is far more efficient. However, there are so-called “fast-summation” methods, like the one by Van Dommelen and Ründensteiner, that can do the integrals very fast too, especially if the region of heat addition is limited.)

The physical meaning of the Green‘s function varies with setting. In heat transfer, it is the solution for a point heat source, in electrostatics a point charge, in gravitation a point mass, in potential flows a point source of fluid, in two-dimensional vortex flows a point vortex, etcetera.

#### 2.2.2 Loose derivation

To verify the two-dimensional Green’s function given in the previous section, the solution to the Poisson equation must be found in which is a delta function spike at some point .

However, dealing with infinite functions like delta functions is a very abstract and fishy problem. Therefore an approach like in the first section will be used. It will assumed that we are really trying to solve the Poisson equation for an arbitrary function . (And so we are, really.) We then mentally cut up this function into spikes. That idea is sketched in two-dimensions in figure fig:2dspike.

The problem for such a narrow, but finite spike can be solved with some physical intuition. The solution will again be called . The total solution will then be the sum of the solutions for all the spikes. Of course, each solution will not just depend on the position at which you want to know the temperature. It will also depends on where the spike is, as indicated by its center point .

Figure 2.7 shows a two-dimensional top view equivalent to figure 2.6. In other words, it shows just the plate, not the function . The dimensions of the little rectangle to which the heat is added by the considered spike will be indicated by . The amount of heat added is , since variations in over the small rectangle can be ignored. Since the Green’s function is the solution for unit added heat flux, the final Green’s function will be obtained by dividing the solution by the amount of heat . (Formally speaking, you would then still need to take the limit , but that becomes trivial under the approximations to be made.)

The mathematical problem being solved is:

For convenience, for now use a polar coordinate system centered around the point of heat addition, as indicated in figure 2.7. Further, since the rectangle is assumed to be very small, almost a single point, you can reasonably assume that the temperature distribution depends only on the distance from the point where the heat is added, not on .

Under those assumptions, it is easiest to simply integrate the mathematical problem above over the inside of a circle of radius :

But , so you can apply the divergence theorem here:

In this two-dimensional problem the “surface” is the perimeter of the circle. And is the radial polar unit vector , which makes the total derivative equal to the derivative with respect to radius, . So you have:

Take the to the other side and integrate to get :

The Green’s function is the solution for unit heat added:

This derivation used a polar coordinate system centered around the heat addition point . To get the expression for whatever the origin of the coordinate system is, substitute

That gives the Green’s function as stated in the overview.

2.2.2 Review Questions
1

Do an analysis similar to either this subsection, or the next one, to derive the Green’s function of the Poisson equation in three dimensional infinite space.

#### 2.2.3 Rigorous derivation

For students who do not like the above derivation with infinitesimal regions, and the assumption that their temperature distribution only depends on distance, here is a mathematically solid derivation.

It will be assumed that a suitable heat distribution is given with at least continuous low-order derivatives. Also that it disappears sufficiently quickly at large distances that you do not have to worry about that region. Then it is to be shown that if you do the Green’s function integration

you get a function satisfying the Poisson equation:

The first thing to check is that you do at least get some function by doing the integration. That is not automatic, since the integrand is infinite when . And integration over an infinite region is not proper either. What you must do is exclude the inside of a very small circle around and the outside of a very large circle from the integration. Then you define to be the limit of the integral when the radius of the small circle becomes zero, and the radius of the the large circle becomes infinite. (Assuming that those limits exist.) See figure 2.8.

A local polar coordinate system will be used centered at the point at which the temperature is desired. Then is the distance that the heat addition point is away from the point at which the temperature is desired. Note that the variables in the integration are and ; is just a fixed point in this entire story. The two-dimensional Green’s function is

in these terms.

The effect of the excluded area outside the large circle may be taken to be vanishingly small if the radius of the large circle, call it , is large, since it was assumed that function vanishes sufficiently quickly at large distances. The effect of the excluded area inside the small circle is can be estimated as

where is the maximum value of within the circle and is the radius of the small circle. Now is finite since is continuous, and the integral becomes vanishingly small when . So the effect of the small circle too can be ignored if it is small enough. In particular, the limit when the radius of the small circle tends to zero and the radius of the large circle tends to infinity exists. So the temperature exists.

But does it satisfy the Poisson equation ? Now what you cannot do here is simply differentiate the Green’s function (the logarithm) within the integral

a couple of times with respect to and . If you differentiate the logarithm, it becomes a more singular function, and you get into trouble.

Instead you need to go back to the basic definition of the partial derivatives. As an example, take

The integral for ;

can be manipulated by defining to become

You can now again drop the bar on since it is just an integration variable whose name makes no difference.

Plug this, and the expression for itself, into the limit above to get:

In the limit , the term in the numerator becomes the partial derivative ! So -derivatives transform to -derivatives on inside the integral. And the equivalent thing happens to -derivatives. Since the derivatives of are assumed to be continuous just like , there are no problems with these integrals.

The net result is that

using the notations

Now you need to show that the integral in the right hand side is equal to to finish the proof that . To shorten the notations, the Green’s function will again be written as , and the question is whether

To show this, exclude again the inside of a small circle of radius and the outside of a large circle of radius around from the integration, as in figure 2.8. In that case, you do not have to worry about infinite quantities, and you can again take the limit and later to get the final result.

If you add a second term,

you can use Green’s second integral identity from section 1.2 in the book. The added term is zero, since away from the singular point . (Remember, is the solution of the Poisson problem where the forcing is a delta function, zero everywhere except at the singular point. So away from the point, and since the components of appear in in exactly the same way as those of , if , then so is .)

Green’s second identity now says that the expression that should equal is the “surface” integral

where is the derivative normal to the surface with respect to the components of . The “surface” is in this two-dimensional case the perimeter of the region of integration. Part of it is the large circle , but you can ignore that since it is assumed that vanishes sufficiently rapidly at large distances. The other part is the perimeter of the little circle that was excluded from the integration. The integral over this little circle cannot be ignored. On it, the normal derivative is minus the radial derivative , (minus since the normal vector must stick out of the region of integration, which here is into the excluded little circle.) So you have on the little circle

The perimeter of the little circle is . So you can estimate the two integrals above as

The first term is certainly no larger than times the maximum value of the gradient of on the spherical surface and disappears when you take the limit . The second term however is the average value of on the perimeter, and that becomes simply when the circle becomes so small that the variations in can be ignored.

So finally, you conclude that is indeed .