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

To calculate wall shear stress in of210 for les simulation

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

Like Tree18Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   December 11, 2012, 04:16
Default
  #21
Member
 
Dejan Morar
Join Date: Nov 2010
Posts: 78
Rep Power: 8
morard is on a distinguished road
Hi Nati,

yes, there are 3 components, in x, y and z direction. If you want to separate them, you can use foamCalc:

foamCalc components wallShearStress

this gives you: wallShearStressx, wallShearStressy and wallShearStressz.

For the lift and drag forces I can't help you much. Try to modify one of the post processing utilities to fit your case. It might help if you take a look into the folder /applications/test. You can find a lot of examples there.

Regards,
Dejan
nati and chengyu like this.
morard is offline   Reply With Quote

Old   December 11, 2012, 04:52
Default
  #22
New Member
 
Join Date: Nov 2011
Posts: 11
Rep Power: 6
nati is on a distinguished road
Hi Morard,

thank you for your help.

Excuse me, Do you know, how to calculate turbulent dissipation energy in each point?

Regards,
Nati


Quote:
Originally Posted by morard View Post
Hi Nati,

yes, there are 3 components, in x, y and z direction. If you want to separate them, you can use foamCalc:

foamCalc components wallShearStress

this gives you: wallShearStressx, wallShearStressy and wallShearStressz.

For the lift and drag forces I can't help you much. Try to modify one of the post processing utilities to fit your case. It might help if you take a look into the folder /applications/test. You can find a lot of examples there.

Regards,
Dejan
nati is offline   Reply With Quote

Old   December 12, 2012, 19:58
Default
  #23
New Member
 
