Аналитические зависимости для скорости процесса

Изучение кинетических закономерностей включает следующие этапы экспериментальное исследование зависимости скорости реакции от условий проведения процесса обоснование механизма процесса на основании полученных экспериментальных данных и вывод аналитических зависимостей скорости процесса от внещних условий на основании экспериментальных данных и принятого механизма процесса обоснование выведенной зависимости путем специально поставленных экспериментов. [c.26]

Рассмотрим процесс заполнения жидкостью всасывающей магистрали для нахождения аналитической зависимости скорости откачки вакуумных устройств в условиях эксплуатации. Из законов термодинамики следует, что при создании разрежения в постоянном объеме (в данном случае в объеме всасывающей магистрали), при скоростях течения газа, меньших скорости звука, расширение его происходит изотермически, т.е. по закону [c.35]

Прямая пропорциональная зависимость скорости процесса от концентрации катализатора в теории промежуточных соединений была аналитически получена для случая, когда концентрация катализатора в реакционной смеси незначительна по сравнению с концентрацией реагирующих веществ Л, В и т. д. [c.84]

Посмотрим теперь, какую информацию мы хотели бы получить из экспериментального уравнения скорости. Прежде всего подход к уравнению будет различным в зависимости от целей — должен ли быть оптимизирован процесс или выяснен механизм реакции. Если в первом случае требуется просто найти аналитическую зависимость скорости от концентрации, то во втором случае уравнение скорости должно дать представление о реакции на молекулярном уровне. Это подразумевает выяснение отдельных элементарных стадий или элементарных реакций, которые происходят одновременно и последовательно внутри общей реакции. Элементарные реакции бывают реакциями первого порядка, например внутримолекулярные перегруппировки или реакции разложения, а чаще всего реакциями второго порядка, когда молекулярное взаимодействие между двумя частицами прямо приводит к продукту. Большинство химических реакций можно представить как те или иные комбинации элементарных реакций. Обычно различают следующие простые комбинации 1) параллельные реакции (конкурирующие, одновременные) 2) последовательные реакции или последовательность реакций 3) обратимые реакции. [c.12]

Так как в химико-технологическом процессе всегда участвует несколько веществ (два минимум), то уравнение (У.16) соответственно усложняется, поскольку появляется несколько коэффициентов диффузии, плотностей и т. п. Поэтому аналитический расчет коэффициента массопередачи практически невозможен и для его определения в каждом конкретном случае нужна постановка специального эксперимента. Сложная зависимость ( .16) может быть упрощена исключением ряда переменных, если известно, в какой области идет процесс — диффузионной, кинетической или переходной. Лимитирующую стадию можно определить, изучая влияние параметров технологического режима на общую скорость процесса и. Если и возрастает с повышением температуры в соответствии с законом Аррениуса (рис. 46) и температурный коэффициент > 1,5, то, как прав 1ло, лимитирующая стадия — химическая реакция, и процесс идет в кинетической области. Если же и растет с увеличением скоростей потоков реагирующих фаз, то лимитирующая стадия — это массообмен между фазами, и процесс идет во внешнедиффузионной области. На рис. 46 показано влияние температуры и скорости газового потока на кинетику процесса в системе Т—Г для обжига, горения, газификации. Из рис. 46 видно, что в области низких температур скорость процесса резко повышается с ростом температуры, так как определяющей стадией служит химическая реакция. В области высоких температур скорость химических реакций настолько возрастает, что процесс переходит в диффузионную область и общую скорость процесса лимитирует степень турбулизации газового потока (пропорциональная скорости газа гш ). Такой вид кривых зависимости скорости процесса или выхода продукта от температуры и скоростей реагирующих фаз (или от степени их перемешивания) характерен и для других гетерогенных систем. [c.109]

В случае линейной формы задания последних членов в правых частях уравнений (3.23), (3.24) (например, для реакций первого порядка в изотермических условиях) задача (3.23)—(3.26) допускает аналитическое решение стандартными методами. При этом удобнее пользоваться постановкой задачи, которая вытекает из диагонализированной формы уравнений (3.19), (3.20) в результате применения к ним интегрального преобразования (3.22). В случае более сложной формы последних членов в правых частях уравнений (3.23), (3.24) (например, при нелинейной зависимости скоростей реакций от состава фаз или когда процесс протекает в неизотермических условиях) решение краевой задачи (3.23)—(3.26) целесообразно искать численными методами. [c.145]

