@Paolo
Its up to you to accept the facts or ignore

All I can say..
Good Luck

Hi Andrew,

I just thought that in principle it is possible to modify everything in OpenFOAM. I do not know, maybe I am wrong.
Anyway, many thanks for the advice!

Cheers,
Paolo

Dear Paolo,

it has been a while since I dealt with this problem. Here, what I found out a while ago for the LUDS scheme

The weights for the surface interpolation are computed in the method linearUpwind.C
and used in the method surfaceInterpolationSchemes.C in
src/finiteVolume/interpolation/surfaceInterpolation/surfaceInterpolationScheme

In linearUpwind.C the method upwind.H is called from
src/finiteVolume/interpolation/surfaceInterpolation/limitedSchemes/upwind In upwind.H the weights are computed according to the UDS scheme, setting  = 1 if the flux is greater
than zero and  = 0 else.

The explicit correction via the gradient
is added in the method linearUpwind.C in
src/finiteVolume/interpolation/surfaceInterpolation/schemes/linearUpwind

I hope this helps?

Kind Regards
Anne

Hi Anne,

Thank you very much for your answer, for your time, and sorry for the belated reply.

Your reply helped me a lot. However, I am still trying to understand how to get the scheme I want, but seems that it will be as much harder as expected. At least, by considering what you said, now it seems that it would be possible... I just need time to understand the structure of the software.

Thanks once again!

Kind Regards,
Paolo

laplacian scheme
 

can anyone suggest which one scheme is to be choosed in case of enthalpy updating method for phase change material.