
[Sponsors] 
May 3, 2012, 06:49 
Macros to compute forces

#1 
New Member
Join Date: Apr 2012
Posts: 17
Rep Power: 6 
Hello,
I'm writing a UDF to describe the motion of a particle in a tube fluid flow by computation of forces on particle surface. I found a FACE MACROS that compute the pressure on faces of the dynamic mesh zone, but I also need of the computation of viscous forces on the particle surface . Is there anyone who knows a MACROS to reach this aim? Thank you Last edited by Faby; May 3, 2012 at 06:52. Reason: I change idea 

May 3, 2012, 10:29 

#2 
Member
Ryne
Join Date: Jan 2010
Posts: 32
Rep Power: 8 
For viscous forces, I have been using
F_STORAGE_R_N3V(f,t1,SV_WALL_SHEAR)[*]; Where the star is the direction vector (0 for x,1 for y). I found this on this site and it was said that this is not found in the UDF manual and that someone found it online. Good Luck 

May 3, 2012, 12:45 

#3  
New Member
Join Date: Apr 2012
Posts: 17
Rep Power: 6 
Quote:
But...what is the correct way to use it? I have to declare face_t f; thread *t1; and SV_WALL_SHEAR ? Is a vector? Do I have to declare it? Besides, to obtain the resultant shear force, for instance, in x direction, do I have to write something like this: Fx = F_STORAGE_R_N3V(f,t1,SV_WALL_SHEAR) [0] ? Or F_STORAGE_R_N3V(f,t1,SV_WALL_SHEAR) computes the value of shear stress? Such as F_P(f,t) provides the value of total pressure and not of pressure force. I'm sorry for so many questions 

May 3, 2012, 16:15 

#4 
Member
Ryne
Join Date: Jan 2010
Posts: 32
Rep Power: 8 
Again, I found this somewhere on this forum and so I'm not completely sure about everything with it, however it is my understanding that this
" Fx = F_STORAGE_R_N3V(f,t1,SV_WALL_SHEAR) [0] ? " is correct. I have manually compared the results of this computation to that of the FLUENT output values for viscous forces and they seemed to be relatively in agreement. As for the pressure, I multiply the pressure by the directional cell projected area to get the force. 

May 4, 2012, 03:21 

#5  
New Member
Join Date: Apr 2012
Posts: 17
Rep Power: 6 
Quote:


May 4, 2012, 04:50 

#6 
New Member
Join Date: Apr 2012
Posts: 17
Rep Power: 6 
if F_STORAGE_R_N3V(f,t1,SV_WALL_SHEAR) [0] computes the value of the stress shear in direction x, to obtain the corresponding force, you should multiply the stress shear for the projected parallel area to get the force, isn't it?


Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
forces in interFoam  Svensson  OpenFOAM PostProcessing  2  May 9, 2012 13:31 
FORCES don't run!  C12Carbon  OpenFOAM  0  September 10, 2011 07:34 
Compute CEs from moments and forces  nico  Main CFD Forum  1  August 6, 2007 19:19 
What are UDF macros C_T_S , C_T_AP,...????  Asghari  FLUENT  0  January 28, 2007 11:54 
Valve Forces in CFdesign  Mike Clapp  Main CFD Forum  3  March 8, 2001 15:09 