I computed the flow around a sphere and also calculated the lift and drag forces/ coefficients. At present I got to high values for Cd (0.61) at Re = 1*10⁴, due to the k-epsilon turbulence model I guess.
But anyway I checked the code for the calculation in the 'forceCoeffs.C' file. Most of it is clear to me, but I did not underśtand a few lines:
Best regards Nico
fm contains the computed forces and moments data (see forces.C). fm.first() contains forces data, fm.second() contains moments data. Forces and moments are stored in pressure and viscous contributions, so that:
fm.first().first() --> normal force
fm.first().second() --> tangential force
fm.second().first() --> moment due to the normal force
fm.second().second() --> moment due to the tangential force
The scope of the two lines you reported is to sum these contributions in order to get the total force and total moment.
|All times are GMT -4. The time now is 06:01.|