CFD Online Discussion Forums

CFD Online Discussion Forums (
-   STAR-CCM+ (
-   -   Plotting Shear Stress (

michael.collins January 21, 2016 08:16

Plotting Shear Stress

I'm modelling a cup and bob viscometer. This consists of two concentric cylinders - a hollow stationary outer cylinder with an inner rotating cylinder. A gap exists between the cylinders which is filled with fluid (let's say water). For the rotational speeds I'm looking at, the flow is laminar, incompressible and viscous. I have the no-slip condition applied to all surfaces so this should create Couette flow (i.e. constant shear stress) at all points across the gap.

I have plotted the velocity profile across the gap and am seeing a linear velocity gradient (consistent with Couette flow and a constant shear stress). I want to plot the shear stress across the gap to verify this constant shear stress. However, when I use the Wall Shear Stress function to plot this I always just get zero shear stress at every point in the flow. Does anybody know why this is happening? Am I using the wrong function?


fluid23 January 21, 2016 12:39

According to the help files, wall shear stress is only defined in cells that exist on a wall face (aka a cell with a wall boundary). Therefore I would venture to guess that is why you see zero everywhere (although you should see something at the walls).

If you want to verify the shear stress in the whole domain you will need to write a field function. It should be fairly straight forward if you work in a cylindrical coordinate system. However, if you have already verified that your velocity gradient is correct then your shear stress field function will also be correct since du/dy (or in cylindrical coordinates dVtheta/dR) is really what drives the stress calculation.

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