Qiang Zhou
Join Date: May 2010
Location: Tongji University, Shanghai, China
Posts: 27
Rep Power: 8
michael1023 is on a distinguished road
What`s the difference between devReff and devBeff? In my case, the wallShearStressLES which include devReff can also be complied. Thank you.
michael1023 is offline   Reply With Quote

Old   December 12, 2012, 20:05
Default
  #24
New Member
 
Qiang Zhou
Join Date: May 2010
Location: Tongji University, Shanghai, China
Posts: 27
Rep Power: 8
michael1023 is on a distinguished road
Quote:
Originally Posted by wikstrom View Post
Reff is the effective (viscous+turbulent) stress tensor (calculated by the turbulence model). For laminar flow (or y-plus low) this would be more or less nu*grad(U).
Hi Wikstrom. What is the difference between the Reff and friction forces on a surface by a turbulence model, such as simulating a flow past over a circular cylinder by LES? Thank you.
michael1023 is offline   Reply With Quote

Old   December 13, 2012, 07:37
Default
  #25
Member
 
Niklas Wikstrom
Join Date: Mar 2009
Posts: 85
Rep Power: 9
wikstrom is on a distinguished road
Look at the wallShearStress.C (applications/utilities/postProcessing/wall/wallShearStress), e.g. in the the calcIncompressible() function.

Reff is a tensor, whereas the wall shear stress is the vector produced by the scalar product of 'wall face unit normal' and the Reff tensor. To get the friction force, then, you would need to surface integrate the wall shear stress.
wikstrom is offline   Reply With Quote

Old   December 18, 2012, 02:28
Default
  #26
New Member
 
Qiang Zhou
Join Date: May 2010
Location: Tongji University, Shanghai, China
Posts: 27
Rep Power: 8
michael1023 is on a distinguished road
Quote:
Originally Posted by wikstrom View Post
Look at the wallShearStress.C (applications/utilities/postProcessing/wall/wallShearStress), e.g. in the the calcIncompressible() function.

Reff is a tensor, whereas the wall shear stress is the vector produced by the scalar product of 'wall face unit normal' and the Reff tensor. To get the friction force, then, you would need to surface integrate the wall shear stress.
Thanks for your quite reply. I have read the code of wallShearStress.C. From what your said, I think that the wall shear stress has the same meaning with the friction force coefficient for three directions (x,y,z). Is it right?
In addition, Reff is a tensor for RAS model and how to calculate the tensor in LES model?
michael1023 is offline   Reply With Quote

Old   December 18, 2012, 04:15
Default
  #27
Member
 
sqing
Join Date: Sep 2012
Location: Dalian
Posts: 77
Rep Power: 5
Sunxing is on a distinguished road
Quote:
Originally Posted by wikstrom View Post
Reff is the effective (viscous+turbulent) stress tensor (calculated by the turbulence model). For laminar flow (or y-plus low) this would be more or less nu*grad(U).
Hi Niklas,

I think that for incompressible flow, the wallShearStress just give us the tau_w/rho, i.e nu*grad(U) , not mu*grad(U). Am i right?


regards,
Sunxing
Sunxing is offline   Reply With Quote

Old   June 24, 2013, 10:43
Smile it error for "essential value entry not provided" and my boundary is symmtryPlane
  #28
New Member
 
Du An
Join Date: Jun 2013
Posts: 2
Rep Power: 0
duan2013xjtu is on a distinguished road
Hi, friends! When I run 'wallShearStressLES', it exit:
--> FOAM FATAL IO ERROR:
essential value entry not provided

file: /home/cylinder/3.800000e+02/phi::boundaryField::sides1_half0 from line 425805 to line 425805.

From function fvsPatchField<Type>::fvsPatchField
(
const fvPatch& p,
const DimensionedField<Type, surfaceMesh>& iF,
const dictionary& dict
)

in file lnInclude/fvsPatchField.C at line 110.

FOAM exiting

my boundary is symmtryPlane. Do I need to change the symmetryPlane to patch?

best wishes,
Linda


Quote:
Originally Posted by morard View Post
Hi Nati,

yes, there are 3 components, in x, y and z direction. If you want to separate them, you can use foamCalc:

foamCalc components wallShearStress

this gives you: wallShearStressx, wallShearStressy and wallShearStressz.

For the lift and drag forces I can't help you much. Try to modify one of the post processing utilities to fit your case. It might help if you take a look into the folder /applications/test. You can find a lot of examples there.

Regards,
Dejan
duan2013xjtu is offline   Reply With Quote

Old   December 27, 2013, 21:21
Default Code
  #29
New Member
 
Oscar Ochoa
Join Date: Nov 2013
Posts: 4
Rep Power: 4
Gatoscar is on a distinguished road
Hello everyone:

Despite the fact that this is an old thread, I attach the code that I modified following the instructions posted here. I verified the results with one case of incompressible fluid in a pipe with dynamicSmagorinsky LES turbulence model. This is for OF2.1.1.

I hope this will be helpful for someone.

Regards
Attached Files
File Type: zip wallShearStressLES.zip (2.8 KB, 133 views)
wyldckat, tiam and Shahr like this.
Gatoscar is offline   Reply With Quote

Old   December 28, 2013, 13:44
Default
  #30
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,551
Blog Entries: 39
Rep Power: 97
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Greetings Oscar,

Many thanks for sharing! Can I convince you to add this to openfoamwiki.net?

You can follow the instructions at this page: http://openfoamwiki.net/index.php/Extend-bazaar
Or follow the examples on this page: http://openfoamwiki.net/index.php/Main_ContribUtilities

Best regards,
Bruno
Gatoscar likes this.
__________________
wyldckat is offline   Reply With Quote

Old   January 10, 2014, 17:27
Default
  #31
New Member
 
Oscar Ochoa
Join Date: Nov 2013
Posts: 4
Rep Power: 4
Gatoscar is on a distinguished road
Hello Bruno

Sorry for the late response. I looked the two links you sent, and I'm convinced to upload the code of the wallShearStressLES, however, I'm gonna make more tests to the code to be completly sure that it works in the right way.

Thanks for the information about the foam-extend and OpenFOAM community.

Regards
wyldckat and vishwakarma like this.
Gatoscar is offline   Reply With Quote

Old   May 11, 2014, 12:28
Default
  #32
Senior Member
 
ArathoN
Join Date: Jul 2011
Posts: 137
Rep Power: 8
ArathoN is on a distinguished road
I tried building this on OF 2.3 but it gives me this error:

/opt/openfoam230/src/turbulenceModels/compressible/LES/LESModel/LESModel.H:150:34: note: no known conversion for argument 4 from ‘Foam::basicThermo’ to ‘const Foam::fluidThermo&’
make: *** [Make/linux64GccDPOpt/wallShearStressLES.o] Errore 1

How can i solve it?
ArathoN is offline   Reply With Quote

Old   May 15, 2014, 08:14
Default
  #33
New Member
 
Hans Barósz
Join Date: May 2014
Posts: 22
Rep Power: 4
HanSolo123 is on a distinguished road
Same problem here. Please let me know when you solved it.
HanSolo123 is offline   Reply With Quote

Old   June 2, 2014, 12:26
Default
  #34
New Member
 
Join Date: Feb 2013
Posts: 18
Rep Power: 5
Shahr is on a distinguished road
Quote:
Originally Posted by ArathoN View Post
I tried building this on OF 2.3 but it gives me this error:

/opt/openfoam230/src/turbulenceModels/compressible/LES/LESModel/LESModel.H:150:34: note: no known conversion for argument 4 from ‘Foam::basicThermo’ to ‘const Foam::fluidThermo&’
make: *** [Make/linux64GccDPOpt/wallShearStressLES.o] Errore 1

How can i solve it?
Well, obviously you should replace the basicThermo to fluidThermo in the source code and header file. It will be fixed.
Shahr is offline   Reply With Quote

Old   June 22, 2014, 13:58
Default
  #35
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,551
Blog Entries: 39
Rep Power: 97
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Greetings to all!

I've done a quick adaptation of the source code by Oscar Ochoa to OpenFOAM 2.2.x and 2.3.x.
It's currently provided here: https://github.com/wyldckat/wallShearStressLES
Download and installation instructions are given on that page.

Best regards,
Bruno
caduqued and Shahr like this.
__________________
wyldckat is offline   Reply With Quote

Old   October 13, 2014, 21:46
Default
  #36
Senior Member
 
Huang Xianbei
Join Date: Sep 2013
Location: CAU,China
Posts: 277
Rep Power: 5
huangxianbei is on a distinguished road
Quote:
Originally Posted by wyldckat View Post
Greetings to all!

I've done a quick adaptation of the source code by Oscar Ochoa to OpenFOAM 2.2.x and 2.3.x.
It's currently provided here: https://github.com/wyldckat/wallShearStressLES
Download and installation instructions are given on that page.

Best regards,
Bruno
Hi,Bruno:
I have a question about the difference between wallShearStress calculated by the code and by the wallGradU. As in the definition, skin-friction coefficient is Cf = tau_wall/(0.5*rho*U^2), and tau_wall = rho*nu*GradU, while in wallShearStress, the wallShearStress is calculated by n_wallnormal &Beff, this doesn't contain rho, right? I compared the wallShearStress calculated by the 2 ways, however, I get quite different values.

Here, calculated by wallShearStressLES and wallGradU are both contained in the .tar
result.tar.gz



The value is of almost 100 times deviation! (The fluid is gas, with rho =1.25kg/m^3 and nu = 1.45e-5 m^2/s) What's wrong with it?

Best regards

Xianbei
huangxianbei is offline   Reply With Quote

Old   October 18, 2014, 15:26
Default
  #37
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,551
Blog Entries: 39
Rep Power: 97
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Hi Xianbei,

If you had looked a bit more about this, specially look for posts where I mentioned "wallShearStress"... with Google, I can find with this:
Code:
site:www.cfd-online.com/Forums "wyldckat" "wallShearStress"
Mmm... a lot of threads... OK, the one where I pretty much answer your questions is in this post: wall shear stress post #6

Best regards,
Bruno
__________________
wyldckat is offline   Reply With Quote

Old   October 18, 2014, 21:45
Default
  #38
Senior Member
 
Huang Xianbei
Join Date: Sep 2013
Location: CAU,China
Posts: 277
Rep Power: 5
huangxianbei is on a distinguished road
Quote:
Originally Posted by wyldckat View Post
Hi Xianbei,

If you had looked a bit more about this, specially look for posts where I mentioned "wallShearStress"... with Google, I can find with this:
Code:
site:www.cfd-online.com/Forums "wyldckat" "wallShearStress"
Mmm... a lot of threads... OK, the one where I pretty much answer your questions is in this post: wall shear stress post #6

Best regards,
Bruno
Hi,Bruno:
Thank you for your reply. From the link, I get to know the two ways to calculate tau, I use the second formula as you quoted, also, if I multiply wallShearStress with rho, they are still not the same I'll check it again.

Xianbei
huangxianbei is offline   Reply With Quote

Old   October 19, 2014, 04:31
Default
  #39
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,551
Blog Entries: 39
Rep Power: 97
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Hi Xianbei,

I think that something very wrong is going on in your simulation. My guess is that the y+ values are way above 300, which results in a bad calculation by wallShearStress and/or wallGradU.

If you can provide more details about your case, it will make it considerably easier to help diagnose the problem. The main details that come t mind as being important:
  1. Is the simulation done with compressible or incompressible modelling?
  2. Regarding turbulence model: Laminar, RAS or LES?
  3. What are the inlet flow speeds and the flow speeds near the wall (the cell closest to the wall) and the sizes of the cells closest to the wall?

Best regards,
Bruno
__________________
wyldckat is offline   Reply With Quote

Old   October 19, 2014, 06:51
Default
  #40
Senior Member
 
Huang Xianbei
Join Date: Sep 2013
Location: CAU,China
Posts: 277
Rep Power: 5
huangxianbei is on a distinguished road
Quote:
Originally Posted by wyldckat View Post
Hi Xianbei,

I think that something very wrong is going on in your simulation. My guess is that the y+ values are way above 300, which results in a bad calculation by wallShearStress and/or wallGradU.

If you can provide more details about your case, it will make it considerably easier to help diagnose the problem. The main details that come t mind as being important:
  1. Is the simulation done with compressible or incompressible modelling?
  2. Regarding turbulence model: Laminar, RAS or LES?
  3. What are the inlet flow speeds and the flow speeds near the wall (the cell closest to the wall) and the sizes of the cells closest to the wall?

Best regards,
Bruno
Hi,Bruno
My case is a 90 bend duct with an aspect ratio of 6. Smagorinsky model is used with mean y+ <4, it's incompressible.
The inlet condition is obtain by interpolate the experimental velocity, and then mapped to the grid points.
The mesh's max aspect ratio is <700, which is OK when doing the checkMesh
Here is the comparison of the two results at inletcompare.ods.tar.gz
The first 7 rows of data are obtained using sample line, the 1st is distance, and following are the wallGradU vector and wallShearStress vector, each has 3 components. The x direction is the studied direction
Xianbei
huangxianbei 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
Measuring wall shear stress in bend pipe liguifan OpenFOAM 70 December 9, 2015 04:49
What is wall shear stress and shear strain rate? warlocklw Main CFD Forum 2 May 21, 2012 06:16
wall shear stress lingo FLUENT 2 June 2, 2003 03:40
Specifying shear stress on wall in 5.5.1 michael rasmussen CFX 1 May 19, 2003 09:49
plotting wall shear stress Ellen FLUENT 1 April 5, 2002 22:50


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