Коэффициенты диффузии в многокомпонентных системах при испарении бинарных азеотропных смесей в инертный газ тема автореферата и диссертации по физике, 01.04.14 ВАК РФ

Введение

Глава 1. Коэффициенты диффузии многокомпонентных систем, методы их расчета.

1.1. Термодинамика необратимых процессов.

1.2. Расчет многокомпонентных коэффициентов диффузии по кинетической теории газов.

1.2.1 Метод Максвелла - Стефана.

1.2.2 Метод Уилка.

1.2.3 Метод Тура (метод линеаризации уравнений).

Глава 2. Экспериментальные методы исследования многокомпонентных коэффициентов диффузии.

2.1. Метод Стефана.:.

2.2. Метод Грэма - Максвелла.

2.3. Метод Лошмидта.

2.4. Хроматографический метод.

Глава 3. Определение равновесных составов в системах бинарная азеотропная жидкость - газ при нормальных условиях.

3.1. Выбор исследуемых веществ и их краткая характеристика.

3.2. Экспериментальная установка для определения равновесных составов многокомпонентных парогазовых систем.

3.3. Исследование изменения состава бинарной азеотропной смеси в процессе испарения и по пути диффузии.

Глава 4 Экспериментальная установка и методика измерения многокомпонентных коэффициентов диффузии паров бинарных азеотропных смесей в инертный газ.

4.1. Выбор метода исследования.

34.2. Описание экспериментальной установки.

4.3. Методика проведения эксперимента.

4.4. Расчетная формула и выбор системы отсчета.

4.4.1. Оценка погрешности опытов.

Глава 5. Результаты измерений многокомпонентных коэффициентов диффузии паров бинарных азеотропных смесей в инертный

5.1. Результаты измерений.

5.2. Сравнение опытных значений полученных многокомпо- 71 нентных коэффициентов диффузии с теоретическими.

Выводы.

Список используемой литературы.

К точности выполнения расчетов, связанных с процессами тепломассообмена, предъявляются все более высокие требования. Это вызвано рядом задач, решаемых при создании современной тепломассообменной аппаратуры, широким применением в расчетах ЭВМ, а также разработкой надежных методов инженерного расчета высокоэффективных тепломассобменных аппаратов и установок.

Эти задачи можно успешно решить, если будет получено более точное описание явлений и внутренних закономерностей, определяющих процессы массопереноса. В значительной мере это касается и описания лимитирующей стадии процесса переноса массы - процесса диффузии. Существующая в настоящее время практика расчета коэффициентов переноса массы в таких массообменных процессах, как ректификация, абсорбция, экстракция, сушка, [1, 2]базируется на использовании приближенных данных о коэффициентах диффузии.

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

Несмотря на то, что в общей теории диффузионных процессов к настоящему времени достигнуты известные успехи [3, 4, 5], теоретические и экспериментальные исследования в области диффузионного массопереноса в многокомпонентных системах далеки от завершенности. Это, с одной стороны, связано со сложностью математического описания процессов переноса в таких системах, а с другой - трудностями экспериментального характера, что обусловлено рядом обстоятельств. Во-первых, методики проведения опытов, 5 на основе которых можно было бы проводить расчет МККД, существуют, в основном, для газовых и жидких систем, во-вторых, для парогазовых систем отсутствуют расчетные уравнения, применимые к конкретному методу исследования.

В настоящее время наиболее полно решена задача динамики частиц с потенциалом взаимодействия "твердых сфер" [5-8], которая, несмотря на упрощенность потенциала, может служить первым приближением для описания свойств реальных газов и жидкостей.

Для описания реальных систем, очевидно, необходимо использование различного рода априорных, феноменологических предположений о динамике частиц, таких, например, как принцип независимости диффузионных потоков [9] или принцип молекулярного хаоса [5, 10, 11]. Подобные предположения позволяют значительно упростить математическое описание взаимодействующих частиц и получить конечные выражения для кинетических коэффициентов через динамические параметры системы.

Как показывает анализ работ, посвященных переносу массы, исследование диффузии в многокомпонентных системах представляет собой сложную задачу [12-24].

На основе строгой кинетической теории разработаны уравнения, описывающие перенос массы в многокомпонентных системах. Эти зависимости проверены только для распределения концентрации компонентов для газов близких по своему состоянию к идеальным [22, 23], в которых роль перекрестных коэффициентов в общем переносе массы мала. Для неидеальных систем и систем под давлением мы располагаем пока линейными соотношениями термодинамики необратимых процессов, для использования которых необходимо знание феноменологических коэффициентов.

В периодической литературе в настоящее время отсутствуют работы, в которых были бы приведены значения МККД для парогазовых систем. 6

Этим объясняется необходимость глубокого экспериментального и теоретического исследования многокомпонентной диффузии, результаты которого могут быть использованы для удовлетворения запросов практики и уточнения вопросов теории.

Целью данной работы является:

1. Определение равновесных составов в системах бинарная жидкость -пар и оценка распределения концентрации паров бинарной азеотропной смеси по пути диффузии.

2.Разработка методики экспериментального исследования и теоретически обоснованного метода расчета многокомпонентных коэффициентов диффузии паров бинарных азеотропных смесей в газ.

3.Проведение комплексных экспериментальных исследований для определения многокомпонентных коэффициентов диффузии для интервала температур 295. .ЗЗОК и нормальном давлении.

4.Численное моделирование прямых и перекрестных коэффициентов диффузии при испарении бинарных азеотропных смесей в инертный газ при температуре 295. .ЗЗОК и нормальном давлении.

Научная новизна заключается в следующем:

1.Определены равновесные составы в системе жидкость-пар для бинарных азеотропных смесей: метанол-бензол, этанол-бензол, н-пропанол-бензол.

2.Проведена оценка распределения концентраций паров бинарных азеотропных смесей метанол-бензол, этанол-бензол и н-пропанол-бензол по пути диффузии и установлена нелинейность распределения концентрации компонентов бинарного азеотропа в газ.

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

4.Впервые определены прямые (Dn, D22) и перекрестные (Di2, D21) коэффициенты диффузии при испарении бинарных азеотропных смесей в инертный газ.

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

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

Автор защищает методики определения равновесных составов систем жидкость-пар бинарных азеотропных смесей, результаты экспериментального исследования и метод расчета прямых и перекрестных коэффициентов диффузии в многокомпонентных парогазовых системах.

