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.