CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Fluent UDF and Scheme Programming

Macros to compute forces

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 3, 2012, 06:49
Default Macros to compute forces
  #1
New Member
 
Join Date: Apr 2012
Posts: 17
Rep Power: 5
Faby is on a distinguished road
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
Faby is offline   Reply With Quote

Old   May 3, 2012, 10:29
Default
  #2
Member
 
Ryne
Join Date: Jan 2010
Posts: 32
Rep Power: 7
Rhyno466 is on a distinguished road
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
Rhyno466 is offline   Reply With Quote

Old   May 3, 2012, 12:45
Default
  #3
New Member
 
Join Date: Apr 2012
Posts: 17
Rep Power: 5
Faby is on a distinguished road
Quote:
Originally Posted by Rhyno466 View Post
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
Wow!!!
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
Faby is offline   Reply With Quote

Old   May 3, 2012, 16:15
Default
  #4
Member
 
Ryne
Join Date: Jan 2010
Posts: 32
Rep Power: 7
Rhyno466 is on a distinguished road
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.
Rhyno466 is offline   Reply With Quote

Old   May 4, 2012, 03:21
Default
  #5
New Member
 
Join Date: Apr 2012
Posts: 17
Rep Power: 5
Faby is on a distinguished road
Quote:
Originally Posted by Rhyno466 View Post
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.
Thank you very much for help! I hope it will work for me!
Faby is offline   Reply With Quote

Old   May 4, 2012, 04:50
Default
  #6
New Member
 
Join Date: Apr 2012
Posts: 17
Rep Power: 5
Faby is on a distinguished road
Quote:
Originally Posted by Rhyno466 View Post

As for the pressure, I multiply the pressure by the directional cell projected area to get the force.
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?
Faby is offline   Reply With Quote

Reply

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
forces in interFoam Svensson OpenFOAM Post-Processing 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


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