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

fast penta diagonal solver

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 22, 2007, 14:20
Default fast penta diagonal solver
  #1
sudhakar
Guest
 
Posts: n/a
hi all,

It is obvious that when we discretise Navier-Stokes equations with central differencing scheme, the resulting matrix will be of penta-diagonal. There are number of special solvers to handle this like Stones algorithm and its variants, Krylov solver etc. Which has the fastest convergence rate among all these solver types from your experience?

thanks for your reply...
  Reply With Quote

Old   August 22, 2007, 16:38
Default Re: fast penta diagonal solver
  #2
Harish
Guest
 
Posts: n/a
I assume that you are talking about implicit algorithms.There are two ways to handle the penta diagonal solver.One is the direct inversion of the pentadiagonal system.The other possible way is to use the approximate factorization to split it into two tridiagonal systems which are much faster to invert than the pentadiagonal system.
  Reply With Quote

Old   August 22, 2007, 16:42
Default Re: fast penta diagonal solver
  #3
ag
Guest
 
Posts: n/a
Or you could use Gauss-Seidel or another iterative approach.
  Reply With Quote

Old   August 22, 2007, 18:24
Default Re: fast penta diagonal solver
  #4
otd
Guest
 
Posts: n/a
Also uncoupled and two space dimensions?
  Reply With Quote

Old   August 23, 2007, 01:23
Default Re: fast penta diagonal solver
  #5
sudhakar
Guest
 
Posts: n/a
hi otd, yes. I am dealing with uncoupled 2D flows only.. I know there are number of iterative schemes to solve penta diagonal matrix.. I need to know which of these has fastest convergence rate.. thanks
  Reply With Quote

Old   August 23, 2007, 09:11
Default Re: fast penta diagonal solver
  #6
Arpiruk
Guest
 
Posts: n/a
Do you mean, for a given precision e, which algorithm can solve the pentadiagonal fastest?

I think multigrid should be the fastest for this type of problem. ILU/GS will eventually stagnate when the low wave number component of residual is dominant. I've seen somewhere that multigrid+SIP(strongly implicit) is pretty good.

One question, is it positive definite? If so preconditioning CG should be very effective.
  Reply With Quote

Old   August 23, 2007, 09:18
Default Re: fast penta diagonal solver
  #7
Tom
Guest
 
Posts: n/a
Actually the best, and quickest, method of solving this is by considering the pentadiagonal system as a tridiagonal system whose entries are 2x2 matrices. The resulting block tridiagonal system can be solved using a variant of the standard tridiagonal solver.
  Reply With Quote

Old   August 23, 2007, 09:36
Default Re: fast penta diagonal solver
  #8
sudhakar
Guest
 
Posts: n/a
thank you friends.. I will look at CG and multi-gird methods.. If you get to know more details regarding this, kindly post a reply..
  Reply With Quote

Old   August 23, 2007, 13:20
Default Re: fast penta diagonal solver
  #9
Neazen
Guest
 
Posts: n/a
Before talking about improvements, you need to look for appropriate methods depending on the matrix properties. You can see in this link a very good introductory in this matter :

http://netlib2.cs.utk.edu/linalg/htm...Templates.html

When a method is working, no matter how long it takes, then you can look for preconditionners.
  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
thobois class engineTopoChangerMesh error Peter_600 OpenFOAM 4 August 2, 2014 10:52
smoothSolver diverges - solution in using PBiCG solver? makaveli_lcf OpenFOAM Running, Solving & CFD 3 September 11, 2013 13:44
Working directory via command line Luiz CFX 4 March 6, 2011 21:02
Using a user-defined solver in OF ozzythewise OpenFOAM Running, Solving & CFD 3 February 8, 2011 16:28
Getting too many iterations by velocity solving (aborting). Changing U - Solver? suitup OpenFOAM Running, Solving & CFD 0 January 20, 2010 08:45


All times are GMT -4. The time now is 04:42.