# how can I input a wall slip law? Which Boundary Condition ?

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

 April 20, 2012, 09:28 how can I input a wall slip law? Which Boundary Condition ? #1 New Member   Join Date: Mar 2012 Posts: 26 Rep Power: 5 Hi, I want to input the wall slip law below with the limit yield stress above which there is a tangantial slip velocity : if if Can you advice me which boundary condition do I have to use ? ***A wall boundary condition with a moving wall and (?) shear rate or ***A velocity boundary condition with the velocity as a component or *** An other way My problem is a circlar pipe with a newtonian or non newtonian fluid which flow inside it. Thank you all...

 April 20, 2012, 10:36 #2 New Member   dengjiajia Join Date: Sep 2011 Posts: 14 Rep Power: 5 what is mean of your t? it is time ?

April 20, 2012, 10:57
#3
New Member

Join Date: Mar 2012
Posts: 26
Rep Power: 5
Quote:
 Originally Posted by dunga82 what is mean of your t? it is time ?
is the yield stress.
When then yield stress is above the limit on the wall, there is a wall slip velocity given by

where is the dynamic viscosity and is the friction coefficient.

(I have omitted in my previous message )

So what do you think about the boundary condition ?

 April 22, 2012, 01:44 #4 New Member   dengjiajia Join Date: Sep 2011 Posts: 14 Rep Power: 5 Wall Motion Based on Velocity Components For problems that include linear or non-linear translational motion of the wall boundary you can enable the Components option and specify the X-Velocity, Y-Velocity, and Z-Velocity of the wall. You can define non-linear translational motion using a boundary profile or a user-defined function for the X-Velocity, Y-Velocity, and/or Z-Velocity of the wall I think this condition is that you want.

April 22, 2012, 02:55
#5
Senior Member

Lucky Tran
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 606
Rep Power: 12
Quote:
 Originally Posted by tsi07 Hi, I want to input the wall slip law below with the limit yield stress above which there is a tangantial slip velocity : if if Can you advice me which boundary condition do I have to use ? ***A wall boundary condition with a moving wall and (?) shear rate or ***A velocity boundary condition with the velocity as a component or *** An other way My problem is a circlar pipe with a newtonian or non newtonian fluid which flow inside it. Thank you all...
You will need a UDF. Your basic wall properties will not handle what you are trying to do. Maybe you can dig one up. You will basically need to program Fluent (using the UDF) to assign the slip velocity to the slip-wall boundary condition each iteration according to the criteria you described.

 April 22, 2012, 06:41 #6 New Member   Join Date: Mar 2012 Posts: 26 Rep Power: 5 see below for the message.

April 22, 2012, 06:45
#7
New Member

Join Date: Mar 2012
Posts: 26
Rep Power: 5
Quote:
 Originally Posted by LuckyTran You will need a UDF. Your basic wall properties will not handle what you are trying to do. Maybe you can dig one up. You will basically need to program Fluent (using the UDF) to assign the slip velocity to the slip-wall boundary condition each iteration according to the criteria you described.

Here is my code below. If I understand well I have to put two component udf: wall_shear_x and wall_velocity_x.

Code:
#define Cf 1            /* friction coefficient */
#define TAU0 0.2            /* limit yield stress for slipping */

{
face_t f;
cell_t c0;
double TAUxy, mu;
double dudy;
{
dudy = C_DUDY(c0,t0);
mu = C_MU_L(c0,t0); /* viscosity */
TAUxy=mu*dudy;
}
}

{
face_t f;
cell_t c0;
double VC[ND_ND];
double TAUxy, VC_MAG, vg, mu;
double dudy;
{
dudy = C_DUDY(c0,t0);
mu = C_MU_L(c0,t0);                        /* viscosity */
VC[0]=-mu*dudy/Cf;
VC[1]=0;                /* stress vector's y componant*/
VC_MAG=NV_MAG(VC);                    /* stress vector's magnitude */
vg=VC[0];
if ( VC_MAG < TAU0)
{
}
}
}
}
And I check wall motion and apply udf wall_velocity_x in X-VELOCITY
and I check specified shear and I put udf wall_shear_x in X-VELOCITY.

