[vjp]

Hello everyone,


I run a CHT simulation on OF v8 with the chtMultiRegionFoam solver.
I use the postProcess utilities to calculate the wallHeatFlux at a fluid/solid interface:

chtMultiRegionFoam -postProcess -func wallHeatFlux -latestTime -region fluid

When I look at the results on ParaView, there are some discontinuities on the map . For other locations there are not such phenomena (see the attached file, there is a difference between central zone and the upper and lower zone).

I noticed that this comportement is also the same with alphat. Alphat is probably used to calcule the wallHeatFlux.

The boundary conditions are not very far away of this surface where I calculate the wallHeatFlux but they did not match with the repartition of the discontinuities.

So I do not know if it is the boundary conditions or the alpaht file that cause this or anything else.

Perhaps someone has faced the same problem ?
Do you have any idea ?


Thanks in advance,