# Pressure varies with cell size!!

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

 December 13, 2009, 08:35 Pressure varies with cell size!! #1 New Member   Join Date: Mar 2009 Posts: 27 Rep Power: 8 Hi, I was using icoFoam and interFoam and noticed that the pressure range (icoFoam) or pressure value (interFoam) varies with cell size and does not seem to converge with refinement. The smaller the cell size, the larger the pressure. If you use the cavity icoFoam case and refine the mesh you see that pressure range continuously increases. Does somebody know what is wrong or experienced this problem?

 December 14, 2009, 05:42 #2 New Member   Join Date: Mar 2009 Posts: 27 Rep Power: 8 Is this a problem due to Rhie-Chow interpolation? Last edited by jugghead; December 15, 2009 at 11:07.

 December 15, 2009, 11:08 #3 New Member   Join Date: Mar 2009 Posts: 27 Rep Power: 8 Anyone?

 December 15, 2009, 21:52 #4 Senior Member   Kevin Smith Join Date: Mar 2009 Posts: 103 Rep Power: 8 I think this is a problem with the boundary conditions on pressure - take a look at the boundary condition on 'p'. May be an issue with the cells at the intersection of the moving wall and the side walls - you can see the pressure peaks in those corners. I've had difficulties with pressure boundary conditions before in openfoam, they can be tricky at times.

 December 16, 2009, 05:22 #5 New Member   Join Date: Mar 2009 Posts: 27 Rep Power: 8 As I said in 1st post, this also happens with interFoam. Even if the pressure is normalized or isn't absolute pressure whatever it should not depend on cell size. It seems to vary linearly with the inverse of cell distance. For example if I have a 1D mesh and have 20 cells the pressure is p, if I have 40 cells the pressure is 2*p and so on. To repliclate this problem just use interFoam with 1 inlet and 1 outlet a simple plate for example 1m x 0.4m and place a alpha=1 at inlet. Probe the pressure near inlet and refine the mesh. Put a velocity of 0.1 m/s at inlet. Start with mesh as 25x10 cells and then put 50x20 and then 75x30, etc. and see what happens to pressure.

 December 16, 2009, 05:29 #6 Senior Member   Anonymous Join Date: Mar 2009 Posts: 110 Rep Power: 8 Jugghead, can you try 100x100 and then 300x300 and report what pressure is doing at these stages compared to 75x30?

 December 16, 2009, 05:51 #7 New Member   Join Date: Mar 2009 Posts: 27 Rep Power: 8 Attention that this is not turbulent flow. So it should converge pretty easily without resorting to high density meshes. Low Reynolds. You can use a very high viscosity and this happens.

 March 26, 2015, 09:02 #8 New Member   Roman Voronov Join Date: May 2010 Location: Russia Posts: 12 Rep Power: 7 Hello! I try to simulate processes in multisection cyclone unit with OpenFOAM. Previously I did an experiment and experimental pressure values are in 2-3 times lower than simulated values. Velocity values are similar to experimental. Tell me please, is there another way to obtain normal values, consistent with experimental values, except splitting the mesh? May be it's possible to set some solver's parameters or even change settings in solver's source code? I tried to simulate my case with this solvers/settings: simpleFoam & pimpleFoam solvers, limitedLinear, linearUpwind & QUICK schemes for div(phi,U), GAMG linear solver for pressure and smoothSolver for other fields with under-relaxation factors 0.3 for pressure and 0.7 for other fields. With all of that, velocity values are consistent with experimental, but pressure values are very big. Thank you!

 Tags cell, pressure, size

 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 shrina OpenFOAM Running, Solving & CFD 10 October 3, 2013 14:38 Giuki FLUENT 1 July 19, 2011 11:35 paka OpenFOAM 15 July 21, 2008 05:03 Antech Main CFD Forum 0 April 25, 2006 02:15 Art Stretton Phoenics 5 April 2, 2002 05:59

All times are GMT -4. The time now is 04:06.