Davlat universiteti қозон миллий тадқИҚотлар технология университети

V BOB. KIMYOVIY REAKSIYALAR KINETIKASINI MATEMATIK MODELLASH

1. Kimyoviy kinetikaning asosiy tushunchalari

Kimyoviy reaksiyalarning tezliklari haqidagi ta‟limotga kimyoviy kinetika deb ataladi. Kimyoviy kinetika fan sifatida XIX asrning 50-70 yillarida shakllana boshladi.

Kimyoviy reaksiyaning tezligi - vaqt birligida (gomogen reaksiyalar reaksiya uchun) yoki sirt birligida (geterogen reaksiyalar uchun) kimyoviy reaksiyada ishtirok etayotgan biror modda konsentratsiyasining o„zgarishi. Shunga ko„ra reaksiya tezligi quyidagicha yoziladi:

• 
  gomogen reaksiya uchun
  

• 
  geterogen reaksiya uchun
  

bunda reaksion faza hajmi, m³; modda miqdori, mol;

reagent konsentratsiyasi, mol/m³; vaqt, s; katalizator sirti, m². (5.58) ifodani quyidagi ko„rinishda yozish mumkin:


85


( )


Doimiy hajmda boruvchi reaksiya uchun (5.4) tenglamadagi ikkinchi had nulga teng bo„ladi, u holda reaksiya tezligi:

Reaksiya hajmining har bir nuqtasida reaksiyaga kirishuvchi moddalarning konsentratsiyalari reaksiya vaqtida vaqt o„tishi bilan doimiy ravishda o„zgarib turadigan reaktorlar uchun kimyoviy reaksiya tezligi hajm birligida modda miqdorining vaqt birligi ichida o„zgarishidir:

yoki geterogen katalitik reaksiya uchun sirt birligi ichida:

bunda reaksion aralashmadagi chi miqdori, mol; reaksion aralashma hajmi yoki katalizator qavati (reaktor hajmi), m³;

katalizatorning solishtirma sirti, m²/m³.
Uzluksiz ishlaydigan reaktor uchun to„liq siqib chiqarish unda urnatilgan rejimda modda konsentratsiyasi apparat uzunligi bo„ylab uzluksiz o„zgarganda kimyoviy reaksiya tezligi reaktor orqali o„tadigan modda mol miqdorining reaksiyaga kirishadigan hajm birligi vaqt birligining o„zgarishidir:



 
bunda reaksion aralashmadagi chi komponentning mol sarfi, mol/s;  reaksion aralashmani berishning hajmiy tezligi, m³/s;  kontakt vaqti, s.
To„liq aralashtiruvchi uzluksiz harakatdagi reaktor uchun o„rnatiladigan rejimda:



bunda reaksion aralashmadagi chi komponentning boshlang„ich va oxirgi miqdori, mol.
86
Amalda odatda reaksiya tezligi mol konsentratsiyada o„lchanadi:

• 
  davriy harakatdagi reaktor uchun:
  

• 
  uzluksiz harakatdagi reaktor uchun:
  

_
 
Agar reaksiya hajmning o„zgarishi bilan bormasi, u holda ideal siqib chiqarish reaktori uchun:









(_)




_



Ideal almashuvchi uzluksiz harakatdagi reaktor uchun:

 



_



bunda reagentning mahsulotga aylanish darajasi:

 o„rtacha bo„lish vaqti, s:

 _
1862-1867 yillarda Norvegiya olimlari Guldberg va Vaage
massalar ta‟siri qonunining dastlabki ta‟rifini berdilar. Massalar ta‟sirning kinetik qonuniga ko„ra, berilgan temperaturada elementar reaksiya tezligi ularning stixiometrik koeffitsiyentlariga teng darajalarda olingan reaksiyaga kirishuvchi moddalar konsentratsiyalarining ko„paytmasiga proporsional.
U holda quyidagi kimyoviy reaksiya o„tganda

 
 

Massalar ta‟siri qonuniga ko„ra reaksiya tezligi quyidagicha yoziladi:
 


87

bunda – chi komponentning konsentratsiyasi; tezlik konstantasi;   moddalarning stexiometrik koeffitsiyentlari.
(5.11) tenglama elementar reaksiyalar uchun o„rinli.
Alohida olingan komponentlar bo„yicha reaksiya tezliklari
bilan umumiy tezlik o„rtasida quyidagi stexiometrik nisbat mavjud:











