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

linear system

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 1, 2002, 06:47
Default linear system
  #1
Annie
Guest
 
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)

Thanks
  Reply With Quote

Old   April 1, 2002, 14:50
Default Re: linear system
  #2
Junseok Kim
Guest
 
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, 19:45
Default Re: linear system
  #3
Sathe
Guest
 
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, 20:00
Default Re: linear system
  #4
ryoga
Guest
 
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, 18:16
Default Re: linear system
  #5
Tony
Guest
 
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

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
CFX11 + Fortran compiler ? Mohan CFX 20 March 30, 2011 18: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 01:23
Solving Linear System of Equations in UDF mohsen zendehbad FLUENT 5 May 8, 2010 10:39
Linear system solver elias FLUENT 0 November 19, 2008 10:30


All times are GMT -4. The time now is 09:56.