При использовании проточного метода с неподвижным слоем катализатора в реакторе обычно допускают, что движение газа в слое катализатора отвечает режиму идеального вытеснения, т. е. пренебрегают радиальными градиентами давления, температуры, концентрации. Соответственно среднюю скорость процесса по высоте слоя Н или по времени контакта т (поскольку т пропорционально Н) определяют интегрированием кинетических уравнений (VI. 1) и (VI. 3). Аналитическое решение кинетических уравнений, как правило, возможно лишь с применением вычислительных машин. При их отсутствии прибегают к графическому дифференцированию зависимости х = /(т), что вносит погрешности. [c.284]

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

Для некоторых частных случаев возможны аналитические или численные решения уравнения (П.45), если известна зависимость скорости реакции г (сд, /) от параметров, определяющих процесс. Однако перенос газа только за счет молекулярной диффузии и условие неподвижности жидкости не являются типичным для промышленных реакторов, в которых одновременно протекают процессы диффузии, конвенции и химическая реакция. Для анализа реальных процессов принимают упрощенные модели, достаточно достоверно отражающие рассматриваемое явление и не требующие большого числа трудно определимых параметров. Наиболее простой и наглядно представляющей процесс переноса вещества из газа в жидкость является пленочная модель. [c.31]

В простейших случаях удается получить аналитическое решение исходных уравнений в явном виде иногда для этого приходится вводить в модель упрощающие допущения. При этом решение описывает объект, но из-за сделанных допущений оно не всегда удовлетворительно. Многие процессы описываются сложнейшими математическими зависимостями. Папример, уравнения химической кинетики, учитывающие зависимость скорости реакции от температуры, нелинейны и часто неинтегрируемы. Многие задачи не имеют аналитического реи ения, т. е. результаты нельзя получить в виде формулы или набора формул. В принципе эти задачи можно решить, например, численным интегрированием, но практически это невыполнимо из-за огромного объема вычислений. [c.322]

Оптимальный размер (диаметр) зерен катализатора da определяют исходя из кинетических н гидродинамических соображений. По данным кинетических опытов или аналитическим расчетом находят максимальный размер зерна, при котором внутридиффузионное торможение существенно не влияет на скорость процесса. В зависимости от пористости зерна и от соотношения скоростей диффузии и реакции для данной реагирующей смеси допусти.мые размеры зерна катализатора колеблются от 0,1 до 3 мм. Однако по результатам последующих гидродинамических расчетов, прежде всего гидравлического сопротивления слоя, максимально допустимый размер зерна может оказаться неоптимальным. Линейная скорость газа в слое, а следовательно, и высота слоя (при заданной объемной скорости) повышаются с увеличением размера зерна. Поэтому для уменьшения гидравлического сопротивления слоя необходимо применять более мелкие зерна, чем максимально допустимые по кинетике процесса. Но для равномерного распределения газа необходима минимальная высота слоя зерен (50—100 мм), и если по расчетам (например, для активного катализатора) получается слой меньшей высоты, то следует увеличивать и высоту, и, соответственно, размеры зерна. Окончательно размеры зерна уточняются после определения гидравлического сопротивления, и линейной скорости газа в реакторе. [c.252]

Зависимость скорости срыва от давления, несомненно, обусловлена частично зависимостью числа Не от плотности. В настоящее время еще не совсем ясно, зависят ли от давления А и 1/т. Другой важной причиной зависимости от давления является суммарный порядок скорости химических процессов, которые в недостаточной степени представлены в уравнениях (14) и (15), прежде всего потому, что Х 11 подробно вычислено только для реакции первого порядка. В настоящее время предпринимаются попытки аналитически определить для химических реакций, порядок которых отличается от первого порядка. На основании представленной формулировки срыва мало что можно сказать о влиянии давления, пока не будет лучше выяснено влияние порядка скорости химической реакции. [c.187]

В политропических устройствах с непрерывным теплообменом р отличие от ступенчато-адиабатических вариантов характер распределения температур по пути следования реагирующего потока определяется в основном теплотехническими факторами. Температура в этих системах является переменной, не зависящей от глубины превращения, и изменение ее принципиально может быть задано тем или иным уравнением распределения температур по длине аппарата или во времени. Определение средних скоростей политропических процессов, в итоге сводящихся к установлению зависимостей г ол от или пол от со а, встречает ряд трудностей вследствие сложности формулирования аналитических зависимостей i = f( o) или 1 = Р(т). [c.105]