Murakkab kimyoviy reaksiyaga massalar ta‟siri qonunini qo„llash uchun uni elementar bosqichlar ko„rinishida ifodalash va bu qonunni har bir bosqichga alohida-alohida tatbiq etish lozim.
Kimyoviy kinetika XIX asrning 80-yillarida Vant-Goff va Arrhenius asarlarida to„liq shakllantirildi; reaksiya tartibining ma‟nosi aniqlandi, faollanish energiyasi, mono-, bi- va polimolekulyar reaksiyalar tushunchalari kiritildi.
Vant-Goff va keyinchalik Arrhenius o„z g„oyalarini rivojlantirgan holda, "temperatura reaksiyaning sababi emas, temperatura reaksiya tezligining o„zgarishiga sabab bo„ladi" deb tasdiqlashdi. Unga ko„ra:

bunda eksponsialoldi ko„paytuvchisi; faollanish energiyasi, j/(mol K); universal gaz doimiysi, j/(mol K); temperatura, K.

Tashqi sharoitlar (temperatura, bosim, tarkib, reaksiya boradigan muhit) berilganda reaksiya tezligi reaksiyaga kirishuvchi moddalarning funksiyasidir, ya‟ni:

(5.4) tenglamani (5.14) ga solishtirib, quyidagi kinetik tenglamani olamiz:

Kimyoviy reaksiya davomida modda konsentratsiyasining vaqt bo„yicha o„zgarishini ifodalovchi tenglamaga kinetik tenglama, egri chiziqqa kinetik egri chiziq deyiladi.

88

2. Gomogen kimyoviy reaksiyalar kinetikasini modellash

Kinetik tenglamalar reaksiya tezligini bog„liq bo„lgan parametrlarga, masalan: konsentratsiya, harorat, bosim, katalizator faolliklariga bog„laydi.

Kimyoviy reaksiyalarning kinetik mexanizmlarini matematik modellash yo„li bilan o„rganish reaksiyaga kirishuvchi moddalar konsentratsiyalarining vaqt bo„yicha o„zgarishini aniqlashdan iborat.
Quyidagi kimyoviy reaksiyalar sodir bo„lsin:
Massalar ta‟siri qonuniga ko„ra kimyoviy reaksiyalarning tezliklari yoziladi va kinetik model tuziladi:






}
bunda moddalar konsentratsiyalari, mol/l.
(5.16) tenglamalar birinchi tartibli oddiy differensial tenglamalar sistemasi bo„lib, oddiy differensial tenglamalarni yechishning sonli usullari (Eyler, Runge-Kutti) yordamida yechish mumkin.
5.1-, 5.2-rasmlarda boshlang„ich moddaning o„zgarish darajasiga va moddalar konsentratsiyasiga haroratning ta‟sirini matematik model (5.16) yordamida o„rganish natijalari keltirilgan.
Olingan natijalar maqsadli mahsulotni olish uchun jarayonning maqbul vaqti haqida xulosa chiqarishga imkon beradi.
Matematik model (5.16) ham xomashyo tarkibining reaksiya mahsulotlari unumiga ta‟sirini o„rganishga imkon beradi.
Ma‟lumki, kimyoviy reaksiyaning tezligi temperaturaga bog„liq, shuning uchun jarayonning borishiga temperaturaning ta‟sirini o„rganish
89
uchun Arrhenius tenglamasi (5.13) yordamida kinetik modelda (5.16) kimyoviy reaksiyaning tezlik konstantasini ifodalash kerak.

5.1-rasm. Vaqt o„tishi bilan reaksiyala ishtirok etgan moddalar konsentratsiyasining o„zgarishi

5.2-rasm. Reagentning mahsulotga aylanish darajasining temperaturaga bog„liqligi

90

3. Geterogen kimyoviy reaksiyalar kinetikasini modellash

Geterogen kimyoviy reaksiyalar kinetikasining asoslari Langmuir, Temkin va boshqalar ishlarida keltirilgan. Ushbu ishlarda gomogen kinetika qarashi bilan o„xshash ideal adsorbsion qavat tushunchasi shakllantirilgan. Bu model quyidagi taxminlardan foydalanishni nazarda tutadi:

