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

Pressure Poisson equation

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 1, 2016, 17:23
Default Pressure Poisson equation
  #1
New Member
 
Join Date: Apr 2014
Posts: 1
Rep Power: 0
Zach_uky is on a distinguished road
Hi, I know there are a lot of thread about pressure poisson equation, I ve read all of them, but I haven't found a solution for my problem. Next let me tell you the details of my problem.

First, I have compressible flow solver, after several years of work, it can solve turbulent problem well. But my adviser think it is too slow to solve an extra energy equation, want me write a incompressible solver based on this frame.

I works several weeks now. This code is a finite volume code. all of the value is in the cell center. I choose projection method for the incompressible solver. First solve NS equation without pressure to get the velocity u* which does not satify divergence free filed u. Then solve the pressure poisson equation and project the velocity to divergence free field. The pressure poisson equation is
Δp=(▽• u*)/t
which could be discretized as
∑A*n• (lr)(p_l-p_r)/(lr)^2=V*(▽• u*)/t
where p_l and p_r is the left and right side pressure, t is time step length, V is volume, A is the face area, n is the normal of the face, (lr) is the vector from left cell center to right cell center. we also can set
AoverH= A*n• (lr)(p_l-p_r)/(lr)^2

Then the projection part is
u=u*-t▽p

Then go to next time step. When I solve the pressure poisson equation, I use
divergence in the cell center as the right hand side(rhs) of my solver, then put the jacobian AoverH into the martix(C).
Then the problem is
C * P=rhs
For the boundary condition, I use neumann condition for inlet, wall, symmetry, but 0 pressure for outlet.
My solver is a domain splitting MPI code. I use GMRES in petsc to solver this problem, the grid I use is a 100000 grid mesh, it is a quarter pipe.

my problem 1: for the neumann boundary face, whether AoverH should be zero? When I set AoverH to zero at the boundary, the iteration number is huge about 1000, the result is right. When AoverH is not zero, the iteration is no more than 20. after converge, the result is wrong.

Problem 2: if AoverH on the neumann boundary should be zero, how to decrease the iteration number.

problem 3: I have not use stagger grid, p and u are in the cell center, is that ok?

Can any one give me some guidance? Thank you in advance
Zach_uky is offline   Reply With Quote

Reply

Tags
cfd, incompressible flow, ppe
Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
question regarding LES of pipe flow - pimpleFoam Dan1788 OpenFOAM Running, Solving & CFD 37 December 26, 2017 14:42
writing an acoustic code based on FW-H equation, have a question regarding pressure mrshb4 OpenFOAM Running, Solving & CFD 3 November 8, 2014 20:12
About theory:how to get the pressure equation when solving u-p simultaneously? lzw2003 Main CFD Forum 8 April 8, 2012 22:13
error message cuteapathy CFX 14 March 20, 2012 06:45
FEM pressure poisson equation (Implicit) cfd-beginner Main CFD Forum 0 August 9, 2005 13:32


All times are GMT -4. The time now is 10:57.