Scientific journal
Modern problems of science and education
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 1,006

METHYLCHLOROSILANE VAPOR AND FLUID HEAT CAPACITY CALCULATION IN EXPERIMENTAL DATA ABSENCE

Pererva O.V. 1 Endovin Yu.P. 1 Sokolov N.M. 1
1 FGUP GNIIChTEOS
Optimum calculation methods were determined allowing to calculate temperature dependence of methylchlorosilane vapor and liquid phase heat capacity in the absence of experimental data with an accuracy sufficient for engineering evaluation. It is shown that for the prediction of methylchlorosilanes liquid heat capacity equation Rowlinson—Bondi provides satisfactory accuracy. When calculating methylchlorosilanes liquid heat capacity, calculation error does not exceed 10% for the reduced temperature below 0.7. To calculate methylchlorosilane vapor phase heat capacity by the additive method, adjusted values of bond components contributions were proposed. Methylchlorosilane vapor heat capacity calculation error with the use of the proposed contributions values is below 5% within the temperature range from 300 to 1000К. Temperature dependence equations coefficients of fluid and vapor heat capacity for high boiling products of methylchlorosilane direct synthesis were found for the first time.
temperature dependence.
vapor heat capacity
fluid heat capacity
methylchlorosilanes
chlorosilanes
Для расчета массообменного и теплообменного оборудования производства метилхлорсиланов требуется знание теплофизических свойств веществ. Для ряда низкокипящих метилхлорсиланов данные по теплоемкости жидкости и пара известны и доступны, в том числе в виде баз данных общеупотребимых пакетов моделирующих программ (например, Aspen или ChemCAD [5]). Для высококипящих метилхлорсиланов экспериментальные данные по теплоемкости жидкости и пара отсутствуют. В настоящее время возрастающее значение приобретает глубокая переработка продуктов, получаемых по известным технологиям синтеза метилхлорсиланов. Если ранее многие побочные продукты рассматривались как отходы, подвергаемые уничтожению, и данных по теплоемкости для этих веществ не требовалось, то сейчас в связи с необходимостью выделения и переработки этих продуктов требуется знание всех теплофизических свойств, в том числе теплоемкости жидкости и пара.

Наиболее надежно использовать экспериментальные методы определения теплоемкости (например, калориметрические измерения). Однако задача осложняется тем, что большинство веществ, теплоемкость которых не известна, трудно получить в чистом виде. Кроме того, таких соединений может быть слишком много, что делает весьма трудоемкой задачу их экспериментального получения, очистки и определения свойств.

Задачей данной работы является на примере низкокипящих метилхлорсиланов определить применимость некоторых известных эмпирических расчетных методов для определения теплоемкости жидкости и пара и определить коэффициенты соответствующих уравнений для расчета свойств высококипящих метилхлорсиланов. В качестве источника данных по теплоемкости низкокипящих метилхлорсиланов использовалась база данных программы ChemCAD [5].

Расчет теплоемкости газа. Достаточно высокую точность расчета теплоемкости газа показали аддитивно-групповые методы [4]. Для соединений, групповые вклады для которых определены на основе экспериментальных значений теплоемкости при различных температурах, погрешность расчета теплоемкости такими методами в большинстве случаев оказывается ниже 5%. Групповые вклады известных методов расчета теплоемкости газа табулированы для различных органических соединений. Возможен учет в структуре соединений таких атомов, как O, Cl, N, S. Данные по атому Si для расчета свойств кремнийорганических соединений отсутствуют.

Для кремнийорганических соединений известен упрощенный вариант аддитивно-группового метода [2]. Согласно этому методу теплоемкость газа метилхлорсиланов при заданной температуре определяется как аддитивная сумма вкладов заместителей при одном атоме кремния. Табулированы вклады заместителей в интервале температур от 300 К до 1000 К. В качестве заместителей рассмотрены такие группы, как Н, Me, Et, Ph и Cl. Проверка точности этого метода для основных метилхлорсиланов и хлорсиланов показала, что погрешность расчета составляет не более 5%. На основании данных работы [2] можно рассчитать теплоемкость различных производных моносилана. Однако при наличии в молекуле двух и более атомов кремния данный метод оказывается неприменимым, так как отсутствуют опубликованные значения вкладов для таких структурных элементов, как Si–Si, Si–O–Si и Si–CH2–Si.

С целью расширения возможностей метода, изложенного в работе [2], нами был выполнен пересчет групповых вкладов на вклады составляющих связей. При этом теплоемкость газа при заданной температуре будет определяться как аддитивная сумма вкладов составляющих связей для данной температуры. Для пересчета использовались данные работ [2] и [3] и известные температурные зависимости теплоемкости газа для кремнийорганических соединений [5]. Рассчитанные значения вкладов составляющих связей приведены в таблице 1.

