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

wallShearStress vs. force functionObject

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

Like Tree2Likes
  • 1 Post By pad
  • 1 Post By mkiewat

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 14, 2013, 03:43
Question wallShearStress vs. force functionObject
  #1
pad
New Member
 
PA
Join Date: Mar 2009
Posts: 25
Rep Power: 17
pad is on a distinguished road
Hello,

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

Code:
forAll(wallShearStress.boundaryField(), patchI)
    {
        wallShearStress.boundaryField()[patchI] =
        (
           -mesh.Sf().boundaryField()[patchI]
           /mesh.magSf().boundaryField()[patchI]
        ) & Reff.boundaryField()[patchI];
    }
while forces.C uses

Code:
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.
__________________
Regards,
Per
pad is offline   Reply With Quote

Old   December 8, 2015, 13:30
Default
  #2
New Member
 
LU
Join Date: Oct 2015
Posts: 11
Rep Power: 11
louvel is on a distinguished road
Hallo

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
Default
  #3
New Member
 
Marco
Join Date: Nov 2014
Location: Germany
Posts: 14
Rep Power: 12
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.
cheng1988sjtu likes this.
mkiewat is offline   Reply With Quote

Reply

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

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Hydrodynamic force? arun7328 STAR-CCM+ 24 July 25, 2020 05:04
ActuatorDiskExplicitForce in OF2.1. Help be_inspired OpenFOAM Programming & Development 10 September 14, 2018 12:12
why does 'sample' do this? wallShearStress question CHARLES OpenFOAM Post-Processing 0 August 7, 2013 20:30
Force can not converge colopolo CFX 13 October 4, 2011 23:03
DEFINE_CG_MOTION and pressure force Teo Fumagalli FLUENT 0 April 11, 2008 11:25


All times are GMT -4. The time now is 06:27.