CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

Pollutant dispersion in the environment

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree1Likes

LinkBack Thread Tools Display Modes
Old   October 7, 2005, 16:47
Default I'm not sure. It's not my fiel
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,911
Rep Power: 28
alberto will become famous soon enoughalberto will become famous soon enough
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):

Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status:
OpenQBMM - An open-source implementation of quadrature-based moment methods.

To obtain more accurate answers, please specify the version of OpenFOAM you are using.
alberto is offline   Reply With Quote

Old   January 29, 2008, 21:27
Default Hi all, So to revisit the q
New Member
Join Date: Mar 2009
Posts: 16
Rep Power: 10
evan is on a distinguished road
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

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



Or is there a better way about going about this? Mattijs' suggestion from before is currently over my head .

evan is offline   Reply With Quote

Old   May 23, 2008, 05:09
Default I have just needed to do this,
Andrew King
Join Date: Mar 2009
Location: Perth, Western Australia, Australia
Posts: 82
Rep Power: 10
andersking is on a distinguished road
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(
"S" + Foam::name(si),


then to solve the (laminar) scalar transport

+ fvm::div(phi, S[si])
- fvm::laplacian(DS[si], S[si])


(corrections welcome)

rob3rt 0ng likes this.
Dr Andrew King
Fluid Dynamics Research Group
Curtin University
andersking is offline   Reply With Quote

Old   September 21, 2010, 03:49
Default Linking species properties with species equation and spacies table
Hrushikesh Khadamkar
Join Date: Jul 2010
Location: Mumbai
Posts: 66
Rep Power: 9
Hrushi is on a distinguished road
Originally Posted by alberto View Post
OK. But hook() just adds an element at a time to the ptrList.

What I'd like to do is to let the user set the number of species, their name and then insert the diffusivities.
The speciesTable manages the names of the species, and I was thinking to store the ordered list of diffusivities in a list like


in the speciesProperties dictionary.

Now, if I use the ptrList, the problem of reading the list remains.

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?

Hrushi is offline   Reply With Quote

Old   November 29, 2014, 14:31
Join Date: Nov 2012
Posts: 30
Rep Power: 7
azna is on a distinguished road

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 ?
azna is offline   Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
pollutant dispersion Andy CFX 1 August 26, 2007 18:41
Simulation of pollutant dispersion in a city George Siemens 1 June 27, 2007 04:26
pollutant source J.Z Phoenics 2 April 3, 2007 18:33
pollutant source Riccardo Buccolieri FLUENT 0 January 22, 2005 04:38
pollutant source Riccardo Buccolieri FLUENT 0 January 20, 2005 05:06

All times are GMT -4. The time now is 23:22.