Личное участие. Все основные результаты работы получены лично автором под научным руководством профессора Я.Д. Золотоносова.

Апробация работы. Основные положения работы были представлены на ежегодных научно-технических конференциях КХТИ им. Кирова в период 1972 ~ 1975г., 1979., 1981г., 1983г., 1986, 1988, 1991; на V Всесоюзной конференции по теплофизическим свойствам веществ, г. Киев, 1974г.; на Российском национальном симпозиуме по энергетике, РНСЭ, г. Казань, 2001г., на XXVI Сибирском теплофизическом семинаре, г. Новосибирск, СО РАН, 8 июнь 2002г.; на Всероссийской школе - семинаре академика В.Е. Алемасова «Проблемы тепломассообмена и гидродинамики в энергомашиностроении», г. Казань, октябрь 2002г.

Публикации. По теме диссертации имеется 7 публикаций. Объем работы. Диссертация изложена на 103 страницах и состоит из введения, пяти глав, выводов, списка литературы из 100 наименований и приложений. Иллюстрационный материал содержит 17 рисунков и 6 таблиц.

Выводы

1. На основании анализа литературных данных установлено, что из всего многообразия разработанных на сегодня способов определения коэффициентов диффузии ? наиболее информативным и надежным для измерения многокомпонентных коэффициентов диффузии в газовой и парогазовой средах является стационарный метод Стефана.

2. На базе статического метода измерения равновесных концентраций бинарных паров,впервые,для интервала температур 295.330К и нормального давления получены экспериментальные и расчетные значения равновесных составов многокомпонентных парогазовых систем: метиловый спирт - бензол-азот, этиловый спирт - бензол-азот, н-пропиловый спирт - бензол-азот.

3. Разработаны методики и проведены экспериментальные исследования для определения области стабильности состава азеотропа и характера изменения концентрации паров испаряющейся бинарной азеотропной смеси по пути диффузии.

Установлено, что при испарении первых 10 весовых процентов бинарной азеотропной смеси, состав жидкой и паровой фаз практически не меняется, а концентрация паров бинарной азеотропной смеси по пути диффузии меняется нелинейно.

4. Посредством стационарного метода Стефана разработана методика, создана экспериментальная установка и проведены измерения многокомпонентных коэффициентов диффузии паров бинарных азеотропных смесей в инертный газ.

Показано, что рассматриваемые системы являются постоянно испаряющимися, а предложенная методика может быть использована для измерения МККД.

81

5. На базе феноменологических уравнений диффузии описан процесс многокомпонентной диффузии паров бинарной азеотропной смеси в инертный газ, определен характер изменения концентрации компонентов по пути диффузии и впервые рассчитаны многокомпонентные коэффициенты диффузии: прямые - Dn, D22 и перекрестные - D12, D21, представленные в виде таблиц.

6. Проведено сравнение прямых коэффициентов диффузии - Дни Q22с бинарными коэффициентами, известными из литературных источников и рассчитанными нами. Показано, что сходимость прямых коэффициентов матрицы МККД с бинарными коэффициентами диффузии не превышает 12%.

82

1. Рид Р., Праусниц Дж., Шервуд Т. Свойства газов и жидкостей. JL: Химия, ЛО, 1982-592с.

2. Шервуд Т., Пигфорд Р., Уилки Ч. Массопередача. М.: Химия, 1982 - 696с.

3. Резибуа П., ДеЛенер М. Классическая кинетическая теория жидкостей и газов. М.: Мир, 1980 - 423с.

4. Балеску Р. Равновесная и неравновесная статическая механика. М.: Мир, т.1, 1978-405с.

5. Гиршфельдер Дж., Кертисс И., Берд Р. Молекулярная теория газов и жидкостей. М.: ИЛ, 1961 - 930с.

6. Van Beijeren, Ernst М.Н. The modified Enskog equation for mixtures//Physica. 1973. v.70.-p. 225-242.

7. Haro M.V., Cohen E.G.D. The Enscog theory for multicomponent mixtures. I Lineor transport theory//J.Chem. Phys. - 1983. v.78, №5. - p. 2747-2759.

8. Van Beijeren. Equlibzium distribution of hard spheres systems and revised Enskod theory//Phys. Lett. - 1983. v.51, №7 - p. 1503-1504.

9. Ферцигер Дж., Капер Г. Математическая теория процессов переноса в газах. М.: Мир, 1976. - 554с.

10. Чепмен С., Каулинг Т. Математическая теория неоднородных газов. М.: ИЛ, 1960-510с.

11. Больцман Л. Лекции по теории газов. М.: ГИТТЛ, 1953 - 556с.

12. Wilke C.R.//Chem. Eng. Progr. 1950. v.46. - p.95.

13. Fairbanks D.F., Wilke C.R.//Ind. Eng. Chem. 1950. v.42. - p.471.

14. Gilliland E.R., Sherwood T.K.//Absorbtion and Exstraction//New York 1937. -p.ll.

15. Toor H.L.//A.I.Ch.E.Journal. 1957. v.3, №2. -p.198.

16. Benedict M., Boas A. // Chem. Eng. Progr. 1951. v.47. p.51, p.l 1183

17. Cichelli M.T., Weathezford, Bowman J.R. 11 Chem. Eng. Progr. 1951. v.47.p.63, p. 123

18. Han Hsien Wen, Bird R.B. I I A.I.Ch.E. Journal. - 1960. v.6. p.516, p.551

19. Shain S.A. //A.I.Ch. E. Journal. -1961.v.7.p.l7

20. Cussler E.L., Lightfoot E.M. // A.I.Ch. E. Journal. 1963. v.9. p.703, p.783

21. Stewart W.E., Prober L. // Ind.Eng.Chem.Fund. 1964. v.3. p.224

22. Toor H.L. // A J. Ch. E. Journal. 1964. v. 10.p.448, p.460

23. Toor H.L. //Altech.- 1962-1963. v.12. p.21-43,p.51-61

24. Берд P., Стьюарт В., Лайтфут E. Явления переноса. М. Химия, 1974-688с.

25. Булатов Н.К., Лундин А.Б. Термодинамика необратимых физико-химических процессов. М.: Химия, 1984-334с.

