|
[Sponsors] |
August 14, 2013, 03:43 |
wallShearStress vs. force functionObject
|
#1 |
New Member
PA
Join Date: Mar 2009
Posts: 25
Rep Power: 17 |
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]; } 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); }
__________________
Regards, Per |
|
December 8, 2015, 13:30 |
|
#2 |
New Member
LU
Join Date: Oct 2015
Posts: 11
Rep Power: 11 |
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 |
|
March 8, 2016, 12:24 |
|
#3 |
New Member
Marco
Join Date: Nov 2014
Location: Germany
Posts: 14
Rep Power: 12 |
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.
|
|
Tags |
force, forces.c, functionobject, wall shear stress, wallshearstress.c |
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |