Hello,

I have a problem with my pressure drag. I calculate the total drag just by "force_x()@WALL" and the friction drag on my body by "areaInt(Wall Shear)@WALL". That both seems to work really well, at least the numbers make sense. However when calulating the pressure drag by "areaInt(Absolute Pressure)@WALL" or "areaInt(Pressure)@WALL" the numbers do not make sense at all. I know that actually I can get the pressure drag just by subtracting but I have intended to make sure that all three calculations fit together. Only in order to double-check.

Has anyone ever had s similar problem or just knows what I have done wrong?

Regards Andy

Here's how i've done it in the past:

X-Viscous force = areaInt(Wall Shear X)@body

Y-Viscous force = areaInt(Wall Shear Y)@body

Z-Viscous force = areaInt(Wall Shear Z)@body

X-Pressure force = areaInt(Pressure*Normal X)@body

Y-Pressure force = areaInt(Pressure*Normal Y)@body

Z-Pressure force = areaInt(Pressure*Normal Z)@body

Out of interest i have checked the CEL computed values with those given in the .out file and guess what, there is a maximum of 10% difference in the values espically in the X component of pressure force (with the CEL giving the smallest value)!

Hope this helps.

James

Dear Andy If you intend to find pressure drop between two locations in the flow path then you can simply calculate the average pressure (use total pressure if two cross-sections are not same) at desired locations and substract as you know. But if you want to find the force on the wall by using areaInt(Absolute Pressure)@WALL" or "areaInt(Pressure)@WALL you will only get the normal force (pressure) acting on the wall not the drag force,as the pressure always acts normal to wall whereas drag is shear force. Bye Ram Dayal

Thanks a lot, you guys have helped me a lot!

@James: I have compared the CEL-value with the given value from the out-file. And I got the similar problem. My CEL-Value for the x-pressure component (calculated with your expression) is about 6% too high. However when adding the out-file-value for the x-pressure to the CEL-value for the x-viscous force the sum is pretty exactly the same as for the CEL-value of the total x-drag [force_x()@body]. That is really strange. And it is even more weird that they show you the values in the out-file but do not provide them automatically in the POST-file. I think there is definitely some space for improvement in user friendliness... ;-) Anyway, thanks a lot and have a nice day!

Yeah, Andy I found exactly the same thing. The viscous drag is pretty much identical in the .out and the CEL, but the pressure drag is quite a bit different! Perhaps there is some rounding error occurring in the CEL or .out routine.

Dear Andy,

Are your boundaries flat? areaInt()@Locator does not do area projection. Pressure is a scalar.

However, the force_x/y/z calls now exactly what is being integrated.. Perhaps that is the difference?

Opaque