Таблица 1

Значения вкладов составляющих связей для расчета теплоемкости идеального газа метилхлорсиланов

тип связи

вклад кал/моль.К при

300 К

400 К

500 К

600 К

800 К

1000 К

С–С

2,06

2,585

2,975

3,235

3,585

3,730

С–Н

1,72

2,1825

2,6373

3,0575

3,7425

4,305

Si–H

1,9605

2,6377

3,283

3,8325

4,5862

5,0117

Si–Cl

5,43

5,8383

6,0565

6,1775

6,3258

6,4463

Si–C

3,587

3,7785

3,9083

3,997

4,1363

4,0815

Si–O

3,004

3,2861

3,5531

3,7558

3,8434

3,9801

Si–Si

3,706

3,5232

3,467

3,46

3,3112

2,8562

Погрешность расчета теплоемкости газа для метилхлорсиланов при использовании полученных вкладов составляющих связей не превышает 5% при температуре 300 К. При температурах более 400 К погрешность расчета теплоемкости газа менее 3%. Исключением является диметилхлорсилан (при 300 К погрешность расчета составляет 7,4%; при 400 К погрешность расчета составляет 4,9%).

Данные по теплоемкости газа метилтрихлорсилана (MeSiCl3), приведенные в базе данных программы ChemCAD [5], кардинально отличаются от всех данных, как опубликованных в литературе, так и рассчитанных по аддитивной схеме.  Поскольку источник данных в [5] не указан, использовать теплоемкость метилтрихлорсилана из базы [5] представляется нерациональным без дополнительной проверки. Для расчета температурной зависимости  теплоемкости газа метилтрихлорсилана целесообразно применять данные работы [3].

С помощью приведенных в таблице 1 значений вкладов составляющих связей нами была впервые рассчитана температурная зависимость теплоемкости газа Cp(T) некоторых высококипящих кремнийорганических соединений, образующихся в ходе прямого синтеза метилхлорсиланов. Полученные данные аппроксимировались следующими уравнениями:

, кал/моль.К                                                          (1)

, Дж/кмоль.К                                  (2)

где            Т – температура, К;

                 A, B, C и D – коэффициенты уравнения

Для инженерных расчетов уравнения (1) и (2) являются равнозначными с точки зрения точности получаемого результата. Выбор одного из этих уравнений для расчета теплоемкости газа определяется возможностями используемых программных средств.

Коэффициенты уравнений (1) и (2) приведены в таблицах 2 и 3. Интервал применимости от 300 К до 1000 К.

 

Таблица 2

Коэффициенты уравнения (1) для расчета теплоемкости идеального газа высококипящих продуктов прямого синтеза метилхлорсиланов (кал/моль.К)

вещество

A

B

С

D

Me3Si–SiClMe2

23,83785

0,10655

-3,18082e-005

-4,31235e-009

Cl3Si–O–SiCl3

21,05719

0,08549

-1,02807e-004

4,29051e-008

Me2ClSi–O–SiClMe2

21,00980

0,12264

-7,17417e-005

1,69133e-008

MeCl2Si–O–SiCl2Me

21,03350

0,10406

-8,72744e-005

2,99092e-008

Me3Si–SiCl2Me

23,84970

0,09727

-3,95746e-005

2,18559e-009

Me2ClSi–SiCl2Me

23,86155

0,08798

-4,73410e-005

8,68353e-009

(SiCl3)2

23,89709

0,06011

-7,06400e-005

2,81774e-008

Me2ClSi–SiClMe2

23,84970

0,09727

-3,95746e-005

2,18559e-009

Me2ClSi–O–SiCl2Me

21,02165

0,11335

-7,95081e-005

2,34113e-008

MeCl2Si–SiCl2Me

23,87339

0,07869

-5,51073e-005

1,51815e-008

Me3Si–CH2–SiClMe2

24,72109

0,13160

-5,01095e-005

1,98157e-009

Me3Si–CH2–SiCl2Me

24,73294

0,12231

-5,78759e-005

8,47951e-009

Me2ClSi–CH2–SiClMe2

24,73294

0,12231

-5,78759e-005

8,47951e-009

Me2ClSi–CH2–SiCl2Me

24,74479

0,11302

-6,56422e-005

1,49775e-008

MeCl2Si–CH2–SiCl2Me

24,75663

0,10373

-7,34086e-005

2,14754e-008

(SiCl3)2SiCl2

35,26003

0,07492

-8,56537e-005

3,27786e-008

 

Таблица 3

Коэффициенты уравнения (2) для расчета теплоемкости идеального газа высококипящих продуктов прямого синтеза метилхлорсиланов (Дж/кмоль.К)

вещество

A

B

С

D

E

Me3Si–SiClMe2

160729,37

322114,11

1226,34

158492,04

527,81

Cl3Si–O–SiCl3

146106,25

52339,68

1544,11

74037,27

685,85

Me2ClSi–O–SiClMe2

158812,53

285567,24

1223,75

156225,70

533,62

MeCl2Si–O–SiCl2Me

167568,34

172897,66

1694,42

131917,40

749,94

Me3Si–SiCl2Me

148726,47

288311,09

1407,02

187048,78

600,90

Me2ClSi–SiCl2Me

159273,14

200867,65

1297,90

112260,18

599,28

(SiCl3)2

141538,57

30021,38

1455,42

52214,24

697,37

Me2ClSi–SiClMe2

148726,47

288311,09

1407,02

187048,78

600,90

Me2ClSi–O–SiCl2Me

177221,64

212627,01

1589,22

152329,87

744,22

MeCl2Si–SiCl2Me

153009,47

146507,99

1326,81

94088,08

617,96

Me3Si–CH2–SiClMe2

181704,34

375840,23

1207,44

174576,48

506,53

Me3Si–CH2–SiCl2Me

170221,50

329170,90

1283,66

187173,22

553,89

Me2ClSi–CH2–SiClMe2

170221,50

329170,90

1283,66

187173,22

553,89

Me2ClSi–CH2–SiCl2Me

194674,00

256654,01

1642,18

171334,59

756,38

MeCl2Si–CH2–SiCl2Me

174692,13

198369,78

1374,28

127491,78

634,05

(SiCl3)2SiCl2

199698,29

38151,03

1450,54

67049,02

700,64

 

Расчет теплоемкости жидкости. Наиболее простыми способами расчета теплоемкости жидкости являются аддитивно-групповые методы и методы, основанные на принципе соответственных состояний [4]. Аддитивно-групповые методы весьма точны при условии, что для рассматриваемого вещества известны все групповые составляющие. Для кремнийорганических соединений определение групповых составляющих ранее никем не проводилось, что делает аддитивные методы неприменимыми.

Алексеев с сотрудниками [1] предложили неаддитивный метод расчета теплоемкости жидких органохлорсиланов. Считая, что зависимость теплоемкости жидкости большинства кремнийорганических соединений при приведенных температурах менее 0,9 близка к линейной, в работе [1] было предложено следующее уравнение:

                                                                      (3)

где             – теплоемкость жидкости, кДж/кг.К;

                 Т – температура, К;

                 М – молекулярная масса, г/моль;

                 – коэффициент, равный –0,1 для хлорсиланов и 0,74 для оксихлорсиланов;

                 – коэффициент, равный 18,0 для хлорсиланов и 15,5 для оксихлорсиланов.

Проверка показала, что данное уравнение имеет низкую точность расчета (табл. 4). Погрешность расчета теплоемкости жидкости составляет в среднем около 20% при приведенной температуре 0,55 и для большинства рассмотренных веществ увеличивается с ростом температуры. Примечательно, что в таблицах свойств, опубликованных авторами [1], для метилхлорсиланов с одним атомом кремния значения теплоемкости жидкости незначительно отличаются от данных [5] и не соответствуют результатам расчета по уравнению (3). Вероятно, уравнение (3) применимо для полиорганосилоксанов, но для более низкомолекулярных соединений (например, метилхлорсиланов) возможность его применения вызывает сомнения.

Таблица 4

Погрешность(*) расчета теплоемкости жидких кремнийорганических соединений

соединение

приведенная температура

0,55

0,7

0,85

0,55

0,7

0,85

погрешность по Алексееву

(уравнение (3), [1])

погрешность по Роулинсону—Бонди (уравнение (4), [4])

дихлорсилан

–26,9

–29,0

–24,8

+16,9

+20,7

+22,3

трихлорсилан

–48,5

–58,4

–58,2

+1,4

+7,1

+9,9

тетрахлорсилан

–48,8

–75,7

–100,9

+1,4

–0,6

–8,0

метилтрихлорсилан

–22,7

–36,9

–45,6

+7,3

+8,0

+6,4

метилдихлорсилан

–23,7

–19,1

–16,5

–1,3

+10,6

+13,8

диметилдихлорсилан

–13,4

–14,0

–10,4

+1,4

+9,4

+14,4

диметилхлорсилан

–1,5

–6,7

–11,0

+8,1

+7,6

+2,7

триметилхлорсилан

+1,8

+4,0

+9,0

–5,2

+3,4

+9,6

тетраметилсилан

+20,1

+15,4

+12,1

+0,1

–1,3

–6,0

гексаметилдисилоксан

–1,3

–4,7

–6,3

+2,6

+2,1

+0,3

(*) Погрешность = [(эксп. – расч.)/ эксп.] × 100%