26. Бережной А.Н., Бикбулатов А.Ш. Установка для определения коэффициентов диффузии паров в сжатых газах //Тр. Казан, хим.-технолог. ин-та. 1969. - вып.43. - с. 47-49.

27. Де Гроот С.Р. Термодинамика необратимых процессов. М.: Гостехиздат, 1956.-370 с.

28. Де Гроот С.Р., Мазур П. Неравновесная термодинамика. М.: Мир, 1964. -457 с.

29. Зубарев Д.Н. Неравновесная статистическая термодинамика. М.: Наука, 1971.-416с.

30. Гуров К.П. Феноменологическая термодинамика необратимых процессов. -М.: Наука, 1978.-128 с.

31. Dunlop P.I., Gosting L.J. //Jann. phys. chem. 1959, v.63, p.86

32. Baldwin P., Dunlop P., Gosting L. //J. Am. Chem. Soc. 1955, v.77 p.5235.

33. Fugita, Gosting L. //J. Am. Chem. Soc. 1956, v.78 p.1099.

34. Ирисов A.C. Испаряемость топлив для поршневых двигателей и методы ее исследования. -М.: Гостоптехиздат, 1955, 306 с.

35. Onda К., Sada Е., Nagasaka М. // Kaguku Kogaku. 1971, v. 35, № 9, p. 62.84

36. Теплопередача /А.Гриб, Б.Уокил // Труды амер. об-ва инж. механиков. 1969. № 2- с.62.

37. Горюнова Н.А., Кувшинский Е.Б. // ИТЖ. 1948. т. XVIII, вып. II. - с. 162.

38. Клибанова Ц.М., Померанцев В.В., Франк-Каменецкий Д.А. // ЖТФ. -1942. т. 12, вып. 1. -с.14.

39. Stefan / J. Sitzungsberg. Akad. Wies. Wien. Abt. 2 1871, Bd. 63. 1873, 68, 385.

40. Stefan / J. Ann. Phys. 1890. № 41, p.725.

41. Кричевский И.Р., Хазанова H.E. // Газовая промышленность. 1962. № 2. -c.29.

42. Maxwell G. /с. Scientific Papers, Dover Publications, N.Y. 1952. № 2. p. 57.

43. Keyes J.J., Pigford R.L. / Chem. End. Sci. 1957, v.6, № 3 - 6. p. 215.

44. Loschmidt / J. Akad. Wiss. Wien. 1870, Bd. 61, Abt. 2, № 3, s. 367 380; Bd. 62, Abt. 2, № 7, s. 468-478.

45. Obermayer.// Akad. Wiss. Wien. 1880, Bd. 81, Abt. 2, № 5, s. 1102 1127; 1882, Bd. 85, Abt. 2, № 1 - 5, s. 147 - 168; 1883, Bd. 87, Abt. 2, № 1, s. 188 -263.

46. Lonius / A.Ann. Phys. 1909. Ziff. № 29. p. 664. Ziff. № 13. p. 1414.

47. Giddings J.C. /J. Chem. Phys. 1959. v. 31. p.1462.

48. Giddings J.C., Seager S.L. /J. Chem. Phys. 1961. v. 35. p.2242.

49. Giddings J.C., Seager S.L. /Ind. End. Chem. Fund. 1962. v. 1. p,277.

50. Huber. / J.F.K.Metodes Phys. d'Analise. 1968. v.4 № 3. p. 285.

51. Физика химических растворов / Ким С.А., Сазонов M.JL, Лозгачев В.И., Жуховицкий А.А.// 1972. с.250

52. Лозгачев В.И., Канчеева О.А. // 1972, ЖФХ, t.XVI. № 5. с. 1242

53. Поливанов М.А., Дьяконов С.Г. Расчет коэффициентов бинарной диффузии в плотных газах. // Труды КХТИ им.Кирова С.М. Казань., 1974. Вып.5385

54. Walker R.E., de Naas N., Westenbergs. // A.A.J.Chem.Phys. 1960. v.32. p.1314

55. Maxwell G.C. On the dynamical theory of gases // Phys. Trans. Roy. Soc. -1867. vol.157, №= l,p.49-88

56. Maxwell G.C. On the dynamical theory of gases //Philos.Mag. 1868. vol.35, S = 236, p.185-217

57. Экспериментальные методы определения коэффициентов взаимной диффузии в газовой фазе: Обзфинформ. / Бережной А.Н., Семенов А.В.; Госстандарт; ВНИЦ MB; Казан.хим-технол. ин-т. М.: Изд-во стндартов, 1986.-39с.

58. Winkelmann. A. Ann. Phys. und Chem.,XXII, 1, 1884, XXII, 154,1884, XXIII, 203, 1884, XXIV, 105, 1885, XXXIII, 485, 1888.

59. Хорсли Jl. Таблицы азеотропных смесей. М: ИЛ, 1951-292с.

60. Огородников С.К., Лестева Т.М., Кочан В.Б. Азеотропные смеси, справочник. Л.О., Химия, 1971 .-848с.

61. Коган В.Б., Фридман В.М., Кафаров В.В. Равновесие между жидкостью и паром (справочник). М.Л.: Наука, 1966.

62. Бретщнайдер С.Т. Свойство газов и жидкостей. М.Л.: Химия, 1966. -422 с.

63. Вревский М.С. Работы по теории растворов, Изд. АН СССР, 1953

64. Коновалов Д.П. Об упругости пара растворов. Wied. Ann.(3), 14. 35, 1881-с.219.; С.-Петербург, ЖРФХО, ч. химич., I, 16

65. Сторонкин А.В. Термодинамика гетерогенных систем. Изд. ЛГУ. 1967

66. Hellund E.J. Phys. Rev., 57, 4, 737, 1940.

67. Sazan A., Sing Y. Canad J. Chem., 44, 2222, 1966.

68. Srivastava B.N., Paul RPhysica, 28, 646,1962.

69. Cunningham R.S., Geankoplis С J. Ind and Eng. Chem. Fund, 7, 3, 429, 1968.

70. Curyiss C.F.J. Chem. Phys., 49, 7, 2917, 1968.

71. Onsager L.Phys. Rev., 37, 405, 1931.86

72. Cussler E.L.J. Chem. Phys., 44, 7, 2829,1966.

