CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

Measuring wall shear stress in bend pipe

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

Like Tree38Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 16, 2014, 13:56
Default
  #61
New Member
 
Nara Shikamaru
Join Date: Apr 2012
Posts: 22
Rep Power: 14
shikamaru is on a distinguished road
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

\tau_{ij}=\mu (\frac{\partial u_i}{\partial x_j}+\frac{\partial u_j}{\partial x_i})

or, The total stress tensor including the pressure term should be used?

\sigma_{ij}=-p\delta_{ij}+\mu (\frac{\partial u_i}{\partial x_j}+\frac{\partial u_j}{\partial x_i})

Last edited by shikamaru; July 16, 2014 at 20:21. Reason: [Typo in equation]
shikamaru is offline   Reply With Quote

Old   July 16, 2014, 16:51
Default
  #62
Senior Member
 
Amir's Avatar
 
Amir
Join Date: May 2009
Location: Montreal, QC
Posts: 735
Blog Entries: 1
Rep Power: 23
Amir is on a distinguished road
Quote:
Originally Posted by shikamaru View Post
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

\tau_{ij}=\mu (\frac{\partial u_i}{\partial x_j}+\frac{\partial u_i}{\partial x_i})

or, The total stress tensor including the pressure term should be used?

\sigma_{ij}=-p\delta_{ij}+\mu (\frac{\partial u_i}{\partial x_j}+\frac{\partial u_i}{\partial x_i})
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:

\sigma_{ij}=-p\delta_{ij}+\mu (\frac{\partial u_i}{\partial x_j}+\frac{\partial u_j}{\partial x_i})

Bests,
__________________
Amir
Amir is offline   Reply With Quote

Old   July 16, 2014, 20:21
Default
  #63
New Member
 
Nara Shikamaru
Join Date: Apr 2012
Posts: 22
Rep Power: 14
shikamaru is on a distinguished road
Thanks for catching the typo, fixed.

Last edited by shikamaru; July 17, 2014 at 00:55.
shikamaru is offline   Reply With Quote

Old   July 17, 2014, 10:16
Default
  #64
New Member
 
Nara Shikamaru
Join Date: Apr 2012
Posts: 22
Rep Power: 14
shikamaru is on a distinguished road
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.
shikamaru is offline   Reply With Quote

Old   July 17, 2014, 10:28
Default
  #65
Senior Member
 
Amir's Avatar
 
Amir
Join Date: May 2009
Location: Montreal, QC
Posts: 735
Blog Entries: 1
Rep Power: 23
Amir is on a distinguished road
Quote:
Originally Posted by shikamaru View Post
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
Amir is offline   Reply With Quote

Old   August 9, 2014, 10:07
Default To Share
  #66
New Member
 
Jianzhi Li
Join Date: Jul 2013
Location: Shanghai, China
Posts: 27
Rep Power: 13
epi_c is on a distinguished road
Send a message via Skype™ to 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.
Attached Files
File Type: gz myWallShearStress.tar.gz (2.2 KB, 84 views)
Mojtaba.a likes this.
epi_c is offline   Reply With Quote

Old   August 17, 2014, 11:56
Default
  #67
New Member
 
Anastasios Stampoultzoglou
Join Date: May 2014
Posts: 21
Rep Power: 12
Tasos is on a distinguished road
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
Tasos is offline   Reply With Quote

Old   September 9, 2014, 12:34
Default
  #68
New Member
 
Daniel
Join Date: Nov 2012
Location: Glasgow, Scotland
Posts: 1
Rep Power: 0
derespin is on a distinguished road
Quote:
Originally Posted by epi_c View Post
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 \frac{kg}{ms^2} and not \frac{m^2}{s^2}. The results are not affected by this, though.

Regards
derespin is offline   Reply With Quote

Old   September 13, 2014, 06:25
Default
  #69
Retired Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,981
Blog Entries: 45
Rep Power: 128
wyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to all
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 View Post
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 View Post
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 View Post
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<volSymmTensorField> devReff() const;
In addition, I'll quote from another post of mine on another thread:
Quote:
Originally Posted by wyldckat View Post
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 View Post
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
JR22, Tasos and Fschi like this.
__________________
wyldckat is offline   Reply With Quote

Old   September 13, 2014, 09:55
Default
  #70
New Member
 
Anastasios Stampoultzoglou
Join Date: May 2014
Posts: 21
Rep Power: 12
Tasos is on a distinguished road
Mr Bruno thank you very very much, your help is priceless.
Tasos is offline   Reply With Quote

Old   December 9, 2015, 04:49
Default
  #71
New Member
 
LU
Join Date: Oct 2015
Posts: 11
Rep Power: 11
louvel is on a distinguished road
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!
louvel is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Water subcooled boiling Attesz CFX 7 January 5, 2013 04:32
2D wall shear stress yossarian88 STAR-CCM+ 0 February 21, 2011 14:12
relationship between wall shear stress and TKE winter FLUENT 0 December 11, 2007 18:11
Macro to set Wall Shear Stress Satish FLUENT 4 November 26, 2003 15:46
WALL SHEAR STRESS Min-Hua Wang CFX 0 April 3, 2002 18:40


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