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

Is there an alternative method for TDMA (Thomas algorithm)

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 18, 2010, 02:29
Default Is there an alternative method for TDMA (Thomas algorithm)
  #1
Member
 
Join Date: Dec 2009
Posts: 49
Rep Power: 16
katakgoreng is on a distinguished road
I'm using TDMA to solve for the SIMPLE algorithm. For small grid, the TDMA converged quite fast but for large grid the convergence rate is killing me. Is there an alternative to TDMA that is simple to code and give faster convergence apart from Gauss elimination (consume alot of memory)?
katakgoreng is offline   Reply With Quote

Old   January 18, 2010, 13:03
Default
  #2
Senior Member
 
Join Date: Nov 2009
Posts: 411
Rep Power: 19
DoHander is on a distinguished road
Try an iterative solver - check lapack library for example, you can find some hints about what solver you can use for a linear system in the book of Peric.

BTW TDMA is just a specialized Gauss elimination for tridiagonal system.

Do
DoHander is offline   Reply With Quote

Old   January 18, 2010, 13:29
Default
  #3
Member
 
Join Date: Dec 2009
Posts: 49
Rep Power: 16
katakgoreng is on a distinguished road
Thanks alot..really appreciate your response..
katakgoreng is offline   Reply With Quote

Old   January 19, 2010, 00:33
Default
  #4
Senior Member
 
Julien de Charentenay
Join Date: Jun 2009
Location: Australia
Posts: 231
Rep Power: 17
julien.decharentenay is on a distinguished road
Send a message via Skype™ to julien.decharentenay
I am slightly puzzled by the use of the term "convergence"... The Thomas algorithm is a direct solver (i.e not iterative), as pointed out.

I think that the Thomas algorithm is part of an algorithm that itself takes a while to converge (ADI, under-relaxation, treatment of non-orthogonal terms, other). You may want to have a check of the algorithm as well.

Hope this help.
Julien
julien.decharentenay is offline   Reply With Quote

Old   January 19, 2010, 07:23
Default Krylov with Schur-complement preconditioning
  #5
jed
Member
 
Jed Brown
Join Date: Mar 2009
Posts: 56
Rep Power: 19
jed is on a distinguished road
Classical SIMPLE is really not a very robust algorithm, a much better approach is to discretize implicitly in time and solve the linear systems with a Krylov method (like GMRES) preconditioned with a Schur-complement scheme. For example,

Code:
@article{elman2008tcp,
  title={{A taxonomy and comparison of parallel block multi-level preconditioners for the incompressible Navier-Stokes equations}},
  author={Elman, H.C. and Howle, V.E. and Shadid, J. and Shuttleworth, R. and Tuminaro, R.},
  journal={Journal of Computational Physics},
  volume={227},
  number={1},
  pages={1790--1808},
  year={2008},
  publisher={Academic Press}
}
As for software, Lapack offers only dense/banded direct solvers (contrary to DoHander's suggestion). I recommend using PETSc because it offers a great deal of composability so that you can try very exotic schemes with low effort (usually without touching your source code).
jed is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Apply Multigrid method to SIMPLE Algorithm yang Main CFD Forum 1 February 25, 2006 12:28
Multigrid method for SIMPLE Algorithm yang Main CFD Forum 2 February 23, 2006 18:11
SIMPLE ALGORITHM approach-Multigrid method Skumar Main CFD Forum 0 November 14, 2005 15:16
SIMPLE Algorithm & Conjugate Gradient Method Abhijit Tilak Main CFD Forum 7 March 20, 2001 01:01
Thomas Algorithm Testing Out of Ideas Main CFD Forum 3 February 8, 2001 23:45


All times are GMT -4. The time now is 14:46.