[vkoppejan]

Hi everyone,

I'm playing around with the new reactingTwoPhaseEulerFoam solver.

I'd like to adapt to include laminar diffusivities per species.

After looking around I think I now have found several options for this

1. Specify a Lewis number for each specie.
2. Rewrite the species transp. eqn to include a mueff / Sc number, and define Sc for each specie.
3. Rewrite the species transport eqn to laminar form (as in scalarTransportFoam) and add define diffusivities for each specie.

One of the problems I have is that I can't seem to find the species transport eqn.

In the YEqns.H file it is defined as:

forAll(Y1, i)
{
tmp<fvScalarMatrix> Y1iEqn(phase1.YiEqn(Y1[i]));

if (Y1iEqn.valid())
{
	Y1iEqn =
	(
		Y1iEqn
	 ==
	   *massTransfer[Y1[i].name()]
	  + fvOptions(alpha1, rho1, Y1[i])
	);

	Y1iEqn->relax();
	Y1iEqn->solve(mesh.solver("Yi"));
}
}

So it appears to get the transport terms from somewhere else. Does anybody now where this definition is?

Btw, I know the Le number is correctly used in all the masstransfer models (apart from a few typos), but I really want to include species diffusivity in the Laplacian term.

Thanks for your help.

Victor

[demichie]

Hi Victor,
you can find the equation here:

applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phaseModel/MultiComponentPhaseModel/MultiComponentPhaseModel.C

        fvm::ddt(alpha, rho, Yi)
      + fvm::div(alphaRhoPhi, Yi, "div(" + alphaRhoPhi.name() + ",Yi)")
      - fvm::Sp(this->continuityError(), Yi)

      - fvm::laplacian
        (
            fvc::interpolate(alpha)
           *fvc::interpolate(this->turbulence().nut()*rho/Sc_),
            Yi
        )
     ==
        this->R(Yi)

      + fvc::ddt(residualAlpha_*rho, Yi)
      - fvm::ddt(residualAlpha_*rho, Yi)

[rarnaunot]

Hi Foamers,

I took a look at the transport equation and I'm confused about some terms of it.

Originally Posted by demichie

Hi Victor,
you can find the equation here:

applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phaseModel/MultiComponentPhaseModel/MultiComponentPhaseModel.C

        fvm::ddt(alpha, rho, Yi)
      + fvm::div(alphaRhoPhi, Yi, "div(" + alphaRhoPhi.name() + ",Yi)")
      - fvm::Sp(this->continuityError(), Yi)

      - fvm::laplacian
        (
            fvc::interpolate(alpha)
           *fvc::interpolate(this->turbulence().nut()*rho/Sc_),
            Yi
        )
     ==
        this->R(Yi)

      + fvc::ddt(residualAlpha_*rho, Yi)
      - fvm::ddt(residualAlpha_*rho, Yi)

What's the meaning of these three terms:

1.

      
      - fvm::Sp(this->continuityError(), Yi)

2.

        

      + fvc::ddt(residualAlpha_*rho, Yi)

3.

        

        - fvm::ddt(residualAlpha_*rho, Yi)

Guess that term 2 and 3 are linked.

Thanks,