Module

for

Dirichlet Problem for the Disk

 

Chapter 12  Fourier Series and the Laplace Transform

 

12.2  The Dirichlet Problem for the Unit Disk

        We have emphasized practical applications of harmonic functions to steady state temperatures, electrostatics, and fluid flow in Chapter 11.  

In Section 11.2 we introduced the N-Value Dirichlet problem, and showed how to solve Laplace's equation  [Graphics:Images/DirichletProblemDiskMod_gr_1.gif],  

in the upper-half plane, for a harmonic function  [Graphics:Images/DirichletProblemDiskMod_gr_2.gif] that has certain specified boundary values on the real axis.  

Now we will develop the solution to the Dirichlet problem in the closed unit disk  [Graphics:Images/DirichletProblemDiskMod_gr_3.gif].  

 

Dirichlet Problem for the Unit Disk

        Given a real valued function  [Graphics:Images/DirichletProblemDiskMod_gr_4.gif]  that is both piecewise continuous and a bounded function.  
        
Let  [Graphics:Images/DirichletProblemDiskMod_gr_5.gif]  be considered as boundary values on the unit circle [Graphics:Images/DirichletProblemDiskMod_gr_6.gif],  in the sense that  

(12.10)           [Graphics:Images/DirichletProblemDiskMod_gr_7.gif],      at points  [Graphics:Images/DirichletProblemDiskMod_gr_8.gif]  on the unit circle.

The Dirichlet problem for the closed unit disk  [Graphics:Images/DirichletProblemDiskMod_gr_9.gif]  is to extend  [Graphics:Images/DirichletProblemDiskMod_gr_10.gif]  to be

                      [Graphics:Images/DirichletProblemDiskMod_gr_11.gif]     for     [Graphics:Images/DirichletProblemDiskMod_gr_12.gif],   

where  [Graphics:Images/DirichletProblemDiskMod_gr_13.gif]  is harmonic in the unit disk and take on the boundary values (12.10)   at points where  [Graphics:Images/DirichletProblemDiskMod_gr_14.gif]  is continuous.   

        Our first method of solution uses the Fourier Series representation for  [Graphics:Images/DirichletProblemDiskMod_gr_15.gif]  that was developed in Section 12.1.  

 

Theorem 12.7 (Extended Fourier Series in the unit disk).   Let  [Graphics:Images/DirichletProblemDiskMod_gr_16.gif]  be the boundary values on the unit circle [Graphics:Images/DirichletProblemDiskMod_gr_17.gif],  

(12.10.1)      [Graphics:Images/DirichletProblemDiskMod_gr_18.gif].  

If  [Graphics:Images/DirichletProblemDiskMod_gr_19.gif]  has the Fourier series representation  

                     [Graphics:Images/DirichletProblemDiskMod_gr_20.gif],  

then

(12.11.1)      [Graphics:Images/DirichletProblemDiskMod_gr_21.gif],  

solves the Dirichlet problem and is a harmonic function in the unit disk  [Graphics:Images/DirichletProblemDiskMod_gr_22.gif].  

        We can motivate Theorem 12.7 if we are allowed to make the following claim:  

The series in (12.11.1) takes on the boundary values in Equation (12.10.1) at points on  [Graphics:Images/DirichletProblemDiskMod_gr_23.gif]  where the radial limits exist, that is

                     [Graphics:Images/DirichletProblemDiskMod_gr_24.gif]  

Furthermore, Exercise 12 in Section 3.3 shows that each of the terms,  [Graphics:Images/DirichletProblemDiskMod_gr_25.gif]   and   [Graphics:Images/DirichletProblemDiskMod_gr_26.gif]   are harmonic,

and so it is reasonable to conclude that the infinite series representing  [Graphics:Images/DirichletProblemDiskMod_gr_27.gif]  in Equation (12.11.1) will be harmonic.  

Remark.  The radial limits will exist except at the finite number of points where  [Graphics:Images/DirichletProblemDiskMod_gr_28.gif]  is discontinuous.  

There are details regarding the convergence of the series and the existence of radial limits that are left for advanced study.

Proof.

 

Example 12.3.  Find the function  [Graphics:Images/DirichletProblemDiskMod_gr_29.gif]  that is harmonic in the unit disk  [Graphics:Images/DirichletProblemDiskMod_gr_30.gif]  

and takes on the boundary values   [Graphics:Images/DirichletProblemDiskMod_gr_31.gif].  

