دوره 16، شماره 2 - ( خرداد و تیر 1396 )
جلد 16 شماره 2 صفحات 97-88 |
برگشت به فهرست نسخه ها
Download citation:
BibTeX | RIS | EndNote | Medlars | ProCite | Reference Manager | RefWorks
Send citation to:
BibTeX | RIS | EndNote | Medlars | ProCite | Reference Manager | RefWorks
Send citation to:
Danaee S, Fazaelipoor M H, Salmani M H. Using of Moisture Term in the Dynamic Modeling in Treatment Process of Hexane-Contaminated Air. TB 2017; 16 (2) :88-97
URL: http://tbj.ssu.ac.ir/article-1-2331-fa.html
URL: http://tbj.ssu.ac.ir/article-1-2331-fa.html
دانایی سروش، فضایلی پور محمد حسن، سلمانی محمدحسین. به‌کارگیری ترم رطوبت در مدل‌سازی دینامیک فرآیند تصفیه هوای آلوده به هگزان . طلوع بهداشت. 1396; 16 (2) :88-97
دانشگاه علوم پزشکی شهید صدوقی یزد ، mhsn06@yahoo.com
متن کامل [PDF 251 kb]
(1853 دریافت)
| چکیده (HTML) (3967 مشاهده)
متن کامل: (1156 مشاهده)
به‌کارگیری ترم رطوبت در مدل‌سازی دینامیک فرآیند تصفیه هوای آلوده به هگزان
نویسندگان:سروش دانایی1،محمدحسن فضایلی پور2، محمدحسین سلمانی3
1.باشگاه پژوهشگران جوان و نخبگان، واحد شهریار، دانشگاه آزاد اسلامی، تهران
2.دانشیار مهندسی شیمی، گروه مهندسی شیمی و پلیمر پردیس فنی مهندسی، دانشگاه یزد
3.نویسنده مسئول: دکتری شیمی تجزیه، دانشکده بهداشت، دانشگاه علوم پزشکی و خدمات بهداشتی درمانی شهید صدوقی یزد تلفن تماس: 09131582878 Email:mhsn06@ssu.ac.ir
چکیده
مقدمه: یک روش مؤثر در تصفیه بیولوژیکی هوای آلوده به ترکیبات فرار آلی روش بیوفیلتراسیون است که با عبور هوا از درون بستر پرشده، ماده آلاینده از فاز گاز به فاز بیوفیلم منتقل‌شده و توسط میکروارگانیسمها تجزیه میگردد. فعالیت میکروب‌های تجزیهکننده به‌شدت تحت تأثیر رطوبت بستر است که مهم‌ترین فاکتور برای کنترل عملکرد بهینه بیوفیلتراسیون به‌حساب میآید. هدف از این مطالعه، به‌کارگیری ترم رطوبت برای توصیف عملکرد مدلی کاربردی در بیوفیلتر ناپایا برای حذف هگزان از جریان هوای آلوده است.
روش بررسی: ابتدا یک مدل دوبعدی از بیوفیلتراسیون هوای آلوده به هگزان ارائه شد که در آن ترمهای رطوبت بستر و پراکندگی محوری در موازنه انتقال جرم در نظر گرفته شدند. اثرات استفاده از سوپر جاذب آب نظیر افزایش مداوم تخلخل بستر و سطح ویژه آن درون معادلات وارد گردید. سپس حل عددی معادلات با روش حجم محدود انجام شد و در انتها، آزمایشهای بیوفیلتراسیون با دبی L/min 7/0 و غلظت ورودی g/m32 جهت تعیین صحت پیشبینیهای مدل انجام گرفت.
یافته‌ها: مدل قادر بود الگویی از افت محتوی رطوبت بستر بعد از قطع آبدهی بستر را ارائه کند. همچنین دادههای آزمایشگاهی افزایش غلظت خروجی آلاینده به دلیل کاهش رطوبت بستر به‌خوبی توسط مدل پیشبینی شدند (خطای استاندارد= 45/0%) و تخمینی از زمان نهایی ازکارافتادن سیستم ارائه شد.
نتیجه‌گیری: نتایج مدل‌سازی نشان داد که در نظر گرفتن ترم رطوبت در شرایط ناپایا میتواند حیطه عملکرد مدل را توسعه دهد و فرایند بیوفیلتراسیون را با تمام پیچیدگیهای آن توصیف نماید.
واژه‌های کلیدی: مدل‌سازی دینامیک، بیوفیلتراسیون، رطوبت بستر، سوپر جاذب آب، هگزان
References
1-Mudliar S, Giri B, Padoley K, Satpute D, Dixit R, Bhatt P. Bioreactors for treatment of VOCs and odours – A review. Environmental Management 2010;91(10):39–54.
2-Khan FI, Ghoshal AK. Removal of volatile organic compounds from polluted air. Loss Prevention in the Process Industries 2000;13(5):27-45.
3-Mudliar S, Giri B, Padoley K, Satpute D, Dixit R, Bhatt P. Bioreactors for treatment of VOCs and odours–a review. Journal of Environmental Management 2010; 91(5):1039-1054.
4-Spigno G, De Faveri DM. Modeling of a Vapor-Phase Fungi Bioreactor for the Abatement of Hexane: Fluid Dynamics and Kinetic Aspects. Biotechnology and Bioengineering 2005;89(3):19-28
5-Arriaga S, Revah S. Mathematical Modeling and Simulation of Hexane Degradation in Fungal and Bacterial Biofilters: Effective Diffusivity and Partition Aspects. Canadian Journal of Civil Engineering 2009;36(19):19-25.
6-Morales M, Hernaändez S, Cornabeä T, Revah S, Auria R. Effect of Drying on Biofilter Performance: Modeling and Experimental Approach. Environment and Science Technology 2003;37(9):85-92.
7-Fazaelipoor MS. Modeling temperature variations and moisture requirements in waste air biofilters under steady state conditions. Environmental Technology 2012; 33(5):7-13.
8-Akbari M.The use of water super absorbent inbiofiltration of organic-contaminated air [MSc. Thesis]. Kerman, Shahid Bahonar university, 2009. [Persian]
9-De Visscher A, Li GQ. Toluene removal biofilter modeling: Optimization and case study. Process Safety and Environmental Protection 2008; 86(4):277-282.
10-Fazaelipoor, MH. Analysis of a dual liquid phase biofilter for the removal of hydrophobic organic compounds from airstreams. Chemical Engineering Journal 2009;147(2): 110-116.
نویسندگان:سروش دانایی1،محمدحسن فضایلی پور2، محمدحسین سلمانی3
1.باشگاه پژوهشگران جوان و نخبگان، واحد شهریار، دانشگاه آزاد اسلامی، تهران
2.دانشیار مهندسی شیمی، گروه مهندسی شیمی و پلیمر پردیس فنی مهندسی، دانشگاه یزد
3.نویسنده مسئول: دکتری شیمی تجزیه، دانشکده بهداشت، دانشگاه علوم پزشکی و خدمات بهداشتی درمانی شهید صدوقی یزد تلفن تماس: 09131582878 Email:mhsn06@ssu.ac.ir
چکیده
مقدمه: یک روش مؤثر در تصفیه بیولوژیکی هوای آلوده به ترکیبات فرار آلی روش بیوفیلتراسیون است که با عبور هوا از درون بستر پرشده، ماده آلاینده از فاز گاز به فاز بیوفیلم منتقل‌شده و توسط میکروارگانیسمها تجزیه میگردد. فعالیت میکروب‌های تجزیهکننده به‌شدت تحت تأثیر رطوبت بستر است که مهم‌ترین فاکتور برای کنترل عملکرد بهینه بیوفیلتراسیون به‌حساب میآید. هدف از این مطالعه، به‌کارگیری ترم رطوبت برای توصیف عملکرد مدلی کاربردی در بیوفیلتر ناپایا برای حذف هگزان از جریان هوای آلوده است.
روش بررسی: ابتدا یک مدل دوبعدی از بیوفیلتراسیون هوای آلوده به هگزان ارائه شد که در آن ترمهای رطوبت بستر و پراکندگی محوری در موازنه انتقال جرم در نظر گرفته شدند. اثرات استفاده از سوپر جاذب آب نظیر افزایش مداوم تخلخل بستر و سطح ویژه آن درون معادلات وارد گردید. سپس حل عددی معادلات با روش حجم محدود انجام شد و در انتها، آزمایشهای بیوفیلتراسیون با دبی L/min 7/0 و غلظت ورودی g/m32 جهت تعیین صحت پیشبینیهای مدل انجام گرفت.
یافته‌ها: مدل قادر بود الگویی از افت محتوی رطوبت بستر بعد از قطع آبدهی بستر را ارائه کند. همچنین دادههای آزمایشگاهی افزایش غلظت خروجی آلاینده به دلیل کاهش رطوبت بستر به‌خوبی توسط مدل پیشبینی شدند (خطای استاندارد= 45/0%) و تخمینی از زمان نهایی ازکارافتادن سیستم ارائه شد.
نتیجه‌گیری: نتایج مدل‌سازی نشان داد که در نظر گرفتن ترم رطوبت در شرایط ناپایا میتواند حیطه عملکرد مدل را توسعه دهد و فرایند بیوفیلتراسیون را با تمام پیچیدگیهای آن توصیف نماید.
واژه‌های کلیدی: مدل‌سازی دینامیک، بیوفیلتراسیون، رطوبت بستر، سوپر جاذب آب، هگزان
مقدمه
با توسعه صنعتی در طی سالیان گذشته هوا، آب‌وخاک در سطح وسیعی آلوده‌شده‌اند. از جمله آلودگیهای زیست محیطی میتوان به ترکیبات آلی فرار و بدبو اشاره نمود که هوای اطراف مراکز صنعتی را آلوده میکنند (1). بااین‌حال، به‌تناسب افزایش سطح آگاهی عمومی در مسائل زیست‌محیطی، مطالبات عمومی گسترش‌یافته است به‌طوری‌که در سالیان اخیر، قوانین سختگیرانهتری در جهت کنترل انتشار این ترکیبات وضع گردیده است. صنایع به‌طورمعمول از دو روش کلی کنترل انتشار آلودگی و اصلاح تجهیزات و فرآیند از طریق حذف یا بازیابی به کاهش ترکیبات آلی فرار از جریان‌های گازی می‌پردازند (2). از میان روشهایی که بر مبنای حذف ترکیبات آلاینده توسعه‌یافته‌اند، بیوفیلتراسیون روشی بیولوژیکی و دوستدار محیطزیست است. در فرآیند بیوفیلتراسیون هوا با عبور از یک ستون که حاوی بستر میکروب‌ها است، ماده آلاینده خود را به فاز بیوفیلم منتقل می‌کند. بیوفیلم تشکیل‌شده روی ذرات بستر (پرکن) حاوی مجموعهای از میکروارگانیسمها است که مسئولیت تخریب آلاینده را در یک واکنش بیولوژیکی بر عهده‌دارند (3). این پدیده به همراه سرعت نفوذ ماده آلاینده به بیوفیلم، محدودیت‌های فرآیند بیوفیلتراسیون را تشکیل می‌دهند که در مدل‌سازی این فرآیند باید مدنظر قرار گیرند. برای توسعه مدل این پژوهش، ابتدا معادلات انتقال جرم شامل آلاینده و رطوبت درون هر بخش دیفرانسیلی از سیستم برقرار می‌شود و پس از گسسته سازی با کد نویسی در محیط نرم‌افزار MATLAB حل می‌شود. همان‌طور که در شکل 1 نشان داده‌شده است، هر بخش شامل دو فاز گاز و بیوفیلم است. مختصات z برای فاز گاز و x برای فاز بیوفیلم در نظر گرفته‌شده است.
شکل 1: نمایی از مکانیسم‌های انتقال جرم درون یک بخش دیفرانسیلی
هرچند مدلهای زیادی برای فرایند بیوفیلتراسیون ارائه‌شده‌اند در اینجا چند مدل مهمی که در سالهای اخیر برای فرایند بیوفیلتراسیون ارائه‌شده است مطرح میشوند.
اسپیگنو و فاوری (2005) مدلی پایا برای یک بیوفیلتر حذف هگزان در حضور قارچ اسپرژیلوس نایجر ارائه کردند که در آن اثرات پراکندگی محوری در نظر گرفته شد. نتایج مدل‌سازی تطابق خوبی بین مقادیر پیش‌بینی‌شده توسط مدل پایا و مقادیر تجربی را نشان داد. به دلیل طبیعت میکروارگانیسم به‌کاررفته، مجموعه آزمایش‌ها نتوانست تکرارپذیری نتایج در تعیین برخی پارامترهای مدل نظیر ثابت اشباعیت، نرخ رشد ویژه و بازده را تائید کند (4). آریاگا و همکارانش (2009) پارامترهای بیوفیلتراسیون هگزان را با استفاده از میکروارگانیسم‌های باکتریایی و قارچی در شرایط آزمایشگاهی تعیین کردند و در ادامه کار این پارامترها را برای تعیین مدل ریاضی در شرایط پایا به کار گرفتند. مدل ارائه‌شده با در نظر گرفتن واکنش مرتبه اول به‌عنوان مرحله محدودکننده، با خطای حدود یک درصد نتایج آزمایشگاهی را پیشبینی کرد. این در حالی است که خطای مدل‌سازی با واکنش مرتبه صفر با محدودکنندگی واکنش خطایی بیشتر از 18 درصد داشت (5). مورالس و همکاران (2003) مدلی برای توصیف کاهش آب در یک بیوفیلتر حذف تولوئن ارائه کردند تا اثر آن را روی عملکرد سیستم پیشبینی کنند. تبخیر درون بیوفیلتر نتیجه زیر اشباع بودن هوای ورودی و همچنین واکنش بیولوژیکی درون بیوفیلم در نظر گرفته شد. ترم تخریب در این مدل شامل دما، محتوی آب و اثرات غلظت به‌صورت توابعی جدا از هم در نظر گرفته شد. نتایج نشان داد که در نقطه بحرانی، عملکرد بیوفیلتر در حذف آلاینده تا 50 درصد کاهش نشان میدهد. مدل، تغییرات در غلظت آلاینده، رطوبت نسبی هوا، دما و محتوی آب محیط را به‌خوبی توصیف کرد (6). فضایلی پور (2012) مدلی برای بیوفیلتراسیون غیر هم‌دما ارائه نمود. در این مدل تبخیر آب و اثرات گرمای نهان تبخیر در نظر گرفته شد و تمامی متغیرها به فرم وابسته به دما به کار گرفته شدند. نتایج مدل‌سازی نشان داد دمای هوا و رطوبت در طول بیوفیلتر در جهت جریان هوا افزایش می‌یابد (7).
یکی از علل به وجود آمدن شرایط نا پایا در بیوفیلتر، قطع ناگهانی آب دهی بستر است. تبخیر مداوم از بستر و انجام واکنش گرمازا سبب خشکی بستر می‌شود. حساسیت بالای میکروارگانیسمها به محتوی آب بستر بیوفیلتر ضرورت پیشبینی شرایط بیوفیلتر با استفاده از مدل را بیشتر مشخص می‌کند. اکبری و همکاران (2011) در تحقیقی بر روی بهبود سیستم تأمین رطوبت درون بیوفیلتر به بررسی استفاده از سوپر جاذب‌های آب پرداختند. سوپر جاذب‌های آب ترکیبات پلیمری هستند که قادر به ذخیره‌سازی آب درون خود می‌باشند و در محیطی با رطوبت کمتر قادر به آزادسازی آب ذخیره‌شده می‌باشند. نتایج نشان داد حضور سوپر جاذب روند کاهش بازده حذف را آرامتر کرده و عملکرد بیوفیلتر را در برابر قطع آب ناگهانی و حتی شوک غلظت مقاومتر می‌کند (8). هدف از این پژوهش، توصیف بهتر فرآیند بیوفیلتراسیون نا‌پایای هوای آلوده به هگزان از طریق ارائه مدلی کاربردی می باشد.
روش بررسی
فرضیات ساده کننده مدل به شرح زیر است: 1- هوا به‌صورت گاز ایده آل و کلیه مشخصات فیزیکی ثابت فرض می‌شود. 2-تغییرات دمایی در طول بیوفیلتر ناچیز است و تعادل گرمایی مابین فازها برقرار است. 3-ذرات متخلخل پرکن و سوپر جاذب به شکل کروی (با قطر متوسط یکسان) در نظر گرفته می‌شود. 4- از تغییرات غلظت در جهت شعاع صرف‌نظر شده و پراکندگی محوری در فاز گاز، با به کار گرفتن ضریب پراکندگی در معادلات منظور می گردد. 5- بیوفیلم به‌صورت یک فاز همگن در نظر گرفته می‌شود که واکنش بیولوژیکی با رابطه Monod درون آن انجام می شود. همچنین اثر ترم رطوبت در واکنش مستقل از غلظت و از طریق ضریبی اعمال میگردد. 6- به دلیل ضخامت کم بیوفیلم نسبت به ذرات پرکن (در این پژوهش پرلیت و سوپر جاذب ذرات پرکن درون بستر هستند)، مختصات دکارتی برای معادلات در نظر گرفته می‌شود. 7- مقاومت در فاز گاز و فصل مشترک دو فاز ناچیز است، درنتیجه غلظت آلاینده در فصل مشترک گاز/ بیوفیلم در تعادل بوده و از قانون هنری تبعیت می‌کند. 8- تبادل رطوبت مابین ذرات سوپر جاذب و پرکن تنها از طریق هوا صورت می‌پذیرد. 9- فرآیند خروج آب از سوپر جاذب به‌صورت تبخیر از سطح یک قطره با قطر هم‌اندازه آن مدل می‌شود. همچنین نرخ فرآیند تبخیر و انتقال به هوا از طریق ضریب انتقال جرم توصیف می‌شود. 10- تخلخل بستر در طول فرآیند متغیر است. تغییرات این ترم به سبب خروج آب از خلل و فرج ذرات پرکن و نیز کوچک شدن قطر ذرات سوپر جاذب اتفاق می‌افتد. این پدیده علاوه بر افزایش تخلخل، کاهش سطح سوپر جاذب را نیز به همراه دارد (10،9).
بر اساس این فرضیات، معادلات مدل با برقراری موازنه‌های جرم و رطوبت بنا شدند (7). جهت بهبود کارایی مدل، چندین تغییر در این معادلات صورت گرفت که ازجمله آنها میتوان به وارد نمودن ترمهای تغییرات بازمان در معادلات آلاینده و رطوبت، تغییر در تخلخل بستر و در نظر گرفتن پراکندگی محوری اشاره نمود. در کنار این معادلات، دو معادله توصیف‌کننده میزان آب موجود در پرکن و سوپر جاذب در هر بخش به کار گرفته‌شده است که به‌صورت معادله 1 و 2 معرفی شدند.
(1)
(2)
که در آنها (m-1)ap و (m-1)as به ترتیب معرف سطح ویژه انتقال جرم آب در پرکن و سوپر جاذب، (kgH2O/m3air)mp و (kgH2O/m3 air)ms به ترتیب بیانگر ذخیره آب پرلیت و سوپر جاذب در بیوفیلتر و (kg H2O/m3 air) mp0 و (kg H2O/m3 air) ms0 ذخیره آب پرلیت و سوپر جاذب در ابتدای شرایط نا پایا هستند. همان‌طور که در بخش فرضیات بیان شد،
سطح ویژه سوپر جاذب در هر تکرار زمانی متغیر بوده حال‌آنکه سطح ویژه پرکن ثابت فرض می‌شود.
همچنین با توجه به دادههای معمول معرفی‌شده برای دو محصول پرلیت (به‌عنوان ماده پرکن) و سوپر جاذب، محاسبات میزان جذب آب این ذرات برای تعیین ذخیره آب اولیه ستون در ابتدای شرایط نا پایا انجام شد. پارامترهای به‌کاررفته در این پژوهش، توسط برازش دادههای منابع مختلف (5،4) و یا محاسبه آنها به دست آمد. دادههای اولیه نیز از منابع ( 8،5) استخراج شدندکه تمامی این داده ها مربوط به بیوفیلترهای تصفیه هگزان تحت شرایط پایا می باشند. از این داده ها در تعیین صحت پیشبینیهای مدل در دوحالت‌پایا و ناپایا استفاده شد.
همچنین آزمایشهای تجربی بیوفیلتراسیون هوای آلوده به هگزان نیز در دو فاز انجام شد. در فاز اول آزمایش‌ها، تاثیر استفاده از سوپرجاذب (پرلیت) در نرخ تبخیر درون یک بیوفیلتر بررسی شد. بدین منظور نرخ تبخیر از دو بستر متفاوت، یکی پرلیت و دیگری مخلوط پرلیت و سوپر جاذب موردبررسی قرار گرفت. برای این کار، دو ستون کوچک 5 سانتی‌متری از بسترها پر شد.
در یک دوره 7 روزه با وزن کردن در زمان‌های مختلف (معمولاً یک‌روزه) درصد کاهش آب نسبت به حالت اولیه اندازه‌گیری شد. دبی جریان L/min 7/0 برای هردو ستون ثابت بود. نتایج این کار تخمین خوبی از میزان تبخیر در بیوفیلتر ارائه داد.
فاز دوم، آزمایشهای بیوفیلتراسیون است که مسیر عملکرد بیوفیلتر به سه بازه زمانی تقسیم شد. اولین بازه زمانی مربوط به افزایش تدریجی رشد میکروبی و تشکیل بیوفیلم درون بستر است. در این بازه به‌طورکلی منحنی بازده حذف شیب صعودی دارد.
بازه دوم مربوط به زمانی است که بیوفیلترها به حداکثر رشد میکروبی رسیده‌اند که منحنی بازده حذف در بیشینه خود به‌صورت افقی ادامه پیدا می‌کند. در بازه سوم که حذف آب زنی (محلول معدنی محیط کشت) بستر انجام می‌شود.
ستون مورداستفاده در این پژوهش، ستون شیشهای استوانهای شکل با قطر 05/0 متر و ارتفاع بستر 5/0 متر بود. از این ستون دو عدد کاملاً مشابه ساخته شد. بستر یکی از دو بیوفیلتر از پرلیت تشکیل‌شده و در بستر دیگری مخلوطی از پرلیت و سوپر جاذب آب به‌کاررفته است. برای این مخلوط مقدار 3/2 درصد خشک سوپر جاذب (با توجه به پژوهش انجام‌شده در تعیین بهترین نسبت) انتخاب شد (8). این کار بستری یکنواخت با نسبت پرلیت به سوپر جاذب 3 :1 تشکیل داد. بسترها به مدت 24 ساعت در محلول معدنی محیط کشت میکروبی قرار گرفتند. درنتیجه با جذب محلول، سوپر جاذب‌ها متورم شده و به قطر برابر با پرلیت رسیدند. بسترهای تهیه‌شده درنهایت به درون ستونها ریخته شده و مجموعه بیوفیلتر آماده تلقیح میکروبی شد. آب زنی بستر تا زمان رسیدن به حالت‌پایا به‌طور مداوم و بافاصله یک روز انجام شد. با برقراری شرایط پایا در غلظت خروجی بیوفیلترها، آب زنی بستر به‌طور کامل قطع گردید.
یافته‌ها
ابتدا نتایج آزمایش‌ها بیوفیلتراسیون ارائه میگردد و در ادامه، نتایج پیشبینی دادههای بیوفیلتراسیون پایا (استخراج‌شده از مقالات) تحلیل میگردد. پس از تائید مدل پایا، به بررسی مدل نا پایا و مقایسه آن با دادههای تجربی پرداخته می‌شود.
بر اساس آزمایشهای اولیه انجام‌شده، در شکل 2 تفاوت میان سرعت تبخیر سوپر جاذب و پرلیت به‌خوبی مشخص است. با انجام آزمون T برای تعیین اختلاف جفت داده‌ها بر روی نتایج حاصل از آزمایش، بافاصله اطمینان 95 درصد اختلاف بین درصد کاهش آب دو بستر در محدوده 672/5- 781/1 قرار داشت.
شکل 2: نمودار تفاوت درصد کاهش آب پرلیت و سوپر جاذب
شکل 3 نتایج دادههای آزمایشگاهی حاصل از دو بیوفیلتر را در دبی L/min7/0 و غلظت ورودی g/m3 2 نشان میدهد. زمان‌ماند ظاهری (EBCT) در بیوفیلترها برابر 40/1 دقیقه محاسبه شد. همان‌طور که گفته شد فراهم نمودن محلول معدنی تا روز 25 ادامه داشته است. پس از تائید شرایط پایا در دو بیوفیلتر، حذف آبدهی بستر، منجر به شرایط نا پایا در دو بیوفیلتر شده است.
شکل 4 تغییرات بازده حذف دو بیوفیلتر را در مدت‌زمان 40 روز نشان میدهد. با شروع روند کاهش رطوبت در بستر بیوفیلترها، درصد حذف نیز کاهش می‌یابد که در شکل، اثر حضور سوپر جاذب در بهبود عملکرد بیوفیلتر به‌خوبی مشخص است.
شکل 3: تغییرات غلظت خروجی (دادههای آزمایشگاهی)
دبی L/min7/0 و غلظت ورودی g/m32
شکل 4: تغییرات بازده حذف بیوفیلترهای آزمایشگاهی با زمان
شکل 5 پیشبینی دینامیکی تغییرات غلظت خروجی در بیوفیلتر با بستر پرلیت و سوپر جاذب (دادههای آزمایشگاهی شکل 3) توسط مدل با در نظر گرفتن ترم تغییرات رطوبت بستر را ارائه میدهد. تفاوت در پیشبینیهای مدل و دادههای تجربی در شروع شرایط نا پایا بسیار کمتر است (1/0% خطا) اما رفته‌رفته در عمل به دلیل تغییر در برخی پارامترها نظیر سطح ویژه بیوفیلم با کاهش رطوبت بستر این تفاوت بیشتر نمایان می‌شود. میزان خطای میانگین برای پیشبینی این داده‌ها 45/0 درصد است.
شکل 5: مقایسه مدل و دادههای آزمایشگاهی بیوفیلتراسیون
نا پایا با بستر حاوی سوپر جاذب و پرلیت
بحث و نتیجه‌گیری
با توجه به نتایج همان‌طور که در فرضیات مدل ذکر شد برای وارد شدن اثر مقاومت درونی سوپر جاذب درون معادلات، از ضریب مشخصی در معادلات مدل استفاده شد. این ضریب سبب اندکی کاهش نرخ تبخیر نسبت به حالت اشباع سطح می‌شود. با توجه به نتایج شکل 2، این ضریب 95/0 تعیین شد (میانگین اختلاف دو نمودار 5 درصد در نظر گرفته شد). همچنین در بررسی شکل 3 مشخص می‌شود تحت شرایط نا پایا، غلظت در خروجی به‌تدریج زیاد می‌شود. حضور سوپر جاذب در ابتدا مقاومت خوبی را برای حفظ شرایط پایا تا سه روز اول نشان میدهد. بعد از روز سوم نیز سوپر جاذب روند افزایش غلظت را کندتر کرده و مانع تغییرات شدید در عملکرد بیوفیلتر می‌شود که با نتایج به‌دست‌آمده در تحقیقات گذشته مطابقت دارد (8).
بررسی شکل 5 نشان می دهد که تغییرات غلظت خروجی با نرخ ثابتی رو به افزایش است. این امر به دلیل افزایش تدریجی تخلخل بستر و اثر این افزایش روی کمیت‌های مختلف است. در این شکل، تطابق خوبی بین دادههای پیش‌بینی‌شده توسط مدل و دادههای تجربی حاصل از آزمایش مشاهده می‌شود. هرچند ثابت فرض کردن مقادیر پارامترهای مدل در طول زمان نا پایا، منجر به اندکی تفاوت در نتایج شده اما روند افزایش غلظت در طول زمان به‌خوبی مدل شده است.
پس از قطع جریان آب یک بیوفیلتر، به دلیل انجام واکنش بیولوژیکیِ تخریب آلاینده که از نوع واکنش‌های گرمازا است، سرعت تبخیر درون بستر زیاد می شود. رفته‌رفته با تبخیر رطوبت و خشکی بستر، از فعالیت میکروارگانیسم‌های مسئول تخریب ماده آلاینده کاسته شده و غلظت خروجی شروع به افزایش می‌کند. هوای تازه با رطوبت کمتر از اشباع، به‌آرامی قسمتهایی از بستر را خشک نموده و از فعالیت می‌اندازد. این مسئله در مدل با استفاده از gε به‌صورت محلی درون بخش دیفرانسیلی به طرز جالبی در نظر گرفته شد. به‌طوری‌که هر بخش مقدار آب مشخصی از کل آب در دسترس بیوفیلتر را در اختیار دارد و با نرخ انتقال جرم مشخصی آب را از دست میدهد. این روند کاهش، بر روی gε محلی اثر می‌گذارد. برخی ترم‌ها در معادلات انتقال جرم نظیر ضرایب معادلات غلظت آلاینده در فازها، سطح ویژه سوپر جاذب که آن‌هم وابسته به شعاع ذرات است (سوپر جاذب با از دست رفتن آب ذخیره‌شده درون آن کوچک می‌شود) و ضریب انتقال جرم در معادله موازنه رطوبت وابسته به تغییرات gε هستند. تغییر gε باعث تغییر در این ضرایب می‌شود و اثر نهایی آن، افزایش غلظت‌ها در آن بخش است. به دلیل آنکه هر بخش، ورودی خود را از بخش قبلی می‌گیرد اثر افزایش غلظت در بخش‌های اولیه با انتقال به بخش‌های دیگر، در خروجی مشاهده می‌شود. پس از به اتمام رسیدن آب درون یک بخش دیفرانسیلی، روند کاهش آب در بخش بعدی آغازشده و به‌این‌ترتیب تا به پایان رسیدن تمام محتوای آب در بستر ادامه می‌یابد. روند افزایش، درنهایت به ازکارافتادن کامل سیستم منجر می‌شود. در این شرایط، غلظت ورودی آلاینده بدون هیچ تغییری در خروجی بیوفیلتر مشاهده می‌شود. معیار مشخص نمودن این شرایط جهت قطع اجرای کدهای برنامه MATLAB، عدم‌تغییر و برابری غلظت‌ها در 10 تکرار آخر است. در حل معادلات نا پایا اگرچه خصوصیات فیزیکی مواد با تغییر رطوبت ثابت فرض شد اما با تعریف دو ضریب در معادلات، اثر کاهش رطوبت اعمال شد. مقادیر این ضرایب با نسبت میزان آب موجود در آن بخش به میزان آب اولیه تعیین شد. چون باگذشت زمان نرخ انتقال آب تحت اثر مویینگی در خلل و فرج پرلیت و مقاومت داخلی سوپر جاذب قرار می‌گیرد، این ضرایب نرخ انتقال آب از سوپر جاذب و پرلیت و میزان نفوذ آلاینده و واکنش درون بیوفیلم را کاهش می‌دهند. این امر بدان دلیل است که رفته‌رفته با کاهش محتوی رطوبت درون بستر، میکروارگانیسمها از بین رفته و فعالیت خود را از دست می‌دهند. درنتیجه عملکرد بیوفیلتر به‌طور مداوم کاهش می‌یابد تا درنهایت عمل حذف آلاینده درون بیوفیلتر به‌کلی قطع میگردد. مدل این پژوهش توانست به‌طور مناسبی شرایط نا پایای حاکم بر بیوفیلتر و روند افزایش غلظت در آن را پیشبینی کند (خطای 45/0 درصد). هرچند دادههای آزمایش‌ها تا رسیدن به غلظت ورودی موجود نیست، مدل این روند را پیشبینی می‌کند و درعین‌حال تخمینی از تعداد روزها (23 روز) تا ازکارافتادن کامل بیوفیلتر را ارائه می‌کند. تحت شرایط عملیاتی بیوفیلتر، از زمان قطع آب دهی تا روز 23 شرایط نا پایا حاکم است و پس‌ازآن غلظت خروجی به غلظت ورودی می‌رسد و بیوفیلتر به‌طور کامل از کار می‌افتد. اثر کاهش رطوبت در بخش قبلی و نتیجه آن کاهش رطوبت فاز گاز در همان بخش، به‌آرامی بر روی بخش بعدی انتقال‌یافته تا درنهایت بخش مذکور با از دست دادن کامل آب در دسترس خود از مسیر محاسبات خارج شود. در این میان، شعاع، سطح و حجم ذرات سوپر جاذب در آن بخش کم شده و مقدار ضریب انتقال جرم افزایش می‌یابد که محاسبات آنها در مدل به‌طور کامل انجام می‌گیرد.
بر این اساس، مدل‌سازی تغییرات رطوبت در بیوفیلتر و اثر آن بر
غلظت خروجی به‌خوبی توانست عملکرد نا پایای بیوفیلتر را توصیف کند. رطوبت در نرخ انتقال آلاینده به بیوفیلم مؤثر است، به‌طوری‌که کاهش رطوبت بستر علاوه بر کاهش نرخ واکنش در بیوفیلم میکروبی، جذب و نفوذ آلاینده در بیوفیلم را نیز تحت تأثیر قرار میدهد. نتایج مدل‌سازی نشان داد در نظر گرفتن ترم رطوبت در شرایط نا پایا میتواند حیطه عملکرد مدل را توسعه دهد و فرایند بیوفیلتراسیون هوای آلوده را با تمام پیچیدگیهای آن توصیف نماید.
تقدیر و تشکر
نویسنده این مقاله، آقای سروش دانایی (دانشجوی مقطع دکتری) از همکاری خانم مهندس ندا بدلی در پیاده سازی مدل، تشکر و قدردانی می نماید.
تضاد منافع
نویسندگان این مقاله اعلام می دارند که هیچ تضاد منافعی وجود ندارد.
با توسعه صنعتی در طی سالیان گذشته هوا، آب‌وخاک در سطح وسیعی آلوده‌شده‌اند. از جمله آلودگیهای زیست محیطی میتوان به ترکیبات آلی فرار و بدبو اشاره نمود که هوای اطراف مراکز صنعتی را آلوده میکنند (1). بااین‌حال، به‌تناسب افزایش سطح آگاهی عمومی در مسائل زیست‌محیطی، مطالبات عمومی گسترش‌یافته است به‌طوری‌که در سالیان اخیر، قوانین سختگیرانهتری در جهت کنترل انتشار این ترکیبات وضع گردیده است. صنایع به‌طورمعمول از دو روش کلی کنترل انتشار آلودگی و اصلاح تجهیزات و فرآیند از طریق حذف یا بازیابی به کاهش ترکیبات آلی فرار از جریان‌های گازی می‌پردازند (2). از میان روشهایی که بر مبنای حذف ترکیبات آلاینده توسعه‌یافته‌اند، بیوفیلتراسیون روشی بیولوژیکی و دوستدار محیطزیست است. در فرآیند بیوفیلتراسیون هوا با عبور از یک ستون که حاوی بستر میکروب‌ها است، ماده آلاینده خود را به فاز بیوفیلم منتقل می‌کند. بیوفیلم تشکیل‌شده روی ذرات بستر (پرکن) حاوی مجموعهای از میکروارگانیسمها است که مسئولیت تخریب آلاینده را در یک واکنش بیولوژیکی بر عهده‌دارند (3). این پدیده به همراه سرعت نفوذ ماده آلاینده به بیوفیلم، محدودیت‌های فرآیند بیوفیلتراسیون را تشکیل می‌دهند که در مدل‌سازی این فرآیند باید مدنظر قرار گیرند. برای توسعه مدل این پژوهش، ابتدا معادلات انتقال جرم شامل آلاینده و رطوبت درون هر بخش دیفرانسیلی از سیستم برقرار می‌شود و پس از گسسته سازی با کد نویسی در محیط نرم‌افزار MATLAB حل می‌شود. همان‌طور که در شکل 1 نشان داده‌شده است، هر بخش شامل دو فاز گاز و بیوفیلم است. مختصات z برای فاز گاز و x برای فاز بیوفیلم در نظر گرفته‌شده است.
شکل 1: نمایی از مکانیسم‌های انتقال جرم درون یک بخش دیفرانسیلی
هرچند مدلهای زیادی برای فرایند بیوفیلتراسیون ارائه‌شده‌اند در اینجا چند مدل مهمی که در سالهای اخیر برای فرایند بیوفیلتراسیون ارائه‌شده است مطرح میشوند.
اسپیگنو و فاوری (2005) مدلی پایا برای یک بیوفیلتر حذف هگزان در حضور قارچ اسپرژیلوس نایجر ارائه کردند که در آن اثرات پراکندگی محوری در نظر گرفته شد. نتایج مدل‌سازی تطابق خوبی بین مقادیر پیش‌بینی‌شده توسط مدل پایا و مقادیر تجربی را نشان داد. به دلیل طبیعت میکروارگانیسم به‌کاررفته، مجموعه آزمایش‌ها نتوانست تکرارپذیری نتایج در تعیین برخی پارامترهای مدل نظیر ثابت اشباعیت، نرخ رشد ویژه و بازده را تائید کند (4). آریاگا و همکارانش (2009) پارامترهای بیوفیلتراسیون هگزان را با استفاده از میکروارگانیسم‌های باکتریایی و قارچی در شرایط آزمایشگاهی تعیین کردند و در ادامه کار این پارامترها را برای تعیین مدل ریاضی در شرایط پایا به کار گرفتند. مدل ارائه‌شده با در نظر گرفتن واکنش مرتبه اول به‌عنوان مرحله محدودکننده، با خطای حدود یک درصد نتایج آزمایشگاهی را پیشبینی کرد. این در حالی است که خطای مدل‌سازی با واکنش مرتبه صفر با محدودکنندگی واکنش خطایی بیشتر از 18 درصد داشت (5). مورالس و همکاران (2003) مدلی برای توصیف کاهش آب در یک بیوفیلتر حذف تولوئن ارائه کردند تا اثر آن را روی عملکرد سیستم پیشبینی کنند. تبخیر درون بیوفیلتر نتیجه زیر اشباع بودن هوای ورودی و همچنین واکنش بیولوژیکی درون بیوفیلم در نظر گرفته شد. ترم تخریب در این مدل شامل دما، محتوی آب و اثرات غلظت به‌صورت توابعی جدا از هم در نظر گرفته شد. نتایج نشان داد که در نقطه بحرانی، عملکرد بیوفیلتر در حذف آلاینده تا 50 درصد کاهش نشان میدهد. مدل، تغییرات در غلظت آلاینده، رطوبت نسبی هوا، دما و محتوی آب محیط را به‌خوبی توصیف کرد (6). فضایلی پور (2012) مدلی برای بیوفیلتراسیون غیر هم‌دما ارائه نمود. در این مدل تبخیر آب و اثرات گرمای نهان تبخیر در نظر گرفته شد و تمامی متغیرها به فرم وابسته به دما به کار گرفته شدند. نتایج مدل‌سازی نشان داد دمای هوا و رطوبت در طول بیوفیلتر در جهت جریان هوا افزایش می‌یابد (7).
یکی از علل به وجود آمدن شرایط نا پایا در بیوفیلتر، قطع ناگهانی آب دهی بستر است. تبخیر مداوم از بستر و انجام واکنش گرمازا سبب خشکی بستر می‌شود. حساسیت بالای میکروارگانیسمها به محتوی آب بستر بیوفیلتر ضرورت پیشبینی شرایط بیوفیلتر با استفاده از مدل را بیشتر مشخص می‌کند. اکبری و همکاران (2011) در تحقیقی بر روی بهبود سیستم تأمین رطوبت درون بیوفیلتر به بررسی استفاده از سوپر جاذب‌های آب پرداختند. سوپر جاذب‌های آب ترکیبات پلیمری هستند که قادر به ذخیره‌سازی آب درون خود می‌باشند و در محیطی با رطوبت کمتر قادر به آزادسازی آب ذخیره‌شده می‌باشند. نتایج نشان داد حضور سوپر جاذب روند کاهش بازده حذف را آرامتر کرده و عملکرد بیوفیلتر را در برابر قطع آب ناگهانی و حتی شوک غلظت مقاومتر می‌کند (8). هدف از این پژوهش، توصیف بهتر فرآیند بیوفیلتراسیون نا‌پایای هوای آلوده به هگزان از طریق ارائه مدلی کاربردی می باشد.
روش بررسی
فرضیات ساده کننده مدل به شرح زیر است: 1- هوا به‌صورت گاز ایده آل و کلیه مشخصات فیزیکی ثابت فرض می‌شود. 2-تغییرات دمایی در طول بیوفیلتر ناچیز است و تعادل گرمایی مابین فازها برقرار است. 3-ذرات متخلخل پرکن و سوپر جاذب به شکل کروی (با قطر متوسط یکسان) در نظر گرفته می‌شود. 4- از تغییرات غلظت در جهت شعاع صرف‌نظر شده و پراکندگی محوری در فاز گاز، با به کار گرفتن ضریب پراکندگی در معادلات منظور می گردد. 5- بیوفیلم به‌صورت یک فاز همگن در نظر گرفته می‌شود که واکنش بیولوژیکی با رابطه Monod درون آن انجام می شود. همچنین اثر ترم رطوبت در واکنش مستقل از غلظت و از طریق ضریبی اعمال میگردد. 6- به دلیل ضخامت کم بیوفیلم نسبت به ذرات پرکن (در این پژوهش پرلیت و سوپر جاذب ذرات پرکن درون بستر هستند)، مختصات دکارتی برای معادلات در نظر گرفته می‌شود. 7- مقاومت در فاز گاز و فصل مشترک دو فاز ناچیز است، درنتیجه غلظت آلاینده در فصل مشترک گاز/ بیوفیلم در تعادل بوده و از قانون هنری تبعیت می‌کند. 8- تبادل رطوبت مابین ذرات سوپر جاذب و پرکن تنها از طریق هوا صورت می‌پذیرد. 9- فرآیند خروج آب از سوپر جاذب به‌صورت تبخیر از سطح یک قطره با قطر هم‌اندازه آن مدل می‌شود. همچنین نرخ فرآیند تبخیر و انتقال به هوا از طریق ضریب انتقال جرم توصیف می‌شود. 10- تخلخل بستر در طول فرآیند متغیر است. تغییرات این ترم به سبب خروج آب از خلل و فرج ذرات پرکن و نیز کوچک شدن قطر ذرات سوپر جاذب اتفاق می‌افتد. این پدیده علاوه بر افزایش تخلخل، کاهش سطح سوپر جاذب را نیز به همراه دارد (10،9).
بر اساس این فرضیات، معادلات مدل با برقراری موازنه‌های جرم و رطوبت بنا شدند (7). جهت بهبود کارایی مدل، چندین تغییر در این معادلات صورت گرفت که ازجمله آنها میتوان به وارد نمودن ترمهای تغییرات بازمان در معادلات آلاینده و رطوبت، تغییر در تخلخل بستر و در نظر گرفتن پراکندگی محوری اشاره نمود. در کنار این معادلات، دو معادله توصیف‌کننده میزان آب موجود در پرکن و سوپر جاذب در هر بخش به کار گرفته‌شده است که به‌صورت معادله 1 و 2 معرفی شدند.
(1)
(2)
که در آنها (m-1)ap و (m-1)as به ترتیب معرف سطح ویژه انتقال جرم آب در پرکن و سوپر جاذب، (kgH2O/m3air)mp و (kgH2O/m3 air)ms به ترتیب بیانگر ذخیره آب پرلیت و سوپر جاذب در بیوفیلتر و (kg H2O/m3 air) mp0 و (kg H2O/m3 air) ms0 ذخیره آب پرلیت و سوپر جاذب در ابتدای شرایط نا پایا هستند. همان‌طور که در بخش فرضیات بیان شد،
سطح ویژه سوپر جاذب در هر تکرار زمانی متغیر بوده حال‌آنکه سطح ویژه پرکن ثابت فرض می‌شود.
همچنین با توجه به دادههای معمول معرفی‌شده برای دو محصول پرلیت (به‌عنوان ماده پرکن) و سوپر جاذب، محاسبات میزان جذب آب این ذرات برای تعیین ذخیره آب اولیه ستون در ابتدای شرایط نا پایا انجام شد. پارامترهای به‌کاررفته در این پژوهش، توسط برازش دادههای منابع مختلف (5،4) و یا محاسبه آنها به دست آمد. دادههای اولیه نیز از منابع ( 8،5) استخراج شدندکه تمامی این داده ها مربوط به بیوفیلترهای تصفیه هگزان تحت شرایط پایا می باشند. از این داده ها در تعیین صحت پیشبینیهای مدل در دوحالت‌پایا و ناپایا استفاده شد.
همچنین آزمایشهای تجربی بیوفیلتراسیون هوای آلوده به هگزان نیز در دو فاز انجام شد. در فاز اول آزمایش‌ها، تاثیر استفاده از سوپرجاذب (پرلیت) در نرخ تبخیر درون یک بیوفیلتر بررسی شد. بدین منظور نرخ تبخیر از دو بستر متفاوت، یکی پرلیت و دیگری مخلوط پرلیت و سوپر جاذب موردبررسی قرار گرفت. برای این کار، دو ستون کوچک 5 سانتی‌متری از بسترها پر شد.
در یک دوره 7 روزه با وزن کردن در زمان‌های مختلف (معمولاً یک‌روزه) درصد کاهش آب نسبت به حالت اولیه اندازه‌گیری شد. دبی جریان L/min 7/0 برای هردو ستون ثابت بود. نتایج این کار تخمین خوبی از میزان تبخیر در بیوفیلتر ارائه داد.
فاز دوم، آزمایشهای بیوفیلتراسیون است که مسیر عملکرد بیوفیلتر به سه بازه زمانی تقسیم شد. اولین بازه زمانی مربوط به افزایش تدریجی رشد میکروبی و تشکیل بیوفیلم درون بستر است. در این بازه به‌طورکلی منحنی بازده حذف شیب صعودی دارد.
بازه دوم مربوط به زمانی است که بیوفیلترها به حداکثر رشد میکروبی رسیده‌اند که منحنی بازده حذف در بیشینه خود به‌صورت افقی ادامه پیدا می‌کند. در بازه سوم که حذف آب زنی (محلول معدنی محیط کشت) بستر انجام می‌شود.
ستون مورداستفاده در این پژوهش، ستون شیشهای استوانهای شکل با قطر 05/0 متر و ارتفاع بستر 5/0 متر بود. از این ستون دو عدد کاملاً مشابه ساخته شد. بستر یکی از دو بیوفیلتر از پرلیت تشکیل‌شده و در بستر دیگری مخلوطی از پرلیت و سوپر جاذب آب به‌کاررفته است. برای این مخلوط مقدار 3/2 درصد خشک سوپر جاذب (با توجه به پژوهش انجام‌شده در تعیین بهترین نسبت) انتخاب شد (8). این کار بستری یکنواخت با نسبت پرلیت به سوپر جاذب 3 :1 تشکیل داد. بسترها به مدت 24 ساعت در محلول معدنی محیط کشت میکروبی قرار گرفتند. درنتیجه با جذب محلول، سوپر جاذب‌ها متورم شده و به قطر برابر با پرلیت رسیدند. بسترهای تهیه‌شده درنهایت به درون ستونها ریخته شده و مجموعه بیوفیلتر آماده تلقیح میکروبی شد. آب زنی بستر تا زمان رسیدن به حالت‌پایا به‌طور مداوم و بافاصله یک روز انجام شد. با برقراری شرایط پایا در غلظت خروجی بیوفیلترها، آب زنی بستر به‌طور کامل قطع گردید.
یافته‌ها
ابتدا نتایج آزمایش‌ها بیوفیلتراسیون ارائه میگردد و در ادامه، نتایج پیشبینی دادههای بیوفیلتراسیون پایا (استخراج‌شده از مقالات) تحلیل میگردد. پس از تائید مدل پایا، به بررسی مدل نا پایا و مقایسه آن با دادههای تجربی پرداخته می‌شود.
بر اساس آزمایشهای اولیه انجام‌شده، در شکل 2 تفاوت میان سرعت تبخیر سوپر جاذب و پرلیت به‌خوبی مشخص است. با انجام آزمون T برای تعیین اختلاف جفت داده‌ها بر روی نتایج حاصل از آزمایش، بافاصله اطمینان 95 درصد اختلاف بین درصد کاهش آب دو بستر در محدوده 672/5- 781/1 قرار داشت.
شکل 2: نمودار تفاوت درصد کاهش آب پرلیت و سوپر جاذب
شکل 3 نتایج دادههای آزمایشگاهی حاصل از دو بیوفیلتر را در دبی L/min7/0 و غلظت ورودی g/m3 2 نشان میدهد. زمان‌ماند ظاهری (EBCT) در بیوفیلترها برابر 40/1 دقیقه محاسبه شد. همان‌طور که گفته شد فراهم نمودن محلول معدنی تا روز 25 ادامه داشته است. پس از تائید شرایط پایا در دو بیوفیلتر، حذف آبدهی بستر، منجر به شرایط نا پایا در دو بیوفیلتر شده است.
شکل 4 تغییرات بازده حذف دو بیوفیلتر را در مدت‌زمان 40 روز نشان میدهد. با شروع روند کاهش رطوبت در بستر بیوفیلترها، درصد حذف نیز کاهش می‌یابد که در شکل، اثر حضور سوپر جاذب در بهبود عملکرد بیوفیلتر به‌خوبی مشخص است.
شکل 3: تغییرات غلظت خروجی (دادههای آزمایشگاهی)
دبی L/min7/0 و غلظت ورودی g/m32
شکل 4: تغییرات بازده حذف بیوفیلترهای آزمایشگاهی با زمان
شکل 5 پیشبینی دینامیکی تغییرات غلظت خروجی در بیوفیلتر با بستر پرلیت و سوپر جاذب (دادههای آزمایشگاهی شکل 3) توسط مدل با در نظر گرفتن ترم تغییرات رطوبت بستر را ارائه میدهد. تفاوت در پیشبینیهای مدل و دادههای تجربی در شروع شرایط نا پایا بسیار کمتر است (1/0% خطا) اما رفته‌رفته در عمل به دلیل تغییر در برخی پارامترها نظیر سطح ویژه بیوفیلم با کاهش رطوبت بستر این تفاوت بیشتر نمایان می‌شود. میزان خطای میانگین برای پیشبینی این داده‌ها 45/0 درصد است.
شکل 5: مقایسه مدل و دادههای آزمایشگاهی بیوفیلتراسیون
نا پایا با بستر حاوی سوپر جاذب و پرلیت
بحث و نتیجه‌گیری
با توجه به نتایج همان‌طور که در فرضیات مدل ذکر شد برای وارد شدن اثر مقاومت درونی سوپر جاذب درون معادلات، از ضریب مشخصی در معادلات مدل استفاده شد. این ضریب سبب اندکی کاهش نرخ تبخیر نسبت به حالت اشباع سطح می‌شود. با توجه به نتایج شکل 2، این ضریب 95/0 تعیین شد (میانگین اختلاف دو نمودار 5 درصد در نظر گرفته شد). همچنین در بررسی شکل 3 مشخص می‌شود تحت شرایط نا پایا، غلظت در خروجی به‌تدریج زیاد می‌شود. حضور سوپر جاذب در ابتدا مقاومت خوبی را برای حفظ شرایط پایا تا سه روز اول نشان میدهد. بعد از روز سوم نیز سوپر جاذب روند افزایش غلظت را کندتر کرده و مانع تغییرات شدید در عملکرد بیوفیلتر می‌شود که با نتایج به‌دست‌آمده در تحقیقات گذشته مطابقت دارد (8).
بررسی شکل 5 نشان می دهد که تغییرات غلظت خروجی با نرخ ثابتی رو به افزایش است. این امر به دلیل افزایش تدریجی تخلخل بستر و اثر این افزایش روی کمیت‌های مختلف است. در این شکل، تطابق خوبی بین دادههای پیش‌بینی‌شده توسط مدل و دادههای تجربی حاصل از آزمایش مشاهده می‌شود. هرچند ثابت فرض کردن مقادیر پارامترهای مدل در طول زمان نا پایا، منجر به اندکی تفاوت در نتایج شده اما روند افزایش غلظت در طول زمان به‌خوبی مدل شده است.
پس از قطع جریان آب یک بیوفیلتر، به دلیل انجام واکنش بیولوژیکیِ تخریب آلاینده که از نوع واکنش‌های گرمازا است، سرعت تبخیر درون بستر زیاد می شود. رفته‌رفته با تبخیر رطوبت و خشکی بستر، از فعالیت میکروارگانیسم‌های مسئول تخریب ماده آلاینده کاسته شده و غلظت خروجی شروع به افزایش می‌کند. هوای تازه با رطوبت کمتر از اشباع، به‌آرامی قسمتهایی از بستر را خشک نموده و از فعالیت می‌اندازد. این مسئله در مدل با استفاده از gε به‌صورت محلی درون بخش دیفرانسیلی به طرز جالبی در نظر گرفته شد. به‌طوری‌که هر بخش مقدار آب مشخصی از کل آب در دسترس بیوفیلتر را در اختیار دارد و با نرخ انتقال جرم مشخصی آب را از دست میدهد. این روند کاهش، بر روی gε محلی اثر می‌گذارد. برخی ترم‌ها در معادلات انتقال جرم نظیر ضرایب معادلات غلظت آلاینده در فازها، سطح ویژه سوپر جاذب که آن‌هم وابسته به شعاع ذرات است (سوپر جاذب با از دست رفتن آب ذخیره‌شده درون آن کوچک می‌شود) و ضریب انتقال جرم در معادله موازنه رطوبت وابسته به تغییرات gε هستند. تغییر gε باعث تغییر در این ضرایب می‌شود و اثر نهایی آن، افزایش غلظت‌ها در آن بخش است. به دلیل آنکه هر بخش، ورودی خود را از بخش قبلی می‌گیرد اثر افزایش غلظت در بخش‌های اولیه با انتقال به بخش‌های دیگر، در خروجی مشاهده می‌شود. پس از به اتمام رسیدن آب درون یک بخش دیفرانسیلی، روند کاهش آب در بخش بعدی آغازشده و به‌این‌ترتیب تا به پایان رسیدن تمام محتوای آب در بستر ادامه می‌یابد. روند افزایش، درنهایت به ازکارافتادن کامل سیستم منجر می‌شود. در این شرایط، غلظت ورودی آلاینده بدون هیچ تغییری در خروجی بیوفیلتر مشاهده می‌شود. معیار مشخص نمودن این شرایط جهت قطع اجرای کدهای برنامه MATLAB، عدم‌تغییر و برابری غلظت‌ها در 10 تکرار آخر است. در حل معادلات نا پایا اگرچه خصوصیات فیزیکی مواد با تغییر رطوبت ثابت فرض شد اما با تعریف دو ضریب در معادلات، اثر کاهش رطوبت اعمال شد. مقادیر این ضرایب با نسبت میزان آب موجود در آن بخش به میزان آب اولیه تعیین شد. چون باگذشت زمان نرخ انتقال آب تحت اثر مویینگی در خلل و فرج پرلیت و مقاومت داخلی سوپر جاذب قرار می‌گیرد، این ضرایب نرخ انتقال آب از سوپر جاذب و پرلیت و میزان نفوذ آلاینده و واکنش درون بیوفیلم را کاهش می‌دهند. این امر بدان دلیل است که رفته‌رفته با کاهش محتوی رطوبت درون بستر، میکروارگانیسمها از بین رفته و فعالیت خود را از دست می‌دهند. درنتیجه عملکرد بیوفیلتر به‌طور مداوم کاهش می‌یابد تا درنهایت عمل حذف آلاینده درون بیوفیلتر به‌کلی قطع میگردد. مدل این پژوهش توانست به‌طور مناسبی شرایط نا پایای حاکم بر بیوفیلتر و روند افزایش غلظت در آن را پیشبینی کند (خطای 45/0 درصد). هرچند دادههای آزمایش‌ها تا رسیدن به غلظت ورودی موجود نیست، مدل این روند را پیشبینی می‌کند و درعین‌حال تخمینی از تعداد روزها (23 روز) تا ازکارافتادن کامل بیوفیلتر را ارائه می‌کند. تحت شرایط عملیاتی بیوفیلتر، از زمان قطع آب دهی تا روز 23 شرایط نا پایا حاکم است و پس‌ازآن غلظت خروجی به غلظت ورودی می‌رسد و بیوفیلتر به‌طور کامل از کار می‌افتد. اثر کاهش رطوبت در بخش قبلی و نتیجه آن کاهش رطوبت فاز گاز در همان بخش، به‌آرامی بر روی بخش بعدی انتقال‌یافته تا درنهایت بخش مذکور با از دست دادن کامل آب در دسترس خود از مسیر محاسبات خارج شود. در این میان، شعاع، سطح و حجم ذرات سوپر جاذب در آن بخش کم شده و مقدار ضریب انتقال جرم افزایش می‌یابد که محاسبات آنها در مدل به‌طور کامل انجام می‌گیرد.
بر این اساس، مدل‌سازی تغییرات رطوبت در بیوفیلتر و اثر آن بر
غلظت خروجی به‌خوبی توانست عملکرد نا پایای بیوفیلتر را توصیف کند. رطوبت در نرخ انتقال آلاینده به بیوفیلم مؤثر است، به‌طوری‌که کاهش رطوبت بستر علاوه بر کاهش نرخ واکنش در بیوفیلم میکروبی، جذب و نفوذ آلاینده در بیوفیلم را نیز تحت تأثیر قرار میدهد. نتایج مدل‌سازی نشان داد در نظر گرفتن ترم رطوبت در شرایط نا پایا میتواند حیطه عملکرد مدل را توسعه دهد و فرایند بیوفیلتراسیون هوای آلوده را با تمام پیچیدگیهای آن توصیف نماید.
تقدیر و تشکر
نویسنده این مقاله، آقای سروش دانایی (دانشجوی مقطع دکتری) از همکاری خانم مهندس ندا بدلی در پیاده سازی مدل، تشکر و قدردانی می نماید.
تضاد منافع
نویسندگان این مقاله اعلام می دارند که هیچ تضاد منافعی وجود ندارد.
References
1-Mudliar S, Giri B, Padoley K, Satpute D, Dixit R, Bhatt P. Bioreactors for treatment of VOCs and odours – A review. Environmental Management 2010;91(10):39–54.
2-Khan FI, Ghoshal AK. Removal of volatile organic compounds from polluted air. Loss Prevention in the Process Industries 2000;13(5):27-45.
3-Mudliar S, Giri B, Padoley K, Satpute D, Dixit R, Bhatt P. Bioreactors for treatment of VOCs and odours–a review. Journal of Environmental Management 2010; 91(5):1039-1054.
4-Spigno G, De Faveri DM. Modeling of a Vapor-Phase Fungi Bioreactor for the Abatement of Hexane: Fluid Dynamics and Kinetic Aspects. Biotechnology and Bioengineering 2005;89(3):19-28
5-Arriaga S, Revah S. Mathematical Modeling and Simulation of Hexane Degradation in Fungal and Bacterial Biofilters: Effective Diffusivity and Partition Aspects. Canadian Journal of Civil Engineering 2009;36(19):19-25.
6-Morales M, Hernaändez S, Cornabeä T, Revah S, Auria R. Effect of Drying on Biofilter Performance: Modeling and Experimental Approach. Environment and Science Technology 2003;37(9):85-92.
7-Fazaelipoor MS. Modeling temperature variations and moisture requirements in waste air biofilters under steady state conditions. Environmental Technology 2012; 33(5):7-13.
8-Akbari M.The use of water super absorbent inbiofiltration of organic-contaminated air [MSc. Thesis]. Kerman, Shahid Bahonar university, 2009. [Persian]
9-De Visscher A, Li GQ. Toluene removal biofilter modeling: Optimization and case study. Process Safety and Environmental Protection 2008; 86(4):277-282.
10-Fazaelipoor, MH. Analysis of a dual liquid phase biofilter for the removal of hydrophobic organic compounds from airstreams. Chemical Engineering Journal 2009;147(2): 110-116.
ارسال پیام به نویسنده مسئول
| بازنشر اطلاعات | |
 |
این مقاله تحت شرایط Creative Commons Attribution-NonCommercial 4.0 International License قابل بازنشر است. |
