[asharma]

Hello,

I'm trying to replicate wind tunnel experiments done for pollutant dispersion inside isolated street canyon. Most of the work is based on the set up published by Meroney et al (1996). I'm having trouble (physics) incorporating the source for the scaler pollutant inside my domain. I've modified pisoFoam to incorporate scaler transport as

 fvm::ddt(T)
	  + fvm::div(phi, T)
	  - fvm::laplacian(turbulence->nut()+DT, T)
	  ==
          source

I've also incorporated nut() value to account for turbulent diffusion. Presently i'm using funkysetfields ( part of swak4foam) to specify source as a constant value in cells where the pollutant is introduced and thus is fixed throughout the run. However my prof tells me this is wrong, as it doesn't fit the physics of the problem. Rather i should force the value of my scaler T to be constant at those points rather than specify source. Is this correct? If so how do i go about this?....

Regards,
Abhinav

[alberto]

You can use setValues to do what you ask for. Please, search on the forum to see how to use it.

[wayne14]

Originally Posted by asharma

Hello,

I'm trying to replicate wind tunnel experiments done for pollutant dispersion inside isolated street canyon. Most of the work is based on the set up published by Meroney et al (1996). I'm having trouble (physics) incorporating the source for the scaler pollutant inside my domain. I've modified pisoFoam to incorporate scaler transport as

 fvm::ddt(T)
	  + fvm::div(phi, T)
	  - fvm::laplacian(turbulence->nut()+DT, T)
	  ==
          source

I've also incorporated nut() value to account for turbulent diffusion. Presently i'm using funkysetfields ( part of swak4foam) to specify source as a constant value in cells where the pollutant is introduced and thus is fixed throughout the run. However my prof tells me this is wrong, as it doesn't fit the physics of the problem. Rather i should force the value of my scaler T to be constant at those points rather than specify source. Is this correct? If so how do i go about this?....

Regards,
Abhinav

Hi, I am doing something similar. Could you explain why using funkysetfields doesn't fit the physics?

Thank you!