CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

Problem on a sample program

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   July 11, 2011, 00:59
Default Problem on a sample program
Join Date: Jul 2011
Posts: 66
Rep Power: 8
houkensjtu is on a distinguished road
Maybe you know the sample program:

It seems a very good one but when I tried to raise the grid number from 10*10 to 12*12,it leads to a divergence.The main error happened in TDMA:
//Forward substitution
for (int i=1; i < n; i++) {
gam(i) = c(i-1) / bet;
bet = b(i) - a(i) * gam(i);

if (bet == 0.0)
throw("Solver Error 2 in triag: Zero in denominator \n");
u(i) = ( r(i) - a(i) * u(i-1) ) / bet;
I can't figure out why it happen and I need your help!
houkensjtu is offline   Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Update boundary conditions calculated by an external program CedricVH OpenFOAM 2 January 15, 2010 12:55
Can some one show me a sample program of 3D k-epsilon? tss Main CFD Forum 2 January 8, 2010 18:52
could someone proose a cfd program for tha above problem kordou Main CFD Forum 0 September 3, 2009 13:57
heat conducting in a solid domain Rogerio Fernandes Brito Phoenics 0 March 18, 2008 18:23
heat conducting in a solid domain Rogerio Fernandes Brito Siemens 0 March 18, 2008 18:23

All times are GMT -4. The time now is 05:35.