A two-dimensional shallow layer model has been developed to predict dense gas dispersion, under realistic conditions, including complex features such as two-phase releases, obstacles and inclined ground. The model attempts to predict the time and space evolution of the cloud formed after a release of a two-phase pollutant into the atmosphere. The air-pollutant mixture is assumed ideal. The cloud evolution is described mathematically through the Cartesian, two-dimensional, shallow layer conservation equations for mixture mass, mixture momentum in two horizontal directions, total pollutant mass fraction (vapor and liquid) and mixture internal energy. Liquid mass fraction is obtained assuming phase equilibrium. Account is taken in the conservation equations for liquid slip and eventual liquid rainout through the ground. Entrainment of ambient air is modeled via an entrainment velocity model, which takes into account the effects of ground friction, ground heat transfer and relative motion between cloud and surrounding atmosphere. The model additionally accounts for thin obstacles effects in three ways. First a stepwise description of the obstacle is generated, following the grid cell faces, taking into account the corresponding area blockage. Then obstacle drag on the passing cloud is modeled by adding flow resistance terms in the momentum equations. Finally the effect of extra vorticity generation and entrainment enhancement behind obstacles is modeled by adding locally into the entrainment formula without obstacles, a characteristic velocity scale defined from the obstacle pressure drop and the local cloud height. The present model predictions have been compared against theoretical results for constant volume and constant flux gravity currents. It was found that deviations of the predicted cloud footprint area change with time from the theoretical were acceptably small, if one models the frictional forces between cloud and ambient air, neglecting the Richardson dependence. The present model has also been validated in widely different experimental conditions such as the Thorney Island instantaneous isothermal releases 8 (unobstructed) and 21 (with semicircular fence), the EEC-55 two-phase propane experiment (with and without linear fence), the Desert Tortoise 4 two-phase ammonia experiment and the Hamburg DAT-638 instantaneous inclined plate experiment and the model predictions were found in reasonable agreement with the experimental data. (C) 2003 Elsevier Science B.V. All rights reserved.