در این پروژه به مطالعه عددی پخش و صعود ستون دود و توزیع غلظت ماده آلاینده در سطح زمین در حالت سه بعدی و غیر دائم پرداخته شده است. هم چنین نحوه گسترش ستون دود با زمان و مکان نیز بررسی شده است. بدلیل اختلاف دمای حدود 100 درجه سانتیگراد بین گازهای خروجی از دودکش و هوا نیروی شناوری از طریق وارد کردن شتاب گرانش به مساله در نظر گرفته شده است. مدلهای توربولانسی مورد استفاده و هستند. شرایط اتمسفر به صورت تغییر سرعت باد و دما با ارتفاع در ورودی جریان باد اعمال شده اند. از برهم کنش بین گازهای خروجی و هوا صرف نظر شده است. برای کاهش زمان محاسبات و بدلیل تقارن هندسه مساله، دامنه حل نصف شده است. نتایج بدست آمده برای پخش ستون دود ناشی از حل عددی با مدل نیمه تجربی گوسی پخش ستون دود مقایسه و تفاوت اندکی مشاهده شد. نتایج صعود ستون دود نیز با روابط نیمه تجربی مقایسه گردیده اند که تطابق خوبی با هم دارند. شبیه سازی عددی صعود ستون دود را اندکی کمتر از مقادیر نیمه تجربی پیش بینی می کند. بنابراین مقدار غلظت آلودگی در سطح زمین در شبیه سازی عددی اندکی بیشتر از مدل نیمه تجربی گوسی بدست آمد. میزان صعود دود در حالت اتمسفر پایدار و ناپایدار به ترتیب کمترین و بیشترین می باشد. بنابراین غلظت آلودگی در سطح زمین در حالت اتمسفر پایدار بیشترین مقدار و در حالت اتمسفر ناپایدار کمترین مقدار است.