[Dicanic]

Hi all,



I have an interesting case where I have a temperature dependent heat source. The heat source is dependent on the average temperature of a cellzone or set of cells and the power source change is significant and depends on the last time (dt) and power as well.



I tried using fvOptions to simply take the average and write it out to a text file to read in on the next timestep, However, when I run fvOptions in parallel, I realize it is only doing this for the average of the partitioned mesh and not the whole set of cells.



So here are a few questions that I have:



Can I reference cells in fvOptions that aren't in the decomposed section?



What order are functions (which I believe I can reference the whole mesh or just cellzones) executed, before or after fvOptions? Also any reference or guide on how to write a function would be very helpful.


Worse comes to worse, can I write a custom solver that would accomplish this and have an input dictionary specify the cellzone/cellset of interest?



I am able to work with any version of OpenFOAM, but my first few attempts were on of8. If there is an easier way to do this on of-2106 or of9, I would like to try it. Any help is appreciated, thank you!