73. By Watts H.Canad.J.Chem., 43, 431, 1965

74. Stewart W.E., Prober L.Ind.Eng.Chem.Fund, 3, 224, 1964

75. Туревский E.H., Александров И.А. ИФЖ, XXIII, 6, 1012, 1972

76. Чермак Д. Диффузия в многокомпонентных системах, (перевод) АН ЧССР, Прага, 1996

77. Варгафтик Н.Б. Справочник по теплофизическим свойствам газов и жидкостей. М.: Наука, 1972, - 598 с.

78. Свентославский В.В. Азеотропия и полиазеотропия. М.: Химия, 1968. -348 с.

79. Новосад Н.Д., Косов Н.Д. Физика, вып.1, 130 Алма-Ата, 1970

80. Tarenko В.Т., Bicbulatov A.Ch., Zolotonosov Y.D. Measuring of the bilateral aseothropic mixture during the process of the evaporation // Proceedings of the Russian National Symposium on Power Engineering. Kazan. RNSPE, 2001. V.l. P. 234-236.

81. Косов Н.Д., Новосад Н.Д., Сб. Вопросы общей и «прикладной физики», Алма-Ата 1969.

82. Степанов В.В. Курс дифференциальных уравнений.- М.: ГТЛ, 1950 254с.

83. Таренко Б.И., Бикбулатов А.Ш., Золотоносов Я.Д. Измерение равновесных составов трехкомпонентных парогазовых смесей. // Известия вузов. Проблемы энергетики. 2002. №7 8. С. 36 - 42.

84. Туревский Е.Н. Александров И.А., Гореченков В.Г. Химия и технология топлив и масел № 5, с. 34 36, 1973.

85. Кириллин П.В., Шейдлин А.Е. Исследование термодинамических свойств веществ. Гоэнергоиздат, 1963.

86. Бриджмен П.В. Физика высоких давлений, ОНТИ, М. Л., 193 5

87. Бриджмен П.В. Новейшие работы в области высоких давлений, И.Л., 1948.

88. Дьяконов С.Г., Казанцев С.А. / Определение коэффициентов диффузии в многокомпонентных смесях методом молекулярной динамики// Теоретические основы химической технологии. М.: 1981, т. XV. С. 20 - 24.

89. Бережной А.Н. /Экспериментальное определение, обобщение и прогнозирование характеристик молекулярного переноса массы в газах: Автореф. дис. . д-ра техн. Наук. М., 1989. - 30с.

90. Суворов А.В. Термодинамическая химия парообразного состояния. М.: Химия. 1970.-238 с.

91. Богомолов В.И., Миначев Х.М. Молекулярная хроматография. М.Л.: Наука. 1964. 138 с.

92. Скляренко С.И., Баранаев М.К. К вопросу о скорости испарения жидкостей // ЖФХ 1935.- Том VI, вып. 9. - С. 1180-1191.

93. Баранаев М.К. Кинетика испарения // Успехи химии -1938.-Том VII, №8,. -С. 12311259.891. DEFINTI-N

94. DECLARE FUNCTION SUM (KD, XF!()) DECLARE SUB MINF (NX, KD, XF(), ST!, EP, FMAX!) DECLARE SUB FPC (II, hh, XF()) DECLARE FUNCTION GODNUP% () DECLARE SUB MENU (IM%, MEN%, TREG$) DECLARE FUNCTION GOYN% ()

95. DECLARE SUB BBODKOR (NNN, AAA(), TEXTCOL$(), NKOL, NSTR, HAG, ZAGOL$) DECLARE SUB SLU (ISLU%, NSLU%, ASLU(), BSLU(), XSLU()) CONST MAXNK = 20 DIM SHARED NEXP, DET, SH(MAXNK), SHH(MAXNK)

96. KEY(15) ON KEY(16) ON KEY(17) ON KEY(18) ON90

97. FOR I = 1 TO NEXP: P(1,1) = C(1,1) * AM(1) * .001: P(2,1) = C(2,1) * AM(2) * .001:

98. P(3,1) = C(3,1) * AM(3) * .001: NEXT

99. BEG: CLS 0: CLS 1: SCREEN 9,1: WIDTH 80, 25:1. COLOR 15, 1

100. FOR I = 1 TO NEXP: DAA(1,1) = TC(I): DAA(I, 2) = h(l): DAA(I, 3) = dh(I): DAA(I, 4) = t(l):

101. DAA(I, 5) = PL(I): DAA(I, 6) = P(l, I): DAA(I, 7) = P(2,1): DAA(I, 8) = P(3,1):1. NEXT1. BEGIN: LNFN = LEN(DDl)

102. CLS 0: CLS 1: SCREEN 9,1: WIDTH 80, 25: COLOR 15,1

103. CATE 2,4: PRINT "Расчёт коэффициентов диффузии в трёхкомпонентной смеси путём аппроксимации» LOCATE3,1: PRINT " «Ввод экспериментальных данных"

104. CALL MENU(1, IREG1, TREG$)

105. ON IREG1 GOTO 1000, 2000, 3000, 4000, 5000 1000 : CLS 0: CLS 1:: SCREEN 9,1: WIDTH 80, 25: CLS : COLOR 15, 1 LOCATE 2, 24: PRINT "Назначение и краткая теория" LOCATE 4, 25:

106. PRINT "D11,D12,D21,D22 для средних температур Т по данным с близкими температурами." PRINT " Расчёт коэффициентов диффузии производится путём решения дифференциальных уравнений,"

107. PRINT " связывающих потоки масс первых двух компонент J1,J2 с градиента -" PRINT " ми их массовых плотностей С1,С2 (при C0=Cl+C2+C3=const) относительно. D11-D22"

108. PRINT" (J 1 * Со-J1 * С1 J2 * С1 )/Co=-D 11 *dC 1 /dx-D 12 *dC2/dx"

109. PRINT" (J2 * Co-J2* C2-J2* C2)/Co=-D21 * dC 1 / dx-D22 * dC2/dx"91

110. PRINT " где Jl=dh*PL*AN(l)/t,J2=dh*PL*N2/t. При решении используются граничные условия."

111. PRINT " С(1,0)С10, С(2,0)=С20, C(l,h)=0, C(2,h)=0 , C0=C(3,h)"

