It has been observed previously that the physical behaviors of Schmidt number (Sc) and Prandtl number (Pr) of an energy-conserving dissipative particle dynamics (eDPD) fluid can be reproduced by the temperature-dependent weight function appearing in the dissipative force term. In this paper, we proposed a simple and systematic method to develop the temperature-dependent weight function in order to better reproduce the physical fluid properties. The method was then used to study a variety of phase-change problems involving solidification. The concept of the "mushy" eDPD particle was introduced in order to better capture the temperature profile in the vicinity of the solid-liquid interface, particularly for the case involving high thermal conductivity ratio. Meanwhile, a way to implement the constant temperature boundary condition at the wall was presented. The numerical solutions of one- and two-dimensional solidification problems were then compared with the analytical solutions and/or experimental results and the agreements were promising.