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

Sparse linear systems in finite volume method

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 12, 2007, 16:34
Default Sparse linear systems in finite volume method
  #1
Lionel S.
Guest
 
Posts: n/a
Hi all,

It seems that finite volume methods makes at each iteration a very big sparse linear system

Ax = b

Does someone know precisely how it is solved ? I suppose that iterative methods as conjugate gradient can be used... What do you think about it ?

It would be very easy if it was possible to make small non sparse matrix from big sparse matrix, as it can be done in finite elements methods (in structure) but I suppose that one cannot do that in finite volumes.

Kindest regards

Lionel
  Reply With Quote

Old   January 12, 2007, 19:17
Default Re: Sparse linear systems in finite volume method
  #2
nico
Guest
 
Posts: n/a
yep, conjugate gradients method are used. The whole matrix is of course not stored. On structured grids, a few diagonals (depending of interpolation scheme) are stored. On unstructured grids, defered correction (upwind interpolation + iterations to get to higher order) is usually used, so only the main diagonal and the contribution of the face-neighbour is stored. Hope this helps.

Nico
  Reply With Quote

Old   January 13, 2007, 06:42
Default Re: Sparse linear systems in finite volume method
  #3
Lionel S.
Guest
 
Posts: n/a
Ok thanks... If ever you know other methods to converge faster than conjugate gradients, please make me know Bye
  Reply With Quote

Old   January 13, 2007, 15:51
Default Re: Sparse linear systems in finite volume method
  #4
raju
Guest
 
Posts: n/a
Sparse linear systems obtained from finite volume discretizations can also be solved using multifrontal solvers. The concept of multifrontal solvers is developed from the frontal methods used in finite element. They can solve the linear system very efficiently. The only disadvantage is its huge memory requirement. For 2D problems, these solvers can be effectively used, but for 3D, memory requirement would become a limitation.

Mulitfrontal solvers can be downloaded free on the internet. Examples of multifrontal solvers are (1) UMFPACK v5.0 (2) MUMPS solver (3) Pardiso etc.
  Reply With Quote

Old   January 14, 2007, 18:48
Default Re: Sparse linear systems in finite volume method
  #5
Lionel S.
Guest
 
Posts: n/a
Thank you for answering...

In fact, I have problems with conjugate gradients because the boundary conditions makes matrix not symmetric, and it doesn't converge.

I can converge if I try to solve (A' * A)x = A' * b

(where A' is A transposed) but the solution is not accruate enough. I could also use a preconditioner such as incomplete LU but inverting these preconditioners is expensive in time and memory.

Do you have an idea ?

  Reply With Quote

Old   January 17, 2007, 10:43
Default Re: Sparse linear systems in finite volume method
  #6
Ben
Guest
 
Posts: n/a
Hi,

Can anyone tell me how does Krylov Space Methods such as GMRES compare to multifrontal solvers?

I'm currently using PETSc which's a Krylov Space package linear solver to solve my 2D FVM NS eqn. The Ns eqn is solved using the fractional step.

I've no experience with multifrontal solvers.

Thank you
  Reply With Quote

Old   January 19, 2007, 11:04
Default Re: Sparse linear systems in finite volume method
  #7
raju
Guest
 
Posts: n/a
I have no experience with fractional step methods. I have worked with steadty state NS implementation using Multifrontal solvers. It is very robust and efficient. When compared with segregated solvers like SIMPLE,SIMPLER etc, multifrontal solvers can be an order of magnitude faster.
  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
Finite Volume Method cfd seeker Main CFD Forum 3 September 8, 2011 04:36
2D Finite Volume Method - Please help WatchFirefly Main CFD Forum 12 September 12, 2010 20:58
Chorin's Projection Method for Finite Volume Scott2 Main CFD Forum 1 August 16, 2010 20:24
Axisymmetrical mesh Rasmus Gjesing (Gjesing) OpenFOAM Native Meshers: blockMesh 10 April 2, 2007 14:00
ALE in finite volume method littlelz Main CFD Forum 5 June 21, 2003 12:50


All times are GMT -4. The time now is 08:05.