Numerical approximation: Forward diffrerence method

Click For Summary

Homework Help Overview

The discussion revolves around deriving difference equations for a numerical approximation problem involving a partial differential equation. The subject area includes numerical methods, specifically focusing on the forward difference, backward difference, and Crank-Nicholson methods in the context of heat equations.

Discussion Character

  • Exploratory, Mathematical reasoning, Assumption checking

Approaches and Questions Raised

  • Participants discuss the use of Taylor expansion to derive expressions for the first and second derivatives. There is uncertainty about the validity of obtaining a second derivative using the forward difference method, with one participant expressing confusion over arriving at zero.

Discussion Status

The discussion is ongoing, with participants exploring different methods to derive the necessary equations. Some guidance has been offered regarding the definition of the forward difference operator, but there is no explicit consensus on the approach to take for the second derivative.

Contextual Notes

Participants note that the problem specifically requires a forward difference representation for the entire problem, which may impose constraints on their usual methods of solving similar equations.

oddiseas
Messages
66
Reaction score
0

Homework Statement



[t]=-U+k[xx] u(x,0)=U(L,0)=0 u(x,0)sin(pix/L)


Write down difference equations for the approximate solution of this problem using the following methods:

1)forward difference
2)backward difference
3)crank nicholson

Homework Equations



I can do part 3, but i am stuck on the first two methods. I can find an exppression for the partial derivative of t, but the second derivative using the forward difference from a taylor approximation is 0 isn't it?

F(x+dx,t)=f(x,t)+[f][x](x,t)dx+[f][xx](x,t)(dx)^2

which if u solve for f''=0

The Attempt at a Solution

 
Physics news on Phys.org
I'm not sure about the Taylor expansion (well, I do think that the relevant equations could be derived that way too, but that's not the usual method, I think).

I'd define the forward difference operator,
[tex]\partial U_i = \frac{U_{i+1}-U_i}{h}[/tex]
and use it successively. Similarly for the backward difference. Does this help?
 
This is the definition for the first partial derivative with respect to x , using the forward difference method.I am trying to figure out if it is possible to find an expression for the second derivative with respect to x using the forward difference method. Because when i try i get zero. Usually for a heat equation i use forward difference in time and central difference for spatial, but this question specifically asks for a forward difference representation for the entire problem, so if anyone can help, because i am stuck.
 
Yeah, I know, that's why I mentioned using the operator I defined in my post several times in a row (i.e. what is [itex]\partial\partial U_i[/itex]). The result you get that way deviates from central differences only by a small amount.
 

Similar threads

  • · Replies 105 ·
4
Replies
105
Views
11K
  • · Replies 4 ·
Replies
4
Views
1K
  • · Replies 2 ·
Replies
2
Views
2K
Replies
6
Views
3K
  • · Replies 11 ·
Replies
11
Views
2K
Replies
4
Views
2K
Replies
1
Views
2K
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
Replies
3
Views
2K