ПРОГНОЗИРОВАНИЕ “ВСЛЕПУЮ” |
Стр. |
На этой странице рассмотрено предсказание комплекса свойств малоизученных веществ при полном отсутствии экспериментальной информации о свойствах исследуемого веществаа. |
|
В настоящее время прогноз свойств молекулярных неорганических жидкостей требует сведений хотя бы об одной экспериментальной характеристике исследуемого вещества – его температуре кипения. Поскольку методы предсказания самой температуры кипения в литературе отсутствуют, невозможен и прогноз физико-химических свойств “вслепую”, т.е. при полном отсутствии исходной информации об исследуемом веществе. Ниже приводится итерация, которая позволяет вычислить “вслепую” как температуру кипения исследуемого вещества, так и целый ряд его свойств.Предварительно рассчитывается парахор П .Сначала задается первое приближение для температуры кипения (Т) :T=200 и для молярного объема жидкости в точке кипения (V): V=27•ПЗатем рассчитываются два определяющих критерия термодинамического подобия: y = 0.1• ln(T) – 0.122• ln(V) +0.006 (1) f = ln(T) – 1.1• ln(П) – 4 (2) Далее следует определение температуры кипения во втором приближении и вычисление ряда физико-химических свойств: T кип=ехр[6.8 y + 1.1• ln(П) + 4.21]T к=1.5• TкипV кип= (– 6.2• f + 31.2)П0.94b=1.27• Vкип – 6 Pк=1.039Tк/b Vк=V1 (2.69 + 0.05• f) Δ Hкип=T1 [73 + 2.09• ln(T1•М/Vкип)]Z к=Pк•Vк/(8.314 •Tк)ρ 293 =(M/Vкип)[1 + 0.46(1 – 293/Ткип)]В конце итерации следует условие проверки заданной точности (0.1 К) abs(T кип – T1)<0.1Если оно выполняется, расчет заканчивается, если нет, то предыдущее значение Т заменяется на новое: T=Tкип и расчет повторяется, начиная с уравнения (1).Прогнозирование свойств “вслепую” проведем на примере оксихлорида ванадия VCl3OДля расчета парахора по стехиометрической формуле заглядываем в табл. 8 SITE.П (VCl3O) = П(V) + 3П(Cl) + П(О) – 0.338•m = 1.55 + 3•1.1 + 0.64 - 0.338•4 = 4.14Далее проводим расчет по итерационной схеме, которая представлена ниже на языке Бейсик. 407 T=200 408 V=27*P 410 F=.1*LOG(T)-.122*LOG(V)+.006 414 F1=LOG(T)-1.1*LOG(P)-4 415 T1=EXP(6.8*F+.21+1.1*LOG(P)+ 426 TK=1.5*T1 430 V1=(-6.2*F1+31.2)*P^.94 434 B=1.27*V1-6 437 PK=1.039*TK/B 440 VK=V1*(2.69+.05*F1) 442 H=T1*(73+2.09*LOG(T1*M/V1)) 446 ZK=PK*VK/(8.314*TK) 448 RR=(M/V1)*(1+.46*(1-293/T1)) 450 IF ABS(T1-T)<.1 GOTO 500 454 T=T1: GOTO 410 500 Вывод результатовРасшифровка обозначений программы .Т - первое приближение для температуры кипения, V - первое приближение для молярного объема жидкости в точке кипения, F - первый определяющий критерий термодинамического подобия, F1 - второй определяющий критерий термодинамического подобия. Далее в скобках указаны относительные ошибки (в процентах) результатов итерации для ряда свойств оксихлорида ванадия. T1 - температура кипения (4) , V1 - молярный объем жидкости в точке кипения (4), TK - критическая температура (-10), PK - критическое давление (-14), VK - критический объем (2), H - теплота испарения жидкости в точке кипения (-6), В - постоянная Ван-дер-Ваальса (4), ZK - критический фактор сжимаемости (-2), RR - плотность жидкости при Т=293 К (-5).Точность метода определена на массиве из 80 соединений элементов 3 – 8 групп Периодической системы. В их число вошли: гидриды, фториды, хлориды, бромиды, йодиды, смешанные галогениды, оксигалогениды. Средняя (по абсолютному значению) погрешность прогноза свойств составила (приводится в скобках в %): Ткип (8), Vкип (3), ΔHкип (10), Тк (8), Vк (3), Pк (10 ), b (4), zк (2), ρ293 (5).Обозначения (в скобках дается обозначение в тексте программы, составленной на языке Бейсик): М (М) – молярная масса, г/моль Ткип (Т1) – нормальная температура кипения, К y (F) – первый определяющий критерий термодинамического подобияf (F1) – второй определяющий критерий термодинамического подобия П (Р) – парахор, Дж 1/4см 5/2моль -1Vкип (V1) – объем жидкости в точке кипения, см 3/моль Тк (ТК) – критическая температура, К V к (VК) – критический объем, см 3/мольРк (РК)– критическое давление, МПаR - универсальная газовая постоянная (8.314 Дж/(моль.К)Δ Нкип (Н) – теплота испарения точке кипения, кДж/мольB (В) –постоянная Ван-дер-Ваальсаρ 293 (RR)– плотность жидкости при Т=293 К, г/см 3
|
home