112. PRINT" (Jl*Co-Jl*Cl-J2*Cl)/Co=-Dll*dCl/dx-D12*dC2/dx"

113. PRINT" (J2 * Co-J2* C2-J2* C2)/Co=-D21 * dC 1 / dx-D22* dC2/dx"

114. Для окончания ввода или просмотра данных нажмите <F10>."

115. PRINT PRINT PRINT PRINT PRINT PRINT PRINT PRINT PRINT PRINT PRINT PRINT N1."

116. PRINT " При необходимости результаты расчётов можно напечатать на принтере." COLOR 4, 7: LOCATE 20, 1: PRINT " РАБОТА с программой производится с помощью меню и функциональных клавиш."

117. CATE 21, 1: PRINT" ИНФОРМАЦИЯ о возможных ситуациях и действиях сообщается по ходу работы."

118. CATE 22,1: PRINT " Для аварийного выхода в меню можно нажать Ctrl+S+Enter или Esc+Enter."

119. COLOR 15, 1: ON GODNUP% GOTO 3000,1000 3000:CLS 0: CLS 1: SCREEN 9,1: WIDTH 80,25: COLOR 15, 1 LOCATE 2,29: PRINT " РАБОТА С БАЗОЙ ДАННЫХ " LOCATE 3,18: PRINT " Выберите режим работы, используя.и (Enter ) "

120. РАСЧЁТ КОЭФФИЦИЕНТОВ ДИФФУЗИИ"

121. PRINT: CALL MENU(2, IREG1, TREG$)

122. ON IREG1 GOTO 3100, 3200, BEGIN 3100 : CLS 0: CLS 1: SCREEN 9, 1: WIDTH 80,25: CLS : COLOR 15, 1 LOCATE 1, 24: PRINT " ПЕРВИЧНЫЙ ВВОД ДАННЫХ"

123. CATE 3, 7: INPUT "Введите имя файла Ваших данных (до пяти символов согласно DOS)»; NFILES

124. CATE 5, 7: INPUT "Введите число экспериментов"; NEXP NFILES = LEFT$(NFILE$, 5) OPEN "R", #1, NDIR$ + NFILES + ".TAP", LNFN: PUT#1, 1, DD(1) CLOSE #1: KEY(10) ON1. ON KEY(10) GOSUB KONBB31

125. CALL BBODKOR(l, DAA(), TEXTCOL$(), 8, NEXP, 10, "Ввод результатов экспериментов")1. GOTO KONB311. KONBB31: KEY(10) OFF1. GOTO KONB31:1. RETURN

126. KONB31: OPEN "R", #1, NDIR$ + NFILES + ".TAP", LNFN: FOR 1=1 TO NEXP

127. CATE 5, 2: PRINT "В памяти; NOIPS; имеются следующие файлы типа ????.ТАР":1. PRINT1. FILES NDIRS + "?????.ТАР"3210 : LOCATE 12,2: INPUT ;" Введите имя файла, который хотите просмотреть (первые символы до точки)"; NFILES NFILES = LEFT$(NFILE$, 5)

128. OPEN "R", #1, NDIRS + NFILES + ".TAP", LNFN: NEXP = LOF(l) / LNFN FOR 1=1 TO NEXP GET #1,1, DD(I)

129. DAA(I, 1) = DD(I).TC: DAA(I, 2) = DD(I).h: DAA(I, 3) = DD(I).dh: DAA(I, 4) = DD(I).t: DAA(I, 5) = DD(I).PL: DAA(I, 6) = DD(I).P1: DAA(I, 7) = DD(I).P2: DAA(I, 8) = DD(I).P3: NEXT CLOSE #1: 3220 : KEY(10) ON ON KEY(10) GOSUB KONBB33

