# Measuring wall shear stress in bend pipe

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 July 16, 2014, 12:56 #61 New Member   Nara Shikamaru Join Date: Apr 2012 Posts: 22 Rep Power: 12 Thanks Amir for the formulation. Having taken solid mech courses I was also skeptical about the traction vector that was being calculated. I see that you are not active anymore in this forum, but big thanks for the helpful posts. I am was working on generating the traction on the interface between two incompressible immiscible fluids (VOF formulation). I am almost done with that. The only part I have doubt about is the stress tensor term being used in this formulation - Should the viscous stress tensor be used for calculating the Traction vector? This is what I would do if I was doing solids or, The total stress tensor including the pressure term should be used? Last edited by shikamaru; July 16, 2014 at 19:21. Reason: [Typo in equation]

July 16, 2014, 15:51
#62
Senior Member

Amir
Join Date: May 2009
Location: Montreal, QC
Posts: 735
Blog Entries: 1
Rep Power: 21
Quote:
 Originally Posted by shikamaru Thanks Amir for the formulation. Having taken solid mech courses I was also skeptical about the traction vector that was being calculated. I see that you are not active anymore in this forum, but big thanks for the helpful posts. I am was working on generating the traction on the interface between two incompressible immiscible fluids (VOF formulation). I am almost done with that. The only part I have doubt about is the stress tensor term being used in this formulation - Should the viscous stress tensor be used for calculating the Traction vector? This is what I would do if I was doing solids or, The total stress tensor including the pressure term should be used?
Hi,

Reasonably, the total traction force would be originated from the total stress tensor consisted of both pressure and viscous terms; however, it has a typo and needs to be modified like this:

Bests,
__________________
Amir

 July 16, 2014, 19:21 #63 New Member   Nara Shikamaru Join Date: Apr 2012 Posts: 22 Rep Power: 12 Thanks for catching the typo, fixed. Last edited by shikamaru; July 16, 2014 at 23:55.

 July 17, 2014, 09:16 #64 New Member   Nara Shikamaru Join Date: Apr 2012 Posts: 22 Rep Power: 12 Since pressure always acts normal to the surface, if anyone only wants the shear stress the pressure term is not needed. However, in my case I am calculating both normal and shear stress on a curved surface. So, I need the pressure term in my calculation. Please correct me if I am wrong.

July 17, 2014, 09:28
#65
Senior Member

Amir
Join Date: May 2009
Location: Montreal, QC
Posts: 735
Blog Entries: 1
Rep Power: 21
Quote:
 Originally Posted by shikamaru Since pressure always acts normal to the surface, if anyone only wants the shear stress the pressure term is not needed. However, in my case I am calculating both normal and shear stress on a curved surface. So, I need the pressure term in my calculation. Please correct me if I am wrong.
Well, this was a concept we were talking about all the time. The viscous term includes both normal and in-plane components, but the pressure term only contains normal component. So if anyone would be in need of shear stress, the pressure term must be definitely omitted and the viscous term should be modified....

Bests,
__________________
Amir

August 9, 2014, 09:07
To Share
#66
New Member

Jianzhi Li
Join Date: Jul 2013
Location: Shanghai, China
Posts: 27
Rep Power: 11
I wrote a post processing utility myWallShearStress according to this thread.

And thanks @Amir a lot.

My utility is general, suitable for incompressible and compressible flows using RAS or LES turbulence models.
Attached Files
 myWallShearStress.tar.gz (2.2 KB, 82 views)

 August 17, 2014, 10:56 #67 New Member   Anastasios Stampoultzoglou Join Date: May 2014 Posts: 21 Rep Power: 10 Hi all, I am a new user on openfoam . I read your posts and it was very useful. I would like to ask you which method the utility "wallShearStress" use (for example log law method or momentum method or Reynolds stress method, etc. ) Thank you very much

September 9, 2014, 11:34
#68
New Member

