- #1
Crispus
- 1
- 0
Crank-Nicolson method for solving hyperbolic PDE?
Hi. I'm not really sure if this is the right part of the forum to ask since it's not really a home-work "problem". Anyway, the question seemed too trivial to ask in the general math forum.
What I'm wondering is wether the Crank-Nicolson method can be used with this PDE (wave equation with a source term i believe):
d2u/dt2 = c^2 * d2u/dx2 + f(t)
f(t) is known for all t.
c is constant.
I have solved the problem with forward euler but the time step has to be really small to have stability. There are actually eight equations to solve so it takes really long time. Solving one of them (using matlab) took 6 min on my computer (2.6GHz).
So I'd like to use an implicit method. And the only one I know of is the Crank-Nicolson method.
The only info I have found about the Crank-Nicolson method in textbooks or on the internet only covers the heat-flow equation.
Can I use normal central difference method for approximation of the d2u/dt2 and then proceed as normal or can this cause stability problems?
(As normal = http://sepwww.stanford.edu/sep/prof/bei/fdm/paper_html/node15.html)
Hi. I'm not really sure if this is the right part of the forum to ask since it's not really a home-work "problem". Anyway, the question seemed too trivial to ask in the general math forum.
What I'm wondering is wether the Crank-Nicolson method can be used with this PDE (wave equation with a source term i believe):
d2u/dt2 = c^2 * d2u/dx2 + f(t)
f(t) is known for all t.
c is constant.
I have solved the problem with forward euler but the time step has to be really small to have stability. There are actually eight equations to solve so it takes really long time. Solving one of them (using matlab) took 6 min on my computer (2.6GHz).
So I'd like to use an implicit method. And the only one I know of is the Crank-Nicolson method.
The only info I have found about the Crank-Nicolson method in textbooks or on the internet only covers the heat-flow equation.
Can I use normal central difference method for approximation of the d2u/dt2 and then proceed as normal or can this cause stability problems?
(As normal = http://sepwww.stanford.edu/sep/prof/bei/fdm/paper_html/node15.html)
Last edited: