Полак_и_др__Вычислительные_методы_в_химической_кинетике (972296), страница 60
Текст из файла (страница 60)
00 Š— 061.00Š— 041.00 Š— 021. 00 Е+00 Количественньв ограничения. В описываемой версии могут рассчитываться кинетические схемы, содержащие не более 200 реакций и БО компонент. Чтобы рассчитывать более сложные схемы, необходимо произвести соответствующие изменения в операторах описания размерности и общем блоке К/ИЕТв основной программе и подпрограммах. При задании уравнений химических реакций предполагается, что все стехиометрические козффициенты равны единице, таким образом, в кипе. тическую схему могу~ входить реакции 1 -»)г, где К = 1, ..., 6, 2 -' К, где /г = 1, ..., Б, 3 К,гдед=1,...,4.
Описанная программа решения прямой кинетической задачи представлена в версии, предназначенной для ЭВМ БЭСМ-6. Все вычисления проводятся с одинарной точностью, так как БЭСМ-6 имеет сравнительно большую 246 лину машинного слова. Для проведения расчетов на ЭВМ серии ЕС целасобраано использовать версию программы с удвоенной точностью. В предпюленном адесь вариепв программа предназначена для расчета прямой 1отермической кинетической аадачи. Однако путем незначительных измезний в программе к системе уравнений химической кинетики могут быть обновлены уравнения для одной нли нескольких температур, плотности, веления и тл.
Чтобы решать расширенную систему уравнений, наобходио изменить весло уравнений в обрацании к подпрограмме ВТ!ЕЕС, предуиотреть заалние начальных условий для новых переменных, запрогрем. ,ировать новые уравнения в подпрограмме О!ЕЕ(3(Ч, куда следует в зтом 1учае перенести обращение к ССАЙ С нз основной программы. Так как рограммирование вычисления якобиана расширенной системы может казаться затруднительным, то следует перейти к режиму интегрирова(ия автомати вским моленным вы молением якобиана (МЕ 22 или 23) . Кроме зтого, укороченный пакет ВТ!ЕЕ [б7). с монитором ЯТ!ЕЕС ожет ислользоватьса автономно при решении задачи Коши для произвольых жестких систем обынновениыхдифференциальных уравнений.
Ниже приводится текст рассмотренной программы. РВСОВХМ КЗ(," '..., "1-чз РХМЗИВХОИ Х(50) "' с(ммоилшпи~а(гаа),РЗ(гоо),хи(гоо)>за(гоо) С4ЬЬ РЬХИЗ РВХит 2200 ВИАР 1,И М,МЬ,ИРМ,ИР,Ы,ЗРЯ,ЖИ амь 1о (мь РВХИЗ 2201,И РВХЩ 2202,М Рихит 2203эмь РВХИР 2204,ИХМ РВХ)гг 2200,ИР РВХИЗ 2205>В РВХИЙ 2206,ЗРВ РВХЯХ 2207,7И хР(му.оз.го.аир.мт.ьз.гз)сото 17 РВХИР 10 ВЗОР 15 с(аггхиси .7)б г хч;и— 2' Х(Х)~0 0 саьь ВЗ01(и х) саьь Всаьз(1() РВХИР 14 саьь Рс(мс(х,и,амь>хси) РВХИЗ 17 саьь ВВ01 (ихм, 11о сдьь ЗВР(ихм,тм)' саьь ссаьс сы~.рьХЩ. ХХиь1 ИЗ~ИМ(ИХМ) саьь итхРРС(З,ЗР,х,и,из,зРВ,ИР) саьь Рьхиз 1 РСВМИ(ИХ) ЭЗ0.1) 14 иыа(ат(!х.
(Виыльниз исицхх(РРациц' ) .241 !Э Томит(!х. »М(минтн нргми))к н которнж пжчатапэся Риижния ) !в Р(имат( ° Йиниии) В!дана МР() 2200 РОИМАТ(31Х(»к И Н 3 т М К А - 8 1») 2201 РОИМАТ(1Х» ЧИСЛО И(8)П(2(3~Т»» 13) 2202 РОЮ(ат(1х»»числО Риаикий»»»13) 2203 РОМ(ат(1х, 'лсГА)а!Ом НАОМХАВИ(ГО мномигжгя ',13) 2204 Ран(АР(11~'число 'гочик нзчвти»' 13) 2203 РО)ШТ(1Х» 'НАЧАЛЬЙ)И МАГ ИНТИГ)4РСВАНИЯ»»»Ж1 2. 4) 2206 РО)агат(1Х»»ОЭНОСИТИЛЬНАЯ ТОЧНОСТЬ ИНТИГРИРСНАЙИЯ ', И! 2 4 ) 220) РОИМАТ(» ТШЭРРАТХРА»» Р8 1,' Е') ггав РО)ют( ммн)л интигрироню)я °, 13) ВТОР ЗМП .АИЬ,ХСИ) ),Р(10) ПО 282»Э»! 10 Е!»(» 1)»2!)+! Иг»2»МХИО(И»Э»10) мрчггиа(!о, Й-1 о»(8-1 ) ) ВО гвэ 1~1',мю 283 Р(1)»х((д 1)»10+1)/га(ь риэит ггаа,(1см(1] ° 1 Е1,К2) ОАЬЬ Рз!24(ин Р) 12(и.ьз.а 10))етаии 282 С(МТХМЦЗ ггаа РОиюат(10(гх, гаа,гх) ) нвтции жив ВОВВ01)тэик сс ОСИНОМ/е)кий/ О),РЖ(200),ти(200),жА(200) »и(200),ЬИ(1800),Р(200),ЭМ(100) »тк»тии, Аиь» и, м,мь, 1ти» 1си(1 а о) ТЖ)а»! ° 987З Э»ТК ЭОЭО 1!и 30 С(1)»ЗХР(РЗ(1)-ЗА(1)/ТМО»ТК»»ти(1) ИЖТЦИИ МЮ ВЦНИОЦТЭИЗ ВХРР!ИФ(ИИИ Ттт ХО Р) с(гаваи/ЖЭЗЖТ/с(200),РЗ(200)»'хм(200),жА(200) ,м(200),ьи(1800),'Р(гоо),эм(100) ' ОЭИЗИВХОИ ХО(ЭО)»РЬО) ПО 1 1»1»и 1Р(ХО(1) ЬТ.О 0)ХО(1)»0.0 1 Р(1)»о,о К»0 ПО 2 1»1,М т с(1) д! ЬИ(Е+1) ПО Э а 1»а1 ЕО»К+2 М КО»ЬВ(ЕО) Э Т»Т"Хо(КО) И(1)»Т ПО а Ю»1»01 ЕО»к+2+»Э КО»ЬИ(ЕО) 4 Р(ко)»Р(ЕО) т Юг ЦКК+2) ПО 3 4»1»)г ка к ю1»2 а ко ьи(ко).
3 Р(ко)»Р(ЕО)+Т 2 Н Е+9 ИЗТЦ)Е ИИП нцвиоцтгиз рипииу(ии,тт,ха,ихпхм) оо) ти(гоо),ио(гоо) о),532(!оо) зм,гси(1оо) вумиингон то(игвгм) соннаигнтсан4Увг(зо, 5 0 ) с(ммоиинннт/с(2 ос ),ги(2 °,м(2 ос),иК(ооо) д(го ,'зм, тгм,'онь, и, м,йьд ко 00121,М во!а=1 И 1 ви(г,а)=о.о во г 3-1,м 12 ьн(к+1) гг(п.им.1)осто з а ! ~ьн(к+3 ) вг(а! л) вг(а),а!)-С(х) аз ьи(к+2) во 4 а 1,аз ао*к+34а ао икао> 4 вг(ао а!)=вг(ао,а!)+0(х) осто 5 3 тг(гг ии 2)аота 6 а! ьн(к+и аг ьн(к+4) в01 м1 1,2 гг(И1 ио 1)0020 21 х1 аг т с(г)"то(а!) оото 22 21 21~а! т~с( 1)" го [аг ) 22 во о нг 1,2 гг а! хт(иг.ио.г)гг«аг о вг(гг г! )-вг(гг,г! )-т аоэЬН(К+2) во 9 У 1ва4 ,УО К+4+а ао ьн(уо) 5 вг(ао и )-вн(ао,г!)+т 1 сситху(аи ООХО™5 6 а)АЛЬМ(К+3) аг ьн(к!4) аз~ьн(к+5) ВОИН! !3 гт(и!.Мо. ! усато 23 тг(и! Ио 2)осто 24 11 аз т с(и то(а!)~хо(аг) Осто 25 23 21ЧТ! т с(г) го(аг) го(,уз) оста 25 24 11~ьУ2 т с(и то(л ) то(аз) 25 во 12 М2~1~3 22 Л гн(и2 но 2)зги п(нг.но.з)гг аз 12 вг(гг и ) вг(гг,г!)-т -14 ьн(к+2) во 13 а 1,а4 .уо к+з»а ао икао) 13 вг(ао 11)~вн(ао,з1)+т 11 с(н(тггмуи со(вне(/зтсаи/и, т,н,Михи,ншх,еРз,иг,нтьоа „тзтАНТ,МАЕОНВ сомма и/зтсом2/ХНАХ со)вкм/зтс(е(5/и™анан с(ез(см/зтс(е(6/Рм со)в(он/зтсом5/РНАте с(евин(/зтс(и(6/навез> иа(ВХВ Вп(тивхои х(50,6),хйох(50)>иаан(50) >Рм(2500),РЗАРХ(150),хо(50) т~тт н-нн ЕРЗ>«ЕРЗЗ МР=МРР ВО 1 Х~1>и х ~ 1)=хо(х) 1 СОИТХИСВ хнивв(нР).ьт.Авз(н))н нР ЕРЗСОМ®1,Е 12 ИХВПИ50 тзтАВТ» о НИАХ~НР НМ(М=1.
Е-10 МАХОНВ= 5 г соитппж хт(т>н.ие.т)аого 11 Рнхнт 7>т>н зтоР с(итхи((в хз о. 005 Х-(М Э Хв ХВ>АВВ(Х(Х~1)) ХЗ=ВРЗСОМ Х(~$РВ ВО 4 Х~1>н НИАХ(х) АМАХ1(Х(Х,1),ХЗ,(.Е-10) 4 СОИТХИСЕ САЬЬ ВТХРР(Х,ИППН) ХР(хтьоа.ие. 0)аото 5 хт(ХРМЗН(и т.н,х озтАНТ).еа.1)сото 5 ХР(т.ат.нР аалто сото 2 5 дЗТАЕТ -1 Н=Н/10.0 оста 2 5 тт=т ВОАХ1>И 6 ХО(х)=Х(х>1) 1 РОВНАТ(1Х, >ВАЛЬНЕЙШЕ МНТЕ'Р((РОВАНМЕ> »> 1Х 'НЕВОВИОННО - Т»Н Т'>2Х» 5-»Е10.0» Н '«Е18»0) Ммо занноатхне зтхРР(х,йтвхи) с(иион /зтсои1/и, т,'н.нмхм,ниох,егз, МР, кгьха,азтхнт,мххсен . со)вне/зтсою/ХИАХ( 1) СОИМ(И/Зтоа 5/Иаай(1) С(ИМОМ/ЗТСОМА/РМ(1) с(иион/втсом5/РВАте(1) с(имои/зтссмо/нсзеВ» масзе0 саима(/зтсхм)/иь,ма' Вп(еивхои х(ихвхм,1),еь(15),та(6) ВАТА Амохве/1.
ое-12/' ВАХА ВРьтен/1. ое-1а/ ИАххтг/5/,шхРАХ/5/ ВАТА Вмххмх/1. он+а~ Йиюа/1о./,ВИХРАХ/г./, хвеьАТ/10/ ВАТА п(сонг/. 25/,Вненю/.1/, па«знь/о. о/, йстезт/. 5/ ВАТА ЗХАЗ1/1.5/, Вь(зг/1. г/, вхАЗ)/1» А/ ьапсоь НТАь)А, с(итен хтьоа о тсо>В т 1ВОСЗ) ТЖВ ИЕ 0))ЕУА)»3А» ТВОЕ )»РАСТОН ХР(УЗТАВТ)100»100»200 100 НСЗЭВ=О, НСВБЕВ»0 Н1»ИХВХМ»1 И2 Н1+НХВЭМ НЗС. И»ИХВХН САЬЬ ВХРРСИ(И,Т,Х,РВАЧЕ) ВО !!О 1 1,И 110 ПХ,2)=РВАЧЕ( 1)"Н Н(~ 1 Ь 2 , ХВООВ»ь+1 НМАХ=Я(ХХНХ НРЗБ-"БСИТ( АИС ХНЕ) ТВВИВ»1, ОЬВЬО=!, Вс»0. НОЬВ»Н МЕХН=ЗИУ10 МХТЖН=МР-1 0 "МНЭН ЕУА3»ГА РАГ БЕ« хР(мхтей ие 0)ЕЧАХ»ГА тнсж СОНУЖХМ,РАХЗЕ сАГь сонет(нихн,нс,жь,тс,иАХВВВ) ЬИАХ»ИАХВЕВ» 1 ЖВ(»=(533(1)»ЖРЗ)»»2 Е» (ТО (2 )»ЖРБ)»»2 НОР»(ТО~Э ~»ИРЗ~»»2 СОТО 200 !Оо хи(н.ес нас)сото 200 Вв»Н/НОЬВ Н НОЬВ САХь ВВБСАь(Х.НХВхм, вн.
ВНАХ.ВО. ь, 200 хР((АВБ(нс-1.).От.встхвт).АЙВЗА Т»Т»Н САГ»ь РВЕВХС( НХВХМ» 1 >НО) 220 ВС 230 1=1,И 2ЭО ЕВВОВ(1) 0 0 1ТЖН 0 САЬЬ В1РРОИ(И,Т,Х,РВАЧЕ(И1 ) ) ХР(.НОТ ЕУАГ»ГА)ООТО 440 ОСТО(240 260 3 10),М1ТЖВ 240 САьь 232)йнГ(й,т,х,ихвхм) В=-Жь(1)»Н ВО 250 1»1»Н ВО 250 О=Х,НБО»НХВХМ 250 РН(0) РИ(0)»в 0070 300 260.БОМ*О 0 Во 270 1 1,И 270 ЗГИ ЗОН+РВАЧЕ( Х»НТВХМ)»»2 ВО=АНБ (Н)»ЗОНТ ( БОМ)» 1, Е+Э»АНОХЗЕ »11»-ИХВХМ ВО 290 3»1»И 01-"З1+КХВХМ ХО=ЮГ,1) В»ЖРЗГ»ХМ)(Х(»7) Н=АНАХ1 а, НО) Па !) 1(2 1)~н РАС1он -Вью!) нгв САЬЬ ВХРРЦН(И~Т~Х,РВАЧЕ) ВО 200 1=1,И 200 РИ(1421) (РВАЧЕ(Х)-РБАУЕ(1+ИХВХМ) 290 П 3»1) 12 300 ВО 305 1 1,ИЩИ! 247 ЗО5 РИ(1)«РИ(1)«!..О ЖЧАХ>ХА..
РАХаж. соичкн«РА)ак КС 1.О САХ>Х> МАХВЕС (И> РЧ>ИХВХМ> РЗА7В(И2) ) сото або З10 )~КР(1) О.1 соичкк=. итак. РО 520 1 1,И ЭгО УИ(1) Х(1, 1)«К (РВАЧИ(1«ИХРХИ) К-Х(1,2)) СА)ь вхууйи(Е,Т,РИ,УВАРЕ) ВО ЗЗО 1=1 И НО«РИ(1)-Х(1>1) РЧ(1)«1.О РАстон=ЙО-еь(1) «н (РВАЧЕ(1)-Рзачж(ихпхм+1) ) РНАРВ(1)«ко(к ХР(АЗЗ(ВО).30.0.0)СОТО 550 ХУ(НАСТОИ.ЕО.О 0)0020 збо Ри(1)-ногуостоЙ УВАРЕ(1)=РЗАЧВ( 1)«РИ(1) ЗЗО СОИТХИОН ЕЧАь А«.РАХЗЕ. кс !.
хко.о РО Зао 1«1 и ИНИОН(1) «КЙКОН( 1) +РВАЧЕ(1) В В+(РВАЧЕ(1)ГХМАХ(1) ) «г ЭАО УЗАЧК(1) Х(1,1)«ВП(1)«ВННОК(1) Сото 550 аао ху(мхтжк.нк.о)сото(або,або,аоо),мхтжн )~о. о ПО 450 1 1>М Н«К«РВАЧЕ(1«ИХРХМ)-Х(1 2) в-в ((н-кинок(1) ухмАХ(1) ) 2 УВАРЕ(1)«х(1,!)«Жп( 1) к азо жинок(1)-н Сото ззо або во а)о 1«1,И 4)0 РВАЧЕ(1«ИХРХМ)«УЗА73(Х+ИХРХМ)«Н-Х(1>2)-КВНОН(1) САПП МАЗБ(2 [И>РИ>ИХРХМ>УВАРЕ(И2)>РЗА73(И1),УЗОРЕ) в=о.о ВО 400 1«1,И жинок( 1) «ЕЙВОИ(1) +РВА73(1) )Фв«(ХНАМ(1)/ХМАХ(1) )««2 аою РВАЧЕ(1)«Х(1, ! )«КП( ! )«ИНИОН(1) Оото 550 490 ВО 500 1«1 И 500 РВАЧЕ(1) Рй(1)«(РВАЧЕ(1«ИХРХМ)%-Х(1>2)-ЕННОН(1) ) в о.о ВО 5!О Х«1>И кинок(1) =жйк(н((1)+узок( х ) Ы«(РВАЧЕ(1)УХМАХ(1) )««2 510 УВАРЕ(1)=Х(1 ° 1)«ЕВ(1)«ХЗНОК(1) ООТО 550 540 УЗА73(1)=Х(1 .1 )«3)>(1)«КННОИ(1) 550 1Р(1ТЕН иж,о)тккив«ОМАХ!(0 9«ТВХМР Р/Р1) ХУ((Р«АИХИ1(1 0>2.0«ХВЕНР) ) Ьж ВИР)СОТО 590 01 Р ХТЕК«ХТХМ«1 хР(1тен ВО ИАххтк)ООХО 560 САХ>Х> ВХУУПИ(И>т>РЗАЧЖ>РВАЧЕ(И1) ) Оотс 440 560 ху(>Ест с(мчкн)оотс 570 ЕЧАХ«)А«тквж ОСТО 220 5)0 Ы 42 «)ВС(УАХ 0,4 5>7 3333333> 670 НРЗИР Н Иа(]ЗЕВ«ИЯ КРЬАО«О РО 680 !«1>Ь РО 680 1 1,И 6ВО Х(1>д)«Х(Х,О)«ж(д)«ЕВВОВ(1) ХР( ХР(Х!В Е0.1 )СОТО 700 ХРО(И]«1РООВ-1 ХР((хвоов.ат.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
