CFD Online Discussion Forums - twoPhaseEulerFoam: "g0" missing in "HrenyaSinclairConductivity.C"

Hello!

In "$FOAM_SOLVERS/multiphase/twoPhaseEulerFoam/kineticTheoryModels/conductivityModel/HrenyaSinclair/HrenyaSinclairConductivity.C", it seems that "g0" is missing from the second term of "kappa":

Code:

    return rhoa*da*sqrt(Theta)*

    (

        2.0*sqr(alpha)*g0*(1.0 + e)/sqrtPi

      + (9.0/8.0)*sqrtPi*0.25*sqr(1.0 + e)*(2.0*e - 1.0)*sqr(alpha)

       /(49.0/16.0 - 33.0*e/16.0)

      + (15.0/16.0)*sqrtPi*alpha*(0.5*sqr(e) + 0.25*e - 0.75 + lamda)

       /((49.0/16.0 - 33.0*e/16.0)*lamda)

      + (25.0/64.0)*sqrtPi

       /((1.0 + e)*(49.0/16.0 - 33.0*e/16.0)*lamda*g0)

    );

Instead, the second term should be:

Code:

      + (9.0/8.0)*sqrtPi*g0*0.25*sqr(1.0 + e)*(2.0*e - 1.0)*sqr(alpha)

       /(49.0/16.0 - 33.0*e/16.0)

I use 2.2.x, but 1.6-ext has the same issue. If someone can confirm on this then we should probably report it as a bug to be fixed. Thanks!

References:
(1) Hrenya, C. M.; Sinclair, J. L. Effects of Particle-phase Turbulence in Gas-solid Flows. AIChE Journal 1997, 43, 853–869.
(2) Van Wachem, B. Derivation, Implementation, and Validation of Computer Simulation Models for Gas-solid Fluidized Beds, Delft University of Technology, 2000, p. 49.