130. CALL BBODKOR(l, DAA(), TEXTCOL$(), 8, NEXP, 10," ВВОД РЕЗУЛЬТАТОВ ЭКСПЕРИМЕНТОВ"93

131. GOTO KONB33 KONBB33: KEY(IO) OFF GOTO KONB33: RETURN

132. KONB33: OPEN "R", #1, NDIR$ + NFILE$ + ".TAP", LNFN: FOR 1=1 TO NEXP

133. CATE 1, 19:PRINT "ПРОВЕДЕНИЕ РАСЧЁТОВ СИСТЕМЫ ПО ДАННЫМ В ФАЙЛАХ" LOCATE 3,2: PRINT "Вы можете проводить проводить расчёты по системам загруженным ранее."

134. CATE 5, 2: PRINT "В памяти "; NDIR$; "имеются следующие файлы типа ?????.ТАР": PRINT1. FILES NDIR$ + "?????.ТАР"

135. CATE 12, 2: INPUT ; "Введите имя файла, по которому будете считать (первые символы до точки) NFILES NFILE$ = LEFT$(NFILE$, 5)

136. OPEN "R", #1, NDIRS + NFILE$ + ".TAP", LNFN: NEXP = LOF(l) / LNFN FOR 1= 1 TO NEXP GET #1,1, DD(I)

137. DAA(I, 1) = DD(I).TC: DAA(I, 2) = DD(I).h: DAA(I, 3) = DD(I).dh: DAA(1,4) = DD(I).t: DAA(I, 5) = DD(I).PL: DAA(I, 6) = DD(I).P1: DAA(I, 7) = DD(I).P2: DAA(I, 8) = DD(I).P3: NEXT CLOSE #1:

138. FOR I = 1 TO NEXP: TC(I) = DAA(I, 1): h(I) = DAA(I, 2): dh(I) = DAA(I, 3): t(I) = DAA(I, 4):

139. PL(I) = DAA(I, 5): P(l, I) = DAA(I, 6): P(2,1) = DAA(I, 7): P(3,1) = DAA(I, 8):1. NEXT4001 : CLS 0: CLS 1:: SCREEN 9,1: WIDTH 80,25: CLS : COLOR 15,1

140. CATE 1,19: PRINT "ПРОВЕДЕНИЕ РАСЧЁТОВ СИСТЕМЫ ПО ДАННЫМ ФАЙЛА " + NFILES

141. PRINT : LOCATE 4, 2: INPUT "ВВЕДИТЕ во сколько раз начальный шаг меньше, чем коэффициент(4, 3,2 )

142. CATE 5, 12: INPUT " ВВЕДИТЕ N1 процент содержания первого компонента "; N1, N2= 100-N1.

143. FOR I = 1 TO NEXP -1: CLS 0: CLS 1: SCREEN 9,1: WIDTH 80, 25: CLS : COLOR 15,1

144. BSLU(l) = (PJ(1,1) * POB(I) (PJ(1,1) + PJ(2,1)) * P(l, I)) / POB(I)

145. BSLU(2) = (PJ(1,1 + 1) * POB(I + 1) (PJ(1,1 + 1) + PJ(2,1 + 1)) * P(l, I + 1)) / POB(I + 1)

146. ASLU(1, 1) = -dcx(l, I): ASLU(1,2) = -dcx(2,1):

147. ASLU(2, 1) = -dcx(l ,1+1): ASLU(2, 2) = -dcx(2,1 + 1)

148. CALL SLU(ISLU, 2, ASLUQ, BSLUQ, XSLUQ)94

149. D11(1) = XSLU(l): D12(I) = XSLU(2):

150. PRINT" Погрешность: (C(l,0)-C10)/C10, (C(2,0)-C20)/C20 для двух температур " CALL FPC(I, 0, XFO)

151. OHHl(I) = (FP(1) P(l, I)) / P(l, I): OHH2(I) = (FP(2) - P(2,1)) / P(2,1) CALL FPC(I+ 1,0, XF())

152. DH11(1) = XF(1): DH12(I) = XF(2): DH21(I) = XF(3): DH22(I) = XF(4) PRINT " С учётом нелинейности: СКО, ШАГ= "; FMAX; ST PRINT DH11(I); DH12(I); DH21(I); DH22(I)

153. PRINT" Погрешности (C(l,0)-C10)/C10, (C2,0 C20)/c20 для двух температур" CALL FPC(I, 0, XF())

154. OHl(I) = (FP(1) P(l, I)) / P(l, I): OH2(I) = (FP(2) - P(2,1)) / P(2,1) CALL FPC(I+ 1,0, XF())

155. PRINT" Графики для";: COLOR IC1, 0: PRINT " CI";: COLOR IC2,0: PRINT "

156. C2 ";: COLOR IC3, 0: PRINT " C3 "

157. XM = (h(I) + h(I + 1)) / 2: YM = (P(3,1) + P(3,1 + 1)) * 1.2 / 2:

158. VIEW (30, 50)-(600, 300), 2,11. WINDOW (0, YM)-(XM, 0)

159. NE (0, 0)-(0, XM), 0,, &HAAAA

160. NE (0, 0)-(0, YM), 0,, &HAAAA RR1. XM/100

161. CIRCLE (0, P(l, I)), RR, IC1 CIRCLE (0, P(l, I + 1)), RR, IC1 CIRCLE (0, P(2,1)), RR, IC2 CIRCLE (0, P(2,1 + 1)), RR, IC2951. NX= 1000: HX = XM/NX

162. CIRCLE (NX * HX, P(3,1)), RR, IC3

163. CIRCLE (NX * HX, P(3,1 + 1)), RR, IC31. FOR KK = 0 TO NX1. X = KK * HX:1. CALL FPC(I, X, XFO)

164. FP(3) = POB(I) (FP(1) + FP(2))

165. PSET (X, FP(1)), IC1: PSET (X, FP(2)), IC2: PSET (X, FP(3)), IC3: NEXT KK

166. CATE 5, 6: PRINT "T="; TC(I); TC(I +1) LOCATE 33, 63: PRINT "X^"; h(I); М(см)" LOCATE 34, 63: PRINT "X="; h(I + 1); "(см)"

167. CLS 0: CLS 1: SCREEN 9,1: WIDTH 80, 25: CLS : COLOR 15, 1

168. CATE 2, 12: PRINT " Матрица коэффициентов многокомпонентной диффузии (мЛ2/с)" PRINT " ДЛЯ СРЕДНИХ ТЕМПЕРАТУР "

169. PRINT" Т Dll D12 D21 D22 СКО" FOR 1=1 ТО NEXP -1:

170. PRINT USING" тм ##.#М#аааа ##.####лллл ##.М##аааа ##.####лллл ##.##ллллм;

171. TC(I) + TC(I + 1)) / 2; DH11(1); DH12(I); DH21(I); DH22(I)'; SH(I) NEXT

172. PRINT " Ошибки % " FORI= 1 TO NEXP-1:

173. PRINT USING " ##.## Ш.Ш#АААА ##.####ЛЛАЛ ##М##ЛЛАА ##.Ж#ааАа (TC(I) +

174. TC(I + 1)) / 2; OHl(I) * 100; OH2(I) * 100; OH3(I) * 100; OH4(I) * 1001. NEXT

175. CATE 25, 15: PRINT " Вы будете печатать результаты? (Y/N)";1. ON GOYN% GOTO PR4,4003:

176. PR4: CLS : LS = 1: LC = 10: NC 8: PRINT1. PRINT " N1="; AN(1);"%"

177. PRINT ""; TEXTCOL$(l); TEXTCOL$(2); TEXTCOL$(3); TEXTCOL$(4); TEXTCOL$(5); TEXTCOL$(6); TEXTCOL$(7); TEXTCOL$(8) FOR 1ST = 1 TO NEXP:

178. PRINT USING "##.##AAAA"; DAA(IST, 1); DAA(IST, 2); DAA(IST, 3); DAA(IST, 4); DAA(IST,5.; DAA(IST, 6); DAA(IST, 7); DAA(IST, 8):1. NEXT

179. PRINT:PRINT: PRINT" T Dll D12 D21 D22(mA2/c)" FOR 1=1 TO NEXP -1:

180. PRINT USING " ##M ##.#mAAAA №.##№AAAA ##M##AAAA ##M##AAAA "; (TC(I) +

181. TC(I + 1)) / 2; DH11(1); DH12(I); DH21(I); DH22(I)1. NEXT

182. RINT:LPRINT:LPRINT " N1="; AN(1);"%"

183. RINT ""; TEXTCOLS(l); TEXTCOL$(2); TEXTCOL$(3); TEXTCOL$(4); TEXTCOL$(5); TEXTCOL$(6); TEXTCOL$(7); TEXTCOL$(8)961. FOR 1ST = 1 TO NEXP:

184. RINT USING "##.##AAAA"; DAA(IST, 1); DAA(IST, 2); DAA(IST, 3); DAA(IST, 4);

185. DAA(IST, 5); DAA(IST, 6); DAA(IST, 7); DAA(IST, 8):1. NEXT

186. RINT: LPRINT" T Dll D12 D21 D22(mA2/c)" FOR I = 1 TO NEXP-1:

187. RINT USING " ##.## ##.####ААЛА ##.####ЛАЛА ##.####AAAA ##.####AAAA(TC(I) +

188. TC(I + 1)) / 2; DH11(1); DH12(I); DH21(I); DH22(I)1. NEXT1.RINT CHR$(12)4003 : LOCATE 25, 14: PRINT " Вы будете менять N1 ? (Y/N) ";

