This paper deals with the interaction of twin elliptic jets with a cooler oncoming crossflow. The jet nozzles are placed tandem with the main flows, three diameters one from another and initially inclined with a 60° angle. The exploration of the resulting flowfield was carried out numerically by means of the finite volume method together with the RSM (Reynolds Stress Model) second order turbulent closure model and non uniform grid system that was particularly refined near the nozzles. After validation with reference to PIV (particle image velocimetry) experimental data, the model was upgraded by discharging a non reactive fume from the jet nozzles and introducing a variable temperature gradient between the interacting flows. We focused mainly in the present work on the determinant role of the temperature gradient on the dispersion of the discharged pollutants within the domain. Such a study is likely to optimize (control, reduce, eliminate, etc.) one of the most alarming nowadays’ environmental concern: the atmospheric pollution.