Is it right ? I will test tomorrow but I wonder if I follow exactly the wall-slip model below with this boundary condition :

if
if

where is the dynamic viscosity and is the friction coefficient.

and with
where is the stress tensor and is the strain tensor.

What do you think about that ? Is the boundary condition appropriate to the model ??

April 22, 2012, 11:35
#8
Senior Member

Lucky Tran
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 606
Rep Power: 12
Quote:
 Originally Posted by tsi07 Thank you for your answer both. What do you think about that ? Is the boundary condition appropriate to the model ?? Thank for your answer...
You are going to apply the slip condition to the entire wall. So your criteria for slip or not, should be based on overall wall shear. Is that how you are checking?

In other words, you calculate wall shear for each individual cell next to the walls and then integrate/add to get the total force, and then use that total wall shear to compare to your slipping wall criteria. Put in other words, you cannot change the wall-slip velocity locally for each cell since the entire wall must have the same slip velocity.

Also, is it necessary to calculate TAUxy? Fluent should already have a variable called x-wall-shear or something that you should be able to access directly and not calculate through UDF.

I am not the best at UDF so I wanted to ask you to make sure.

 April 23, 2012, 04:10 #9 New Member   Join Date: Mar 2012 Posts: 26 Rep Power: 5 I don't find anywhere how to access directly the wall shear stress. So, in my case I have to calculate the element of the strain tensor in each cell to calculate the wall shear stress and to be able to compare this value with the limit stress . This comparison has to be made in each cell to know if there is a wall slip or not. My real purpose is to have a numerical wall slip model for non newtonian flow. That's why I have to calculate in each cell. With a newtonian flow it's not very usefull to make a loop on each cell. You say that it's impossible to have different slip velocity for each cell with a wall motion boundary condition (?). So how can I make this model ? Thank you.

April 23, 2012, 07:49
#10
Senior Member

Lucky Tran
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 606
Rep Power: 12
Quote:
 Originally Posted by tsi07 I don't find anywhere how to access directly the wall shear stress. So, in my case I have to calculate the element of the strain tensor in each cell to calculate the wall shear stress and to be able to compare this value with the limit stress . This comparison has to be made in each cell to know if there is a wall slip or not. My real purpose is to have a numerical wall slip model for non newtonian flow. That's why I have to calculate in each cell. With a newtonian flow it's not very usefull to make a loop on each cell. You say that it's impossible to have different slip velocity for each cell with a wall motion boundary condition (?). So how can I make this model ? Thank you.
Newtonian or non-newtonian that is a simple fluid property. Whether you have newtonian or non-newtonian the wall simply does not move the way you are trying to make it move. How can you have a section of wall moving in one location and non-moving in another? This is a boundary-condition problem, not a fluid property problem. The fact that you are trying to simulate a non-newtonian flow is immaterial to the fact that you are trying to make the wall move. For non-newtonian flow you just need to patch the viscosity based on the localized strain-rate. Recall again that you called it "wall" slip, not localized wall motion or portions of wall slip, it is simply [entire] wall slip. There is no need to model this type of behavior since it simply does not occur, the wall is bound by physical kinematic constraints.

The steady flow, the wall shear stress is again stored in a simple variable wall-shear-x y and z. Since it can be outputted to file / printed to screen, that variable should be accessible. Anyway, you can calculate if you really want to spend that effort.

 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 Saturn CFX 34 October 16, 2014 05:27 Alex F. FLUENT 14 August 6, 2014 05:08 HMR CFX 3 March 6, 2011 21:10 mehrdadeng CFX 10 February 25, 2011 06:25 Sharad Yeri FLUENT 0 May 18, 2008 12:32

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