# 3D Poisson's equation

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 June 8, 2005, 14:22 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

 June 8, 2005, 16:22 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.

 June 8, 2005, 16:45 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?

 June 9, 2005, 00:10 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

 June 9, 2005, 00:32 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/

 June 9, 2005, 04:38 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.

 June 9, 2005, 04:56 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 ?

 June 9, 2005, 06:20 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 ?

 June 9, 2005, 06:47 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/

 June 9, 2005, 16:02 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.

 June 9, 2005, 16:06 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).

 June 9, 2005, 16:21 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.

 June 10, 2005, 13:25 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.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Mihail CFX 7 September 7, 2014 06:27 Sas CFX 15 July 13, 2010 08:56 saii CFX 2 September 18, 2009 08:07 Anil CFX 8 September 25, 2007 12:57 dean schrage Main CFD Forum 16 April 11, 2000 13:46

All times are GMT -4. The time now is 13:55.

 Contact Us - CFD Online - Privacy Statement - Top