[JoeFriend]

This last month I have tried to extend RichardsFoam2 functionality by adding the heat transport, but with no success.



I have added the convective-diffusive equation in the solver as it follows

 fvm::ddt(rhoCpe,T) == fvm::laplacian(Ke,T) - rhoCpf*fvm::div(phi,T)

where rhoCpe is the effective volumetric heat capacity (solid + gas + water), Ke is the effective thermal conductivity and rhoCpf is the effective volumetric heat capacity of the fluid phases (gas + water).


For some reason the convective term is making my test cases lose a lot of energy, which makes no sense to me when I try to validate them.


I have attached a file with the solver's main code. I would appreciate if anyone could help me