The air temperature variation of a closed room, well insulated, during the initial time of operation of air-conditioning systems up to temperature stabilization, is simulated by a two-dimensional integral model as a quasi-steady-state phenomenon. The model equipped with a conservation equation for tracer concentration or relative temperature, including the stratification parameter, is well qualified. The flow leaving the air conditioning device forms an inclined buoyant jet which bends over and meets the room floor, where it spreads sideways forming a layer with jet temperature. A sequence of layers, which affect the jet temperature through entrainment, are produced by a novel bottom-up technique. The layer air temperatures are calculated through the bulk dilution of a near bottom jet cross-section, which feeds each new layer. The model simulated a real case and predicted the transient variation of room air and buoyant jet temperatures up to stabilisation. It also predicted the time needed for stabilisation, the cooling rates of the room and jet air temperatures, the Brunt-Väisälä frequency occurring during the temperature transitions, and more. The results are promising as they agree with observations. Thus, the model could be used to evaluate the effectiveness of relevant HVAC systems operating in such rooms.