Adding temperature equation to twoPhaseEulerFoam
 

I added temperature equation for both the phases. Then I ran "bed" case which by default has a vertical channel with inlet at bottom. It ran properly. I wanted to see flow profile and temperature variation in a horizontal channel. For that 1. I changed "g" vector value so that now it is normal to right wall. 2. I made necessary changes in blockMeshDict to get a horizontal channel. Both time I am getting a diverging temperature field. I am not understanding why code is not working for horizontal channel but for vertical channel? Can someone help me? I am working on OF2.2.2