آلودگی خاک به مواد نفتی از چالش های مهم عصر کنونی است. به منظور بررسی کمی ویژگی های هیدرولیکی خاک های آلوده به مواد هیدروکربنی، منحنی های رطوبتی بنزین، گازوئیل، نفت سفید، نفت خام و آب در چهار نوع محیط متخلخل خاک شامل بافت های sandy loam، silt loam، silty clay و شن خلیج فارس در شرایط دو فازی هواnapl- بدست آمد. همچنین، منحنی های رطوبتی این خاک ها در شرایط سه فازی آب-هوا-napl در سه نسبت هیدروکربن-آب شامل 25-75، 50-50 و 75-25 تعیین شد. مواد هیدروکربنی مورد مطالعه شامل بنزین، گازوئیل، نفت سفید و نفت خام بود. هدایت هیدرولیکی اشباع هر یک از خاک ها برای بنزین، گازوئیل، نفت سفید و نفت خام و آب نیز به روش بار ثابت بدست آمد. سپس پارامترهای مربوط به منحنی های رطوبتی در شرایط دوفازی و سه فازی بر اساس مدل های ون گنوختن، بروکس و کوری و کمپبل بدست آمد. همچنین منحنی های هدایت هیدرولیکی غیر اشباع حالت های دو فازی آلاینده های بنزین، گازوئیل، نفت سفید و نفت خام و آب توسط مدل های ون گنوختن-معلم، بروکس و کوری-معلم، بروکس و کوری-بوردین و کمپبل بدست آمد. نتایج نشان داد در هر دو حالت دو فازی و سه فازی پارامتر مکش ورود هوا در توابع هیدرولیکی نسبت به پارامتر دامنه توزیع خلل و فرج حساسیت بیشتری به تغییرات ویژگی های فیزیکی سیالات دارد. در اغلب موارد، مکش ورود هوا برای نفت خام نسبت به آب افزایش و برای سه آلاینده ی دیگرکاهش یافت. همچنین پارامتر بیان کننده ی توزیع خلل و فرج، در اغلب موارد برای بنزین بیشترین و برای نفت خام کمترین مقدار را داشت. مقدار این پارامتر برای گازوئیل و نفت سفید نیز نسبت به آب کمتر بود. هدایت هیدرولیکی اشباع به ترتیب برای بنزین، نفت سفید و گازوئیل بیشتر و برای نفت خام کمتر از آب بود. اما در مکش های نزدیک به اشباع هدایت هیدولیکی نفت سفید، بنزین و گازوئیل کمتر از آب بود. هدایت هیدولیکی نفت خام اغلب در مکش های نزدیک به یک اتمسفر بیش از سه napl دیگر بود. تقریبا در تمامی سیستم های سه فازی، منحنی های مشخصه napl-آب بالاتر از منحنی آب و ترکیب 75-25 napl-آب بالاتر از دو ترکیب دیگرقرار گرفت. در بیشترمحیط ها، ترکیبات ماده آلاینده نفت خام، منحنی رطوبتی بالاتری نسبت به بقیه ترکیبات نشان داد. در اغلب سیستم های سه فازی، پارامتر بیان کننده توزیع تخلخل کاهش و مکش ورود هوا افزایش یافت. منحنی های مشخصه سه فازی نفت خام اغلب بیشترین تغییر را نشان داد. پس از نفت خام، گازوئیل و سپس نفت چراغ و بنزین قرار گرفتند. با این حال، رفتار محیط های سه فازی حاوی بنزین و نفت سفید شبیه به یکدیگر بود. در بسیاری از موارد، با کاهش مقدار ماده هیدروکربنی در ترکیبات از 75 به 25 درصد، مقادیر پارامتر شاخص توزیع تخلخل افزایش و مکش ورود هوا کاهش نشان داد و به مقادیر این پارامتر ها در منحنی مشخصه آب نزدیک شدند. نتایج این پژوهش نشان می دهد که خاک های خشک توانایی بهتری برای عبور آلاینده های هیدروکربنی از محیط غیر اشباع به آب های زیرزمینی دارند. همچنین بنزین به دلیل حلالیت در آب و هدایت هیدرولیکی زیاد محیط های متخلخل برای این سیال، به مراتب آلاینده ی خطرناکتری برای سلامتی آب های زیرزمینی محسوب می شود.

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