Take the 2-minute tour ×
Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. It's 100% free, no registration required.

I'm trying to derive the LTE for CN applied to the linear heat equation; $u_t = u_{xx}$.

The problem is that I end up with terms of the form $\frac{{\Delta t}^k}{{\Delta x}^2}$ when using a two dimensional Taylor expansion around $(x,t)$ for the term:

${\Delta x}^2 {\delta^2_x} = (u_{i+1}^{n+1} - 2 u_{i}^{n+1} + u_{i-1}^{n+1})$

What am I doing wrong?

(trying to prove LTE for the Crandall-Douglas scheme)

share|improve this question
    
I have the same problem! I did the expansion around x_i, t_n+1/2 and those terms cancel out and we have LTE h^2+k^2 but doing the expansion around any other point doesn't give the same result. What is wrong? Then when there is any method how do we say what is LTE if no point is specified? should not method have the same LTE around any point inside of the domain? thanks! –  user9104 Apr 4 '11 at 19:04
add comment

2 Answers 2

I don't know what is the $k$ in your $\frac{\Delta t^k}{\Delta x^2}$. And I cannot tell what's wrong in your result since you didn't provide enough details.

For the LTE (Local Truncation Error) of the C-N (Crank-Nicolson scheme) for the 1-$D$ heat equation, there is a complete discussion in Chapter 2.10 in the following book

Numerical Solution of partial Differential Equations (2nd Edition) by Morton and Mayers

where Crank-Nicolson scheme is called $\theta$ method where $\theta=\frac{1}{2}$.

share|improve this answer
add comment

your taylor expansions should be around (x_i,t_n+1)

share|improve this answer
add comment

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.