Parameterize the transition boundaries of multi-phase heat transfer

Right now the band of heat transfer smoothing is fixed at +/- 0.02. I suggest allowing for a parameter that the user can access.