The purpose of this section is to introduce the Green’s function ideas.
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
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 19.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
19.1. Figure 19.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 19.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
The slope does change from one side of the spike region to the other.
In fact, integration of the equation above produces
If we approximate the spike area as , where
is the center point of the spike, we get
Now first consider an idealized spike problem. In this idealized spike problem, the spike is given a unit area. Then the limit is taken that the width of the spike becomes zero. (In that limit, the height of the spike must go to infinity to keep the area constant.) The limiting infinitely narrow, infinitely high, spike is called the “Dirac delta function” .
The corresponding solution is called the
Green’s function
. It is equal to
Back to the original spike of small area . We can
approximate its solution in terms of the Green’s function above as
The total solution to the original problem is obtained by
summing the contribution of all the spikes:
ExampleQuestion: Find Green’s function approximations to the solution of the Poisson problem
Use various spike widths . Verify that you do seem to get the exact solution when as claimed above.Solution:
Figure 19.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 19.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 3, so no spikes are needed outside that range. In the top left graph, the interval 3 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 3 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 19.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.
Solve the Poisson equation
Include your code.
Show that
Figure 19.5 shows the definition of the one-dimensional delta 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 delta 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.
For example, the Green's function of the Poisson equation in one
dimension formally satisfies
Another thing to note is that Green’s functions in infinite
domains are usually not unique. The most general Green’s function
for the Poisson equation in one dimension is
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.