katalizator sirti barcha bo„limlarining bir xilligi va ximisorbsiya energiyasining sirtni turli adsorbentlar bilan to„ldirish darajasiga bog„liqmasligi;
katalizatorning o„zgarmasligi va uning xossalarining reaksion aralashmaning tarkibi va uning katalizatorga ta‟siridan o„zgarmasligi;. energiyaning teng taqsimlanishi.
Qattiq sirtlardagi elementar jarayonlar uchun massalar ta‟siri kinetik qonunining formal analogi ta‟sir etuvchi sirtlar qonuni hisoblanadi.
Kimyoviy reaksiya tezligining daslabki ta‟rifiga ko„ra kimyoviy reaksiya tezligi reaksiyaga kirishayotgan moddalar sirt konsentratsiyalarining stexiometrik nisbatlar darajasiga olingan qo„paytmaga to„g„ri proporsional:
<sup>     

</sup> 
∏

bunda tezlik doimiysi; adsorbsiyalangan zarracha bilan
band bo„lgan sirt ulushi; erkin sirt ulushi;  bosqichlarning stexiometrik koeffitsiyentlari;  kimyoviy reaksiya o„tganda mollar sonining o„zgarishi.
Quyidagi elementar reaksiya o„tadigan bo„lsin:
   
Bunda barcha moddalar adsorbsiyalangan modda bilan ta‟sirlashgan holatda reaksiyaga kirishadi.
Adsorbsiyalangan chi modda bilan bo„lgan sirt ulushini bilan belgilaymiz. U holda sirt ta‟siri qonuniga binoan (5.18) qaytmas reaksiya uchun tezlik ifodasi quyidagicha yoziladi:
∏ ^{   }


91
Agar barcha moddalar adsorbsiyalangan holatdan ta‟sirlashmasa, faqat gaz fazadan to„g„ridan-to„g„ri ta‟sirlashsa, u holda sirt ta‟sirlari qonunini umumiy holda quyidagicha ifodalash mumkin:
∏ ^ ∏ ^

bunda gaz fazadan reaksiyaga kirishayotgan moddaning parsial bosimi (konsentratsiyasi); katalizator sirtida va gaz fazadan reaksiyaga kirishadigan moddalar sonlari.

Masalan, Z katalizatorning faol markazida vodorodning
adsorbsiyasi adsorbirlangan sirtda ZH₂ birikmasini hosil qilish bilan borgan bo„lsin:

u holda sirt ta‟sirlari qonuniga asosan berilgan elementar reaksiyaning tezligini quyidagicha ifodalash mumkin:

Kinetik bog„liklikni aniqlovchi asosiy omil sifatida dastlab siqib chiqarish omili – katalizator sirtidagi joy uchun reaksion aralashma komponentlarining “kurashi” qaraldi.
Bunda haqiqiy kimyoviy o„zgarishga qaraganda adsorbsiya va desorbsiya bosqichlarining yuqori tezlikda borishi haqida qo„shimcha taxmin qabul qilingan.
Keyingi tadqiqotlar ushbu taxminlarning sezilarli cheklovlarga egaligini ko„rsatdi. Shunga qaramay Xinshelvud, Shvab, Xougen, Uotson va boshqalar ma‟lum bir parametr o„zgarishlari oralig„ida kinetik tajribani qoniqarli tasvirlovchi tenglamalarni olishdi.
Bu taxminlarga mos kinetik tenglamaning umumiy ko„rinish quyidagicha:

^
( ∏ )

( ∑ )

bunda tezlik konstantasi; gaz fazadagi -chi reagentning
konsentratsiyasi; -chi komponent adsorbsiya bosqichining muvozanat konstantasi;  stexiometrik koeffitsiyentlar.


92

Misol

Toluolning girokrekingi murakkab geterogen kimyoviy reaksiyani qaraymiz.
Kimyoviy reaksiyaning mexanizmi:

1. 
   

→

←

2. 
   

3. 
   

→

←



bunda katalizator sirtidagi faol markazlar; va boshqalar adsorbirlangan oraliq mahsulotlar.
Sirt ta‟siri qonuni bo„yicha reaksiyaning elementar bosqichlari uchun reaksiya tezligini yozamiz:





Ushbu kimyoviy jarayonning matematik modeli kuzatilayotgan moddalar va oraliq birikmalar konsentratsiyalarining vaqt bo„yicha o„zgarishini ifodalovchi differensial tenglamalar sistemasidan iborat, ya‟ni:


93

(5.22) differensial tenglamalar sistemasini yechishda Eyler va Range-Kuttilarning sonli usullaridan foydalanish mumkin.

4. Geterogen kimyoviy reaksiyalarning kinetik modellarni tuzish usullari