CFD中文网

return neg(Re - 1000)*24*(1.0 + 0.15*pow(Re, 0.687)) + pos0(Re - 1000)*0.44*max(Re, residualRe_);

https://www.openfoam.com/documentation/guides/v2206/api/reactingEuler_2multiphaseSystem_2interfacialModels_2dragModels_2SchillerNaumann_2SchillerNaumann_8C_source.html

结合SchillerNaumann阻力模型的公式，你看一下，公式我笔记里面有

算不准大概率是网格问题，不仅仅是数量，还有质量，另外，系数值的获取是如何计算，自己手动计算还是直接通过软件监测？若通过软件监测，则需注意参考值的设定。