[vabishek]

Hello FOAMers,

I am trying to add a new viscosity model that is dependent on temperature. I was able to successfully compile the model. It works fine with geometries that have just one domain/region. However, when I try extending it to multi-regions, the model fails. In this case, the viscosity gets updated before the first step (based on my model) but it doesn't get updated after that i.e. during each time-step in the simulation.

I have a fluid region surrounded by four solid regions. I have attached a pdf to illustrate this.

I ran a test problem with just one region and printed the viscosity values at each time-step, and it does get updated correctly. I am struggling to see what needs to be modified in order to get this model working with multi-region geometry.

I have attached the code for the viscosity model.

Thanks!