2. Kineetika rakendused simulatsioonis
2.3. Keemiliste reaktsioonide rakendamine simulatsioonis
Erinevad võimalused, kuidas näha mingi kindla ainete kogumiku reaktsioone. Näidiseks on valitud g = GRI 30 (53 liikide kogumik). Selleks, et näha konkreetselt viiendat reaktsiooni tuleb sisestada käsk:
>> e8 = reactionEqn(g,5)
Vastuseks tuleb:
e8 =
'H2O2 + O <=> HO2 + OH' Samuti võib valida reaktsioone ka mingis kindlas vahemikus, näidis (GRI 30 reaktsioonid vahemikus 1-8).

Joonisel on näha GRI30 kogumiku elementide reaktsioonid vahemikus 1- 8 , juhul kui on soov näha reaktsioone näiteks vahemikus 1-20, siis koodi esimeses reas tuleb kirjutada mitte 1:8, vaid 1:20.
Selleks, et saaks näha kõikide ainete omavahelist reaktsioone tuleb sisestada käsk:
>> eqs = reactionEqn(g)
Programm väljastab tulemuse:

Programmi poolt väljastatud tulemus (esimesed 15 reaktsiooni, kokku näitab 325).
Vedelike/gaaside ülekandeomadused
Selles osas on toodud välja erinevad funktsioonid, mida kasutatakse transportimise omaduste välja arvutamiseks. On oluline teha kindlaks millist elementide kogumiku hakatakse katsetama. Praegu kasutatakse kahte põhilist mudelit ning mõlemad põhinevad kineetilise teooria väljenditel. Esimeses mudelis on mitmekomponendiline koostis ning teine mudel on lihtsustatud, kuna seal kasutatakse ainult 1-3 liikide segudega seotud väljendeid, et keskmistatud koostis.
Esimene mudel:
>> g1 = GRI30('Multi')
Teine mudel:
>> g2 = GRI30('Mix')
Omadused, mida saab leida mõlemal mudelil:
Viskoossus
>> visc = [viscosity(g1), viscosity(g2)]
visc =
9.0002e-06 9.0002e-06
Soojusjuhtivus
>> lambda = [thermalConductivity(g1), thermalConductivity(g2)]
lambda =
1.8691e-01 1.8679e-01
Binaarsed difusioonikoefitsiendid
>> bdiff1 = binDiffCoeffs(g1)
>> bdiff2 = binDiffCoeffs(g2)
Näidis, kuidas antud kood näeb välja MatLabi käsuaknas:

Esimese mudeli binaarsed difusioonikoefitsiendid.
Segukeskmised difusioonikoefitsiendid
>> dmix2 = mixDiffCoeffs(g1)
>> dmix1 = mixDiffCoeffs(g2)
MatLabi käsuakna vaade:

Esimese mudeli segukeskmised difusioonikoefitsiendid.
Omadused, mida saab leida ainult esimesel mudelil:
Mitmekomponendilised difusioonikoefitsiendid
>> dmulti = multiDiffCoeffs(g1)

Esimese mudeli mitmekomponendilised difusioonikoeditsiendid.
Käsk multiDiffCoeffs(g1) tagastab mitmekomponendilise difusioonikoefitsientide maatriksi, mille iga element Dij, kirjeldab liigi i difusioonivoolu tundlikkust liigi j kontsentratsioonigradientide suhtes. Maatriks on suurusega liikide arv × liikide arv (nt 53×53 GRI30 korral) ja annab täieliku Maxwell–Stefani mudeli koefitsiendid. Seega tabel ei näita üksikuid difusioonikoefitsiente, vaid kõigi liikide omavahelisi difusioonisuhteid.
Termilise difusiooni koefitsiendid
>> dt = thermalDiffCoeffs(g1)

Esimese mudeli termilise difusiooni koefitsiendid.