CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Post-Processing

wallShearStress vs. force functionObject

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

Like Tree1Likes
  • 1 Post By pad

LinkBack Thread Tools Display Modes
Old   August 14, 2013, 02:43
Question wallShearStress vs. force functionObject
New Member
Join Date: Mar 2009
Posts: 23
Rep Power: 10
pad is on a distinguished road

I have been simulating an impeller and extracting the power by using both wallShearStress utility (for moment calc. in paraview) and the force functionObject for gettting the moment directly.

Comparing the data I realized that the moment from viscous forces are in opposite direction, i.e. with wallShearStress utility and paraview I get the following moment from viscous forces

(0.300199 0.191866 1.55875)

while the force function object gives me

(-0.3002453 -0.1916713 -1.556391)

Looking at the code for wallShearStress and forces.c it seems that wallShearStress.C uses

forAll(wallShearStress.boundaryField(), patchI)
        wallShearStress.boundaryField()[patchI] =
        ) & Reff.boundaryField()[patchI];
while forces.C uses

forAllConstIter(labelHashSet, patchSet_, iter)
            label patchi = iter.key();

            vectorField Md
                mesh.C().boundaryField()[patchi] - coordSys_.origin()

            vectorField pf(Sfb[patchi]*(p.boundaryField()[patchi] - pRef));

            fm.first().first() += rho(p)*sum(pf);
            fm.second().first() += rho(p)*sum(Md ^ pf);

            vectorField vf(Sfb[patchi] & devRhoReffb[patchi]);

            fm.first().second() += sum(vf);
            fm.second().second() += sum(Md ^ vf);
Why does wallShearStress.C incorporate the neg. sign while forces.C does not? Which approach is the correct one to use when trying to calculate the pump / impeller power?
louvel likes this.
pad is offline   Reply With Quote

Old   December 8, 2015, 13:30
New Member
Join Date: Oct 2015
Posts: 11
Rep Power: 3
louvel is on a distinguished road

do you already konw the answer?
i am not also have the same questiones
I donnot understand the difference to calculate visocus forces between wall shear stress and force utilities .the code are different but the result should be the same
louvel is offline   Reply With Quote

Old   March 8, 2016, 12:24
New Member
Join Date: Nov 2014
Location: Germany
Posts: 14
Rep Power: 4
mkiewat is on a distinguished road
I would guess that wallShearStress gives you the shear stress forces acting on the fluid while forces returns the total force (pressure AND shear stress forces) acting on a specified patch (the wall of an object). Because actio = reactio, the sign is opposite, the force vectors point in oopposite direction. So if you want to compute the moments acting on your impeller or whatever you should go with the "forces" function object.
mkiewat is offline   Reply With Quote


force, forces.c, functionobject, wall shear stress, wallshearstress.c

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
Hydrodynamic force? arun7328 STAR-CCM+ 17 September 17, 2014 03:38
ActuatorDiskExplicitForce in OF2.1. Help be_inspired OpenFOAM Programming & Development 8 July 3, 2014 10:54
why does 'sample' do this? wallShearStress question CHARLES OpenFOAM Post-Processing 0 August 7, 2013 19:30
Force can not converge colopolo CFX 13 October 4, 2011 22:03
DEFINE_CG_MOTION and pressure force Teo Fumagalli FLUENT 0 April 11, 2008 10:25

All times are GMT -4. The time now is 23:58.