К неаддитивным методам также относятся методы, основанные на принципе соответственных состояний. Погрешность расчета теплоемкости жидкости данными методами обычно составляет порядка 5–10%. Наиболее простым и удобным из данных методов является расчет по модифицированному уравнению Роулинсона—Бонди ([4], формула 5.8.2):

      (4)

где  – теплоемкость жидкости, кал/моль.К

      – теплоемкость газа, кал/моль.К;

     R – универсальная газовая постоянная,  кал/моль.К;

     ω – фактор ацентричности;

     Tr – приведенная температура ;

     Tc – критическая температура, К.

Проверка этого уравнения для метилхлорсиланов с известными теплоемкостями жидкости и газа показала (табл. 4), что метод имеет удовлетворительную точность в интервале приведенных температур ниже 0,7 (для большинства рассмотренных веществ погрешность менее 10%). В интервале  приведенных температур от 0,7 до 0,85 погрешность расчета возрастает. При приведенных температурах выше 0,85–0,9 погрешность превышает 10%, поэтому уравнение Роулинсона–Бонди нами не рекомендуется использовать в этих условиях.

С помощью уравнения (4) нами была впервые рассчитана температурная зависимость теплоемкости жидкости некоторых высококипящих кремнийорганических соединений, образующихся в ходе прямого синтеза метилхлорсиланов. Полученные данные аппроксимировались с помощью уравнения (5).

,  Дж/кмоль.К                              (5)

Коэффициенты уравнения (5) приведены в таблице 5.

 

Таблица 5

Коэффициенты уравнения (5) для расчета теплоемкости жидких высококипящих продуктов прямого синтеза метилхлорсиланов (Дж/кмоль.К)

вещество

интервал применимости, К

коэффициенты уравнения (5)

min

max

A

B

С

D

E

Me3Si–SiClMe2

273

520

7,7282e005

-6,1936e003

27,0782

-0,0500

3,4674e-005

Cl3Si–O–SiCl3

273

520

7,5908e005

-6,1763e003

26,1619

-0,0484

3,3516e-005

Me2ClSi–O–SiClMe2

273

520

7,8023e005

-6,2883e003

27,5287

-0,0510

3,5422e-005

MeCl2Si–O–SiCl2Me

273

525

7,3357e005

-5,7777e003

24,7443

-0,0454

3,1241e-005

Me3Si–SiCl2Me

273

530

7,2303e005

-5,5816e003

24,0008

-0,0437

2,9941e-005

Me2ClSi–SiCl2Me

273

535

7,0708e005

-5,4066e003

22,9646

-0,0417

2,8390e-005

(SiCl3)2

273

540

6,7348e005

-5,1225e003

21,0672

-0,0380

2,5740e-005

Me2ClSi–SiClMe2

273

535

7,0111e005

-5,3042e003

22,7233

-0,0412

2,8002e-005

Me2ClSi–O–SiCl2Me

273

530

7,1354e005

-5,4805e003

23,5772

-0,0430

2,9401e-005

MeCl2Si–SiCl2Me

273

545

6,5853e005

-4,8314e003

20,1171

-0,0360

2,4172e-005

Me3Si–CH2–SiClMe2

273

535

7,2001e005

-5,2438e003

22,8899

-0,0414

2,8069e-005

Me3Si–CH2–SiCl2Me

273

550

6,7470e005

-4,6976e003

20,1609

-0,0360

2,4022e-005

Me2ClSi–CH2–SiClMe2

273

565

6,2603e005

-4,0827e003

17,3689

-0,0304

1,9950e-005

Me2ClSi–CH2–SiCl2Me

273

580

5,8933e005

-3,6559e003

15,2371

-0,0263

1,6958e-005

MeCl2Si–CH2–SiCl2Me

273

585

5,7570e005

-3,5290e003

14,4721

-0,0248

1,5944e-005

(SiCl3)2SiCl2

273

595

6,5122e005

-3,7492e003

14,4549

-0,0244

1,5499e-005

Выводы

Выполнена проверка применимости известных методов расчета теплоемкости жидкости и пара к определению свойств метилхлорсиланов. Показано, что в случае отсутствия экспериментальных данных по температурной зависимости теплоемкости метилхлорсиланов известные методы расчета могут быть использованы для прогнозирования необходимых свойств для соединений данного класса.

Для аддитивного расчета теплоемкости идеального газа метилхлорсиланов рассчитаны ранее не известные значения вкладов составляющих связей.

Впервые определены коэффициенты температурной зависимости теплоемкости жидкости и пара для высококипящих продуктов прямого синтеза метилхлорсиланов.

Рецензенты:

Гартман Т.Н., д.т.н., профессор, РХТУ им. Д.И. Менделеева, зав. кафедрой информатики и компьютерного проектирования, г. Москва;

Щербакова Г.И., д.х.н., ведущий научный сотрудник лаборатории 4, ФГУП «ГНИИХТЭОС», г. Москва.