Solution.

        Using the result of Example 12.1 in Section 12.1, we write the Fourier series for  [Graphics:Images/DirichletProblemDiskMod_gr_32.gif],   

                    [Graphics:Images/DirichletProblemDiskMod_gr_33.gif].  

Using Equation (12.11.1) for the extended Fourier series solution of the Dirichlet problem, we obtain

                    [Graphics:Images/DirichletProblemDiskMod_gr_34.gif].  

    The series representations of    [Graphics:Images/DirichletProblemDiskMod_gr_35.gif]  take on the prescribed boundary values at points where  [Graphics:Images/DirichletProblemDiskMod_gr_36.gif]  is continuous.

The boundary function  [Graphics:Images/DirichletProblemDiskMod_gr_37.gif]  is discontinuous at  [Graphics:Images/DirichletProblemDiskMod_gr_38.gif],  which corresponds to  [Graphics:Images/DirichletProblemDiskMod_gr_39.gif];  which are points where  [Graphics:Images/DirichletProblemDiskMod_gr_40.gif]  was not prescribed.

The approximations  [Graphics:Images/DirichletProblemDiskMod_gr_41.gif]   and   [Graphics:Images/DirichletProblemDiskMod_gr_42.gif],  

and the true solutions   [Graphics:Images/DirichletProblemDiskMod_gr_43.gif]   and   [Graphics:Images/DirichletProblemDiskMod_gr_44.gif]   are shown in the figures below.  

                                

                    Figure 12.13.A. The functions  [Graphics:Images/DirichletProblemDiskMod_gr_47.gif]  and   [Graphics:Images/DirichletProblemDiskMod_gr_48.gif].  

                                

                    Figure 12.13.B. The functions  [Graphics:Images/DirichletProblemDiskMod_gr_51.gif]  and   [Graphics:Images/DirichletProblemDiskMod_gr_52.gif].  

Explore Solution 12.3.

 

 

    Our next result, Theorem 12.8, is the analog to Theorem 11.3, Poisson's integral formula for the upper half-plane, that was introduced in Section 11.3.

It shows that the value of the harmonic function  [Graphics:Images/DirichletProblemDiskMod_gr_90.gif]  inside the unit disk is a special "average of the values  [Graphics:Images/DirichletProblemDiskMod_gr_91.gif]  on the boundary."  

In the integrand the function [Graphics:Images/DirichletProblemDiskMod_gr_92.gif] is multiplied by the Poisson kernel  [Graphics:Images/DirichletProblemDiskMod_gr_93.gif]  which includes the variables  [Graphics:Images/DirichletProblemDiskMod_gr_94.gif].

 

Theorem 12.8 (Poisson Integral Formula in the unit disk).   Let  [Graphics:Images/DirichletProblemDiskMod_gr_95.gif]  be the boundary values on the unit circle [Graphics:Images/DirichletProblemDiskMod_gr_96.gif],  

(12.10.2)        [Graphics:Images/DirichletProblemDiskMod_gr_97.gif].  

If  [Graphics:Images/DirichletProblemDiskMod_gr_98.gif]  is both piecewise continuous and bounded, then  

(12.11.2)        [Graphics:Images/DirichletProblemDiskMod_gr_99.gif],      [Graphics:Images/DirichletProblemDiskMod_gr_100.gif],    

solves the Dirichlet problem and is a harmonic function in the unit disk  [Graphics:Images/DirichletProblemDiskMod_gr_101.gif].  

Observation.  The integrand is taken over the unit circle where the parameter of integration is  [Graphics:Images/DirichletProblemDiskMod_gr_102.gif],  and the numerator includes the function  [Graphics:Images/DirichletProblemDiskMod_gr_103.gif].

Proof.

 

Corollary 12.1  (N-Value Dirichlet solution for the unit disk).   Assume that  [Graphics:Images/DirichletProblemDiskMod_gr_104.gif].  

If  [Graphics:Images/DirichletProblemDiskMod_gr_105.gif]  is the boundary value function on the unit circle  [Graphics:Images/DirichletProblemDiskMod_gr_106.gif],  

(12.10.3)        [Graphics:Images/DirichletProblemDiskMod_gr_107.gif]  

then   

(12.11.3)        [Graphics:Images/DirichletProblemDiskMod_gr_108.gif],  

solves the Dirichlet problem and is a harmonic function in the unit disk  [Graphics:Images/DirichletProblemDiskMod_gr_109.gif].  

Proof.  Use the indefinite integral  

                    
[Graphics:Images/DirichletProblemDiskMod_gr_110.gif].

Then

                    
[Graphics:Images/DirichletProblemDiskMod_gr_111.gif]

