THM/F Correct history variable update

The strain accrued by the solid matrix before onset of freezing is one of the central ingredients of the mass balance equation in the THM+freezing formulation, especially, regarding the way it must be "updated" and "pushed back" within the incremental treatment. It appears that the aforementioned procedures (blocks) implemented in the previous version of code were not entirely correct. This became evident (was detected) when we designed and computed the appropriate benchmark which aimed at modeling a 9% volumetric strain expansion of water during the phase-change. With the amended current version of the code, the simulated results of the benchmark are correct --- the corresponding strains transit from 0 to the reference magnitude of 0.03 during the freezing, as required.

1. Feature description was added to the changelog
2. Tests covering your feature were added?
3. Any new feature or behavior change was documented?