اگر حبابهای گازی و ذرات سیلیکاتی ذوب نشده در مذاب باقیمانده باشد کیفیت محصولات شیشه افت میکند.
لذا هر چه عیوب ناشی از وجود حبابها در محصولات شیشه کمتر باشد، بهتر است.
برای مثال حضور 6 حباب در هر تن شیشه تلویزیون موجب به هدر رفتن 10 درصد محصولات نهائی میگردد و اینگونه عیوب در خیلی از موارد دیگر مثل شیشه اتومبیل، ساختمان، ظروف و غیره میتواند اثرات خیلی زیانآور اقتصادی داشته باشد، لذا عمل حبابزدائی یکی از اساسیترین کارها در صنعت تولید شیشه میباشد.
حرارت در کورههای ذوب شیشه به صورت تابشی از بالای محفظه مذاب تأمین میگردد، اگر حبابها بالا بیایند و در سطح مذاب جمع شوند، تشکیل لایهای کف حبابی روی سطح مذاب میدهند و باعث جلوگیری از انتقال حرارت تابشی به عمق مذاب میشوند و اثرات خیلی مخربی در این حالت در راندمان انرژی، کیفیت شیشه، توید، انتشار آلودگی موجب میشوند.
طی پروسه ساخت شیشه، بچ ترکیبی شیشه وارد کوره مذاب شده و در آن جا به سبب مشعلهای حرارتی تمامی ترکیبات شیشه ذوب، و در اثر جریانات کنوکسیونی عمل یکنواختی مذاب انجام میگیرد.
عمل ذوب مواد خام شیشه یک پروسه فیزیکوشیمیائی خیلی پیچیدهای است که شامل تعداد زیادی از واکنشهای شیمیائی و تبدیلات فازی در محدوده دمائی 800 الی 1200 میشود.
برای مثال در تهیه شیشه و ظروفی واکنشهای شیمیائی پایهای زیر در بچ حاوی سیلیس، کربنات سدیم و کربنات کلسیم انجام میگیرد.
 محدوده دمای 550 C 600-83C 720-900C 600-900C گاز Co2 در نتیجه سه واکنش فوق در این نوع ترکیب شیشه تولید میشود، در مذاب نفوذ میکند.
کسر کوچکی از این گاز در جوانهزنی غیریکنواخت حبابهای گازی درون یا زیربچ جوانهزنی میکند.
برخی از این حبابها از درون بچ عبور و به فضای احتراق وارد میشود در صورتی که بقیه در مذاب گرفتار میشوند و تحت جریانات کنوکسیونی مذاب حمل میگردند.
عوامل حبابزدا که در واکنشهای اکسیداسیون و احیای تعادلی در تولید یا و صرف گازها شرکت میکنند به منظور حذف حبابهای نامطلوب در مذاب شیشه به ترکیب افزوده میگردد.
سه نوع از عوامل حبابزدا مورد استفاده قرار میگیرد: 1- اکسیدهای فلزی با ظرفیتهای متغیر به منظور ایجادگاز اکسیژن: Sb2O3، AS2O3، CeO2.
2- سولفات و سولفتیها که مخلوطی از گاز اکسیژن و So2 تولید میکنند.
3- کلرید، برمید و یدیدها که در دمای تصفیه بخار میشوند.
در محدوده دمائی بالا تعادل واکنشهای شیمیائی مربوط به حبابزدائی منجر به تولید گاز میشود و گاز ایجادی درون شیشه مذاب به درون حبابهای گازی موجود در مذاب نفوذ میکند.
گازهای موجود در درون حبابها با گازهائی که از این طریق وارد حبابها میشوند، مخلوط میگردند که به نوبه خود نفوذ گازها از مذاب به درون حبابهای رشد یابنده را افزایش میدهد.
نفوذ گازها به درون حبابها تا زمانی ادامه مییابد که حباب تا اندازهای بزرگ شود که بتواند به سطح مذاب حرکت کند.
در دماهای پائین، تعادل واکنشهای شیمیائی مربوط به حبابزدائی در جهت مصرف گازها پیش میرود و منجر به این میشود که گازهای موجود در حبابها به درون مذاب نفوذ میکند.
از این طریق حبابهای کوچک که هنوز نتوانستهاند به اندازهای رشد کنند که بتوانند خارج شوند، در مذاب تجزیه شده و از بین میروند.
II.
تئوری 1- مکانیسم نفوذ اکسیژن انتقال اکسیژن از درون مذاب شیشه توسط محققین زیادی مورد مطالعه قرار گرفته است.
نفوذی ماشینی Seld-diffusion یا نفوذ بین نشینی Inter diffusion مورد بررسی قرار میگیرد.
نفوذ جانشینی به این ترتیب است که یک یون یا یک مولکول در یک جهت در محلول توسط رابطه  به پیش میرود به طوری که  میانگین مربع جابجائی عامل نفوذکننده پس از زمان t میباشد که میتواند به طور آزمایشی با مطالعه جایگزینی O16 با ایزوتوپهای O17 و O18 اندازهگیری شود.
احتمالاً یونهای اکسیژن غیرپلزن در شبکه با یونهای نشاندار (ایزوتوپها) تعویضی میگردند.
نفوذ بین نشینی به این ترتیب است که اکسیژن مولکولی یا اتمی در ساختار مذاب از طریق درز و شکاف حرکت میکند و به نفوذ فیزیکی یا مولکولی اکسیژن مشهور است.
Rund و همکارانش نفوذ اکسیژن مولکولی را مطابق زیر مورد مطالعه قرار دادهاند.
اغلب اکسیدهای فلزی با ظرفیت متغیر به عنوان عامل حبابزدا به مذاب شیشه افزوده میشوند: ،  و .
اکسیژن مولکولی از طریق واکنشهای اکسیداسیون و احیا یونهای با ظرفیت متغیر، تولید یا مصرف میشوند و شکل عمومی این واکنشها به صورت زیر نوشته میشود: (1-1) ثابت تعادل Kp برای واکنش فوق یا شرط ثابت بودن فعالیت یون اکسیژن در مذابمطابق زیر نوشته میشود: (1-2) مقدار Kp بستگی به ترکیب شیشه و دما ارد.
مطابق قانون هندسی غلظت اکسیژن حل شده به صورت مولکولی در حالت تعادل، برای حالتی که پیوندهای مولکول O2 شکسته نشده باشد، مطابق رابطه زیر به فشار محیط ارتباط پیدا میکند: (1-3) در این معادله حلالیت اکسیژن مولکولی یا فیزیکی است.
نسبت غلظت ناحیهای اکسیژن محلول به صورت مولکولی به ضریب هنری (L)، فشار اکسیژن تعادلی در ناحیه معین از مذاب را بدست میدهد که میتواند توسط سنسورهای اکسیژن تعیین گردند.
حلالیت فیزیکی گازهای نادر و نیتروژن در شیشههای سودار لایم به ترتیب 10-6 تا 10-5 مول بر مترمکعب در فشار یک پاسکال میباشد و به مقدار جزئی به دما بستگی دارد.
با استفاده از دو رابطه اخیر، ثابت تعادل Kc بر پایه غلظت نمونههای مربوطه میتواند مطابق زیر تعیین شود.
(1-4) مقادیر KP و KC بستگی به دما و میزان قلیائیت سیستم (ترکیب شیشه) دارد.
در مذابهای شیشه تعادل شیمیائی موضعی واکنشهای اکسیداسیون و احیا نسبت به سرعت نفوذ پائین اکسیژن، خیلی سریعتر حاصل میشود و کنترل کننده واکنشهای اکسیداسیون و احیا را نفوذ اجزاء کنترل میکند.
معادله دیفرانسیل انتقال جرم از طریق نفوذ در یک بعد بدون کنوکسیون اما با واکنش شیمیائی در حال تعادل به صورت زیر داده میشود: (1-5) اولین جمله سمت راست دانسیته جریان جرمی در فاصله (dx) را نشان میدهد و تولید به مصرف موضعی اکسیژن درون مذاب توسط واکنشهای اکسیداسیون و احیا در واحد زمان و واحد حجم با P(t) داده میشود.
مطابق واکنش شیمیایی اکسیداسیون و احیا به ازای هر مول از یونهای MK+ تولیدی،  مول از اکسیژن مصرف میشود:  (1-6) از نفوذ یونهای با ظرفیت متغیر در مذاب صرفنظر میشود و از مقایسه اطلاعات آزمایشی نفوذ یون فلزی با مقادیر نفوذ اکسیژن مولکولی در مذابهای شیشه این فرض معقول است.
با توجه به اینکه مجموع غلظت یونهای Mk+ و Mn+ اجزای با ظرفیت متغیر طی پروسه نفوذ ثابت میماند.
(1-7) نتیجه زیر با استفاده از روابط بالا حاصل میگردد: (1-8) با دیفرانسیل گرفتن از رابطه فوق نسبت به زمان و جایگذاری در رابطه (1-6) معادله دیفرانسیل غلظت اکسیژن محلول به صورت فیزیکی بدست میآید و فقط میتواند به روش عددی حل گردد.
(1-9) (1-10) (1-11) DC هم بعد با D میباشد.
تغییرات غلظت اکسیژن نسبت به زمان بستگی به غلظت یونهای با ظرفیت متغیر Cm، ثابت تعادل KC وابسته به دما و غلظت اکسیژن در مذاب دارد.
برای اینکه DC به غلظت موضعی اکسیژن وابسته و نیز تابعی از زمان و مکان میباشد.
مقادیر KC در مذابهای شیشه سودالایم برای تعادل  از دمای 1200 تا 1400، از 1 تا 10  تغییر میکند.
غلظت اکسیژن محلول به صورت فیزیکی [O2] از 1/0 تا 1 دارد.
برای و k=5 و cn=3 نسبت  مقادیری بین 01/0 تا 1 دارد.
برای اینکه Cm، KC، [O2] و (k-n)2 در معادله (11) مقادیر مثبتی میباشند، DC همیشه پائینتر از D میباشد.
برای مذاب عاری از اجزاء با ظرفیت متغیر Cm = 0 و D = DC میباشد.
برای دماهای پائین،  خیلی بیشتر از یک میباشد لذا  و بدین معنی است که واکنش اکسیداسیون و احیا با مشکل پروسه نفوذ O2 در دماهای پائین رات حت تأثیر قرار میدهد.
در دماهای بالا  در مقایسه با 1 خیلی کوچک است، در این حالت مخرج معادله فوق نزدیک به یک میباشد و مقدار DC تقریباً همان D میباشد.
چگالی جریان جرمی اکسیژن از اتمسفر گازی به درون مذاب مطابق زیر میباشد:  x فاصله از سطح مذاب میباشد و x = 0 سطح مذاب شیشه را نشان میدهد و J(O2) چگالی جریان جرمی اکسیژن به درون مذاب میباشد.
دانسیته جریان اکسیژن از صفحه عمود به جهت انتقال جرم با  داده میشود.
مقدار J(O2) برای پروسههای احیائی منفی میباشد.
به خاطر اینکه DC همیشه از D کوچکتر میباشد، مقدار گرادیان غلظت اکسیژن محلول مولکولی طی فرایندهای اکسیداسیون و احیا در مذابهای حاوی یونهای با ظرفیت متغیر نسبت به مذابهای عاری از یونهای فوق بزرگتر است.
همانطوری که در شکل نشان داده میشود تغییرات منحنی غلظت نسبت به فاصله برای مذابهای حاوی یونهای با ظرفیت متغیر سریعتر است.
تفاوت بین دو منحنی به خاطر واکنش شیمیایی اکسیژن محلول فیزیکی با یونهای با ظرفیت متغیر در مذاب ایجاد میگردد.
برای انتقال جرم به درون یک مذاب نیمه محدود برای مقدار ثابت D و بدون واکنشهای شیمیایی در مذاب، چگالی جریان اکسیژن به درون مذاب به صورت زیر تعریف میگردد.
  غلظت اولیه در نمونه و [O2]I غلظت اکسیژن در سطح مشترک مذاب شیشه - گاز برحسب  میباشد.
