[chefkoch89]

Dear Community,

this is my first thread and I'm still new with OpenFOAM.
I'm working on multiphase solver wich is based on InterFoam with some enhancements.
Right now I want to implement a porosity criteria funtion based on Darcys law into the code.

For these computations I need a loop, which is checking, if the temperatur in each cell and the cells within ten times the diameter of that cell are in a critical stage of solidification (this value is given by a temperature range).

My problem is, that I don't know how I can identify one cell and check if the temperatures of the surrounding cells are below the given temperature range. Finally I want to add their volumes to a new value.

I hope you are able to understand the problem. I image that the needed code is quite simple and not a big problem

Thank you!