I'm not sure. It's not my fiel
 

I'm not sure. It's not my field. But you could give a look to the papers on dispersion in the environment where flacs calculations are presented.
They validated their models on large scale data and experiments.

You can find the manual here (you need to subscribe for free):

http://www.gexcon.com/index.php?src=...alidation.html

Regards,
Alberto

Hi all, So to revisit the q
 

Hi all,

So to revisit the question raised by Alberto, does anyone know (or be willing to discuss) how to make an indexed list of constants that is read from a dictionary. So, the diffusivity example Alberto brought up would work great. That is, say you have n diffusion constants listed in your constants dictionary as:

dimensionedScalar D1 (dictionary.lookup("D1"));
.
.
.
dimensionedScalar Dn (dictionary.lookup("DY"));

but you want to index them in your solver (along with other indexed fields) as:

for(label i=0; i<n; i++)
{
volScalarField& Yi = Y[i];

// call D[i] list here

solve
(
fvm::ddt(Yi)
+ fvm::D[i]*fvm::laplacian(Yi)

);

etc.

Or is there a better way about going about this? Mattijs' suggestion from before is currently over my head http://www.cfd-online.com/OpenFOAM_D...part/happy.gif.

Best,
Evan

I have just needed to do this,
 

I have just needed to do this, so I've posted what I have used.

for creating the scalar fields:

Info<< "Reading number of scalars...";

scalar nScalar(readScalar(transportProperties.lookup("nSc alar")));
Info<< nScalar <<endl;
PtrList<dimensionedscalar> DS(nScalar);
PtrList<volscalarfield> S(nScalar);

forAll(S, si)
{
Info<< "Reading diffusivity DS" << si << endl;

dimensionedScalar DSi(transportProperties.lookup("DS" + Foam::name(si)));
DS.set(si,new dimensionedScalar(DSi));

Info<< "Reading field S" << si << endl;

S.set(si,new volScalarField(
IOobject
(
"S" + Foam::name(si),
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
));

}

then to solve the (laminar) scalar transport

forAll(S,si)
{
solve
(
fvm::ddt(S[si])
+ fvm::div(phi, S[si])
- fvm::laplacian(DS[si], S[si])
);

}
}

(corrections welcome)

Regards,
Andrew

Linking species properties with species equation and spacies table
 

Hi Alberto,

I am also trying to develop steady state liquid phase multispecies
solver. Can you elaborate more on how to link individual species properties in the species equation?

Hrushikesh

Hi,

I am working on a multiphase flow ( air-water) project in Fluent. I would like to assume that at time=0 (initial condition), nitrogen exists in the water ( without any inlet). I was wondering how can I define water+nitrogen, for e.g. 50mg/l nitrogen exists in water.?

is there any section I can define concentration of nitrogen is water ?