CFD Online Discussion Forums - Using RAS turbulence modelling in multiphaseEulerFoam (OF 2.1.1)

Hello everyone!
From what I've gathered, multiphaseEulerFoam is a LES-based solver. Since I intend to include population balance equations to this solver, I would like to use a less computationally intensive turbulence modelling (RAS). I'm trying to modify the solver so that the turbulence model may be switched from the case directory using turbulenceProperties. I made the following changes:
1. In the file multiphaseEulerFoam.C

Code:

#include "LESModel.H"

Code:

#include "turbulenceModel.H"

and,

Code:

sgsModel->correct();

Code:

turbulence->correct();

2. In the file createFields.H

Code:

autoPtr<incompressible::LESModel> sgsModel

    (

        incompressible::LESModel::New(U, phi, fluid)

    );

Code:

autoPtr<incompressible::turbulenceModel> turbulence

    (

        incompressible::turbulenceModel::New(U, phi, fluid)

    );

3. In the file UEqns.H

Code:

volScalarField nuEff(sgsModel->nut() + iter().nu());

Code:

volScalarField nuEff(turbulence->nut() + iter().nu());

4. In the file Make/options added the following lines:

Code:

-I$(LIB_SRC)/turbulenceModels/incompressible/turbulenceModel \

Code:

-lincompressibleTurbulenceModel \

Is this all that should be changed in the solver? The solver compiles fine after making these changes in OF 2.3.x. In OF 2.1.1, however, I get a lot of errors.