On Oct 20, 2009, at 4:34 AM, Ferry Schrijer wrote:
My apologies, it seems that I made a mistake. The solution is correct. The reason that I was doing this test is because when I use a convective heat transfer boundary condition (q = h*[T_w - T_aw], where h and T_aw are given), I do not get an accurate solution.
T_w is the value of phi at the boundary? If so, please see the FAQ: "How do I apply an outlet or inlet boundary condition?". If you still have questions, please ask.
Anyway, I will keep trying, sorry to have bothered you.
No bother.
timeStepDuration = 0.9 * dx**2 /(2*alpha)
There is no need to restrict your time step this way. You have defined a fully implicit diffusion equation, so there is no stability restriction on the time step.