189. ON GOYN% GOTO 4001, BEGIN:5000 : CLS 0: CLS 1 :LOCATE 2, 18: PRINT ""

190. CATE 5, 24: PRINT " Конец работы до свидания "1. CLOSE1. END1. KONES:

191. KEY(15) OFF: KEY(16) OFF: KEY(17) OFF: KEY(18) OFF: KEY(19) OFF: KEY(20) OFF: CLOSE : GOTO BEG RETURN

192. SUB BBODKOR (NBB, AAA(), TEXTCOL$(), NCOL, NSTR, HAG, ZAGOL$) LSTR = 1: IF NBB = 3 THEN LSTR = 5 LCOL = HAG

193. DIM TEXTSTR(4, 20) AS STRING * 24' LSTR 'DIM TEXTCOL$(NBB, NSTR) AS STRING * LSTR FOR I = 1 TO NSTR: TEXTSTR(1,1) = STR$(I): NEXT" SCREEN 0: COLOR 15, 1 ISTR = 1: ICOL = 1

194. BEGBB: CLS:LOCATE 1, 23: PRINT ZAGOL$

195. CATE 2, 1: PRINT " используя <стрелки> числа ,<Enter>- Введите или просмотрите данные в произволь-"

196. CATE 3,1: PRINT " ном порядке. Для окончания работы и записи данных НАЖМИТЕ <F10>."1. FOR ICO = 1 ТО NCOL:

197. CATE 6, LSTR + (ICO -1) * LCOL: PRINT " TEXTCOLS(ICO): NEXT

198. FOR ICO = 1 TO NCOL: FOR 1ST =1 TO NSTR:

199. CATE 7 1 + 1ST, LSTR + (ICO -1) * LCOL: PRINT USING "##.### AAAA"; AAA(IST, ICO): NEXT NEXT

200. BB4: LOCATE 4, 1, 0, 4, 7: PRINT SPACE$(78)

201. CATE 4, 20: PRINT "Полное значение текущей переменной"; AAA(ISTR, ICOL) LOCATE 7 -1 + ISTR, LSTR + (ICOL -1) * LCOL, 1, 5, 8: LOCATE,, 1 BB6: C$ = INKEY$: IF C$ = "" THEN GOTO BB6 IF LEN(C$) = 2 THEN C$ = RIGHT$(C$, 1) С = ASC(C$)

202. С = 72 AND ISTR > 1 THEN ISTR = ISTR 1: GOTO BB4

203. CATE 4, 1: PRINT SPACE$(78):LOCATE 4, 1, 1, 5, 8: LOCATE ,, 1 INPUT " Введите значение переменной "; AAA(ISTR, ICOL): GOTO BEGBB: END IF GOTOBB6 ENDSUB

204. SUB FPC (II, XX, XF()) hh = h(II) * .01 :X = XX * .01

205. F1=(B2* A12-B1 * A22)/A: F2 = (B1 * A21 -B2* All)/A DET = (A11 + A22) A 2 / 4 + (A12 * A21 A11 * A22) SELECT CASE DET CASE IS < 0

206. AL1 = (A11 + A22) / 2: AL2 = SQR(-DET)

207. BSLU(l) = -F1 / EXP(AL1 * hh): BSLU(2) = -F2 * A12 / EXP(AL1 * hh)

208. ASLU(1, 1) = COS(AL2 * hh): ASLU(1, 2) = SIN(AL2 * hh)

209. ASLU(2,1) = (AL1 A11) * COS(AL2 * hh) - AL2 * SIN(AL2 * hh)

210. ASLU(2,2) = AL2 * COS(AL2 * hh) + (AL1 All)* SIN(AL2 * hh)

211. CALL SLU(ISLU, 2, ASLU(), BSLU(), XSLU())

212. CK1 = XSLU(l): CK2 = XSLU(2)

213. CA1 = (CK1 * (AL1 A11) + CK2 * AL2) / A12:

214. CA2 = (CK2 * (AL1 All) - CK1 * AL2) / A12:

