Bouyant Pimple Foam: p_rgh and p increases
 

Dear OpenFoam users,

I am simulating heat transfer through an enclosure filled with non-Newtonian fluid. The viscosity and density of the fluid are considered to be a function of temperature. It is noticed that after several hours of simulation the pressure starts increasing to a higher value in the domain. Please suggest the possible reason for this.

Thermophysical properties are modified to relate density, viscosity, and conductivity with the temperature only.

External heat transfer coefficients are given at the enclosing surface.

pref=101325
pRefPoint = top most point of the enclosure