
[Sponsors] 
August 29, 2005, 04:42 
FISHPACK F90 package for NavierStokes solvers

#1 
Guest
Posts: n/a

Hi all, I am trying to implement FISHPACK F90 package to my existing code for the solution of pressure Poisson equation instead of the classical iterative technique.
Does anyone know any successfull implementation? Are there any critical issues which I should pay special attention? Any resource might be helpful (paper, book etc.) Thanks in advance, Ahmet 

August 29, 2005, 12:35 
Re: FISHPACK F90 package for NavierStokes solvers

#2 
Guest
Posts: n/a

Have you looked at the examples they have? Can you describe in more detail where you are stuck?


August 30, 2005, 04:18 
Re: FISHPACK F90 package for NavierStokes solvers

#3 
Guest
Posts: n/a

Hi agg,
By the way I checked my message archive and I realized that this is not the first time you're helping me As far as I remember you adviced me to look for the FISHPACK package. From that moment on I did quite a lot of things. I would like thank you for everything before I start The situation is as follows: I downloaded the FISHPACK F90 package and carefully went over the relevant examples and tutorials. As a starting point I attempted to solve a simple Poisson type equation on cartesian coordinates and validated the solution I obtained with FISHPACK against the classical iterative solution. Having seen that I obtained successful results I proceeded with NavierStokes solution. At this point a little info on the CFD code I am currently working on: The code is a transient code based Method of Lines (MOL) (an efficient semidiscrete technique for the solution of PDEs). In MOL the spatial derivatives are evaluated using high order approximations by employing readily available discretization packages such as 4th order lagrange interpolation polynomials. Then after spatial discretization, you obtain a set of ODEs in time and you integrate these ODEs using a powerful ODE solver such as LSODES. My overall scheme is in a way similar to PISO (pressure implicit with splitting operators), I predict velocities then correct it once, then predcit temperatures, correct temperatures etc.. I implemented the FISHPACK package to my CFD code. I first tested the code on a isothermal problem and obtained excellent results. While doing this, I computed the Poisson source term necessary for FISHPACK using the same discretization packages mentioned above. Then I tested my code on a nonisothermal (walls of the pipe is heated) turbulent (Re=4200) test case. I obtained meaningful results. However, I realized the followings: *The magnitude of the pressures computed by FISHPACK are hihgly sensitive to grid resolution. *At the first 4 time steps I observe negative pressure values (which is not reasonable)then after the 5th time step reasonable pressures are obtained then code starts to produce somewhat nice results till the steady state. *Although the results are reasonable, they are not as good as I want them to be. (10% relative error) I do know that success of a CFD code depends on many factor, so many that we cannot list them all here. The reason I posted this message was, to learn if someone else who had previously worked with FISHPACK encountered similar problems. More specifically; *Can you comment on my discretization scheme? Is it okey to use such a package for discretization while computing the Poisson source term? *While choosing the grid resolutions, should I establish it in such a way that delta r= delta z? You told me describe in detail and perhaps this message was too much in detail Sorry if I took a lot from your precious time and lots of thanks in advance.. Ahmet PS. You are always a helping hand After your previous message the whole course of my study has changed dramatically. Thanks one more time 

September 1, 2005, 03:52 
Re: FISHPACK F90 package for NavierStokes solvers

#4 
Guest
Posts: n/a

Hi agg,
I just recieved an email today that you have posted a message in response to my previous, rather long, email. However when I click the link I see that the link is dead. If possible can you post your reply again? Thanks in advance, Ahmet 

September 2, 2005, 11:10 
Re: FISHPACK F90 package for NavierStokes solvers

#5 
Guest
Posts: n/a

Hi,
I think my message was one of those accidentally deleted by the admin. Here are my thoughts. 1) Since you were able to obtain good solutions for the nonisothermal case, the problem COULD be somewhere else and not the Poisson solver (bouyancy term in momentum equations, just a guess) 2) Fishpack is a second order finite difference solver. Ideally, you should not use the 4th order scheme to compute the source term of your Poisson equation. Remember, d2P is divergence of the pressure gradient. Hence, the divergence and pressure gradient operator should be consistent. You could still use the 4th order Lagrangian scheme for other variables. I had problems with stability for my finite difference solutions and so I moved pressure to center of the cell. Your solution however, seems stable, but the inconsistency in the differentiation could be responsible for the inaccuracies that you are noticing. Here are some references that will be helpful before you decide to make any changes to your code 1) "Numerical heat transfer and fluid flow", Patankar (chapter 6) 2)"Numerical study of the turbulent flow past an airfoil with trailing edge separation", Rhie and Chow, AIAA Journal, Vol. 21, No. 11, p. 15251532 (1983) 3)"Numerical solutions for the pressure Poisson equation with Neumann boundary conditions using a nonstaggered grid, I", Abdallah, JCP (70), p. 182192 (1987) 4)"The discrete continuity equation in primitive variable solutions of incompressible flow", Sotiropoulos and Abdallah, JCP (95), p. 212227 (1991) 5)"Alternate formulations for the pressure equation Laplacian on a collocated grid for solving the unsteady incompressible NavierStokes equations", Danesh Tafti, JCP(116), p. 143153 (1995) Good Luck! 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Extend Project 1.6ext Binary Release for Mac OS X  hjasak  OpenFOAM Announcements from Other Sources  26  November 5, 2013 17:50 
OpenFoam 1.6ext  RPM build errors  preibie  OpenFOAM  12  September 8, 2011 03:12 
OpenFOAM 1.5.x package  CentOS 5.3 x86_64  linnemann  OpenFOAM Installation  7  July 30, 2009 03:14 
Errors running allwmake in OpenFOAM141dev with WM_COMPILE_OPTION%3ddebug  unoder  OpenFOAM Installation  11  January 30, 2008 21:30 
Navier Stokes solvers for growing domains  Naga Kiran Potluri  FLUENT  0  August 11, 2005 10:24 