I'm looking for a simple approach rather than the BCVTB and Matlab approach used in this paper presented at eSim 2014, since it needs to be built into a larger piece of software without adding any new dependencies.
The idea I had was to model a second tank which is only available to be charged once the first tank is at the melting point of the PCM and so acts like it's storing the latent heat of the phase change. Is this an approach anyone has tried? Am I missing any obvious flaws here? Or is there a better approach you can suggest?