Pollutant dispersion over a specified terrain
 

Hi,

I am very new to OpenFOAM and CFD (few weeks). I would like to study pollutant dispersion over a specified terrain.

Currently, I barely managed to set up the initial flow boundary conditions which are mostly based on OpenFOAM tutorial "turbinesiting".

I have no idea as to how to include pollution into the picture. I read a few threads which refer me to the tutorial ScalarTransportFoam but I am still totally lost. Can anyone help explain as to how can dispersion study be built on the existing model (turbinesiting)?

you should solve a transport equation beside p and U to simulate dispersion of pollutant. Therefore you need some programming to add this transport equation to the basic solver that is simpleFoam or pimpleFoam or pisoFoam.

Hi, thanks for the reply. But I would really appreciate it if you could put it into a simpler term. I am really new to this.

As of current, I am using a k-epsilon turbulent model and within the turbinesiting/constant folder, there is already a transportProperties file which defined a Newtonian transportModel (do not know what does this do).

Thank you very much in advance

well, as i said before you need programming :)
now you only solve continuity equation and momentum equations which return U and p.
for simulation of pollutant dispersion, you need to solve extra transport equation, look this.

this equation should be added to your basic solver. so you should modify and extend the existing solver. to learn how to add extra equation to existing solver, read this

WOW!!! that is a very comprehensive read.
I did further reading and came across reactingFOAM, will it do the trick as well. By using it, I could avoid all the programming which are very complicated for me.
Another thing, after defining the solver, how can the pollution be introduced?

And THANK YOU SO MUCH for entertaining my noobish question. feeling totally overwhelmed with this CFD and OpenFOAM.

Dear Nimasam,

I have read your post writing about programming,
May you help to modify for me any simple solver including attached c++ code of the solid motion.

I want it to solve for x_t(the displacement of a point).

I have been trying following the link you share of adding temperature in icoFoam but failed.


I would appreciate you support.