Subsections

### 2.1 Introduction

The purpose of this section is to introduce the Green’s function ideas.

#### 2.1.1 The one-dimensional Poisson equation

This subsection will consider a very simple problem, the Poisson equation in one-dimensional infinite space. The solution will be obtained using a Green’s function approach.

In general, the Poisson equation reads

where is a given function and the unknown to be found. In one-dimensional infinite space that becomes

Of course, this equation is trivial to solve. That makes it such a good example to understand the Green’s function approach.

It may be noted that the solution is not quite unique; adding to any solution, with and constants, produces another solution. Therefore, solving the problem will simply be taken to be finding a solution, whichever one.

Figure 2.1 shows a sketch of an arbitrary given function .

The basic idea of a Green’s function approach is to chop the function into narrow spikes and solve for each spike separately.

Consider an arbitrary example spike, shown in grey in figure 2.1. Figure 2.2 shows this one spike separately. The solution due to this one spike, call it , is shown in red. The total solution can be obtained by summing the solutions for all the spikes together:

To be sure, solving the problems for the spikes exactly is just as difficult as solving the original problem. But if the spikes are narrow, approximations can be made. Before doing so however, consider the exact solution in figure 2.2 more closely. Note that the solution is linear everywhere except in the narrow region of the spike. That is because it satisfies the Poisson equation

Now , as shown in green in figure 2.2, is zero everywhere outside the spike region. And if the second derivative of is zero, then is linear.

The slope changes from one side of the spike region to the other. In fact, integration of the equation above produces

Figure 2.2 took the slopes equal and opposite at both sides

To approximate this solution, the trick is to replace the actual spike by an infinitely narrow one. That idea is sketched in figure 2.3. Note that for a valid approximation, the change in slope must be the same as that for the exact solution. If the change in slope is different, then there is no way that the approximate solution can meaningfully approximate the exact solution at both sides of the original spike. To be precise, a small error is allowed. It can be seen that it is accurate enough to approximate the exact integral over the spike as

where is the center point of the interval and is the width of the spike. The infinitely thin spike must produce a slope change equal to that amount.

To formalize the approximation, it is useful to introduce a few new concepts. First of all,

A delta function is defined as an infinitely narrow, infinitely high spike that integrates to 1.
By convention, in one dimension a delta function located at a point is indicated as .

Now the approximate spike in figure 2.3 must integrate to , not 1. So that spike is

Next

A Green’s function is defined as a solution for a delta function inhomogeneous term.
In one dimension a Green’s function is typically indicated as where is the location of the delta function. Since the delta function in figure 2.3 is multiplied by a factor , so is the Green’s function solution.

All that is left to do now is find the Green’s function. That is not difficult. Recall that the change in slope going from one side of the spike to the other must equal the integral of the spike. For a delta function, the integral is 1. So the slope must change by 1. A solution that does that is

 (2.1)

Indeed, for , the absolute signs do nothing, so the slope is . For , the absolute signs produce a change of sign so the slope is . That makes the total change in slope 1.

In summary, to solve

first chop the function into narrow spikes. Each spike produces an approximate solution

where is the center point of the spike and its width. To get the total solution , sum all these approximate solutions together. The following example explores that idea.

Example

Question: Find Green’s function approximations to the solution of the Poisson problem

Solution:

Figure 2.4 shows some results obtained using matlab. First of all, this problem has an exact solution

where erf is the so-called error function. This exact solution is indicated by the blue dots in figure 2.4.

The question is now, how good is a Green’s function approximation for this problem?

The right hand side in the Poisson problem is negligibly small outside the range , so no spikes are needed outside that range. In the top left graph, the interval was chopped up into two “spikes.” Each spike was approximated by a delta function spike at its center as described above. Then the two Green’s function solutions of these spikes were added to give the red solid line. You can see the locations of the delta functions from the kinks in this solution. Obviously, this Green’s function solution is not accurate.

It gets better if the interval is divided up into 5 narrower spikes, and each one is approximated by a delta function spike. The solution for that is shown in the top right graph of figure 2.4. The next graph shows that for 10 spikes, the Green’s function solution is quite close to the exact solution. However, there are still visible kinks at the locations of the delta function spikes. At 20 spikes, the kinks are virtually invisible.

As the example suggests, you get the exact solution to the Poisson problem by letting the width of the spikes become infinitesimal. In that case, summation becomes integration:

 (2.2)

You can verify this solution by splitting the integral into two and integrating by parts. It is somewhat messy to do so, however.

2.1.1 Review Questions
1

Solve the Poisson equation

numerically using Green’s functions.
2

Show that

is a solution to

You can assume that function becomes zero rapidly at large . (If you want, you can assume it is zero beyond some value of .)

#### 2.1.2 More on delta and Green’s functions

Figure 2.5 shows the definition of the one-dimensional Green’s function. Note that the function value of is zero at all points except at the single point . At that single point however, the function value is infinite.

Of course, infinite function values are invalid mathematics. The Green’s function is not a properly defined function. The best way to deal with that as an engineer is to mentally not make the delta function infinitely narrow. Instead think of a delta function as an extremely narrow, extremely high spike that integrates to 1. Mathematicians have better but more complicated ways of dealing with the problem, {A.1}.

Usually delta functions are used as inhomogeneous terms in differential equation problems. The solutions to these problems are called Green’s functions. Fortunately, it turns out that while the delta function is not well defined, the Green’s function typically is. In the limit that the width of the delta function becomes zero, the Green’s function stays a perfectly good function.

Like the example of the precious section illustrates, Green’s functions in infinite domains are usually not unique. The most general Green’s function for the Poisson equation in one dimension is

where and are arbitrary constants. The final two terms are a solution of the homogeneous equation.

You would typically like the Green’s function to be zero at large distances. But for the one dimensional Green’s function above, (as well as for the two-dimensional equivalent, for that matter), there is no way to do it. There is no way to choose and so that is zero at both . The best you can do is make the derivatives at as small as possible. If is nonzero, the derivative at either or is greater than in magnitude. So you take zero so that neither derivative exceeds in magnitude. There is nothing defensible that you can take for the constant , so you take it also zero.

It may be noted that in wave propagation problems, trying to make the wave function as small as possible typically does not work. Instead you take the Green’s function so that at large distances it describes waves that move away to infinity. Green’s functions that describe waves that come in from infinity are physically undesirable.