CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   WallShearStress Utility a little doubt (http://www.cfd-online.com/Forums/openfoam-solving/59064-wallshearstress-utility-little-doubt.html)

 sek March 6, 2008 18:31

Has anyone had luck with wall

Has anyone had luck with wall shear stress calculation with OpenFOAM with RANS? First, let me say that I'm not using wall functions - my near wall mesh is around y+=1. We're using Spalart-Allmaras model. And we get an incredibly small wall shear stress, a few order of magnitude smaller than what it should be.

The WallShearStress utility computes the wall shear from the inner product of the wall normal and stress tensor (R). For all the turbulence models I've seen in OF, R is computed using "turbulent" viscosity (nut) and the velocity gradient, which is a bit surpprising to me. I suppose nuEff be used instead?

Since I'm not using wall function, I presume that OF does not mess up with nut, which would become very small in viscous sublayer and nutdU.dy is NOT the wall shear stress! can someone enlighten me?

 egp March 6, 2008 19:23

Hi Sung-Eun, wallShearStres

Hi Sung-Eun,

wallShearStress utility is only good for wall-functions.

If you have sublayer resolved grids, nut is essentially zero in the near-wall cells.

As you well know, you calculate wall shear for this type of simulation using the wall normal gradient, and the molecular viscosity. A quick way to do it is to run wallGradU, and multiply by mu in your post-processing software (e.g., Tecplot).

I agree that it would be best if wallShearStress was smart enough to handle laminar flow, sublayer resolution, and wall-functions. Moreover, it would be nice if the wall-functions were adaptive and could handle the transition from sublayer to y+~500 resolutions.

Eric

 pbo March 6, 2008 22:01

Have you ever tried to use the

Have you ever tried to use the turbulentForceAndTorque function object (OF-dev only) to compute the wall shear stress, instead of using the wallShearStress utility?

sources in \$FOAM_SRC/src/postProcessing/incompressible/forceAndTorque/turbulentForceAndTorq ueFunctionObject.C

 sek March 7, 2008 09:31