Диссертация (1172877), страница 14
Текст из файла (страница 14)
Комаров; правообладатели М.А. Грохотов, И.Р.112Бегишев, А.К. Беликов, А.А. Комаров; заявл. 04.10.2017; зарегистр. 22.11.2017;опубл. 22.11.2017, Бюл. № 1. – 1 с.61. ГОСТхарактеристики8.585–2001.Термопары.преобразованияНоминальные[Электронныйресурс]:статическиегосударственныйстандарт // Гарант: инф.-прав. об-ние. – Эл. дан. – М., 2017. – Доступ из лок-нойсети б-ки Академии ГПС МЧС России (дата обращения 22.09.2017)62.
Термогазодинамика пожаров в помещениях [Текст] / М.В. Астапенко,Ю.А. Кошмаров, И.С. Молчадский, А.Н. Шевляков / Под ред. Ю.А. Кошмарова. –М.: Стройиздат, 1988. – 448 с.63. МанометрысусловнымиишкаламивакуумметрытиповМОдеформационныеиВОобразцовые[Электронныйресурс]:паспорт 3.9060.303 ПС // ООО «Приборика»: сайт. – Режим доступа:https://priborika.ru/omishop/01/0101/010108/mo015.pdf (дата обращения 22.09.2017)64.
Грохотов, М.А. Параметры пожаровзрывоопасности хлорпроизводныхметана в смеси с хлором в модели фототеплового воспламенения [Текст] /М.А. Грохотов, И.Р. Бегишев, А.К. Беликов, А.А. Комаров // Пожарнаябезопасность. – 2017. – № 3. – С.
43 – 49.65. Беликов, А.К. Воспламенение смесей хлорметана с хлором поддействием ультрафиолетового излучения [Текст] / А.К. Беликов, И.С. Никитин,И.Р. Бегишев // Пожаровзрывобезопасность. – 2010. – № 10.– С. 9–12.66. Бегишев, И.Р. Формирование очага воспламенения в смеси CH3Cl + Cl2при действии УФ-света [Текст] / И.Р. Бегишев, А.К. Беликов, П.В. Комраков,И.С.
Никитин // Журнал физической химии. – 2016. – № 7. – С. 993–999.67. Комраков,П.В.Взрывобезопасныеусловияфотохлорированиядихлорметана в газовой фазе / П.В. Комраков, А.К. Беликов, И.Р. Бегишев //симпозиум: XII Симпозиум по горению и взрыву. Химическая физика процессовгорения и взрыва. – Черноголовка, 2000. – Т. 3. – С. 48–49.68. Комраков, П.В. Особенности воспламенения и распространения пламенив газовых смесях дихлорметана с хлором при зажигании ультрафиолетовымсветом [Электронный ресурс] / И.Р. Бегишев, А.К.
Беликов, П.В. Комраков //113Пожарноедело:Электронныйжурнал.––2007.Режимдоступа:http://pojar01.ru/11/PROCESS_GOR/ST/ST_KOMRAKOV_/text.html69. Беликов,А.К.Ингибирующийэффектприфлегматизациитрихлорметаном смесей хлористого метила и хлора при их фотовоспламенении[Текст] / А.К. Беликов, И.С. Никитин, И.Р. Бегишев // Химическая физика. –2011. – № 10.
– С. 1–4.70. Никитин, И.С. Флегматизация четырѐххлористым углеродом смесейхлорметана и хлора при их фотовоспламенении [Текст] / И.С. Никитин,А.К. Беликов, И.Р. Бегишев // Пожаровзрывобезопасность. – 2015. – № 2. – С. 32–35.71. Тимофеев,А.Ф.Техникабезопасностиприхранении,транспортировании и применении хлора [Текст] / А.Ф. Тимофеев, А.Ф. Мазанко,Б.Ю.
Ягуд и др. – М.: Химия, 1990. – 336 с.72. Зельдович,Я.Б.Теориядетонации[Текст]/Я.Б.Зельдович,А.С. Компанеец. – М.: Гос. изд-во технико-теоретической литературы, 1955. – 268 с.73. Лойцянский, Л.Г. Механика жидкости и газа [Текст] / Л.Г. Лойцянский.–М.: Наука, 1970. – 904 с.74. Варнатц,Ю.Горение.Физическиеихимическиеаспекты,моделирование, эксперименты, образование загрязняющих веществ [Текст] /Ю. Варнатц, У. Маас, Р.
Диббл / Пер. с англ. Г.Л. Агафонова; под ред.П.А. Власова. – М.: Физматлит, 2003. – 352 с.75. Грохотов, М.А. Расчѐт СРФП пламени при дефлаграционном взрыве[Текст] / М.А. Грохотов [и др.] // Материалы 6-й Международной науч.-практ.конф.: Пожаротушение: проблемы, технологии, инновации. – М.: Академия ГПСМЧС России, 2018. – С. 301–303.76. Зельдович, Я.Б. Теория горения идетонации газов [Текст] /Я.Б. Зельдович.
‒ М.: Изд. АН СССР, 1944. – 71 с.77. Щелкин, К.И. Неустойчивость горения и детонации газов [Текст] /К.И. Щелкин // Успехи физических наук. – 1967. – Т. 87 (вып. 2). – С. 273–302.78. Зельдович, Я.Б. О распределении давления и скорости в продуктахдетонационного взрыва, в частности, при сферическом распространении114детонационной волны [Текст] / Я.Б. Зельдович // Журнал экспериментальной итеоретической физики. – 1942. – Т. 12 (вып. 9).
– С. 389–406.79. Teodorczyk, A. The structure of fast turbulent flames in very rough, obstaclefilled channels / A. Teodorczyk, J.H.S. Lee, R. Knystautas // International Symposium onCombustion, 1991, vol. 23, No. 1, pp. 735–741.80. Wolanski P. Detonative propulsion / P. Wolanski // Proceedings of theCombustion Institute, 2013, vol. 34, No. 1, pp. 125–158.81.
Азатян, В.В. Проблемы горения, взрыва и детонации газов в теориинеизотермических цепных процессов [Текст] / В.В. Азатян // Журнал физическойхимии. – 2014. – Т. 88. – № 5. – С. 759–771.82. Азатян, В.В. Цепные реакции в процессах горения, взрыва и детонациигазов [Текст] / В.В. Азатян. – Черноголовка: Изд-во ИПХФ РАН, 2017. – 448 с.83. Семѐнов, В.Н. Физика быстропротекающих процессов. Горение идетонация газовых смесей [Текст] / В.Н. Семѐнов.
– М.: Институт проблембезопасного развития атомной энергии, 2006. – 66 с.84. Войцеховский, Б.В. Структура фронта детонации в газах [Текст] /Б.В. Войцеховский, В.В. Митрафанов, М.Е. Топчиян. – Новосибирск: Изд-во СОАН СССР, 1963.85. Гвоздева, Л.Г. Преломление детонационных волн при падении награницу раздела газовых смесей [Текст] / Л.Г. Гвоздева // ЖТФ, 1961.
– Т. 31. –№ 6. – С. 731–739.86. Гвоздева,Л.Г.Экспериментальноеисследованиедифракциидетонационных волн в стехиометрической смеси метана с кислородом [Текст] /Л.Г. Гвоздева // ЖПМТФ. – 1961. – № 5. – С. 53–56.87. Bredley D. Mathematical solutions for explosions in spherical vessels /D. Bredley, A. Mitcheson // Combustion and Flame, 1976, vol.
26, No. 2, pp. 201–207.88. Berman, M. A critical review of recent large scale experiments on hydrogenair detonations / M. Berman // Nuclear Science and Engineering, 1986, vol. 93, pp.321–347.11589. Макеев, В.И. Горение и переход к детонации газовых смесей впространстве с частичным загромождением / В.И. Макеев, А.А.
Пономарев,В.В. Строганов, В.Л. Карпов // Пожаровзрывобезопасность. – 1993. - № 2. – С. 3–6.90. Макеев, В.И. Экспериментальное исследование ускорения пламени всвободных сферических объѐмах [Текст] / В.И. Макеев, В.В. Строганов,А.П. Чугуев, Ю.Н.
Чернушкин // В кн. Горение и проблемы тушения пожаров. –М.: ВНИИПО, 1981. – С. 17–20.91.Cook G. Flixborough petrochemical plant [Электронный ресурс] /G. Cook. – Режим доступа: http://www.cookeonfire.com/pdfs/Flixborough.pdf (датаобращения 15.03.2018)92. ПоследнийденьзаводавФликсборо[Электронныйресурс]/И. Сорокина // Новые химические технологии: аналитический портал.
– Режимдоступа:http://www.newchemistry.ru/printletter.php?n_id=302(датаобращения15.03.2018)93.Маршалл, В. Основные опасности химических производств [Текст] /В. Маршалл. – М.: Мир, 1989. – 672 с.94.Пособие по обследованию и проектированию зданий и сооружений,подверженных воздействию взрывных нагрузок [Текст]. – М.: ЦНИИПромзданий,2000.
– 122 с.95.Бузаев,Е.В.Разработкаметодовпрогнозированияпараметроввзрывоопасных зон при аварийных выбросах горючих веществ [Текст]: дис. …канд. техн. наук: 05.26.03 / Бузаев Евгений Владимирович. – М., 2015. – 124 с.96. Комаров,А.А.Особенностигидродинамическихрасчѐтовприпроведении технической экспертизы аварийного взрыва [Текст] / А.А. Комаров,Д.А. Корольченко, Н.В. Громов // Пожаровзрывобезопасность. – 2017. – Т.26. –№ 10.
– С. 36–42.116Приложение А(обязательное)АКТЫ ВНЕДРЕНИЯ117118119120121Приложение Б(обязательное)КОД ПРОГРАММЫ«РАСЧЁТНЫЙ МЕТОДОПРЕДЕЛЕНИЯ ПОКАЗАТЕЛЕЙ ПОЖАРОВЗРЫВООПАСНОСТИ ГАЗОВЫХСМЕСЕЙ.122function varargout = menu(varargin)% MENU MATLAB code for menu.fig%MENU, by itself, creates a new MENU or raises the existing%singleton*.%%H = MENU returns the handle to a new MENU or the handle to%the existing singleton*.%%MENU('CALLBACK',hObject,eventData,handles,...) calls the local%function named CALLBACK in MENU.M with the given input arguments.%%MENU('Property','Value',...) creates a new MENU or raises the%existing singleton*.
Starting from the left, property value pairsare%applied to the GUI before menu_OpeningFcn gets called. An%unrecognized property name or invalid value makes propertyapplication%stop. All inputs are passed to menu_OpeningFcn via varargin.%%*See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one%instance to run (singleton)".%% See also: GUIDE, GUIDATA, GUIHANDLES% Edit the above text to modify the response to help menu% Last Modified by GUIDE v2.5 09-Jul-2017 22:47:37% Begin initialization code - DO NOT EDITgui_Singleton = 1;gui_State = struct('gui_Name',mfilename, ...'gui_Singleton', gui_Singleton, ...'gui_OpeningFcn', @menu_OpeningFcn, ...'gui_OutputFcn', @menu_OutputFcn, ...'gui_LayoutFcn', [] , ...'gui_Callback',[]);if nargin && ischar(varargin{1})gui_State.gui_Callback = str2func(varargin{1});endif nargout[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});elsegui_mainfcn(gui_State, varargin{:});end% End initialization code - DO NOT EDIT% --- Executes just before menu is made visible.function menu_OpeningFcn(hObject, eventdata, handles, varargin)set(handles.pushbutton3, 'Enable', 'off') %делает кнопку не доступнойset(handles.pushbutton4, 'Enable', 'off') %делает кнопку не доступной% This function has no output args, see OutputFcn.% hObjecthandle to figure% eventdata reserved - to be defined in a future version of MATLAB% handlesstructure with handles and user data (see GUIDATA)% varargincommand line arguments to menu (see VARARGIN)123% Choose default command line output for menuhandles.output = hObject;% Update handles structureguidata(hObject, handles);% UIWAIT makes menu wait for user response (see UIRESUME)% uiwait(handles.figure1);% --- Outputs from this function are returned to the command line.function varargout = menu_OutputFcn(hObject, eventdata, handles)% varargout cell array for returning output args (see VARARGOUT);% hObjecthandle to figure% eventdata reserved - to be defined in a future version of MATLAB% handlesstructure with handles and user data (see GUIDATA)% Get default command line output from handles structurevarargout{1} = handles.output;% --- Executes on button press in pushbutton1.function pushbutton1_Callback(hObject, eventdata, handles)set(handles.pushbutton1, 'Enable', 'off') %делает кнопку не доступнойset(handles.pushbutton2, 'Enable', 'on') %делает кнопку доступнойset(handles.pushbutton3, 'Enable', 'off') %делает кнопку не доступнойset(handles.pushbutton4, 'Enable', 'off') %делает кнопку не доступной% hObjecthandle to pushbutton1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handlesstructure with handles and user data (see GUIDATA)clear allcla% --- Executes on button press in pushbutton2.function pushbutton2_Callback(hObject, eventdata, handles)set(handles.pushbutton2, 'Enable', 'off') %делает кнопку не доступной Стартset(handles.pushbutton1, 'Enable', 'on') %делает кнопку доступной Стоп% hObjecthandle to pushbutton2 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handlesstructure with handles and user data (see GUIDATA)global e1 p kt Et R k5global k11 k12 k13 k14EE11 EE12global k21 k22 k23 k24EE21 EE22global k31 k32 k33 k34EE31 EE32global k41 k42 k43 k44EE41 EE42global Q1 Q2 Q3 Q4k11=0.2*25118864.32; k12=0.2*7943282.347; k13=36307805.48; k14=398107170.6;%%%m^3/(mol*s) метанk21=33113112.15; k22=3981071.706; k23=3981071.706; k24=251188643.2;%%%m^3/(mol*s) хлорметанk31=26915348.04; k32=1000000; k33=2511886.432; k34=223872113.9;%%%m^3/(mol*s) дихлорметанk41=6918309.709; k42=549540.8739; k43=50118.72336; k44=63095734.45;%%%m^3/(mol*s) трихлорметанk5=19952.62315; %%%m^6/(mol^2*s)EE11=16035.44; EE12=9629.64; EE21=13816.44; EE22=12560.4; %%% Дж/мольEE31=12560.4; EE32=16747.2; EE41=13983.912; EE42=20934; %%%% Дж/мольe1=2; p=41.6; kt=10^8; Et=196648; R=8.314;Q1=98920; Q2=105060; Q3=99270; Q4=92000;124%%%%%%%%%%%% Построение сетки по X %%%%%%%%%%%%npx=60;% число ячеек по хx1=0;x2=0.04;%in mdx=(x2-x1)/npx; %in m шаг по хx=linspace(x1+dx/2,x2-dx/2,npx);% массив иксов, linspace разбиваеит наравные части%%%%%%%%%%%% ГРАНИЧНЫЕ УСЛОВИЯ X %%%%%%%%%%%%% 0 - свободная, 1 - жесткая, 2 - задается%леваяjxgrl=ones(size(x));% jxgrl=zeros(size(x));% jxgrl(10:30)=0;%jxgrl=293*ones(size(x));% 1 - жесткая,% 0 - свободная,%праваяjxgrp=ones(size(x));% jxgrp=zeros(size(x));%jxgrp=293*ones(size(x));%%%%%%%%%%%% Построение сетки по X %%%%%%%%%%%%npy=27;y1=0;y2=0.018;dy=(y2-y1)/npy;y=linspace(y1+dy/2,y2-dy/2,npy);IIIII0=zeros(size(y));%%%%%%%%%%%%%%if (dx~=dy)msgbox('НЕ РАВНОМЕРНАЯ СЕТКА','Предупреждение', 'warn')end%%%%%%%%%%%%%%%%%%%%%%%%%% ГРАНИЧНЫЕ УСЛОВИЯ Y %%%%%%%%%%%%%левая%iygrl=ones(size(y));iygrl=zeros(size(y));% iygrl=293*ones(size(y));% 0 - свободная,%праваяiygrp=ones(size(y));% iygrp=zeros(size(y));% iygrp=293*ones(size(y));%%%%%%%%%%%%% ИСХОДНЫЙ ПРОФИЛЬ %%%%%%%%%%%%%%%%[X,Y]=meshgrid(x,y);T=ones(size(X));Tpr=zeros(size(X));%температура в КTprPr=zeros(size(X));%температура в КCLOR=zeros(size(X));%концентрация хлора в моль/м3125CH4=zeros(size(X));%концентрация метанаWW1=zeros(size(X));WW11=zeros(size(X));CH3Cl=zeros(size(X));%концентрация хлорметанаWW2=zeros(size(X));WW22=zeros(size(X));CH2Cl2=zeros(size(X));%концентрация дихлорметанаWW3=zeros(size(X));WW33=zeros(size(X));CHCl3=zeros(size(X));%концентрация трихлорметанаWW4=zeros(size(X));WW44=zeros(size(X));CCl4=zeros(size(X));%концентрация тетрахлорметанаWW5=zeros(size(X));Cl=zeros(size(X));%концентрация радикалов хлораDx=zeros(size(X));nn=zeros(size(X));Cv1=zeros(size(X));Cv2=zeros(size(X));Cv3=zeros(size(X));Cv4=zeros(size(X));Cv5=zeros(size(X));Cv6=zeros(size(X));%%%%%%%%%%%%%%%%%%%%%t=str2num (get (handles.edit8, 'String')); %%время расчѐта, секdt=str2num (get (handles.edit10, 'String')); %%шаг по времени, секI0=str2num (get (handles.edit6, 'String')); %%%% интенсивность светаM=41.6; %количество молей смесиSS=t/dt;nstep=fix(SS); %%%количество шагов%%%%%%%%%%%%%%%%%%%%%%%%Концентрация вещесвт в смеси (горючего) задаѐтся в процентах%%%ConCH4=str2num (get (handles.edit1, 'String'));%метанConCH3Cl=str2num (get (handles.edit2, 'String')); %хлорметанConCH2Cl2=str2num (get (handles.edit3, 'String')); %дихлорметанConCHCl3=str2num (get (handles.edit4, 'String')); %трихлорметанConCCl4=str2num (get (handles.edit5, 'String'));%четырехлористый угT0=str2num (get (handles.edit7, 'String')); %начальная температура смеси, КRashGaz=1; %расширение газа 1 - ДА, 0 - НЕТblah=get(handles.uipanel3, 'SelectedObject');whatever=get(blah, 'Tag');if(whatever == 'radiobutton1')IIIII0(1:npy)=1; %включение УФ-света 1 - ДА, 0 - НЕТelseIIIII0(1:npy)=0; %включение УФ-света 1 - ДА, 0 - НЕТend%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%PPPP=(100-ConCH4-ConCH3Cl-ConCH2Cl2-ConCHCl3-ConCCl4); %хлор в процентах%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%CH40=M*ConCH4/100;CH3Cl0=M*ConCH3Cl/100;CH2Cl20=M*ConCH2Cl2/100;CHCl30=M*ConCHCl3/100;CCl40=M*ConCCl4/100;PromCon=(CH40+CH3Cl0+CH2Cl20+CHCl30+CCl40);CLOR0=M-PromCon;T=T0*T;126for j=1:npx%%% Xfor i=1:npy%%%YCLOR(i,j)=CLOR0;% начальные условияCH4(i,j)=CH40;CH3Cl(i,j)=CH3Cl0;CH2Cl2(i,j)=CH2Cl20;CHCl3(i,j)=CHCl30;CCl4(i,j)=CCl40;endendddtend=dt*nstep;nit=0;ttt=zeros(1,nstep);%TTpoint1=ttt;TTpoint2=ttt;%TTpoint3=ttt;TTT=ttt;CLORR=ttt; % хлор средеобъемныйWWW1=ttt; %метанWWW2=ttt; %хлорметанWWW3=ttt; %дихлорметанWWW4=ttt; %трихлорметанWWW5=ttt; %тетрахлорметанCLORRR=ttt; %хлор в точкеDxx=ttt;nnn=ttt;%%первая точка%%%jxpoint1=25;%iypoint1=20;%%вторая точка%%jxpoint2=30;iypoint2=1;%%третья точка%%%jxpoint3=76;%iypoint3=20;nprint=120; %% шаг вывода графикаddtnprint=dt*nprint;for it=1:nstep%%%%%%%nit=nit+1;nnout=nit-nprint*fix(nit/nprint);ttt(it)=it*dt;%TTpoint1(it)= T(iypoint1,jxpoint1);TTpoint2(it)= T(iypoint2,jxpoint2);%TTpoint3(it)= T(iypoint3,jxpoint3);WWW1(it)=CH4(iypoint2,jxpoint2); %концентрация метана в точкеWWW2(it)=CH3Cl(iypoint2,jxpoint2); %концентрация хлорметана в точкеWWW3(it)=CH2Cl2(iypoint2,jxpoint2); %концентрация дихлорметана в точкеWWW4(it)=CHCl3(iypoint2,jxpoint2); %концентрация трихлорметана в точкеWWW5(it)=CCl4(iypoint2,jxpoint2); %концентрация тетрахлорметана в точкеCLORRR(it)=CLOR(iypoint2,jxpoint2); %концентрация хлора в точкеDxx(it)=Dx(iypoint2,jxpoint2);nnn(it)=100*nn(iypoint2,jxpoint2);%%%ТЕМПЕРАТУРА НА ГРАНИЦЕ%%%TTT0xpr=str2num (get (handles.edit9, 'String'));TTT0xlev=str2num (get (handles.edit9, 'String'));127TTT0ypr=str2num (get (handles.edit9, 'String'));TTT0ylev=str2num (get (handles.edit9, 'String'));bbb=sum(sum(T))/npx/npy; %%%% расчѐт среднеобъемной температурыCCv=((35.6-R)*ConCH4+(40.7-R)*ConCH3Cl+(50.8-R)*ConCH2Cl2+(65.3R)*ConCHCl3+(82.6-R)*ConCCl4+(33.9-R)*PPPP)/100;%%теплоемкость смеси вначальный моментfor j=1:npx%%% Xfor i=1:npy%%%Y%%%%%%%%%%%%%%%%%%%%%%%%%% BS 11111j2=j+1;if(j==npx)if (iygrp(i)==0)POTX2=0;elsePOTX2=TTT0xpr-T(i,npx);endelsePOTX2=T(i,j+1)-T(i,j);endN%%правая граница по Х %%j1=j-1;if (j==1)if (iygrl(i)==0)POTX1=0;elsePOTX1=T(i,1)-TTT0xlev;end%%левая граница по Х%%elsePOTX1=T(i,j)-T(i,j-1);end%%%%%%%%%%%%%%%%%%%%i2=i+1;if (i==npy)if (jxgrl(j)==0)POTY2=0;%%левая граница по Y%%elsePOTY2=(TTT0ylev-T(npy,j));endelsePOTY2=(T(i+1,j)-T(i,j));endi1=i-1;if (i==1)if (jxgrp(j)==0)POTY1=0;%%правая граница по Y%%elsePOTY1=(T(1,j)-TTT0ypr);endelsePOTY1=(T(i,j)-T(i-1,j));end%%%%%%%МЕТАН%%%%%%%if(CH40>0)128if (T(i,j)<=1300)AA1=-0.703029; BB1=108.4773; CC1=-42.52157; DD1=5.862788; EE1=0.678565;elseAA1=85.81217; BB1=11.26467; CC1=-2.114146; DD1=0.138190; EE1=-26.42221;endCv1(i,j)=AA1+BB1*(T(i,j)/1000)+CC1*(T(i,j)/1000)^2+DD1*(T(i,j)/1000)^3+EE1/(T(i,j)/1000)^2-R; %теплоѐмкость метана Дж/(моль*К)elseCv1(i,j)=0;end%%%%%%%ХЛОРМЕТАН%%%%%%%if (CH3Cl0>0)if (T(i,j)<=1200)AA2=3.524690; BB2=136.9277; CC2=-82.14196; DD2=20.22797; EE2=0.278032;elseAA2=88.66691; BB2=10.06834; CC2=-1.928737; DD2=0.128191; EE2=-19.76840;endCv2(i,j)=AA2+BB2*(T(i,j)/1000)+CC2*(T(i,j)/1000)^2+DD2*(T(i,j)/1000)^3+EE2/(T(i,j)/1000)^2-R; %теплоѐмкость хлорметана Дж/(моль*К)elseCv2(i,j)=0;end%%%%%%%ДИХЛОРМЕТАН%%%%%%%if(CH2Cl20>0)if (T(i,j)<=1200)AA3=19.17351; BB3=136.8444; CC3=-95.12993; DD3=26.03105; EE3=-0.119405;elseAA3=95.12993; BB3=6.721722; CC3=-1.288196; DD3=0.085646; EE3=-14.93885;endCv3(i,j)=AA3+BB3*(T(i,j)/1000)+CC3*(T(i,j)/1000)^2+DD3*(T(i,j)/1000)^3+EE3/(T(i,j)/1000)^2-R; %теплоѐмкость дихлорметана Дж/(моль*К)elseCv3(i,j)=0;end%%%%%%%ТРИХЛОРМЕТАН%%%%%%%if(CHCl30>0)if (T(i,j)<=1200)AA4=44.24706; BB4=114.6734; CC4=-88.81837; DD4=25.89992; EE4=-0.522808;elseAA4=101.4762; BB4=3.429756; CC4=-0.657348; DD4=0.043707; EE4=-9.348019;endCv4(i,j)=AA4+BB4*(T(i,j)/1000)+CC4*(T(i,j)/1000)^2+DD4*(T(i,j)/1000)^3+EE4/(T(i,j)/1000)^2-R; %теплоѐмкость трихлорметана Дж/(моль*К)elseCv4(i,j)=0;end%%%%%%%ТЕТРАХЛОРМЕТАН%%%%%%%if(CCl40>0)AA5=103.1134; BB5=4.188644; CC5=-1.126475; DD5=0.095677; EE5=-1.919624;Cv5(i,j)=AA5+BB5*(T(i,j)/1000)+CC5*(T(i,j)/1000)^2+DD5*(T(i,j)/1000)^3+EE5/(T(i,j)/1000)^2-R; %теплоѐмкость трихлорметана Дж/(моль*К)elseCv5(i,j)=0;129end%%%%%%%ХЛОР%%%%%%%if(CLOR0>0)if (T(i,j)<=1000)AA6=33.0506; BB6=12.2294; CC6=-12.0651; DD6=4.38533; EE6=-0.159494;elseAA6=42.6773; BB6=-5.00957; CC6=1.904621; DD6=-0.165641; EE6=-2.09848;endCv6(i,j)=AA6+BB6*(T(i,j)/1000)+CC6*(T(i,j)/1000)^2+DD6*(T(i,j)/1000)^3+EE6/(T(i,j)/1000)^2-R; %теплоѐмкость метана Дж/(моль*К)elseCv6(i,j)=0;endDx(i,j)=(3.22*10^-2*(ConCH4+ConCH3Cl+ConCH2Cl2+ConCHCl3+ConCCl4)+0.88*10^-2*PPPP)/100/CCv/p*(T(i,j)/293); %%%коэффициент температуропроводностиTprPr(i,j)=T(i,j)+Dx(i,j)*(POTX2-POTX1)/4/dx^2*dt+Dx(i,j)*(POTY2*(i+1)POTY1*(i-1))/4/dy^2*dt;if(RashGaz==1)aaa=bbb/T(i,j); %%% расчет отношения сред.объемной темпер-ры к темп-ре вточкеelseaaa=1;endif (IIIII0(i)>0)WF=2*e1*CLOR(i,j)*aaa*I0*IIIII0(i)*exp(-e1*CLOR(i,j)*aaa*dx*j)*exp(10.59*dx*j);elseWF=0;endWT=2*kt*M*aaa*CLOR(i,j)*aaa*exp(-Et/R/TprPr(i,j));WW=WF+WT;if (CLOR(i,j)>=0)A=WW*CLOR(i,j)*aaa;if (CH4(i,j)>=0)A1=k11*k14/k12*CH4(i,j)*aaa*exp(-(EE11-EE12)/R/TprPr(i,j));elseA1=0;endif (CH3Cl(i,j)>=0)A2=k21*k24/k22*CH3Cl(i,j)*aaa*exp(-(EE21-EE22)/R/TprPr(i,j));elseA2=0;endif (CH2Cl2(i,j)>=0)A3=k31*k34/k32*CH2Cl2(i,j)*aaa*exp(-(EE31-EE32)/R/TprPr(i,j));elseA3=0;130endif (CHCl3(i,j)>=0)A4=k41*k44/k42*CHCl3(i,j)*aaa*exp(-(EE41-EE42)/R/TprPr(i,j));elseA4=0;endA5=k5*M*aaa;Cl(i,j)=(A/(A1+A2+A3+A4+A5))^0.5;elseCl(i,j)=0;end%%%%%%%%%%%%%%% 1-я стадия %%%%%%%%%%%%%%if (CH4(i,j)>=0)WW1(i,j)=k11*CH4(i,j)*aaa*Cl(i,j)*exp(-EE11/R/TprPr(i,j))*dt;%%%образование хлорметанаif((CH4(i,j)-WW1(i,j))>=0)CH4(i,j)=CH4(i,j)-WW1(i,j); %концентрация метана на следующий шагCH3Cl(i,j)=CH3Cl(i,j)+WW1(i,j); %концентрация хлорметанаelseCH4(i,j)=CH4(i,j)-CH4(i,j); %концентрация метана на следующий шагCH3Cl(i,j)=CH3Cl(i,j)+CH4(i,j); %концентрация хлорметанаWW1(i,j)=CH4(i,j);endelseCH3Cl(i,j)=0;WW1(i,j)=0;end%%%%%%%%%%%%%%% 2-я стадия %%%%%%%%%%%%%%if (CH3Cl(i,j)>=0)WW22(i,j)=CH3Cl(i,j)-WW1(i,j);WW2(i,j)=k21*WW22(i,j)*aaa*Cl(i,j)*exp(-EE21/R/TprPr(i,j))*dt;%%%образование дихлорметанаif((CH3Cl(i,j)-WW2(i,j))>=0)CH3Cl(i,j)=CH3Cl(i,j)-WW2(i,j); %концентрация хлорметана на следующий шагCH2Cl2(i,j)=CH2Cl2(i,j)+WW2(i,j); %концентрация дихлорметанаelseCH3Cl(i,j)=CH3Cl(i,j)-WW22(i,j); %концентрация хлорметана на следующий шагCH2Cl2(i,j)=CH2Cl2(i,j)+WW22(i,j); %концентрация дихлорметанаWW2(i,j)=WW22(i,j);endelseCH3Cl(i,j)=0;WW2(i,j)=0;end%%%%%%%%%%%%%%% 3-я стадия %%%%%%%%%%%%%%if (CH2Cl2(i,j)>=0)WW33(i,j)=CH2Cl2(i,j)-WW2(i,j);WW3(i,j)=k31*WW33(i,j)*aaa*Cl(i,j)*exp(-EE31/R/TprPr(i,j))*dt;%%%образование трихлорметанаif((CH2Cl2(i,j)-WW3(i,j))>=0)CH2Cl2(i,j)=CH2Cl2(i,j)-WW3(i,j); %концентрация дихлорметана на следующийшаг131CHCl3(i,j)=CHCl3(i,j)+WW3(i,j); %концентрация трихлорметанаelseCH2Cl2(i,j)=CH2Cl2(i,j)-WW33(i,j); %концентрация дихлорметана на следующийшагCHCl3(i,j)=CHCl3(i,j)+WW33(i,j); %концентрация трихлорметанаWW3(i,j)=WW33(i,j);endelseCH2Cl2(i,j)=0;WW3(i,j)=0;end%%%%%%%%%%%%%%% 4-я стадия %%%%%%%%%%%%%%if (CHCl3(i,j)>=0)WW44(i,j)=CHCl3(i,j)-WW3(i,j);WW4(i,j)=k41*WW44(i,j)*aaa*Cl(i,j)*exp(-EE41/R/TprPr(i,j))*dt;%%%образование тетрахлорметанаif((CHCl3(i,j)-WW4(i,j))>=0)CHCl3(i,j)=CHCl3(i,j)-WW4(i,j); %концентрация трихлорметана на следующийшагCCl4(i,j)=CCl4(i,j)+WW4(i,j); %концентрация тетрахлорметанаelseCHCl3(i,j)=CHCl3(i,j)-WW44(i,j); %концентрация трихлорметана на следующийшагCCl4(i,j)=CCl4(i,j)+WW44(i,j); %концентрация тетрахлорметанаWW4(i,j)=WW44(i,j);endelseCHCl3(i,j)=0;WW4(i,j)=0;endCLOR(i,j)=CLOR(i,j)-WW1(i,j)-WW2(i,j)-WW3(i,j)-WW4(i,j);хлора%% концентрацияif(CLOR(i,j)>0)QQQQQ=(WW1(i,j)*Q1+WW2(i,j)*Q2+WW3(i,j)*Q3+WW4(i,j)*Q4);QQQQQ=QQQQQ/(Cv1(i,j)+Cv2(i,j)+Cv3(i,j)+Cv4(i,j)+Cv5(i,j))/p;тепла (Kельвин) прибавилосьelseCLOR(i,j)=0;QQQQQ=0;endnn(i,j)=1.90*10^-14*(T(i,j))^4;Tpr(i,j)= TprPr(i,j) + (1-nn(i,j))*QQQQQ;end%%%Yend%%% XT=Tpr;TTT(it)=sum(sum(T))/npx/npy;CLORR(it)=sum(sum(CLOR))/npx/npy;if (nnout==0)surf(X,Y,T)%%сколько132%colormap autumnвероятности уровни)colormap jetвероятности)(для давления)xlabel('x,m')zlabel('T,K')ylabel('y,m')%axis([x1 x2 y1 x2 0 50])axis([x1 x2 y1 y2 0 2000])%'pause Step'pause (0.01)% от красного к желтому(для% от синего до красного (дляendend%%%%%for it=1:nstepsave W1.mat WWW1 -asciisave W2.mat WWW2 -asciisave W3.mat WWW3 -asciisave W4.mat WWW4 -asciisave W5.mat WWW5 -asciisave TT.mat TTpoint2 -asciisave t1.mat ttt -asciiset(handles.pushbutton3, 'Enable', 'on') %делает кнопку доступной График1set(handles.pushbutton4, 'Enable', 'on') %делает кнопку доступной График2function edit6_Callback(hObject, eventdata, handles)% hObjecthandle to edit6 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handlesstructure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit6 as text%str2double(get(hObject,'String')) returns contents of edit6 as adouble% --- Executes during object creation, after setting all properties.function edit6_CreateFcn(hObject, eventdata, handles)% hObjecthandle to edit6 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handlesempty - handles not created until after all CreateFcns called% Hint: edit controls usually have a white background on Windows.%See ISPC and COMPUTER.if ispc && isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');endfunction edit7_Callback(hObject, eventdata, handles)% hObjecthandle to edit7 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handlesstructure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit7 as text%str2double(get(hObject,'String')) returns contents of edit7 as adouble133% --- Executes during object creation, after setting all properties.function edit7_CreateFcn(hObject, eventdata, handles)% hObjecthandle to edit7 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handlesempty - handles not created until after all CreateFcns called% Hint: edit controls usually have a white background on Windows.%See ISPC and COMPUTER.if ispc && isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');endfunction edit8_Callback(hObject, eventdata, handles)% hObjecthandle to edit8 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handlesstructure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit8 as text%str2double(get(hObject,'String')) returns contents of edit8 as adouble% --- Executes during object creation, after setting all properties.function edit8_CreateFcn(hObject, eventdata, handles)% hObjecthandle to edit8 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handlesempty - handles not created until after all CreateFcns called% Hint: edit controls usually have a white background on Windows.%See ISPC and COMPUTER.if ispc && isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');endfunction edit9_Callback(hObject, eventdata, handles)% hObjecthandle to edit9 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handlesstructure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit9 as text%str2double(get(hObject,'String')) returns contents of edit9 as adouble% --- Executes during object creation, after setting all properties.function edit9_CreateFcn(hObject, eventdata, handles)% hObjecthandle to edit9 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handlesempty - handles not created until after all CreateFcns called% Hint: edit controls usually have a white background on Windows.%See ISPC and COMPUTER.if ispc && isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');134endfunction edit1_Callback(hObject, eventdata, handles)% hObjecthandle to edit1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handlesstructure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit1 as text%str2double(get(hObject,'String')) returns contents of edit1 as adouble% --- Executes during object creation, after setting all properties.function edit1_CreateFcn(hObject, eventdata, handles)% hObjecthandle to edit1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handlesempty - handles not created until after all CreateFcns called% Hint: edit controls usually have a white background on Windows.%See ISPC and COMPUTER.if ispc && isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');endfunction edit2_Callback(hObject, eventdata, handles)% hObjecthandle to edit2 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handlesstructure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit2 as text%str2double(get(hObject,'String')) returns contents of edit2 as adouble% --- Executes during object creation, after setting all properties.function edit2_CreateFcn(hObject, eventdata, handles)% hObjecthandle to edit2 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handlesempty - handles not created until after all CreateFcns called% Hint: edit controls usually have a white background on Windows.%See ISPC and COMPUTER.if ispc && isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');endfunction edit3_Callback(hObject, eventdata, handles)% hObjecthandle to edit3 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handlesstructure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit3 as text%str2double(get(hObject,'String')) returns contents of edit3 as adouble% --- Executes during object creation, after setting all properties.function edit3_CreateFcn(hObject, eventdata, handles)135% hObject% eventdata% handleshandle to edit3 (see GCBO)reserved - to be defined in a future version of MATLABempty - handles not created until after all CreateFcns called% Hint: edit controls usually have a white background on Windows.%See ISPC and COMPUTER.if ispc && isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');endfunction edit4_Callback(hObject, eventdata, handles)% hObjecthandle to edit4 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handlesstructure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit4 as text%str2double(get(hObject,'String')) returns contents of edit4 as adouble% --- Executes during object creation, after setting all properties.function edit4_CreateFcn(hObject, eventdata, handles)% hObjecthandle to edit4 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handlesempty - handles not created until after all CreateFcns called% Hint: edit controls usually have a white background on Windows.%See ISPC and COMPUTER.if ispc && isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');endfunction edit5_Callback(hObject, eventdata, handles)% hObjecthandle to edit5 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handlesstructure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit5 as text%str2double(get(hObject,'String')) returns contents of edit5 as adouble% --- Executes during object creation, after setting all properties.function edit5_CreateFcn(hObject, eventdata, handles)% hObjecthandle to edit5 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handlesempty - handles not created until after all CreateFcns called% Hint: edit controls usually have a white background on Windows.%See ISPC and COMPUTER.if ispc && isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');end136function edit10_Callback(hObject, eventdata, handles)% hObjecthandle to edit10 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handlesstructure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit10 as text%str2double(get(hObject,'String')) returns contents of edit10 as adouble% --- Executes during object creation, after setting all properties.function edit10_CreateFcn(hObject, eventdata, handles)% hObjecthandle to edit10 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handlesempty - handles not created until after all CreateFcns called% Hint: edit controls usually have a white background on Windows.%See ISPC and COMPUTER.if ispc && isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');end% --- Executes on button press in pushbutton3.function pushbutton3_Callback(hObject, eventdata, handles)% hObjecthandle to pushbutton3 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handlesstructure with handles and user data (see GUIDATA)load W1.mat -asciiload W2.mat -asciiload W3.mat -asciiload W4.mat -asciiload W5.mat -asciiload t1.mat -asciis1='g'; %зеленаяs2='r'; %краснаяs3='m'; %фиолетоваяs4='b';% синийs5='--'; %пунктирнаяfigureylabel('С, моль/м3')xlabel('t,с')hold onplot(t1,W1,s1)%%метан %%зеленаяhold onplot(t1,W2,s2)%%хлорметан %%краснаяhold onplot(t1,W3,s3)%% дихлорметан %%фиолетоваяhold onplot(t1,W4,s4)%% трихлорметан %%синияhold onplot(t1,W5,s5)hold on%%тетрахлорметан %%пунктирнаяlegend('метан', 'хлорметан', 'дихлорметан', 'трихлорметан','тетрахлорметан')hold on137% --- Executes on button press in pushbutton4.function pushbutton4_Callback(hObject, eventdata, handles)% hObjecthandle to pushbutton4 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handlesstructure with handles and user data (see GUIDATA)load TT.mat -asciiload t1.mat -asciis2='r'; %краснаяfigureylabel('T,K')xlabel('t,с')hold onplot(t1,TT,s2)hold on138Приложение В(обязательное)СВИДЕТЕЛЬСТВАО РЕГИСТРАЦИИ ПРОГРАММЫ ДЛЯ ЭВМ139140.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
