10.34  Numerical Methods: Chem Eng

The solution to the final exam is posted under Materials:General

Boundary Condition lecture notes from Nov. 2 class are posted.

Note that almost all of the lectures by both Prof. Swan and Prof. Green have been videotaped. We have added another link to the second repository of lecture videos.

Prof. Green's short BVP notes and the Matlab program using Galerkin method and local basis set are posted. Galerkin program seems to run nicely, but convergence rate as number of basis function is increased is not so good. As shown in class, the Collocation method using cosine basis set seems to converge rapidly (?) but with the current initial guess fsolve fails at N=30. Prof. Swan suggests changing the initial guess to help fsolve converge well, e.g. use the same Tguess as in the Galerkin method, which corresponds to d = cosmat \ Tguess. It would also be worth checking if it is computing dT/dr at r=R correctly. Give it a try!

Both the draft Collocation program we wrote in class on Oct 28, and a finished version of the Collocation program have been posted under "Lecture Notes". Check them out, and give the finished version a try! Look how the solution changes between N=20 and N=30.