# Numerical problem in buoyantBoussinesqPimpleFoam

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

 July 3, 2012, 08:43 Numerical problem in buoyantBoussinesqPimpleFoam #1 New Member   Vitor Campos Join Date: Jun 2012 Posts: 1 Rep Power: 0 Sponsored Links Hello, I am trying to simulate heat transfer in a small vertical cilindrical pipe. The pipe is established with its axis in the y direction, which is the direction of the acceleration g. Air comes from the inlet which is the lower part of the pipe, flows through it and exits at the upper part. It looks similar to this: outlet | | | | walls | | | | inlet I am using the buoyantBoussinesqPimpleFoam with kOmegaSST turbulence model. The boundary conditions for time 0 are as following: U: inlet fixedValue (0 8 0); outlet zeroGradient; wall fixedValue (0 0 0) (no slip condition) T: inled fixedValue, 300; outlet zeroGradient; wall of the pipe fixedValue 450; p: calculated uniform 0 in all the boundaries; p_rgh: buoyantPressure uniform 0 in all the boundaries; alphat: inlet fixedValue 0; outlet zeroGradient; wall alphatWallFunction uniform 0 kappat: inlet calculated uniform 0; outlet zeroGradient; wall kappatJayatillekeWallFunction value=uniform 0 and Prt 0 0.85; nut: inlet calculated uniform 0; outlet zeroGradient; wall nutkWallFunction uniform 0; k: inlet fixedValue uniform 0.3267; outlet zeroGradient: wall kqRWallFunction uniform 0.3267 omega: inlet fixedValue uniform 3760; outlet zeroGradient: wall omegaWallFunction uniform 3760 The problem appears when I try to run the solver and it asks for a reference point and value for the field p_rgh. I have set the reference value as 0, since my conditions are null pressure at outlet. Therefore I decided to set my reference point at some cell on the outlet boundary. However, no matter which cell or reference point I choose, in the outlet or not, it comes to a numerical error with unphysical velocity and pressure values (for example U = 500 m/s) at that exact point. Can any of you help me finding my error on the simulation?

 July 4, 2012, 05:12 #2 Senior Member   Tom Fahner Join Date: Mar 2009 Location: Delft, Netherlands Posts: 370 Rep Power: 14 Dear Vitor, Try to put a fixedValue uniform 0 for p_rgh on the outlet instead of buoyantpressure there, I think the solver then also does not need the reference cell or reference value, which it should only need in a closed volume. Regards, Tom

 July 7, 2013, 07:55 #3 New Member   Korichi Abdelkader Join Date: Jan 2013 Posts: 11 Rep Power: 6 Dear Tom, I have a problem to simulate the natural convection in channel, I can't found the appropriate boundary conditions for p and p_rgh at the inlet and outlat channel. Thanks Aek.

 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 venkataramana OpenFOAM 3 December 1, 2013 08:30 Wouter Fluent UDF and Scheme Programming 6 June 6, 2012 04:43 fatb0y CFX 3 January 14, 2012 20:42 Kai_kc FLUENT 1 October 27, 2010 05:29 johny Main CFD Forum 7 September 5, 2005 05:53