Domains are simulated using a new model in which the electron transport and relaxation coefficients depend on the mean energies of the electrons in the valleys, which are deduced from the energy balance within the domain. In GaAs, the accumulation layers are narrowed by cooling effects.