CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

Problem on a sample program

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 11, 2011, 00:59
Default Problem on a sample program
  #1
Member
 
HouKen
Join Date: Jul 2011
Posts: 64
Rep Power: 5
houkensjtu is on a distinguished road
Maybe you know the sample program:
http://www.students.ncl.ac.uk/f.j.he...webs/codes.php

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

Reply

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 11:55
Can some one show me a sample program of 3D k-epsilon? tss Main CFD Forum 2 January 8, 2010 17: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 CD-adapco 0 March 18, 2008 18:23


All times are GMT -4. The time now is 03:23.