Hello,<div><br></div><div>I am trying simulate a backward facing step flow at Re = 5000 with a passive scalar using Schmidt number = 0.71. I am having some trouble with the boundary condition for the passive scalars.</div>
<div><br></div><div>At the wall I use zero gradient boundary condition by using insulated boundary ("I  ") as the flag in the .rea file. After several flow thorough times, the passive scalars just fill the domain up and become uniform everywhere. This can be observed from the following 2 movies. </div>
<div><br></div><div>1. This is the contour plot of the passive scalar. The "Min" value of the scalar concentration keeps increasing. </div><div><a href="http://www2.uic.edu/~hkanch1/s3.mpeg">http://www2.uic.edu/~hkanch1/s2.mpeg</a></div>
<div><br></div><div>2. In this movie, the "Min" value of scalar is fixed. It can be observed that the scalar concentration is becoming uniform throughout the domain.</div><div><a href="http://www2.uic.edu/~hkanch1/s3.mpeg">http://www2.uic.edu/~hkanch1/s3.mpeg</a>
</div><div><br></div><div>The simulation works just fine if I set the passive scalar value to "zero" on the wall.</div><div><a href="http://www2.uic.edu/~hkanch1/s4.mpeg">http://www2.uic.edu/~hkanch1/s4.mpeg</a>
</div><div><br></div><div>I don't understand why this is happening. May be something wrong in my set up of the problem. </div><div><br></div><div>Regards,</div><div><br></div><div>Harish.</div>