Remark.  When applying formula (12.11.3) it is necessary to pay attention and to use appropriate branches and to beware of branch cuts.

 

Extra Example 1.  Find the function  [Graphics:Images/DirichletProblemDiskMod_gr_112.gif]  that is harmonic in the unit disk  [Graphics:Images/DirichletProblemDiskMod_gr_113.gif],  

and takes on the boundary values   [Graphics:Images/DirichletProblemDiskMod_gr_114.gif]  

                    [Graphics:Images/DirichletProblemDiskMod_gr_115.gif]            

                    Figure 1.  The graphs of    [Graphics:Images/DirichletProblemDiskMod_gr_117.gif]    and    [Graphics:Images/DirichletProblemDiskMod_gr_118.gif].    

 

Solution using Fourier series.

In Exercise 9 in Section 12.1, we showed that the Fourier series for  [Graphics:Images/DirichletProblemDiskMod_gr_119.gif] can be written as  

                    [Graphics:Images/DirichletProblemDiskMod_gr_120.gif].  

Use Equation (12.11.3) to construct the extended Fourier series solution of the Dirichlet problem.

Therefore, the Fourier series solution is

                    [Graphics:Images/DirichletProblemDiskMod_gr_121.gif].  

Solution using Poisson's integral.

Use the known anti-derivative   [Graphics:Images/DirichletProblemDiskMod_gr_122.gif]   and get  


                    [Graphics:Images/DirichletProblemDiskMod_gr_123.gif]          
Therefore, the Poisson integral solution is

                    [Graphics:Images/DirichletProblemDiskMod_gr_124.gif].

Solution using N-Value Dirichlet formula.

        Set  [Graphics:Images/DirichletProblemDiskMod_gr_125.gif]  and  [Graphics:Images/DirichletProblemDiskMod_gr_126.gif]  and then use the formula  (12.11.3)  and get

                    [Graphics:Images/DirichletProblemDiskMod_gr_127.gif]

Therefore, the N-Value Dirichlet solution is

                    [Graphics:Images/DirichletProblemDiskMod_gr_128.gif].
         
Remark.  This is just a straightforward way of calculating the Poisson integral solution, and is recommended for working some of the exercises.

Explore Extra Solution 1.

 

Extra Example 2.  Find the function  [Graphics:Images/DirichletProblemDiskMod_gr_190.gif]  that is harmonic in the unit disk   [Graphics:Images/DirichletProblemDiskMod_gr_191.gif],  

and takes on the boundary values   [Graphics:Images/DirichletProblemDiskMod_gr_192.gif].  

                    [Graphics:Images/DirichletProblemDiskMod_gr_193.gif]            

                    Figure 2.  The graphs of    [Graphics:Images/DirichletProblemDiskMod_gr_195.gif]    and    [Graphics:Images/DirichletProblemDiskMod_gr_196.gif].    

Explore Extra Solution 2.

 

 

Remarks.  Some integrals of the form   [Graphics:Images/DirichletProblemDiskMod_gr_210.gif]  are readily available.  

For   [Graphics:Images/DirichletProblemDiskMod_gr_211.gif]   the computer algebra systems Mathematica and Maple can be used to obtain the following formulae.


                    [Graphics:Images/DirichletProblemDiskMod_gr_212.gif]
                    
                    
                    [Graphics:Images/DirichletProblemDiskMod_gr_213.gif]  
                    
                    
                    [Graphics:Images/DirichletProblemDiskMod_gr_214.gif]  

The software Maple will also compute the above integrals, but the Maple syntax will use the notation  [Graphics:Images/DirichletProblemDiskMod_gr_215.gif]  and  [Graphics:Images/DirichletProblemDiskMod_gr_216.gif].  

The first integral   [Graphics:Images/DirichletProblemDiskMod_gr_217.gif]   can easily be verified using techniques from calculus.

Caveat.  It is beyond the scope of this book to use the other two integrals involving  [Graphics:Images/DirichletProblemDiskMod_gr_218.gif].  

 

Exercises for Section 12.2.  The Dirichlet Problem for the Unit Disk

 

Library Research Experience for Undergraduates

Fourier Series

Fourier Series and Transform

Dirichlet Problem

Neumann Problem

Poisson Integral

 

 

Download This Mathematica Notebook

Download The Maple Worksheet

 

  

The Next Module is

Vibrations in Mechanical Systems

 

 

 Return to the Complex Analysis Modules

 

 

Return to the Complex Analysis Project

 

 

 

(c) 2010 John H. Mathews, Russell W. Howell