wallGradU and wallShearStress, inconsistent result for RAS model

lzhou · August 19, 2014, 12:08

Hi,

I am doing a simple simulation of water flowing on a plate with RAS k-epsilon model.
My goal is to verify the simulation results with different roughness height on the surface of the plate and to calculate the wall shear stress of each case.

However, the utility wallGradU gives (36.5451 -8.58721e-05 0.000132998) for one point on the bottom wall and the utility wallShearStress gives (-0.000808236 3.08435e-09 -2.97856e-09) at this point. As I am doing a simulation of water so the dynamic (shear) viscosity is 8.9e-4 Pa·s. And the wall shear stress can be calculated by the equation:

So I think at this point ,the wall shear stress should be 8.9e-4*36.5451=0.0325 pa.

I would like to ask why there is the difference between the results of wallGradU and wallShearStress ? If possible, could anyone tell me the equations used in OpenFoam to calculate the wall shear stress ? Thanks !

lzhou · August 19, 2014, 12:22

To add some details. The results I got is for a smooth plate with nukRoughWallFunction and Ks=0 Cs=0.5 as the boundary condition at the bottom wall. And I did not find a place to add the density of the fluid into the simulation so I guess the wall shear stress is calculated with the kinematic viscosity rather than the dynamic viscosity. Or it assumes that the density is 1 kg/m^3. I am a starter of CFD and OpenFoam so can someone help me please. Thanks !

August 19, 2014, 12:08	wallGradU and wallShearStress, inconsistent result for RAS model	#1
lzhou New Member Lu ZHOU Join Date: Jul 2014 Location: Lyon, France Posts: 12 Rep Power: 12	Hi, I am doing a simple simulation of water flowing on a plate with RAS k-epsilon model. My goal is to verify the simulation results with different roughness height on the surface of the plate and to calculate the wall shear stress of each case. However, the utility wallGradU gives (36.5451 -8.58721e-05 0.000132998) for one point on the bottom wall and the utility wallShearStress gives (-0.000808236 3.08435e-09 -2.97856e-09) at this point. As I am doing a simulation of water so the dynamic (shear) viscosity is 8.9e-4 Pa·s. And the wall shear stress can be calculated by the equation: So I think at this point ,the wall shear stress should be 8.9e-4*36.5451=0.0325 pa. I would like to ask why there is the difference between the results of wallGradU and wallShearStress ? If possible, could anyone tell me the equations used in OpenFoam to calculate the wall shear stress ? Thanks !

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
wallGradU or wallShearStress on Time-Averaged Velocity Field	minger	OpenFOAM Post-Processing	6	November 7, 2019 06:30
wallShearStress, wallGradU are inconsistent	dancfd	OpenFOAM Post-Processing	5	July 18, 2018 09:33
calculate and sample wallGradU & wallShearStress every output time	ArathoN	OpenFOAM Post-Processing	8	July 11, 2018 13:14
k-omega-SST model (OF 1.6) - turbulent flat plate	cboss	OpenFOAM Running, Solving & CFD	25	August 9, 2016 10:53
wallShearStress VS wallGradU	bigrock	OpenFOAM Post-Processing	0	July 2, 2013 06:54

August 19, 2014, 12:22		#2
lzhou New Member Lu ZHOU Join Date: Jul 2014 Location: Lyon, France Posts: 12 Rep Power: 12	To add some details. The results I got is for a smooth plate with nukRoughWallFunction and Ks=0 Cs=0.5 as the boundary condition at the bottom wall. And I did not find a place to add the density of the fluid into the simulation so I guess the wall shear stress is calculated with the kinematic viscosity rather than the dynamic viscosity. Or it assumes that the density is 1 kg/m^3. I am a starter of CFD and OpenFoam so can someone help me please. Thanks !