Adjustment of nodal temperatures after phase change
After phase changes, the change in element temperature needs to be distributed over the nodes. Currently, this is done by summing the change in nodal temperature to the element temperature: Element[i].Te += delta Node[i].T; Element.Te[i] += delta Node[i+1].T;
However, I argue that this is only correct for the highest and lowest node. The other nodes are connected two 2 elements and every element has 2 nodes. So the temperature should only be adjusted half: Element[i].Te += 0.5* delta Node[i].T; Element[i].Te += 0.5* delta Node[i+1].T;
It can be easily seen by a though experiment where all the nodes (and thus elements) are initially 274K. Using my proposed approach, all elements and nodes will be 273K after phase changes, whereas in the current approach, only the highest and lowest node are 273K, all the other nodes are 272K.