Module

for

Elliptic P.D.E.'s

 

Background for Elliptic Equations

    
As examples of
elliptic partial differential equations, we consider the Laplace equation, Poisson equation, and Helmholtz equation. Recall that the Laplacian of the function  u(x,y)  is  

            
[Graphics:Images/EllipticPDEMod_gr_1.gif].  

With this notation, we can write the Laplace, Poisson, and Helmholtz equations in the following forms:  

            

[Graphics:Images/EllipticPDEMod_gr_2.gif]

  

[Graphics:Images/EllipticPDEMod_gr_3.gif]

[Graphics:Images/EllipticPDEMod_gr_4.gif]

  

[Graphics:Images/EllipticPDEMod_gr_5.gif]

[Graphics:Images/EllipticPDEMod_gr_6.gif]

  

[Graphics:Images/EllipticPDEMod_gr_7.gif]

  


It is often the case that the boundary values for the function  
u(x,y)  are known at all points on the sides of a rectangular region  R  in the plane. In this case, each of these equations can be solved by the numerical technique known as the finite-difference method.

Proof  Elliptic PDE's  Elliptic PDE's  

 

Computer Programs  Elliptic PDE's  Elliptic PDE's  

Program (Dirichlet Method for Laplace's equation)  To approximate the solution of the Laplace's equation  [Graphics:Images/EllipticPDEMod_gr_8.gif]  over the rectangle  [Graphics:Images/EllipticPDEMod_gr_9.gif]  with  [Graphics:Images/EllipticPDEMod_gr_10.gif],  [Graphics:Images/EllipticPDEMod_gr_11.gif],  for  [Graphics:Images/EllipticPDEMod_gr_12.gif],  and  [Graphics:Images/EllipticPDEMod_gr_13.gif],  for  [Graphics:Images/EllipticPDEMod_gr_14.gif].  It is assumed that  [Graphics:Images/EllipticPDEMod_gr_15.gif]  and integers  n  and  m  exist so that  [Graphics:Images/EllipticPDEMod_gr_16.gif].  

[Graphics:Images/EllipticPDEMod_gr_17.gif]

Example 1.  Solve Laplace's equation over a  9 by 9  grid with boundary conditions
    Top:       180
    Left:         80
    Bottom:   20
    Right:        0
Solution 1.

 

Example 2.  Solve Laplace's equation over a  21 by 21  grid with boundary conditions
    Top:       180
    Left:         80
    Bottom:   20
    Right:        0
Solution 2.

 

 

Program (Neumann Method for Laplace's equation)  To approximate the solution of the Laplace's equation  [Graphics:Images/EllipticPDEMod_gr_44.gif]  over the rectangle  [Graphics:Images/EllipticPDEMod_gr_45.gif]  with  [Graphics:Images/EllipticPDEMod_gr_46.gif],  [Graphics:Images/EllipticPDEMod_gr_47.gif],  for  [Graphics:Images/EllipticPDEMod_gr_48.gif],  and  [Graphics:Images/EllipticPDEMod_gr_49.gif],  for  [Graphics:Images/EllipticPDEMod_gr_50.gif].  It is assumed that  [Graphics:Images/EllipticPDEMod_gr_51.gif]  and integers  n  and  m  exist so that  [Graphics:Images/EllipticPDEMod_gr_52.gif].  

[Graphics:Images/EllipticPDEMod_gr_53.gif]

Example 3.  Solve Laplace's equation with the Neuman boundary condition over a  9 by 9  grid with boundary conditions
    Top:            180
    Left:              80
    Bottom:    [Graphics:Images/EllipticPDEMod_gr_54.gif] = 0
    Right:             0
Solution 3.

 

Improvements in the graphics.  Making better graphics which would include realistic labels is more complicated.  If you don't mind "editing pictures" then we can "patch up" the graph in Example 1 by using graphics commands which "label what ever we want."  

 

Example 4.  Solve Laplace's equation over the square [Graphics:Images/EllipticPDEMod_gr_70.gif] where the boundary conditions are  
    Top:       180  
    Left:         80  
    Bottom:   20  
    Right :        0  
Solution 4.

 

Example 5.  Solve Laplace's equation over the square [Graphics:Images/EllipticPDEMod_gr_79.gif] where the boundary conditions are  
    Top:       180
    Left:         80
    Bottom:   20
    Right:        0
Solution 5.

 

Example 6.  Solve Laplace's equation over the square [Graphics:Images/EllipticPDEMod_gr_92.gif] where the boundary conditions are  

     [Graphics:Images/EllipticPDEMod_gr_93.gif],  [Graphics:Images/EllipticPDEMod_gr_94.gif],  for  [Graphics:Images/EllipticPDEMod_gr_95.gif],  
and  
     [Graphics:Images/EllipticPDEMod_gr_96.gif],  for  [Graphics:Images/EllipticPDEMod_gr_97.gif].
Solution 6.

 

Research Experience for Undergraduates

Elliptic PDE's  Elliptic PDE's  Internet hyperlinks to web sites and a bibliography of articles.  

 

Download this Mathematica Notebook Elliptic P.D.E.'s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004