Pressure Correction Method for Cylindrical Coordinates

Eslam Reda · July 2, 2009, 08:33

Hello every body, I am facing a trouble with the pressure correction method in a simple flow in pipe problem. When I use very low values for density (1 kg/m3), the code converges. While when using heavy densities (1000 kg/m3) the program diverges. Please I need your help in this annoying trouble.

RenardP · July 2, 2009, 14:23

Hello,

What kind of help do you need ? Do you have a specific question ? What solver do you use ?

RenardP

Eslam Reda · July 2, 2009, 23:39

Hello, let's start:

I have a simple laminar water flow in a pipe problem. The inlet and outlet pressures are known (0.035 Pa and 0 Pa respectively), the inlet radial velocity is set to zero while the outlet one is floating and also the inlet and outlet axial velocities are floating.
I use an ordinary SIMPLE algorithm, Fortran program. I tried it as explicit and implicit giving the same problem.
The trouble is the divergence of the reults.
The strange thing is when I set the density to unity, rather than 1000 kg/m3 for water, the program converges and gives the correct results.

I need your help in this problem as it is considered as my start to the more difficult practical problems.

Thanks so much.
Eslam Reda.

RenardP · July 3, 2009, 07:50

OK,

the SIMPLE algorithm is a classical one so I think you have implemented a correct version.

In your place I would check the grid and the implementation of the boundary conditions. Usually a divergent solution is caused by bad implementation of the boundary conditions.

Also if you can find an analytical solution for a simplified version of your tube problem you can test various parts of your implementation (e.g a Poiseuille flow in a duct).

Paul

Eslam Reda · July 4, 2009, 02:15

Hello again,
Surely I checked the code for a simple Poiseuielle flow problem. The program diverges at high densities but gives the expected analytical results for low densities. I wish to get the data, Dx, Dy, under-relaxation factor and Dt, for a solved water flow in a pipe problem.

July 2, 2009, 08:33	Pressure Correction Method for Cylindrical Coordinates	#1
Eslam Reda New Member Eslam Reda Join Date: Jun 2009 Posts: 12 Rep Power: 3	Hello every body, I am facing a trouble with the pressure correction method in a simple flow in pipe problem. When I use very low values for density (1 kg/m3), the code converges. While when using heavy densities (1000 kg/m3) the program diverges. Please I need your help in this annoying trouble.

July 2, 2009, 23:39	Pressure Correction Method for Cylindrical Coordinates	#3
Eslam Reda New Member Eslam Reda Join Date: Jun 2009 Posts: 12 Rep Power: 3	Hello, let's start: I have a simple laminar water flow in a pipe problem. The inlet and outlet pressures are known (0.035 Pa and 0 Pa respectively), the inlet radial velocity is set to zero while the outlet one is floating and also the inlet and outlet axial velocities are floating. I use an ordinary SIMPLE algorithm, Fortran program. I tried it as explicit and implicit giving the same problem. The trouble is the divergence of the reults. The strange thing is when I set the density to unity, rather than 1000 kg/m3 for water, the program converges and gives the correct results. I need your help in this problem as it is considered as my start to the more difficult practical problems. Thanks so much. Eslam Reda.

Thread Tools
Show Printable Version Email this Page
Display Modes
Linear Mode Switch to Hybrid Mode Switch to Threaded Mode

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Pressure correction problem	richard_larson	OpenFOAM Running / Solving / CFD	8	October 30, 2008 09:48
Streamfunction-vorticity method - Pressure find	Perumal	Main CFD Forum	1	September 23, 2008 14:43
"Pressure Inlet" Boundary Setup	Wijaya	FLUENT	11	December 23, 2006 06:21
Pressure Correction method, Finite Volume Method	Seeker01	Main CFD Forum	2	January 13, 2003 02:49
Gas pressure question	Dan Moskal	Main CFD Forum	0	October 24, 2002 22:02

July 2, 2009, 14:23		#2
RenardP Member Join Date: Mar 2009 Posts: 59 Rep Power: 3	Hello, What kind of help do you need ? Do you have a specific question ? What solver do you use ? RenardP

July 3, 2009, 07:50		#4
RenardP Member Join Date: Mar 2009 Posts: 59 Rep Power: 3	OK, the SIMPLE algorithm is a classical one so I think you have implemented a correct version. In your place I would check the grid and the implementation of the boundary conditions. Usually a divergent solution is caused by bad implementation of the boundary conditions. Also if you can find an analytical solution for a simplified version of your tube problem you can test various parts of your implementation (e.g a Poiseuille flow in a duct). Paul

July 4, 2009, 02:15		#5
Eslam Reda New Member Eslam Reda Join Date: Jun 2009 Posts: 12 Rep Power: 3	Hello again, Surely I checked the code for a simple Poiseuielle flow problem. The program diverges at high densities but gives the expected analytical results for low densities. I wish to get the data, Dx, Dy, under-relaxation factor and Dt, for a solved water flow in a pipe problem.