|
[Sponsors] |
July 11, 2011, 00:59 |
Problem on a sample program
|
#1 |
Member
HouKen
Join Date: Jul 2011
Posts: 67
Rep Power: 14 |
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! |
|
|
|
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 17:23 |
heat conducting in a solid domain | Rogerio Fernandes Brito | Siemens | 0 | March 18, 2008 17:23 |