CFD Online Logo CFD Online URL
Home > Forums > General Forums > Main CFD Forum

linear system

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

LinkBack Thread Tools Search this Thread Display Modes
Old   April 1, 2002, 07:47
Default linear system
Posts: n/a
Dear all,

When solving 2d or 3d pressure Poisson equation by iterative methods, we usually have to compose a big matrix ((nx*ny,nx*ny) for 2D). For large nx or ny, say 200, the resultant matrix will be very large. Is there any means to cope this problem?

(I ever heard that since the matrix is usually sparse, there may be some tricks in storing the matrix)

  Reply With Quote

Old   April 1, 2002, 15:50
Default Re: linear system
Junseok Kim
Posts: n/a
You can use pointwise iteration method with multigrid method. Then you don't have to solve the big matrix.

  Reply With Quote

Old   April 1, 2002, 20:45
Default Re: linear system
Posts: n/a
You may use matrix free methods. A good reference for that is "A Globally Convergent Matrix-Free Algorithm For Implicit Time-Marching Schemes Arising In Finite Element Analysis In Fluids" Johan, Hughes, Shakib - Computer Methods in Applied Mechanics and Engineering 87 1991 (281-304)
  Reply With Quote

Old   April 1, 2002, 21:00
Default Re: linear system
Posts: n/a
I don't know if this will answer your question or not, but Argonne National Laboratory has a "suite" of computational tools that you may find useful. It is called PETSc (Portable, Extensible Toolkit for Scientific Computation). It has a nice library of iterative solvers -- including Krylov Subspace Methods -- and preconditioner routines. Furthermore, it has the capability to use Matrix-Free solvers, although that is more of an advanced feature and it may require some programming on your part. In my experience, it works quite well and should easily handle systems as large as 200 x 200. Hope this helps!
  Reply With Quote

Old   April 2, 2002, 19:16
Default Re: linear system
Posts: n/a
Hi there,

The nodal support is very narrow and the matrix must be symmetric. You can take the advantage of the sparsity when assembly the matrix. It is not a problem for nowadays computers to restore 200x200 or bigger matrices.
  Reply With Quote


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
CFX11 + Fortran compiler ? Mohan CFX 20 March 30, 2011 19:56
Singular Linear System - Pressure Poisson Equation Allan Main CFD Forum 11 March 17, 2011 10:39
solution diverges when linear upwind interpolation scheme is used subash OpenFOAM 0 May 29, 2010 02:23
Solving Linear System of Equations in UDF mohsen zendehbad FLUENT 5 May 8, 2010 11:39
Linear system solver elias FLUENT 0 November 19, 2008 10:30

All times are GMT -4. The time now is 17:23.