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

poisson equation

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 15, 2006, 07:02
Default poisson equation
  #1
Beginner
Guest
 
Posts: n/a
Hi all,

I'm dealing with a problem since one month and finally I prefer to ask in order to have your point of view.

I have to solve a poisson equation, and I've tried many numerical methods : direct and iterative ones. Matrix are very large [800][800] is the minimum, and until now I still didn't find a solution. No one of the method I applied converged... Jacobi, Gauss-Seidel

So anyone have any suggestion about this.
  Reply With Quote

Old   October 15, 2006, 19:32
Default Re: poisson equation
  #2
zxaar
Guest
 
Posts: n/a
you might be making a mistake, which is difficult to tell from what you have posted, try reading how to discretise poisson equation from cfd-wiki (on cfd-online.com). I wrote small portion describing that. and yes, it converges well with gauss seidel etc.
  Reply With Quote

Old   October 15, 2006, 20:40
Default Re: poisson equation
  #3
Harish
Guest
 
Posts: n/a
Make sure that the applied boundary conditions are correct.If you use a tridiagonal solver with implicit scheme check for diagonal dominance.What kind of source terms and boundary conditions do you apply ?

A simple mathematical statement might help the cause

  Reply With Quote

Old   October 15, 2006, 20:56
Default Re: poisson equation
  #4
Renato.
Guest
 
Posts: n/a
First of all, make sure that your solver is correctly implemented and that you're using it correctly. It can be easily done by using a short and "controled" system (one that you have the solution) of equations (let's say 5x5 or 3x3). After that, if your method keep diverging you can go to check your Poisson discrete solution. As zxaar said, the wiki page can help you in this phase.

Cheers

Renato.

  Reply With Quote

Old   October 15, 2006, 23:24
Default Re: poisson equation
  #5
Ramp
Guest
 
Posts: n/a
To make sure, your code is working fine, Benchmark your code with simple problems like. cavity flow or backward step facing flow etc...
  Reply With Quote

Old   October 16, 2006, 06:57
Default Re: poisson equation
  #6
Beginner
Guest
 
Posts: n/a
Well, first of all, thanx for considering my post.

I'll try to answer all the question you've asked in order to make the situation clear.

I'm performing incompressible planar bidimensional and laminar Navier Stokes solver. I've derived consistant pressure equation (laplacian) and boundaries conditions (newmann).

I verified linear system Matrix that I reduced to [n][5] and vector construction.

The only problem is the method. I've read somewhere that sometimes a method don't work with in some specific cases.

I'm really discouraged, I've verified so many times all the algorithms.

Concerning Gauss Seidel and Jacobi, I've developped them before this project and I'm sure they are working. Actually I'm trying the SOR by lines and it seems to give logic results. (I didn't have any material providing pressure results to compare with) But its take long time to converge (8 min for a [2500][5] for a 10-4 eps)!!!

I have another question : do you think conjugate gradient method can be used in this case? It is really fast?
  Reply With Quote

Old   October 16, 2006, 16:06
Default Re: poisson equation
  #7
dontknow
Guest
 
Posts: n/a
Hi!

Gauss-Seidel and Jacobi are very old algorithms for iterative solution of linear equation systems. They have a rather bad convergence rate for big sparse matrices (and your equation system is big). Better algorithms are for example the SIP-algorithm or cgstab-family algorithms. If you want to use gauss-seidel or jacobi then only in combination with multigrid methods which are very powerful and have a fast convergence behaviour independent of size of matrices- they have optimal convergence behaviour. Google... try the book of ferziger and peric for a compact introduction of linear equation system solvers. you can also download many codes there... good luck!

  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
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 06:27
Poisson equation combined with transport equation Se-Hee CFX 0 December 27, 2007 02:00
pressure poisson equation on staggerred grid chandra Main CFD Forum 9 May 5, 2005 06:01
Poisson Equation in CFD Maciej Matyka Main CFD Forum 9 November 10, 2004 12:30
Poisson equation vs continuity equation DJ Main CFD Forum 1 August 5, 2004 20:01


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