[Solved] M232 Homework 7-Numerical Solution of Differential Equations II

1. Consider the following method for solving the heat equation u_t = u_xx:

• Determine the order of accuracy of this method (in both space and time).
• Suppose we take k = h² for some fixed > 0 and refine the grid. For what values of (if any) will this method be Lax-Richtmyer stable and hence convergent?

Hint: Consider the MOL interpretation and the stability region of the time-discretization being used.

• Is this a useful method?

2. The m-file m solves the heat equation u_t = u_xx using the Crank-Nicolson method. Run this code, and by changing the number of grid points, confirm that it is second-order accurate. (Observe how the error at some fixed time such as T = 1 behaves as k and h go to zero with a fixed relation between k and h, such as k = 4h.)

You might want to use the function error_table.m to print out this table and estimate the order of accuracy, and error_loglog.m to produce a log-log plot of the error vs. h. See bvp_2.m for an example of how these are used.

3. Modify m to produce a new m-file heat_trbdf2.m that implements the TR-BDF2 method on the same problem. Test it to confirm that it is also second order accurate. Explain how you determined the proper boundary conditions in each stage of this Runge-Kutta method.
4. Modify m to produce a new m-file heat_FE.m that implements the forward Euler explicit method on the same problem. Test it to confirm that it is O(h²) accurate as h 0

provided when k = 24h² is used, which is within the stability limit for = 0.02. Note how many more time steps are required than with Crank-Nicolson or TR-BDF2, especially on finer grids.

5. Modify m to solve the heat equation for 1 < x < 1 with step function initial data

With appropriate Dirichlet boundary conditions, the exact solution is

where erfc is the complementary error function

erfc(

• Test this routine m = 39 and k = 4h. Note that there is an initial rapid transient decay of the high wave numbers that is not captured well with this size time step.
• How small do you need to take the time step to get reasonable results? For a suitablysmall time step, explain why you get much better results by using m = 38 than m = 39. What is the observed order of accuracy as k 0 when k = h with suitably small and m even?