[Kronecker delta]

Hello all,


I am attempting to model melting by calculating enthalpy in the PCM subdomain. The domain is attached below.


Basically, I am using a UDS for enthalpy and in order to describe its change on the PCM side of the boundary I need the value of heat flux transported through the boundary, so I would need to obtain the value of heat flux on the wall side (or temperature gradients normal to boundary). I have tried using BOUNDARY_HEAT_FLUX macro, but to no avail. As I understand it, I would need to get access to cell and face values of temperature of the adjacent (wall) cells, and I am unable to do so, and I would appreciate your help.



Also, if there is a better/simpler way that I have not considered or if I am wrong about something, please notify me, I would be grateful for any kind of advice you might have.