We studied the formation of a thermal front during the expansion of hot plasma into colder plasma. We used a three-dimensional electromagnetic particle-in-cell model that includes inductive effects. In early phases, in the area of the expanding hot plasma, we found several thermal fronts, which are defined as a sudden decrease of the local electron kinetic energy. The fronts formed a cascade. Thermal fronts with higher temperature contrast were located near plasma density depressions, generated during the hot plasma expansion. The formation of the main thermal front was associated with the return-current process induced by hot electron expansion and electrons backscattered at the front. A part of the hot plasma was trapped by the thermal front while another part, mainly with the most energetic electrons, escaped and generated Langmuir and electromagnetic waves in front of the thermal front, as shown by the dispersion diagrams. Considering all of these processes and those described in the literature, we show that anomalous electric resistivity is produced at the location of the thermal front. Thus, the thermal front can contribute to energy dissipation in the current-carrying loops of solar flares. We estimated the values of such anomalous resistivity in the solar atmosphere together with collisional resistivity and electric fields. We propose that the slowly drifting reverse drift bursts, observed at the beginning of some solar flares, could be signatures of the thermal front.