CFD Online Discussion Forums - Activ Scalar transport equation

Hello everyone,

I'm trying to add an activ scalar equation on the solver buoyantSimpleFoam.

The goal :
- 2 species involved (1 standard + 1 pollutant)
- density varying on temperature (large scale of T)
- turbulent flow

So I guess I need to declare 2 differents thermo in createFields.H :

Code:

autoPtr<rhoThermo> pThermo1(rhoThermo::New(mesh));

rhoThermo& thermo1 = pThermo1();

thermo1.validate(args.executable(), "h", "e");



volScalarField rho1

(IOobject

(

...

),

thermo1.rho1()

);



autoPtr<rhoThermo> pThermo2(rhoThermo::New(mesh));

rhoThermo& thermo2 = pThermo2();

thermo2.validate(args.executable(), "h", "e");



volScalarField rho2

(IOobject

(

...

),

thermo2.rho2()

);



// scalar definition

volScalarField alpha2

(IOobject

(

...

),

mesh

);



volScalarField rho

(IOobject

(

...

),

alpha2*rho2+ (1-alpha2)*rho1

);

Does it look like the right way to code this kind of solver ?

Thank you for your help.

NB : For information, I already performed this for a passiv scalar (density of scalar = density of fluid).