thomas algorithm

cvkeshav · March 28, 2009, 08:14

Hey, im working on a tridiagonal solver, aka the Thomas Algorithm to solve boundary value problems for my graduate course in cfd at the university of arizona. I am writing my codes in c++. the problem with my code is that my solution oscillates at every point, one a negative value and the next being a positive value, all through my domain. this is very weird to me and i dont know what to do...
Thanks for your help in advance.
C.V.Keshav

mr_fluent · March 28, 2009, 09:48

take a look at the algo given at our wiki reference. That algorithm works well because i am using it now a days without any problems. In fact my new solver works heavly around it and it performs beautifully.

Check your coefficients carefully by printing them, or use this algo to validate small test problem first.

cvkeshav · March 29, 2009, 06:37

Hey, I found the error in my code. It is working fine now. Thanks.

March 28, 2009, 08:14	thomas algorithm	#1
cvkeshav New Member C.V.Keshav Join Date: Mar 2009 Posts: 2 Rep Power: 0	Hey, im working on a tridiagonal solver, aka the Thomas Algorithm to solve boundary value problems for my graduate course in cfd at the university of arizona. I am writing my codes in c++. the problem with my code is that my solution oscillates at every point, one a negative value and the next being a positive value, all through my domain. this is very weird to me and i dont know what to do... Thanks for your help in advance. C.V.Keshav

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
About Phase Coupled SIMPLE (PC-SIMPLE) algorithm	Yan Kai	Main CFD Forum	0	April 18, 2007 04:48
About Phase Coupled SIMPLE (PC-SIMPLE) algorithm	Yan Kai	FLUENT	0	April 14, 2007 00:17
mach-uniform algorithm for LES/DNS	ilyas	Main CFD Forum	0	February 22, 2007 11:53
Thomas Algorithm Testing	Out of Ideas	Main CFD Forum	3	February 8, 2001 23:45
Why Thomas Algorithm for tridiagonal systems?	Yogesh Talekar	Main CFD Forum	6	July 15, 1999 14:56

March 28, 2009, 09:48		#2
mr_fluent Member MrFluent Join Date: Mar 2009 Posts: 33 Rep Power: 17	take a look at the algo given at our wiki reference. That algorithm works well because i am using it now a days without any problems. In fact my new solver works heavly around it and it performs beautifully. Check your coefficients carefully by printing them, or use this algo to validate small test problem first.

March 29, 2009, 06:37		#3
cvkeshav New Member C.V.Keshav Join Date: Mar 2009 Posts: 2 Rep Power: 0	Hey, I found the error in my code. It is working fine now. Thanks.