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

Quadradiagonal solver

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 8, 2006, 08:53
Default Quadradiagonal solver
  #1
Frank
Guest
 
Posts: n/a
Hello Folks, I am writing an implicit solver that solves the transonic flow over an airfoil. The scheme I'm trying to use involves the solution of a quadradiagonal set of equations.

I'm used to solving tridiagonal sets of equations for a similar application and use the standard thomas algorithm. I've heard there's a similar method for the quad case but I can't seem to find the particular algorithm anywhere.

I realise that I could come up with something from first principles but feel that an efficient, tried and tested method would be better than anything I'd cobble together myself.

If someone knows of a good algorithm to use could you please point me in the right direction?

Many thanks, Frank
  Reply With Quote

Old   May 8, 2006, 11:22
Default Re: Quadradiagonal solver
  #2
ganesh
Guest
 
Posts: n/a
Dear Frank,

You could always make use of relaxation procedures/iterative solvers to get the solution in an efficient manner, rather than inverting a sparse matrix.

Regards,

Ganesh
  Reply With Quote

Old   May 8, 2006, 12:38
Default Re: Quadradiagonal solver
  #3
Frank
Guest
 
Posts: n/a
Hi Ganesh, I'm already using an iterative proceudre. The quadridiagonal equations result from solving one vertical column (in 2D cartesian mesh) for an iterative SLOR algorithm. If I nested one iterative procedure within another I guess this would take a significant amount of computer time.

To be honest I'm not even too sure how I would go about solving a quadradiagonal set of equations with an iterative solver. I'm relatively new to the whole CFD thing and until now have relied upon mainly black box solvers. My previous work required the solution of tridiagonal sets of equations within an SLOR scheme, for which I used the Thomas algorithm, and I've been told that quad-equations can be solved using an extension of this technique. This would make my life much easier as I'm familiar with the technique and it'd mean the structure of my new code would be quite similar to previous ones. Cheers, Frank
  Reply With Quote

Old   May 8, 2006, 14:42
Default Re: Quadradiagonal solver
  #4
ganesh
Guest
 
Posts: n/a
Dear Frank,

Please take a look at the book by Pozrikidis Pozrikidis (1998), pp. 144-146. I think the title is Theorretical and Computational Fluid Dynamics, where it tals about a Modified Thomas algorithm for sparse banded matrices, tri/quad/penta. Also note that Thomas algorithm is a special form of Gaussian elimination, and if you are not really worried about the computational speed and efficiency any iterative solver(like Jacobi, Gauss-Siedel) etc.. in conjunction with LU decomposition can do the job. Do take a look at some of the literature on implicit schemes in CFD, for eg. the point implicit method employed routinely by Lohner et.al. among others.

Hope this helps

Regards,

Ganesh

  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
thobois class engineTopoChangerMesh error Peter_600 OpenFOAM 4 August 2, 2014 09:52
A New Solver for Supersonic Combustion nakul OpenFOAM 6 March 30, 2014 03:00
A New Solver for Supersonic Combustion nakul OpenFOAM Announcements from Other Sources 18 February 19, 2013 08:48
Working directory via command line Luiz CFX 4 March 6, 2011 21:02
why the solver reject it? Anyone with experience? bearcat CFX 6 April 28, 2008 14:08


All times are GMT -4. The time now is 15:32.