To improve the light-technical properties of HID lamps the use of complex mixtures of metal halides (MH) is increasing. To optimise existing HID lamps and to predict properties of future ML lamps, numerical models can be of great use. Results of such a model developed at Philips Lighting were already reported in [1]. This model solves the energy balance for a cylindrical geometry taking into account chemical equilibrium and species transport by diffusion and radiation. An important shortcoming of this model is that natural convection is not included. For metal halide lamps it is known that natural convection is an important driving force behind radial and axial segregation, which has a considerable influence on the light output properties. ...