[sreeyuth]

Dear All,

I am simulating passive scalar transport from the bottom of a porous material which is connected to a free wind flow domain (Setup.png). I give a scalar value of T = 0.1 at the bottom of the porous material and then would like to see how it is diffused and advected by the wind. In scalar_source.png, part of the area which acts as a scalar source can be clearly seen. All the other cells have an initial scalar concentration of 0.

My DT value is 0.000025 m2/s. I give deltaT=0.0001. As is the requirement of scalartransportFoam, I have already solved the velocity field and it is given as an input (mean wind speed is 1.5 m/s). The scalartransportFoam solver has been modified to account for turbulent diffusivity. I successfully used this very same solver to calculate scalar transfer from a flat plate.

However, when I simulate scalar transport from porous media as shown in Setup.png, my solution starts diverging around 2.5 seconds. When I first found out this problem, my scalar value was T=10. Upon decreasing the value to T=0.1, the situation still persists. According to checkMesh, the mesh is OK.

Can anyone tell me where I could be wrong? Please let me know if you need any additional data.

Best regards,
Sreeyuth.