template<class BasicTurbulenceModel>
          twophaseMixingLengthBerzi class
        
        Mixing length model for two-phase flow (only valid for 1D case with Y the wall-normal direction)
References:
"Turbulence Locality and Granularlike Fluid Shear Viscosity in Collisional Suspensions." Berzi, D. & Fraccarollo, L. (2015); Physical Review Letters 115 (19), 194501. "Stresses and Drag in Turbulent Bed Load From Refractive Index-Matched Experiments." Ni, W-J & Capart, H. (2018). Geophysical Research Letters 45 (14), 7000–7009.
The default model coefficients correspond to the following:
    twophaseMixingLengthBerziCoeffs
    {
        Cmu         0.09;
        expoLM      1.0;
        alphaMaxLM  0.55;
        kappaLM     0.225;
    }
This file is in the