Daniel
Join Date: Nov 2012
Location: Glasgow, Scotland
Posts: 1
Rep Power: 0
Quote:
 Originally Posted by epi_c I wrote a post processing utility myWallShearStress according to this thread. And thanks @Amir a lot. My utility is general, suitable for incompressible and compressible flows using RAS or LES turbulence models.
Thanks for the utility! That being said, the units for the compressible stress tensor are wrong, as they should be and not . The results are not affected by this, though.

Regards

September 13, 2014, 05:25
#69
Retired Super Moderator

Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,969
Blog Entries: 45
Rep Power: 127
Greetings to all!

This has become a rather long and somewhat confusing thread, but unfortunately I will only be able to try and clarify Tasos' question:
Quote:
 Originally Posted by Tasos I read your posts and it was very useful. I would like to ask you which method the utility "wallShearStress" use (for example log law method or momentum method or Reynolds stress method, etc. )
The answer is given in some of the previous posts, from which I'll quote the most important details I could quickly spot:
Quote:
 Originally Posted by ngj BTW: In stead of making the stress tensor yourself, you could benefit from the turbulence model by calling devReff(). The only thing to remember is that this tensor is also non-tangential with the wall, so you have to perform the same projection (as you already do) onto the boundary face. The upside of using the turbulence model is that you do not have to implement turbulence model specific shear stress methods, because the turbulence model itself tells, how the shear stress is defined.
And:
Quote:
 Originally Posted by Nucleophobe Edit: I just tried this for my case, and I get the same result doing either of the following: tau = (wallShearStress) * rho tau = (wallGradU) * mu Is that expected? Was the wallShearStress utility fixed? I am solving a laminar incompressible flow with a density rho close to that of water.
@Tasos: you'll have to draw your own specific conclusions, but I'd say from the quotes above, the utility basis closely resembles the "Reynolds stress method" you've considered, namely "devReff()": https://github.com/OpenFOAM/OpenFOAM...inar/laminar.H
Quote:
 Code:  //- Return the effective sub-grid turbulence stress tensor // including the laminar stress virtual tmp devReff() const;
In addition, I'll quote from another post of mine on another thread:
Quote:
Originally Posted by wyldckat
My guess is that you did not take into account the fact that the wallShearStress utility will give you the results in the same units of the pressure field: http://www.cfd-online.com/Wiki/Wall_shear_stress
Which for incompressible flow, in OpenFOAM, it usually is in "kinematic pressure", namely in "m2/s2".

In post #46 of the other thread you pointed out, you'll find this comment:
Quote:
 Originally Posted by Nucleophobe Edit: I just tried this for my case, and I get the same result doing either of the following: tau = (wallShearStress) * rho tau = (wallGradU) * mu
This means that the value are seeing in OpenFOAM does not yet take into account the density of the fluid. You can use the Calculator filter in ParaView to multiply the "wallShearStress" field by the density value you have for your fluid.
I hope this answers your question and hopefully clear up some details asked above.

Best regards,
Bruno
__________________

 September 13, 2014, 08:55 #70 New Member   Anastasios Stampoultzoglou Join Date: May 2014 Posts: 21 Rep Power: 10 Mr Bruno thank you very very much, your help is priceless.

 December 9, 2015, 03:49 #71 New Member   LU Join Date: Oct 2015 Posts: 11 Rep Power: 9 Hello Amir could you tell me why n=-mesh.Sf().boundaryField()[patchi] / mesh.magSf().boundaryField()[patchi] there is a minus "-" here And I also check the code in Openfoam for the forces calculation,there also calculate the viscous force which is the also wall shear stress but in the utility ,the code is: sfb[patchI]&devrhoReffb[patch] do you think there are the same? Best regards!

 Thread Tools Search this Thread Search this Thread: Advanced Search 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 Off Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Attesz CFX 7 January 5, 2013 03:32 yossarian88 STAR-CCM+ 0 February 21, 2011 13:12 winter FLUENT 0 December 11, 2007 17:11 Satish FLUENT 4 November 26, 2003 14:46 Min-Hua Wang CFX 0 April 3, 2002 17:40

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

 Contact Us - CFD Online - Privacy Statement - Top