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

3D Poisson's equation

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 8, 2005, 14:22
Default 3D Poisson's equation
  #1
Davoche
Guest
 
Posts: n/a
Hi everybody, I'm searching for the best (rapid) method to resolve a 3D Poisson's equation (in space) with stretched grid for 2 directions (the last one is uniform). I precise that's it should be applicable to a finite difference resolution. For the moment, I'm thinking about the Multigrid Method or Non Uniform FFT... Thanks
  Reply With Quote

Old   June 8, 2005, 16:22
Default Re: 3D Poisson's equation
  #2
Hubert Janocha
Guest
 
Posts: n/a
I am a fan of multigrid. There are some more interesting points of it, e.g. whole SIMPLE as a smoother on ALL grid levels (idea by Milovan Peric), and the strategy of the "Nested Iteration". I am sure multigrid is more efficient than FFT. In "Numerical Recipes" you can find the linear and the FAS one as subroutines for finite difference discr. , but the FAS isn't really needed. I can give you some papers for multigrid if you want it, but only in german language.
  Reply With Quote

Old   June 8, 2005, 16:45
Default Re: 3D Poisson's equation
  #3
Mani
Guest
 
Posts: n/a
Is there such a thing as a nonuniform FAST Fourier transform? I would like to know more about that. Any good references?
  Reply With Quote

Old   June 9, 2005, 00:10
Default Re: 3D Poisson's equation
  #4
shekharc
Guest
 
Posts: n/a
Could any of you please suggest me a reference from where I can get to know about implementation of FFT to solve a PDE like Poission's equation. I am also facing the same problem of slow convergence rate and dont want to implement the multigrid method because it's a little complex

  Reply With Quote

Old   June 9, 2005, 00:32
Default Re: 3D Poisson's equation
  #5
agg
Guest
 
Posts: n/a
Here's a package that might be useful (http://www.fftw.org). They have both serial and parallel versions.

For finite difference here are 2 packages i know of http://www.scd.ucar.edu/css/software/fishpack/ http://www.scd.ucar.edu/css/software/mudpack/
  Reply With Quote

Old   June 9, 2005, 04:38
Default Re: 3D Poisson's equation
  #6
Davoche
Guest
 
Posts: n/a
Thanks everybody

Hubert> Thanks for your help but I dont speak a word of german.

Mani> For Non Uniform FFT, I suggest you "Fast Fourier transforms for nonequispaced data:A tutorial" D.Potts, G.Steidl, M.Tasche. You can also find some C routine of the method on the Web (don't remember where).

shekharc> Standard FFT => Numerical recipes

Agg> Thank you. I will look this soon.
  Reply With Quote

Old   June 9, 2005, 04:56
Default Re: 3D Poisson's equation
  #7
Davoche
Guest
 
Posts: n/a
Hubert> For the Poisson's equation, I don't need the FAS for nonlinear elliptic operator because Laplacian is linear. Isn't it ?

  Reply With Quote

Old   June 9, 2005, 06:20
Default Re: 3D Poisson's equation
  #8
Davoche
Guest
 
Posts: n/a
Sorry, another question: Is Numerical Recipes method for Multigrid easily converted to 3D and non uniform grid in 2 ( on 3) directions ?
  Reply With Quote

Old   June 9, 2005, 06:47
Default Re: 3D Poisson's equation
  #9
andy
Guest
 
Posts: n/a
Fourier and cyclic reduction approaches are usually 'best' whenever they are applicable.

A multigrid approach is usually the 'best' general approach. Here is a web site full of multigrid algorithms and code:

http://www.mgnet.org/

  Reply With Quote

Old   June 9, 2005, 16:02
Default Re: 3D Poisson's equation
  #10
Hubert Janocha
Guest
 
Posts: n/a
Yes that's right. In certan circumstances you nead FAS if you want to solve a nonlinearity in the iteration step, e.g. heat conduction as a function of temperature itself. or in compressible flows (especially high speed gas flows). The density as a function of temperature, this introduces strong nonlinearities into the solution, causing divergence. The root of this problem is solving the total enthalpy or rothalpy (in general codes), which causes errors in velocities and turbulent kinetic energy, which leads to erros in temperatur.
  Reply With Quote

Old   June 9, 2005, 16:06
Default Re: 3D Poisson's equation
  #11
Hubert Janocha
Guest
 
Posts: n/a
Yes it is very easy to converted it into 3d, because it's just only the quit normal Poison equation without SIMPLE or turbulence or Navier-Stokes are any thing else we found in CFD normaly. It is only a short subroutine (for both linear and FAS).
  Reply With Quote

Old   June 9, 2005, 16:21
Default Re: 3D Poisson's equation
  #12
Hubert Janocha
Guest
 
Posts: n/a
There is at least one very important thing to keep in mind: I don't really know much about FFT and its behauviour for convergence, but every algorithm i know for solving the diagonal matrix has problems with very large meshes. First of all the simple Gaus Seidel relaxation. So when you want to do aero dynamics you need very large grids, and at a special border of mesh size (depended to the solver algorithm) the cells react only of the changing of there neighbour cells, and don't progress in convergence. Never the less you will see a flat line for the residuals. That's a great problem especially with large tet meshes, because the neihgbour coupling consists only of 4 cells instead of 6 in hex meshes. Multigrid converges linear, means the mesh size is (almost) independed to the convergence behaviour. I think in 10 years we will have so powerfull computers that we want bigger and bigger meshes, and we all will use multigrid as a must.
  Reply With Quote

Old   June 10, 2005, 13:25
Default Re: 3D Poisson's equation
  #13
Mani
Guest
 
Posts: n/a
shekharc, what you are looking for is called "spectral method" (more general but including FFT). There is a comprehensive Annual Review in Fluid Mechanics by Hussaini and Zang (1987)... old but still valid.

The principle is simple though, it's basically the same thing you did in your undergraduate analysis class: use the Fourier transform to solve a linear PDE analytically. Only in this case you do it numerically for nonlinear problems.
  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
Constant velocity of the material Sas CFX 15 July 13, 2010 08:56
mass flow in is not equal to mass flow out saii CFX 2 September 18, 2009 08:07
Solve Poissons Equation Anil CFX 8 September 25, 2007 12:57
pressure Poissons equation dean schrage Main CFD Forum 16 April 11, 2000 13:46


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