Для аппарата, реакционный объем которого вычисляется по формуле V =УЛС,, вместо скорости процесса и/ л ставится его выражение, причем X, 2 и Я определяются теми аналитическими зависимостями, которые нами получены (8), а входящая во все эти уравнения величина у (в данном случае это аргумент, относительно которого у нас найдены все функции) имеет значение, полученное в результате экспериментального распределения продуктов реакции при заданной степени преврашения целевого продукта. Аналогичным образом определяется выход и других продуктов реакции при заданной степени превращения. С точки зрения математической обработки расчет реакционного объема аппарата, работающего в режиме смешения, не вызывает трудности. Для сравнения приведем метод расчета аппарата вытеснения на мак-ро- и микроуровнях. Длина такого аппарата должна определяться из уравнения (66). В этом уравнении при [c.90]

Именно на основании данных о зависимости концентрации от времени С — I или зависимости скорости от концентрации V — С удается выяснить формальную кинетику процесса, т. е. зависимость скорости от аналитических концентраций исходных веществ. Зависимость У — I представляет реальный интерес лишь в случае превращений первого порядка, например радиоактивных превращений. [c.73]

Например, исходя из тепловой теории воспламенения получены аналитические зависимости [167], которые используются лишь для качественной оценки влияния различных физических факторов на процесс воспламенения. В частности, получена зависимость пространственной скорости распространения пламени и от удельной поверхности дисперсной фазы [c.101]

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

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

Таким образом, в настоящее время уже накапливается опыт промышленных методов сжигания газа в кипящем слое. Залогом дальнейшего совершенствования этих методов является создание теории горения газа в кипящем слое и прежде всего определение аналитической зависимости для эпюры распределения температуры в слое (что определяет как стойкость самих горелок, так и устойчивость от слипания всего слоя) в зависимости от калорийности газа, характера смесеобразования, выходных скоростей, фракции и удельного веса частиц, числа псевдоожижения и температуры слоя. Выявление этой зависимости позволит выбрать для конкретного технологического процесса оптимальный режим сжигания газа и рациональную конструкцию горелок для обеспечения этого режима. [c.15]

Аналитическое определение примесей в винилацетате не всегда дает достаточно данных для того, чтобы можно было судить о его поведении при полимеризации при заданном количестве катализатора. На практике лучшие результаты можно получить анализом кривых зависимости вязкости от времени полимеризации при определенной температуре и концентрации инициатора. В зависимости от характера кривых можно рассчитать то количество инициатора, которое необходимо для обеспечения заданной скорости процесса полимеризации. [c.277]

Многие физические и химические свойства вещества (например, электропроводность, люминесценция, радиационная стабильность и т. п.) зависят от степени его чистоты часто ничтожные примеси в концентрации 10 2—10 % резко изменяют наблюдаемые свойства вещества, что исключает возможность их применения во многих химических процессах и приборах. Необходимым условием проведения почти любого корректного научного исследования является применение только особо чистых соединений в условиях, исключающих их загрязнение. Поэтому анализ загрязнений в чистых соединениях и определение примесей в товарных продуктах в настоящее время является один из основных направлений развития современной аналитической химии [1]. К чистоте исходных веществ в реакциях полимеризации также предъявляются очень высокие требования. В табл. 1 [4] в качестве примера приведены требования к чистоте углеводородных мономеров. В зависимости от реакционной способности предельная концентрация примесей не должна превышать —10 %. Содержание примесей в растворителях и в других веществах, применяемых для реакции, не должно превышать, по-видимому, предельно допустимых норм для мономеров. Столь высокие требования по чистоте исходных веществ, естественно, предъявляются только в отношении вредных примесей, образующих нежелательные продукты или существенно уменьшающих скорость процесса. [c.327]

По другим представлениям происходит электроокисление молекул и радикалов гидразина. В МЭИ получено много косвенных данных об адсорбции гидразина сдвиг потенциала пика на вольт-амперной кривой с увеличением скорости развертки потенциала, экстремальная зависимость скорости реакции от концентрации гидразина. Кроме того, получены прямые данные об адсорбционных процессах при электроокислении гидразина аналитическим методом, методом анодных импульсов и методом окисления в адсорбционном слое. Методом палладиевой мембраны [Л. 81] было установлено, что при электроокислении гидразина происходит сорбция водорода палладием, причем содержание водорода в металле может расти с увеличением анодной поляризации. Водород был обнаружен в палладии в процессе электроокисления гидразина при достаточно высоких анодных потенциалах (до Ег=, 0 В). [c.128]

