Possible wrong code implementation for ErgunWenYuDragDragForce

Summary

The code implemented for ErgunWenYuDrag inside src/lagrangian/intermediate/submodels/Kinematic/ParcelForces/Drag/ErgunWenYuDrag/ErgunWenYuDragFroce.C seems different than the source cited in the corresponding header file. There seems to be a '(1 - alphac)' missing in the numerator of the return value of the calcCoupled function, whereas there is an extra 'alphac' in the denominator that is absent in the cited papers.

• OpenFOAM version : v2112