CFD Online Logo CFD Online URL
Home > Forums > Phoenics

zero velocity gradient

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

LinkBack Thread Tools Display Modes
Old   February 10, 2006, 10:31
Default zero velocity gradient
Posts: n/a
Hi, What is the correct way of setting zero velocity gradient at the outlet. Thanks.
  Reply With Quote

Old   February 10, 2006, 14:17
Default Re: zero velocity gradient
Posts: n/a
By default all boundaries in PHOENICS are zero flux. i.e d(phi)/dn = 0.

You should set a boundary condition that sucks the same mass out from each cell that is entering it. You could probably do this using inform, or easily with ground. But you will need to calculate\estimate the rho*vel*Area of the fluid as it enters the outlet cell for this you will need the previous cell density and velocity and cell face area. This could be done in a flash using ground.
  Reply With Quote

Old   February 11, 2006, 08:53
Default Re: zero velocity gradient
Posts: n/a
Thank you Mick. Best regards
  Reply With Quote

Old   February 11, 2006, 09:40
Default Re: zero velocity gradient
Posts: n/a
We are using version 3.6. In outlet attributes window, for velocities we have USER SET, IN CELL, DEDUCED options. İn order to assign zero velocity gradient at the outlet which one we should choose.And also for zero temperature gradient which option should be selected.Thank you.
  Reply With Quote

Old   February 16, 2006, 06:01
Default Re: zero velocity gradient
Posts: n/a
I think that you need to set the outlet BC yourself rather then use standard outlet. As I described above, this will be more like an inlet but sucking mass out of the domain.

As I said before you should be sucking out the same mass from the oulet cell that enters it. I would use ground coding, but you might try experimenting with INFORM to see what you can do. You need access to the velocity of the preceeding cell and density and cell face area .

Perhaps some ground coding.

in Q1 : Patch(exit,east,nx,nx,1,ny,1,nz,1,lstep)

coval(exit,P1,fixflu, grnd)

In ground.for in group 13, something like

IF(NPATCH(1:4).eq.'EXIT')THEN l0val=L0F(VAL) L0AE=L0F(EAST) L0DEN=L0F(DEN1) L0VEL=L0F(WEST(U1)) ! Try this to get Vel for previous cell otherwise you will have to have an Icell2 to pick up slab NX-1

do ix=1,nx

do iy=1,ny

icell =iy+(ix-1)*ny


enddo enddo

NOTE; you will probably need to fix a refrence pressure in the domain somewhere to anchor pressure

With a normal outlet Pressure is just set, In this case I think that if you want DT/Dn=0 the you should set ; COVAL(EXIT,TEM1,0.0,SAME) as this will ensure that if any temperature enters through the oulet because of the exit pressure value, then it will enter with the same TEMP as in the last cell.

In terms of the MENU for OUTLET type BC and TEM1 this should be set to 'IN-CELL'.

Hope this helps, has anyone else out there got any better ideas or method on this. What about trying INFORM anyone!.


  Reply With Quote


Thread Tools
Display Modes

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

Similar Threads
Thread Thread Starter Forum Replies Last Post
Velocity u.Gradient X Mgtripple CFX 11 January 28, 2014 11:01
velocity gradient Merinadica FLUENT 0 December 16, 2006 05:53
Velocity gradient Stephen FLUENT 0 April 7, 2003 10:29
Velocity gradient Stephen Main CFD Forum 0 April 7, 2003 10:28
Velocity gradient Salman FLUENT 1 April 1, 2003 00:43

All times are GMT -4. The time now is 20:49.