Forces calulated through pressure
Hello my fellow Foamers,
I have to address something I find confusing:
I have found two ways to calculate forces.
The first one is with the help out of this forum. Someone suggested to hang this piece of code to the controlDict
My second way of calculating the forces is by using the command patchIntegrate. If I use this command with the p field and the wall it should (thats what I thought) give me the forces resulting through pressuredifferences.
My problem is that the two methodes come to different solutions. If I run the simpleFoam/airfoil tutorial in timestep 500 I get these two solutions:
controlDict: (-97.6164 755.598 -6.01592e-14)
patchIntegrate: (-81.0095 627.052 -4.99247e-14)
The two solutions are not completly different, but still the question remains:
Why are they different?
Your forces differ all by the factor 1.205 - which is exactly the value of your reference density rhoInf.
If you check out the 2D-Airfoil tutorial you will notice, that the pressure dimension is set to m²/s², which means you actually deal with p/rho (that's at least my assumption).
And so does of course patchIntegrate.
edit: In the transportProperties dictionary rho is set to 1 kg/m³. Changing that to your reference density will also produce matching results.
You are right.
|All times are GMT -4. The time now is 00:11.|