برای حالتی که واکنشهای اکسیداسیون و احیا وجود دارد که شامل یونهای با ظرفیت متغیر میشود معادله فوق برای مقدار ثابت DC با فاکتور تصحیح  به صورت زیر تعریف میشود: فاکتور  مساوی با عاملی است که به واسطه آن ضخامت لایه مرزی در مقایسه با مذاب بدون یونهای با ظرفیت متغیر کاهش یافته است.
DC در معادله فوق مقدار ثابتی است.
اما DC مطابق معادله (11-1) بستگی به غلظت موضعی اکسیژن محلول فیزیکی دارد و مقدار ثابتی است.
از نتایج محاسبات عددی با استفاده از معادله (10-1) و تعیین چگالی جریانی [O2] به درون مذاب از معادله (12-1)، دانسیته جریانی یا خطای 10% با جایگزینی در معادله (14-1) بوسیله  تخمین زده میشود و مقدار میانگین نسبت در لایه مرزی مطابق زیر است: با جایگزینی [O2]i در معادله (11-1) و  با جایگزینی  در این معادله محاسبه میشود و از اینرو دانسیته جریان اکسیژن به درون یا بیرون از مذاب با تقریب خوبی میتواند بدون عمل محاسبات عددی با معادله زیر حاصل شود.
شیشهای که تحت اتمسفر گازی با ترکیب ثابت ذوب میشود پس از چند ساعت با اتمسفر فوق به تعادل خواهد رسید.
حل فیزیکی گازها درون مذاب تقریباً مستقل از دما میباشد به این معنی که با تغییر دمای مذاب شیشه مجاری از یونهای با ظرفیت متغیر، جایگزینی جزئی اکسیژن توسط اتمسفر رخ میدهد.
حلالیت شیمیائی مولکولهای اکسیژن در مذاب شیشهای که حاوی عوامل حبابزدا میباشد به طور قابل ملاحظهای با دمای مذاب تغییر میکند.
برای اینکه ثوابت تعادل واکنشهای اکسیداسیون و احیا شدیداً وابسته به دما میباشد.
پائین آوردن دما موجب تعادل به سمت حالتهای اکسیداسیون یونهای با ظرفیت متغیر میگردد.
لذا کاهش دما منجر به کاهش غلظت اکسیژن در مذاب میشود و چون که غلظت اکسیژن در سطح مشترک مذاب گاز ثابت میباشد، نفوذ اکسیژن به درون مذاب انجام خواهد گرفت.
همانطوری که قبلاً نیز بیان شده است، زمان نفوذ در مقایسه با زمان واکنش خیلی طولانی است.
تغییر در غلظت اکسیژن محلول فیزیکی در مذاب فوراً پس از تغییر دما مساوی ضرب در تغییر غلظت بخش احیائی اجزاء با ظرفیت متغیر میگردد: 1 مربوط به حالت پس از تغییر به دمای T1 و صفر مربوط به حالت اولیه T0 میباشد.
برای هر دو دمای T0 و T1 تعادل شیمیائی مطابق زیر بدست میآید.
Co2 غلظت اکسیژن، V(t,r) سرعت شعاعی حل شدن حباب،  ضریب نفوذ اکسیژن در مذاب، جمله مربوط به واکنش و r فاصله اندازهگیری شده از مرکز حباب میباشد.
RA شعاع حباب، PaO2 فشار جزئی اکسیژن درون حباب و Lo2 ثابت هنری مربوط به مولکول اکسیژن میباشد.
 حالت اکسیداسیون M با ظرفیت متغیر میباشد و  حالت احیائی یون فوق میباشد.
