
[Sponsors] 
March 3, 2010, 16:01 
[HELP] Slip flow boundary condtion in CFX

#1 
New Member
Mengbai Wu
Join Date: Mar 2010
Posts: 7
Rep Power: 8 
Hey friends:
I am dealing with this gas flow field calculation in CFX. I have a relatively large Knudsen number, which is about 0.03. In this range NS equations still work but the noslip boundary condition needs to be revised accoding to either the first order or the second order slip velocity condition. I am not sure whether CFX can handle this problem. I am new to CFX. Thanks for any possible help. 

March 3, 2010, 17:57 

#2 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,638
Rep Power: 98 
You will need to define a wall shear stress to do this. Hopefully you can do this using CEL.


March 3, 2010, 18:47 

#3 
New Member
Mengbai Wu
Join Date: Mar 2010
Posts: 7
Rep Power: 8 
Thanks very much.
I am using ANSYS CFX 11.0. I don't think it can define wall shear stress directly. I just saw a ANSYS 12.1 technical specification note, which showed that ANAYS 12.1 can do finite slip and wall specified shear stress. CEL is a good idea, I am thinking about do a user subroutine now. Thank you again. 

June 8, 2011, 03:26 

#4 
New Member
siva
Join Date: Sep 2009
Posts: 7
Rep Power: 9 
hai
I too have the same problem. i can see some documentation in ansys 12 where wall shear stress components can be given but this option is not available in ansys 11. have u identify this option in ansys 12. 

September 14, 2012, 06:14 

#5 
New Member
Ramin Mirzazadeh
Join Date: Aug 2012
Posts: 16
Rep Power: 6 
Hello
Having a knudsen number around 0.01, I tried to do a laminar flow analysis in a simple rectangular micro duct, I used CEL for introducing the first order slip boundary model through specified shear stress on walls. I got consistent results comparing with analytical solution for steady state analysis, which shows the CEL is physically right, when I do the same analysis but in transient approach, cfx gives error "Floating point exception: Overflow". I read http://www.cfdonline.com/Wiki/Ansys...publishable.3F, and I tried all of its suggestions (mesh change, time step change,short intial time steps,upwind option, no. of coef. loops change )but does not work. can anyone give any hint on it? is there any problem in CFX for transient analysis of slip boundary? Thanks in advance 

September 14, 2012, 07:48 

#6 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,638
Rep Power: 98 
First of all I would run your simulation with normal no slip boundaries. That makes sure your general setup is OK.


September 14, 2012, 09:29 

#7 
New Member
Ramin Mirzazadeh
Join Date: Aug 2012
Posts: 16
Rep Power: 6 
exactly, I did that already. as soon as I use the CEL code that in steady state works, CFX cant handle anymore


September 14, 2012, 09:39 

#8 
New Member
Ramin Mirzazadeh
Join Date: Aug 2012
Posts: 16
Rep Power: 6 
for further information, I used the following CEL
ShearSlipU = ave(Velocity u)@wall1*1.8E5[N s m^2]/(1.2222*MeanFreePath) MeanFreePath = boltzmann*298[K]/(sqrt(2)*pi*ave(Absolute Pressure )@REGION:B16*MolecularDiam^2) MolecularDiam = 3.7e10[m] 1.8e5 is viscosity and 1.2222 is a constant related to tangential momentum accommodation coefficient. 

September 14, 2012, 10:27 

#9 
Senior Member
Lance
Join Date: Mar 2009
Posts: 595
Rep Power: 12 
Wouldnt ShearSlipU become an extremely small number? like 1e17 or so? What do you estimate it to be?


September 14, 2012, 10:54 

#10 
New Member
Ramin Mirzazadeh
Join Date: Aug 2012
Posts: 16
Rep Power: 6 
ciao Lance,
shear stress at wall, in CFX for steady state analysis yields to 0.2499 Pa and analytically it should be 0.226 Pa 

September 14, 2012, 11:00 

#11 
Senior Member
Lance
Join Date: Mar 2009
Posts: 595
Rep Power: 12 
I see. Well, one reason to the overflow could be that during the startup of your simulation one of the variables (U or abs pressure) fluctuate a bit and causes your expression to overflow.


September 14, 2012, 11:08 

#12 
New Member
Ramin Mirzazadeh
Join Date: Aug 2012
Posts: 16
Rep Power: 6 
I had set initialisation values to Zero, it might be the reason for overflow, I will do analysis with U value half of expected steady state value which is 0.6225 m/s. can this help analysis to converge?
thanks for you idea 

September 15, 2012, 07:12 

#13 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,638
Rep Power: 98 
If you are still having convergence difficulties then I suspect the issue will be the way you are coupling the velocity at the control volume adjacent to the wall to the wall shear stress. This, in effect, sets up a new numerical relationship which needs solution. As the ave() function is evaluated based on the variabel values from the last iteration, this means you are in effect doing an explicit calculation of the the wall condition and this means your allowable time step size will be massively reduced or the simulation goes unstable and diverges.
If this is correct then the fix is simple  a smaller time step size. And the required time step might be massively smaller than you currently run  like 1000 times smaller, not half. 

September 15, 2012, 07:17 

#14 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,638
Rep Power: 98 
Also, why do you use an average function anyway? If you define ShearSlipU = Velocity u*1.8E5[N s m^2]/(1.2222*MeanFreePath) then the function can vary over the wall, which sounds more correct to me.


September 18, 2012, 10:14 

#15 
New Member
Ramin Mirzazadeh
Join Date: Aug 2012
Posts: 16
Rep Power: 6 
changing the initialisation values did not help.
thanks Glenn for the idea, I tried sensitivity analysis and the problem converged. I used 1e7 s timsteps, also I used tight RMS residual target 1e10 with max coeff. loop of 100. I am not professional with CFX and I read cfx manual and this forum to get information about RMS residual target and coeff loops. cfx manual suggests RMS 1e5 and up to 10 loops for most of problems but still not sure what this value should be for Microscale problems or the importance of coeff loops. it would be helpful if you can give some hints about them. Thanks 

September 18, 2012, 19:38 

#16 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,638
Rep Power: 98 
The documentation's guide is good for most classes of flow. As long as your Re>0.01 or so you should be fine. And for single phase flows you should probably aim for 35 coeff loops for time step, 10 is a bit too many.


March 13, 2013, 11:14 

#17 
New Member
kullimanali
Join Date: Dec 2012
Posts: 13
Rep Power: 5 

March 13, 2013, 11:17 

#18 
Senior Member
Lance
Join Date: Mar 2009
Posts: 595
Rep Power: 12 
1) why not try the fluent forum?
2) CEL is not available in fluent 3) you dont calculate wall shear stress with CEL 

March 13, 2013, 11:17 

#19  
New Member
kullimanali
Join Date: Dec 2012
Posts: 13
Rep Power: 5 
Quote:
i am doing same thing,can you explain if u got solution, how to calculate shear stress at wall 

March 13, 2013, 11:38 

#20 
New Member
kullimanali
Join Date: Dec 2012
Posts: 13
Rep Power: 5 

Tags 
slip flow 
Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Flow Direction: normal to boundary!!!  Atit  CFX  1  August 2, 2015 13:42 
3D Compressible Flow Boundary Conditions  Tyler  FLUENT  4  February 5, 2009 20:58 
Mass flow and UMom flow in CFX  Zhihua Xie  CFX  0  September 3, 2007 09:49 
Slip boundary condition what is inside  normunds  OpenFOAM Running, Solving & CFD  2  June 4, 2007 06:45 
maintaining a logarithmic velocity distribution  Morten Andersen  CFX  1  January 8, 2007 12:37 