Pulsar wind nebulae (PWN) are fascinating systems and archetypal sources for high-energy astrophysics in general. Due to their vicinity, brightness, to the fact that they shine at multi-wavelengths, and especially to their long-living emission at gamma rays, modelling their properties is particularly important for the correct interpretation of the visible Galaxy. A complication in this respect is the variety of properties and morphologies they show at different ages. Here, we discuss the differences among the evolutionary phases of PWN, how they have been modeled in the past and what progresses have been recently made. We approach the discussion from a phenomenological, theoretical (especially numerical) and observational point of view, with particular attention to the most recent results and open questions about the physics of such intriguing sources.
