CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > CFX

pressure calculation-junction box

Register Blogs Community New Posts Updated Threads Search

Like Tree2Likes
  • 1 Post By ghorrocks
  • 1 Post By hosseinhgf

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 1, 2010, 19:47
Default pressure calculation-junction box
  #1
Member
 
h-h
Join Date: Sep 2009
Posts: 38
Rep Power: 16
hosseinhgf is on a distinguished road
hello
in journal bearing, half sommerfeld condition must be applied.
i.e we must set the negative pressure to zero.
i read a little of junction box routines help but i don t know what should i do?
how we can call pressure on all cell and return it?
hosseinhgf is offline   Reply With Quote

Old   November 1, 2010, 20:37
Default
  #2
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,703
Rep Power: 143
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
Can you explain what boundary condition you want to apply? I do not know what the Sommerfield condition is.
ghorrocks is offline   Reply With Quote

Old   November 2, 2010, 02:02
Default
  #3
Member
 
h-h
Join Date: Sep 2009
Posts: 38
Rep Power: 16
hosseinhgf is on a distinguished road
in the journal bearing pressure is positive in converging portion and negative where it diverges.but Such pressures are rarely encountered in real bearings.
Mineral oils contain between 8 and 12% dissolved air. This air will start to come out of solution whenever the pressure falls below the saturation pressure. In many situations the saturation pressure is similar to the ambient pressure surrounding the bearing, and in these cases gas liberation will maintain the pressure in the divergent clearance space at close to the ambient level.to predicted this we use half sommerfeld condition.
we setting all values of negative gauge pressures to the value of the operating pressure. Thus, the gauge pressures for the top of journal are set to zero when negative values occur.
it is done before with fluent udf with this macro by "K.P. Gertzos" :
DEFINE_ON_DEMAND macro
thread_loop_c, begin_c_loop
thread_loop_f, begin_f_loop
i think this condition can be applied just for cell on the faces that has negative pressure and it is not necessary to add the all cell in the volume because pressure gradient in r direction is zero.
hosseinhgf is offline   Reply With Quote

Old   November 2, 2010, 05:20
Default
  #4
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,703
Rep Power: 143
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
This boundary is ill-defined and is not physical. Either you specify a wall and let the solver calculate the pressure for you, or you specify the pressure and let the solver calculate the flow velocity and direction. The second does not make sense in case as the walls are, well... walls. So use a wall boundary condition.

To account for the air coming out of solution you will need to put in a physical model for that. You might be able to write a custom density function or maybe use the cavitation model, adapted to dissolved air.

What are you trying to get out of this simulation?
ghorrocks is offline   Reply With Quote

Old   November 2, 2010, 09:49
Default
  #5
Member
 
h-h
Join Date: Sep 2009
Posts: 38
Rep Power: 16
hosseinhgf is on a distinguished road
Thank you
These are from some journals:


-“Cavitation is modeled by setting all calculated negative pressure and their gradients equal to zero throughout the iteration. The iteration is repeated until the pressure satisfies the convergence criterion”


-“When the film pressure drops down to the atmospheric pressure due to heavy external load or high operating rotational speed, rupture of the film or cavitation occurs. Taking into account this situation, the following boundary condition is used which applies the non-sub-atmospheric pressure constraint:
p>=p_ambient at pi<theta<2*pi

When cavitation is included, the vaporization pressure, the surface tension coefficient and the non-condensable gas mass fraction are also needed.”


-“As a constraint, the hydrodynamic pressure needs to be greater than the liquid cavitation pressure everywhere in the flow domain, i.e.
p>=p_cav at pi<theta<2*pi



In all journals that simulated “journal bearing” with cfd, this condition is applied.
Cfd simulation of journal bearing is part of my thesis.
It is done befor with fluent and cfx-Tascflow.
the half Sommerfeld solution results in more realistic predictions of load components, but this simple approach leads to a violation of the continuity of mass flow at the outlet end of the pressure curve. A better boundary condition is the Reynolds boundary condition but in one journal says that Reynolds BC cannot be use in fluent.
simulation based on the assumption of incompressible fluid and laminar flow
in the bearing.
What is cfx-Tascflow?
can junction box routines do it(like udf)?
hosseinhgf is offline   Reply With Quote

Old   November 2, 2010, 16:51
Default
  #6
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,703
Rep Power: 143
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
You have to be careful about these sort of models. The ones which describe simply limiting the pressure to a certain value which are not physically correct and are never going to converge anyway. CFX has a cavitation model buit-in which you should use instead.

I don't understand the film pressure boundary condition. Where is it applied? At the edge of the film wedge?

CFX-Tascflow is an old CFD solver which CFX purchased years ago. The rotating frames of reference, and basic solver approach in CFX5+ was adapted from CFX-Tascflow. Tascflow has been superseded by CFX5+.
elcino likes this.
ghorrocks is offline   Reply With Quote

Old   November 2, 2010, 17:31
Default
  #7
Member
 
h-h
Join Date: Sep 2009
Posts: 38
Rep Power: 16
hosseinhgf is on a distinguished road
the right bearing has "half sommerfeld condition"
http://www.4shared.com/photo/-tAGtC5-/jb_online.html

it applied to all cell in domain with negative pressure.
negative pressure occur in divergent portion.
as your suggestion i use cavitation model. it has not bad result.
i don t know value of oil saturation pressure.
sozen likes this.
hosseinhgf is offline   Reply With Quote

Old   May 5, 2015, 12:40
Default
  #8
New Member
 
Join Date: Oct 2014
Posts: 14
Rep Power: 11
sozen is on a distinguished road
Dear hosseinhgf,
can you explain how did you solve the problem? txh
sozen is offline   Reply With Quote

Reply


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
OpenMP in Junction Box Routine Hannes_Kiel CFX 10 September 21, 2010 13:51
Pressure Rise Error emueller CFX 0 May 5, 2009 11:08
Pressure calculation in Navier-Stokes equations Hadian Main CFD Forum 1 November 10, 2008 02:47
calculation of pressure difference wrong? Ralf Schmidt FLUENT 0 January 5, 2006 10:06
Gas pressure question Dan Moskal Main CFD Forum 0 October 24, 2002 22:02


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