Skip to content

Wrong CdRe in sphereDragForce.C

Hello, While investigating the particleFoam algorithm I notice a possible mistake.

The calculation of CdRe in sphereDragForce.C reads:

" if (Re > 1000.0) { return 0.424Re; } else { return 24.0(1.0 + 1.0/6.0*pow(Re, 2.0/3.0)); } "

But according to the paper: Liu, A.B., Mather, D., and Reitz, R.D. 1993. “Modeling the Effects of Drop Drag and Breakup on Fuel Sprays”, SAE Paper 930072. it should be: if (Re > 1000.0) { return 0.424; } else { return 24.0/Re*(1.0 + 1.0/6.0*pow(Re, 2.0/3.0)); }

Greetings,

Maxim V.C.