Из-за сложности процессов течения при смешении аналитическое описание напряженного состояния смеси пока невозможно, т. е. не может быть определена временная и координатная зависимости скоростей у и напряжений т сдвига. [c.83]

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

Учет всех этих явлений при анализе работы кристаллизационной колонны весьма затруднителен. Можно, однако, с целью получения хотя бы приближенных аналитических зависимостей принять во внимание и рассматривать лишь одну, основную стадию процесса. Такой основной стадией при незначительном эффекте перекристаллизации будет, по-видимому, по аналогии с ректификацией, диффузионный межфазовый массообмен. 3 соответствии с этим допущением примем далее, что скорость межфазового массообмена в кристаллизационной колонне лимитируется скоростью диффузии примеси б твердой фазе. [c.102]

Общая схема расчета параметров фильтрационного процесса при дв женин в пористой среде пластических жидкостей основана на двух осно ных предположениях. Первое касается описания течения флюида на ми роуровне - в элементарном капилляре (поровом канале) - и состоит в то что в качестве закона трения при движении флюидов в капиллярах прин мается некоторая аналитическая зависимость скорости сдвиговой дефо мации Ys от касательных напряжений Tj. При этом закон течения вязко несжимаемого флюида в капилляре имеет вид (У.Л. Уилкинсон, 1964 г.) [c.70]

Уравпеппя (66) и (67) могут применяться для чисто аналитических решений, когда известны функциопальные зависимости скорости от глубины процесса или времени реагирования, а также для графоаналитических расчетов, когда имеются только эмпирические 1 ннетпческие графики. В последнем случае по графикам предварительно находятся величипы т . — то, чие- [c.128]

II этапа с повышением температуры окисления значительно возрастают, что объясняется сдвигом равновесия в сторону диссоциации пачки. Изменение температуры окисления влияет не только на скорость процесса, но также на состав получаемого битума на его свойства. С использованием величины стехиометрических коэффициентов можно легко рассчитать количество кислорода (воздуха), необходимое дая получения асфальтенов заданной концентрации. Учитывая, что между их концентрацией в битуме и температурой размягчения последнего существует определенная зависимость, и используя данные по юшетике, можно выполнить аналитическое решение практических задач, связанных с расчетом процессов и аппаратов для получения битумов разных марок. [c.777]

При промышленном оформлении химических процессов аналитические зависимости w = P x) и х = (р(г) весьма часто бывают еще не достаточно изученными и разработку схем и конструкций реакционных устройств приходится базировать на простейших эмпирических данных, представленных в виде тех или иных кинетических графиков. В этих условиях нахождение Тр и Топ может легко производиться графическим методом. По асимптотам кривых X —т или X—ft), построенных для разных температур, находятся значения Тр, а графическим диференцированием — скорости образования целевых продуктов при различных х и Т. Представляя последние в координатах скорость—>температура в виде семейства кривых для ряда значений х= onst, по точкам их перегиба можно определить максимальные скорости и соответствующие им оптимальные температуры, как это показано, например, на фиг. 68. [c.242]

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

Такая форма записи удобна тец, что (скорости накопления концентраций промежуточных веществ) могут быть явно выражены в виде полиномов от концентраций реагентов. Эта особенность позво-ляет получать оценки величины Тпел, рассматриваемой как функции от независимых переменных хГ, причем такие оценки могут быть получены при отсутствии аналитических зависимостей Xi(i) и значений Х °°. Соотношения (47) и (48) при 1- °° совпадают, но последнее использовалось в [43] без допущения о близости процесса к стационарному состоянию. Установлено, что т ел в любой момент времени может иревышать т не более чем в п раз в случае двуста-дийпой нелинейной реакции п-то порядка. На примере трехстадийной реакции показано, что т ел может обращаться в бесконечность.. [c.221]

Для выделения адсорбированных компонентов адсорбент нагревают (внутри колонны имеется система нагревательных элементов, не показанных на схеме). Выделяющиеся компоненты откачивают через системы трубок А, В ш С. Применяя соответствующие вращающиеся переключатели, аналогичные тем, которые мы применили ранее для аналитических целей [102], внутри колонны можно создавать движущееся тепловое поле с некоторым градиентом температуры, обеспечивая таким путем перемещение и лучшее разделение компонентов. В зависимости от числа разделяемых компонентов или фракций и необходимой их чистоты компоненты откачивают через те или иные трубки систем А, В и С. Степень разделения смеси скорость процесса регулируют температурой нагрева при помощи реостата, а также путем подбора соответствующего режима работы нереключателя. Нагрев адсорбента можно проводить также путем впуска перегретого пара через трубки 1, 2 ж 3. Эту отпарку можно гочетать с действием нагревательных элементов. [c.150]

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

