Realistic modeling of residential users is important to achieve accurate and reliable results in harmonic studies of public LV grids. This modeling is complex, because multiple factors influence the emission and propagation of harmonics through the network, like time-variation of number and types of connected devices, the supply voltage distortion and the network impedance. This paper discusses the advantages and disadvantages of state of the art modeling methodologies and it describes a novel modeling approach based on measurement of harmonic emission of individual residential users and aggregates of residential users. The developed model can represent the stochastic variation of harmonic currents and is useful e.g. to analyze the impact of new equipment technologies (like EV chargers) on the summation and propagation of harmonic currents in public LV grids.