Boundary Conditions : Total Pressure or Velocity

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 February 12, 2010, 17:36 Boundary Conditions : Total Pressure or Velocity #1 Member   Jérémy Bulle Join Date: Nov 2009 Posts: 93 Rep Power: 7 Hello, I have to simulate a curved pipe. Here are the caracteristics of the flow : Incompressible Re = 40 000 Hydraulic diameter = 0.08m So we found U characteristic = 15.7m/s k = 0.3676 epsilon = 0.3328 (with a turbulent intensity of 6, don't remember the characteristical length used to calculate this) Since I begin in CFD, I've asked help to a professor wich told me that I should impose these BC : Outlet = Ambient pressure Inlet = Total pressure To calculate the total pressure, we used empirical relation for pressure loss and we find pstatic by adding the pressure loss to the outlet pressure and we just have to add rho u² to have total pressure (u characteristic is used for u) But in all the examples I found over the web, openfoam, ... BC used are always : Inlet : velocity Outlet : static pressure I don't understand why I have to fix different BC for my case. My professor told us that this was more physical and it probably will cause less numerical problems but he is used to work in compressible flows and said he was not familiar with incompressible solver. Also in openfoam as well as in Fluent, the BC with the total pressure seem to cause problem (convergence in openfoam) When I fix the BC with the velocity @15.7m/s at the inlet, my total pressure is quite different from the one we calculate with empirical formulas and I don't understand why (We found 102 Pa total pressure and the numerical simulation with velocity inlet give us 500 Pa in fluent) Can someone help me to solve my problem with the boundary conditions with the total pressure? Our professor told us it could be the mesh ... Thanks for the help!!

 February 28, 2011, 17:37 #2 New Member   Mark Beal Join Date: Feb 2011 Posts: 24 Rep Power: 6 Did you find the solution to this? Mark

February 28, 2011, 22:18
#3
Senior Member

Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 418
Rep Power: 15
Quote:
 Originally Posted by Gearb0x Hello, I have to simulate a curved pipe. Here are the caracteristics of the flow : Incompressible Re = 40 000 Hydraulic diameter = 0.08m So we found U characteristic = 15.7m/s k = 0.3676 epsilon = 0.3328 (with a turbulent intensity of 6, don't remember the characteristical length used to calculate this)
OK.

Quote:
 Since I begin in CFD, I've asked help to a professor wich told me that I should impose these BC : Outlet = Ambient pressure Inlet = Total pressure To calculate the total pressure, we used empirical relation for pressure loss and we find pstatic by adding the pressure loss to the outlet pressure and we just have to add rho u² to have total pressure (u characteristic is used for u)
It's Ok.

Quote:
 But in all the examples I found over the web, openfoam, ... BC used are always : Inlet : velocity Outlet : static pressure
It's usual because you know the volumetric flux or mass flux of your pipe, then setting an static pressure at the outlet you can obtain the pressure losses. Another related problem is to select the correct diameter to have the desired volumetric flux with no more of X pressure loss.

Quote:
 I don't understand why I have to fix different BC for my case. My professor told us that this was more physical and it probably will cause less numerical problems but he is used to work in compressible flows and said he was not familiar with incompressible solver. Also in openfoam as well as in Fluent, the BC with the total pressure seem to cause problem (convergence in openfoam) When I fix the BC with the velocity @15.7m/s at the inlet, my total pressure is quite different from the one we calculate with empirical formulas and I don't understand why (We found 102 Pa total pressure and the numerical simulation with velocity inlet give us 500 Pa in fluent) Can someone help me to solve my problem with the boundary conditions with the total pressure? Our professor told us it could be the mesh ... Thanks for the help!!
Setting total pressure at inlet allows you to obtain a completely developed profile in all pipe extension, which should match the pressure losses better than a velocity inlet. This is because velocity inlet BC generates a developing zone, near the inlet, which is larger as viscosity increases.
The differences are probably due meshing issues, you need to satisfy y+ parameters in near wall regions in order to use k-epsilon model appropriately (check these Guidelines). Once you have a good mesh then results for pressure-pressure and velocity-pressure should be slightly different for Re outside of laminar region (region of validity of k-epsilon model and short development zone)

Regards.
__________________
Santiago MÁRQUEZ DAMIÁN, Ph.D.
Post-doctoral Fellow
Research Center for Computational Mechanics (CIMEC) - CONICET/FICH-UNL
T.E.: 54-342-4511594 Ext. 1005
Güemes 3450 - (3000) Santa Fe
Santa Fe - Argentina
http://www.cimec.org.ar

 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 prashanthreddyh FLUENT 1 December 2, 2009 12:06 Vijay FLUENT 2 January 14, 2009 13:24 Chebeba CFX 1 March 16, 2008 03:00 Tudor Miron CFX 17 March 19, 2004 20:23 Tom Main CFD Forum 0 November 5, 2002 02:54

All times are GMT -4. The time now is 23:55.

 Contact Us - CFD Online - Top