The written project is due on Friday, Nov. 15 (tentative date).
This should consist of either:
Printout of power point slides, or pdf slides, or other projected images during the presentation, and an appendix (which may be handwritten) with details of calculations skipped in the presentation
A written report, if you do plan to use the chalk board during your presentation. This should be approximately 5 typewritten pages. The equations can be written if necessary.
I expect that most of you will want to use computers in some capacity. Equations can be inserted into MS Word documents. You can also do the write-up in Mathematica or PowerPoint. For the actual math part, you will probably want to use Mathematica or Matlab to plot figures. The easiest way to do a project is to use one of these programs, or some other programming language such as C++, to do some "number crunching." If you would rather do a "pencil-and-paper" calculation that is possible too.
I want you to give me a title and a short description of you project on Monday, October 28. Please talk to me before then if you have any questions about what to do, or what would be appropriate.
The project should be about something that you are interested in, and it should involve PDEs (or boundary value problems of ODEs). It would be great if there is some problem with an application in one of the sciences that has captured your interest. On the other hand, I do not want you to duplicate something you have seen elsewhere. Here are some ideas for projects. I have labeled them as pencil and paper (P&P) or computer (C) projects. Of course, the (P&P) projects will probably need computers to produce figures, and the (C) projects will involve a lot of written analysis before you go to the computer for the number crunching.
(P&P) Investigate one or more of the
unusual coordinate systems
in which the Laplacian is separable.
Use them to solve Laplace's equation, and/or to find eigenvalues of the Laplacian in certain geometries. For example, use bipolar coordinates to solve Laplace's equation on the disk with u = 1 on the upper semi-circle and u = -1 on the lower semi-circle.
Or, use bipolar coordinates to solve Laplace's equation in the region of the plane outside two circle, each held at a different temperature (or electrostatic potential).
Or, use elliptic coordinates to find the eigenvalues of the Laplacian inside and ellipse with u = 0 on the boundary.
(C) Complete problem 2.3.2g, as described in class when I showed this demonstration.
(C) Use the L matrix for approximating the 1-D Laplacian using finite differences to do one of these, numbered L1, L2, etc.
(L1) Solve the wave equation with specific non-constant density rho(x). Pick one, for example rho(x) = 1 + (x/L)^2. Your code will define the function rho(x), so that it can be easily changed and re-run. Use equation (4.2.8) with the boundary conditions(4.4.2). Use the separation of variable like in section 4.4, but solve the eigenvalue problem for phi(x) numerically using the L matrix. Show pictures of the eigenfunctions. If you have time, also solve the PDE with the initial conditions (4.4.3), which some specific choice of f(x) and g(x). This will require you to find the orthogonality condition for the eigenfunctions. Hint: the eigenvalue problem you solve is a regular Sturm-Liouville problem.
(L2) Solve the 1D heat equation with non-constant thermal conductivity and/or other properties. This would involve solving an eigenvalue problem numerically.
(L3) Solve the 1D Schroedinger equation from quantum mechanics numerically (with a potential other than Coulomb's law or Hooke's Law).
(L4) Solve the radial Schroedinger's equation for a central force other than Coulomb's law or Hooke's law. Perhaps the 6-12 rule?
The following few involve eigenvalues and eigenfunctions of the Laplacian in polar coordinates. After assuming cos(m theta) or sin(m theta) dependence, you just need to find a function of r.
(C) Find the eigenvalues of the Laplacian on the disk with fixed temp, zero derivative, and Newton's Law of cooling at the boundary.
(C) Find the eigenvalues of the Laplacian on an annulus, with u = 0 on the boundaries.
The next two projects
(C) Find the eigenvalues of the Laplacian in some 2-D region using finite difference. For example find the eigenfunction/eigenvalues on a cross using a rectangular grid, or on a hexagon or triangle with a triangular grid.
(C) Solve the 2D Schrödinger equation numerically, with a non-central force.
(C) Implement one or more of the numerical schemes described in Chapter 6 of our textbook.
(P&P) Solve this system of two coupled PDEs for heat flow in a porous medium. (See
the work of Peter Vadasz.)
cu du/dt = Ku d2u/dx2 + k(v-u)
cv dv/dt = Kv d2v/dx2 + k(u-v)
u(0, t) = v(0, t) = u(L, t) = v(L, t) = 0 for all t
u(x, 0) = f(x), and v(x, 0) = g(x), where f and g are given functions.
Jim Swift's home page
Dept of Mathematics and Statistics
NAU Home Page