# Pressure Correction Method for Cylindrical Coordinates

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

 July 2, 2009, 08:33 Pressure Correction Method for Cylindrical Coordinates #1 New Member   Eslam Reda Join Date: Jun 2009 Posts: 16 Rep Power: 9 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, 14:23 #2 Member   Join Date: Mar 2009 Posts: 62 Rep Power: 9 Hello, What kind of help do you need ? Do you have a specific question ? What solver do you use ? RenardP

 July 2, 2009, 23:39 Pressure Correction Method for Cylindrical Coordinates #3 New Member   Eslam Reda Join Date: Jun 2009 Posts: 16 Rep Power: 9 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.

 July 3, 2009, 07:50 #4 Member   Join Date: Mar 2009 Posts: 62 Rep Power: 9 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 New Member   Eslam Reda Join Date: Jun 2009 Posts: 16 Rep Power: 9 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.

 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 Wijaya FLUENT 15 May 18, 2016 10:08 richard_larson OpenFOAM Running, Solving & CFD 8 October 30, 2008 09:48 Perumal Main CFD Forum 1 September 23, 2008 14:43 Seeker01 Main CFD Forum 2 January 13, 2003 03:49 Dan Moskal Main CFD Forum 0 October 24, 2002 22:02

All times are GMT -4. The time now is 11:09.