215. FP(1) CK1 * EXP(AL1 * X) * COS(AL2 * X) + CK2 * EXP(AL1 * X) * SIN(AL2 FP(2) = CA1 * EXP(AL1 * X) * COS(AL2 * X) + CA2 * EXP(AL1 * X) * SIN(AL2 CASE IS > 0

216. AL1 = (A11 + A22) / 2 + SQR(DET): AL2 = (A11 + A22) / 2 SQR(DET): CA1 = (AL1 - A11) / A12: CA2 = (AL2 - A11) / A12: BSLU(l) = -F1: BSLU(2) = -F2

217. ASLU(1,1) = EXP(AL1 * hh): ASLU(1, 2) = EXP(AL2 * hh)

218. ASLU(2, 1) CA1 * EXP(AL1 * hh): ASLU(2, 2) = CA2 * EXP(AL2 * hh)

219. CALL SLU(ISLU, 2, ASLUO, BSLU(), XSLU())

220. CK1 = XSLU(l): CK2 = XSLU(2)

221. FP(1) = CK1 * EXP(AL1 * X) + CK2 * EXP(AL2 * X) + F1

222. FP(2) = CK1 * CA1 * EXP(AL1 * X) + CK2 * CA2 * EXP(AL2 * X) + F2

223. PRINT "1"; CK1; CA1; AL1; Fl; "kl"

224. PRINT "2м; CK2; CA2; AL2; F2; "k2"1. CASE 01. AL = (A11 + A22) / 2

225. BSLU(l) = -Fl / EXP(AL * hh): BSLU(2) = -F2 / EXP(AL * hh) * A12

226. ASLU(1, 1) = 1: ASLU(1, 2) = hh

227. ASLU(2, 1) = AL A11: ASLU(2, 2) = 1 + (AL - A11) * hh

228. CALL SLU(ISLU, 2, ASLU(), BSLU(), XSLU())

229. CK1 = XSLU(l): CK2 = XSLU(2)

230. CA1 = (CK1 * (AL A11) + CK2) / A12: CA2 = CK2 * (AL - A11) / A12:

231. FP(1) = (CK1 + CK2 * X) * EXP(AL * X) + Fl

232. FP(2) = (CA1 + CA2 * X) * EXP(AL * X) + F21. END SELECT1. END SUB

233. FUNCTION GODNUP% LOCATE 25, 17:

234. PRINT " Листать далыне-PgDn, листать назад -PgUP"; FL: C$ = INKEY$: IF C$ = "" THEN GOTO FL С = ASC(C$)

235. LEN(C$) = 1 THEN GOTO FL C$ = RIGHT$(C$, 1): С = ASC(C$) IF С = 81 THEN GODNUP% 1: EXIT FUNCTION IF С = 73 THEN GODNUP% = 2: EXIT FUNCTION GOTO FL END FUNCTION1. FUNCTION GOYN%1. STATIC C$ YN: C$ = INKEY$

236. C$ = "Y" OR C$ = "y" OR C$ = "K" OR C$ ="-" THEN GOYN% = 1

237. ELSEIF C$ = "N" OR C$ = "n" OR C$ = "'" OR C$ = "в" THEN GOYN% = 2 ELSE GOTO YN END IF END FUNCTION

238. CATE 7, 13: PRINT "Выберите в меню режим работы, используя CHR$(24);""; CHR$(25);" Ё Enter)": PRINT FOR I = 1 TO N(IM%): LOCATE NSTR + 2 * 1,10: PRINT I;: NEXT FOR I = 1 TO N(IM%): LOCATE NSTR + 2 * 1,14: PRINT TM(IM%, I) 1=1:

239. LOCATE NSTR + 2 * I, 8, 1: PRINT "#" 6 C$ = INKEY$: IF C$ = "" THEN 6 IF LEN(C$) = 2 THEN C$ = RIGHT$(C$, 1) С = ASC(C$)

240. С = 13 THEN MEN% = I: TREGS = TM(IM%, I): EXIT SUB IF С = 72 THEN 1С = 1:1 = I 1: GOTO 5 IF С = 80 THEN 1С = I: I = I + 1: GOTO 5 GOTO 699

241. IF I = N(IM%) + 1 THEN 1=1 IF I = 0 THEN I = N(IM%) LOCATE NSTR + 2* 1С, 8, 1: PRINT SPACE$(1) GOTO 4 STOP END SUB

242. SUB SLU (ISLU, NSLU, ASLU(), BSLU(), XSLU()) SELECT CASE ISLU

243. CASE 1 FOR I = 1 TO NSLU -1 FOR J = I + 1- TO NSLU ASLU(J, I) = -ASLU(J, I) / ASLU(I, I) FOR К = I + 1 TO NSLU:

244. ASLU(J, K) = ASLU(J, K) + ASLU(J, I) * ASLU(I, K) NEXT К

245. BSLU(J) = BSLU(J) + ASLU(J, I) * BSLU(I) NEXT J NEXT I

246. XSLU(NSLU) = BSLU(NSLU) / ASLU(NSLU, NSLU) FOR I = NSLU 1 TO 1 STEP -1: HSLU = BSLU(I) FOR J = I + 1 TO NSLU: HSLU = HSLU - XSLU(J) * ASLU(I, J): NEXT J XSLU(I) = HSLU / ASLU(I, I) NEXT I CASE 2 MS = 0

247. RS = MS * ASLU(I, 0) LS ASLU(K, 0) ASLU(K, 0) = LS * ASLU(I, 0) + MS * ASLU(K, 0) ASLU(I, 0) = RS1. NEXT К NEXT I

248. FOR I = NSLU TO 1 STEP -1: MS = 0

249. FOR К = 0 TO NSLU -1 1: MS = MS + ASLU(0, NSLU - K) * ASLU(I, NSLU - K) NEXT К

250. ASLU(0,1) = (ASLU(1,0) MS) / ASLU(I, I) NEXT I

251. FOR К = 1 TO NSLU: XSLU(K) = ASLU(0, K): NEXT CASE 3100

252. DSLU = ASLU(1, 1) * ASLU(2, 2) ASLU(2,1) * ASLU(1, 2) XSLU(l) = (BSLU(l) * ASLU(2, 2) - BSLU(2) * ASLU(1, 2)) / DSLU XSLU(2) = (BSLU(2) * ASLU(1, 1) - BSLU(l) * ASLU(2, 1)) / DSLU END SELECT END SUB

253. FUNCTION SUM (II, XFO) SU = 01. CALL FPC(II,0, XF())

254. SU = SU + ((FP(1) P(l, II)) / P(l, II))л 2 + ((FP(2) - P(2, II)) / P(2, II))A 21. CALL FPC(II, h(II), XF())

255. FP(3) = POB(II) (FP(1) + FP(2))

256. SU = SU + ((FP(3) P(3, II)) / P(3, II)) A 21. CALL FPC(II+1,0, XF())

257. SU = SU + ((FP(1) P(l, II + 1)) / P(l, II + 1)) A 2 + ((FP(2) - P(2, II + 1)) / P(2, II + 1)) A 2 CALL FPC(II + 1, h(II + 1), XF()) FP(3) = POB(II + 1) - (FP(1) + FP(2)) SU = SU + ((FP(3) - P(3, II + 1)) / P(3, II + 1))A 2 SUM=SU/6 END FUNCTION

258. Завлабораторией оптимизации технологичес""" -----------""кокрасочных производств, к.т.н.

259. Старший научный сотрудник лаборатории аппаратуттп-фртиплпгического оформления процессов синтеза лаковых смол, к.т.н

260. Казанское ордена Октябрьской революции производственное объединение „Органический синтез"420051 Казань Телеграф «Синтез» телетайп № 168 с/с 9225403

261. Начальник производства PI^I^S^* ■ Недбавлюк Б.Е.