تغییرات غلظت O2،  و  به همدیگر وابسته میباشد.
بنابراین سرعت مولی تولید اکسیژن در واحد حجم  مطابق زیر داده میشود.
3- فرمولاسیون کلی رفتار یک حباب درون مذاب شیشه: هدف این مطالعه ارائه آنالیز ریاضی توأم با واقعیت از رفتار حباب در کوره مذاب شیشه با استفاده از تئوری بالانس جمعیتی میباشد.
دستهای کامل از معادلات، از جمله معادلات مربوط به غلظت عوامل حبابزدا، غلظت گازهای حل شده در مذاب شیشه و تابع دانسیته حباب توأم با شرایط مرزی مربوطه را ارائه میدهد.
این مدل انتقال کنوکسیونی سه بعدی عامل حبابزدا، گازها، حبابها، به علاوه رشد آنها در نتیجه نفوذ گازهای چندگانه و جوانهزنی حباب روی دیوارههای نسوز و واکنشهای حبابزدائی را در نظر دارد.
به منظور حل ریاضی مسئله تولید و انتقال حباب در محفظه مذاب شیشه فرضیات زیر در نظر گرفته میشود: 1- اثرات حبابها، گازهای حل شده بر سرعت و دما و خواص ترموفیزیکی مذاب شیشه مورد بررسی قرار نمیگیرد.
2- حبابها به شکل کره در نظر گرفته میشوند.
3- حبابها دارای اینرسی قابل صرفنظر میباشند این فرض نظر به اینکه حبابهای خیلی کوچکی در مذاب شیشه وجود دارند، معقول است.
4- اجزاء بردار سرعت حباب همانند مذاب شیشه در