CFD Online URL
[Sponsors]
Home > Forums > Main CFD Forum

Pressure Correction Method for Cylindrical Coordinates

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 2, 2009, 09:33
Default Pressure Correction Method for Cylindrical Coordinates
  #1
New Member
 
Eslam Reda
Join Date: Jun 2009
Posts: 12
Rep Power: 7
Eslam Reda is on a distinguished road
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.
Eslam Reda is offline   Reply With Quote

Old   July 2, 2009, 15:23
Default
  #2
Member
 
Join Date: Mar 2009
Posts: 62
Rep Power: 7
RenardP is on a distinguished road
Hello,

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

RenardP
RenardP is offline   Reply With Quote

Old   July 3, 2009, 00:39
Default Pressure Correction Method for Cylindrical Coordinates
  #3
New Member
 
Eslam Reda
Join Date: Jun 2009
Posts: 12
Rep Power: 7
Eslam Reda is on a distinguished road
Hello, let's start:
  1. 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.
  2. I use an ordinary SIMPLE algorithm, Fortran program. I tried it as explicit and implicit giving the same problem.
  3. The trouble is the divergence of the reults.
  4. 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.
Eslam Reda is offline   Reply With Quote

Old   July 3, 2009, 08:50
Default
  #4
Member
 
Join Date: Mar 2009
Posts: 62
Rep Power: 7
RenardP is on a distinguished road
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
RenardP is offline   Reply With Quote

Old   July 4, 2009, 03:15
Default
  #5
New Member
 
Eslam Reda
Join Date: Jun 2009
Posts: 12
Rep Power: 7
Eslam Reda is on a distinguished road
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.
Eslam Reda is offline   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
"Pressure Inlet" Boundary Setup Wijaya FLUENT 13 April 11, 2013 09:50
Pressure correction problem richard_larson OpenFOAM Running, Solving & CFD 8 October 30, 2008 10:48
Streamfunction-vorticity method - Pressure find Perumal Main CFD Forum 1 September 23, 2008 15:43
Pressure Correction method, Finite Volume Method Seeker01 Main CFD Forum 2 January 13, 2003 03:49
Gas pressure question Dan Moskal Main CFD Forum 0 October 24, 2002 23:02


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