[ejaz8894]

I was working on the icoReactingMultiphaseInterFoam solver & figured out that there is no mass diffusivity term in the species transport equation. I was trying to prescribe mass diffusivity similar to the YEqn.H file in the reactingFoam solver. The governing equation of species transport is not present in the main solver directory, rather it is declared in the src directory: OpenFOAM-v2012/src/phaseSystemModels/multiphaseInter/phasesSystem/phaseModel/MultiComponentPhaseModel.

Original implementation of the governing equation is:

fvScalarMatrix YiEqn
(
fv::EulerDdtScheme<scalar>(mesh).fvmDdt(Yi)
+ fv::gaussConvectionScheme<scalar>
(
mesh,
phi,
upwind<scalar>(mesh, phi)
).fvmDiv(phi, Yi)
==
Su[i] + fvm::Sp(Sp[i], Yi)
);

YiEqn.solve();

I tried to put modification in order to include mass diffusivity similar to the way it is included in the YEqn.H of reactingFoam solver:

fvScalarMatrix YiEqn
(
fv::EulerDdtScheme<scalar>(mesh).fvmDdt(Yi)
+ fv::gaussConvectionScheme<scalar>
(
mesh,
phi,
upwind<scalar>(mesh, phi)
).fvmDiv(phi, Yi)
- fvm::laplacian(turbulence->muEff(),Yi) //ADDED
==
Su[i] + fvm::Sp(Sp[i], Yi)
);

YiEqn.solve();



After hitting the wmake command, following error is showing:

lnInclude/MultiComponentPhaseModel.C:407:32: error: ‘turbulence’ was not declared in this scope
- fvm::laplacian(turbulence->muEff(),Yi)

Can anyone help me, how can I declare this 'turbulence' function? I would also appreciate if anyone can suggest any better approach to include mass diffusivity in species transport equation.