Однако ни в этой работе, ни в других исследованиях влияния температуры на протекание [28, 30, 32, 160, 161, 162[ органических реакций в начале 40-х годов не был изучен механизм этого процесса. Причинами этого, но нашему мнению, явились а) довольно слабая изученность в это время механизмов органических реакций и особенно их элементарных актов б) канонизация уравнения Аррениуса (17) как аналитической формы телшературной зависимости скоростей реакций, что превратило фактически это уравнение из объекта исследования в метод установления механизмов реакций. [c.72]

Сложность процесса и недостаточная его изученность исключают пока возможность установления строгой аналитической зависимости между всеми величинами, определяющими производительность валковых машин. Ниже приводятся некоторые данные окружных скоростей V (в м1мин), принятые в промышленности для валковых машин непрерывного действия [c.170]

Оптимальный размер (диаметр) зерен катализатора определяется исходя из кинетических и гидродинамических соображений. По данным кинетических опытов или аналитическим расчетам находят максимальный размер зерна, при котором внутридиффузи-онное торможение существенно не влияет на скорость процесса. В зависимости от пористости зерна и от соотношения скоростей диффузии и реакции для данной реагирующей смеси допустимые [c.276]

Таким образом, введение ряда ограничений и использование несксльких эмпирических зависимостей позволило получить для оптимальных условий не очень сложные аналитические зависимости (32.30) и (32.31). При этом оптимизация проводилась лишь по одному воздействию (скорость газа), которое было принято за основное. Остальные факторы, влияющие на процесс, были учтены при формулировании условий. Такая постановка вопроса несомненно, является упрощенной, но к ней прибегают довольно часто, так как в противном случае задача слишком усложняется. [c.192]

Одним из факторов, определяющих осаждение из газовой фазы, является скорость осаждения, которая в общем случае может испытывать кинетические, термодинамические и диффузионные ог-ранирения. Если при этом она подчиняется закону Аррениуса и пропорциональна первой степени концентрации реагирующего компонента, то ее аналитическая зависимость от параметров процесса в общем случае определяется как [1—3] [c.81]

При оценке эффекта разделения в ректификационной колонне еще Е. Сорелем [4, 166] было показано, что разделительная способность колонны, работающей в отборном режиме, всегда меньше, чем ее разделительная способность в безотборном режиме, если соблюдается постоянство всех остальных условий процесса. Знание этой закономерности имеет очень важное значение при решении практических вопросов конструирования и эксплуатации ректификационных установок. В общем виде эту зависимость можно установить методом расчета от тарелки к тарелке или графическими методами. Но, как уже отмечалось выше, в ряде случаев, особенно когда приходится иметь дело с колоннами, содержащими большое число тарелок, для этой цели удобнее пользоваться соответствующей аналитической зависимостью, полученной при допущениях, справедливых для заданных условий. Этим объясняется большое число работ, посвященных установлению аналитической взаимосвязи между скоростью отбора (флегмовым числом) продукта, его составом и числом теоретических тарелок колонны при использовании тех или иных допущений для различных типов колонн. В большинстве из этих работ вопросы применения ректификации для глубокой очистки веществ не рассматриваются, однако разработанные в них методы расчета можно использовать и для решения указанной задачи. В связи с этим большой интерес представляет сопоставление результатов, которые дают наиболее строгие из этих методов. С этой целью рассмотрим вначале метод, предложенный К. Торманом [2], который, пожалуй, является самым наглядным, будучи прямым аналитическим эквивалентом метода расчета от тарелки к тарелке . [c.73]

Уравнение (3.284) позволяет оценить концентрацию примеси в отобранном продукте в любой момент времени при известных параметрах процесса и величине Nd/Nq. Если же, наоборот, величина X]) задана, то с помощью уравнения (3.284) можно определить выход продукта, допустимую скорость отбора или время процесса. Нетрудно видеть, что если жидкостные захваты ректифицирующей части и конденсатора колонны незначительны и ими можно пренебречь, выражения (3.282) — (3.284) упрощаются и приводятся к виду известных из литературы соответствующих аналитических зависимостей, полученных другими способами [555, 560J. Так, например, выражение (3.282) в этом случае преобразуется в соотношение [c.142]

Справочник химика 21

Химия и химическая технология