Гидродинамический ламинарный поток

Механизм процесса переноса массы сводится к молекулярной и турбулентной диффузии. При молекулярной диффузии, происходящей в неподвижной фазе и ламинарном потоке, перенос массы характеризуется коэффициентом диффузии ), который рассчитывают по формулам (631)—для газов и (633)—для жидкости. При турбулентной диффузии перенос вещества осуществляется движущимися частицами среды и определяется гидродинамическим состоянием потока. Механизм переноса вещества через поверхность раздела фаз является кардинальным вопросом теории массопередачи и окончательно не решен. Предполагая, что диффузионные сопротивления в жидкой и газообразной фазах обладают свойством аддитивности, можно записать основное уравнение массопередачи [c.336]

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

Сопротивление при ламинарном потоке в трубах на участке гидродинамической стабилизации. Дтя ламинарного стабилизированного потока в круглой грубе профиль скорости описывается выражением (2.11), из которого получено уравнение (2.13) для величины касательных напряжений в зависимости от радиуса г. При г = Гв абсолютное значение касательных напряжений из (2.13) равно [c.76]

Теплоотдача при переходном режиме. Переходный режим в области значений Не = 2200 10000 существенно отличается от ламинарного режима гидродинамическими свойствами потока и механизмом переноса тепла коэффициент теплоотдачи при переходном режиме значительно выше, чем при ламинарном. Для переходной области непригодны уравнения (276) и (280). Коэффициенты теплоотдачи наиболее надежно определяют непосредственно из опыта, но приближенно их можно вычислить по функции В = /(Не, =-- [c.118]

Теория Поттера (пограничных слоев). Поттер [77] рассматривает молекулярную диффузию в жидкости, двигающейся упорядоченно, принимая за основу выводов гидродинамические отношения, т. е. относительное движение ламинарных потоков, двигающихся в том же направлении. Для такой модели массо-перенос определяется коэффициентом диффузии О в степени п, изменяющейся в зависимости от отношения количеств фаз ЕЩ. Показатель степени. имеет значения в пределах =0,33 0,5 [c.78]

Предположим, что мы имеем ориентированный по оси л капилляр радиуса г и длины I, наполненный жидкостью, к концам которого приложена разность потенциалов Е (рис. 30). Под влиянием электрического поля происходит электроосмотический перенос жидкости с некоторой скоростью причем в результате такого течения жидкости создается некоторая разность давлений Р. Описание движения вязкой, несжимаемой жидкости под влиянием электрического поля и при наличии гидростатического давления может быть сделано с использованием гидродинамических уравнений Навье—Стокса. Для данного случая — ламинарного потока жидкости в направлении оси л — в стационарном состоянии в соединении с уравнением несжимаемости жидкости уравнение Навье—Стокса сводится к следующему выражению [c.54]

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

Рис. 2.22. Гидродинамическая стабилизация потока в трубе при ламинарном (а) и турбулентном (б) течениях.

Некоторые характерные черты присущи теплообмену, связанному с ламинарным потоком через каналы с некруглыми поперечными сечениями. Этот случай изучался аналитически [Л. 96] для стабилизованного теплового и гидродинамического потоков через канал, поперечное сечение которого имеет форму сектора круга и для условия, при котором поток тепла от стенки канала в жидкость постоянен в направлении оси канала. Было найдено, что локальный коэффициент теплообмена значительно изменяется по периферии канала, приближаясь к нулевому значению в углах, и что средний коэффициент теплообмена во многом зависит от граничных условий. Были рассмотрены два граничных условия по окружности канала температура стенки, которая является постоянной по периферии, и локально постоянный тепловой поток. Найдено, что критерий Нуссельта, усредненный по окружности для постоянной температуры стенки, в 7 раз больше его величины для постоянного потока тепла, когда угол вершины сектора был равен 20°. Для угла у вершины в 60° соотношение этих двух чисел Нуссельта равно 2,5. Коэффициент теплообмена в числах Нуссельта определяется как осред-ненный тепловой поток у стенки, деленный на разность между объемной температурой жидкости и средней температурой стенки (осредненной по периферии канала). [c.251]

В ламинарном потоке, стабилизированном гидродинамически, но не стабилизированном в тепловом и диффузионном отношении, [c.260]

Входящее в выражение для ядро коагуляции /С(со, V), называемое иногда константой коагуляции, обусловливает частоту столкновения капель объемами со и У и может быть определено в результате исследования относительного движения двух капель под действием различных сил взаимодействия — гравитационной, гидродинамической, молекулярной. Характер гидродинамической силы зависит от структуры потока. В ламинарном потоке относительное движение капель различного размера происходит за счет гравитационного осаждения и градиента скорости несущей среды. При этом /С(со, V) определяется сечением столкновения капель и находится в результате анализа траекторий движения одной капли относительно другой [см. раздел 13.1]. В турбулентном потоке сближение капель происходит за счет хаотических пульсаций, приводящих к большему по сравнению с ламинарным потоком числу актов столкновения в единицу времени. Существуют три основных механизма, обусловливающих сближение и коагуляцию капель в турбулентном потоке инерционный, за счет различных скоростей движения отличных по размерам капель сдвиговый (градиентный), вызванный наличием сдвигового течения в окрестности рассматриваемой капли при обтекании ее пульсациями различного масштаба турбулентная диффузия, в основе которой лежит предположение об аналогии между процессом диффузии и движением капель под действием случайных турбулентных пульсаций. В разделе 13.6 показано, что применительно к рассматриваемому процессу основной вклад в скорость коагуляции капель в турбулентном потоке дает механизм турбулентной диффузии, и ядро коагуляции с учетом сил гидродинамического и молекулярного взаимодействия капель имеет вид [c.547]

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

Постановка задачи. Выбор системы координат. Рассмотрим трехмерную задачу о стационарной конвективной диффузии к поверхности твердой или жидкой частицы произвольной формы, обтекаемой ламинарным потоком вязкой несжимаемой жидкости. Как и ранее, предполагается, что число Пекле Ре = aUD велико здесь а — характерный размер частицы (в качестве которого обычно выбирается радиус эквивалентной по объему сферы а ), U — характерная скорость потока (на бесконечности), D — коэффициент диффузии. Считается также, что на поверхности частицы и вдали от нее концентрация принимает постоянные значения, равные нулю и Соо, а поле течения жидкости известно иа решения соответствующей гидродинамической задачи об обтекании частицы. [c.126]

Коснемся оценки погрешности описанного интегрального метода. Во-первых, здесь пмеет место определенная аналогия с интегральными методами, неплохо зарекомендовавшими себя в теории гидродинамического пограничного слоя [63, 64, 109]. Во-вторых, как следует из результатов этого параграфа, в общем случае произвольного обтекания при О для п = 2 погрешность метода составляет менее 6%. В-третьих, в задачах о диффузии к реагирующим каплям (пузырям) в ламинарном потоке вязкой жидкости и к частицам в идеальной жидкости (что соответствует значению г = 1 в выражении (7.5)) предлагаемый метод является точным. Это можно показать путем непосредственной подстановки формулы (7.5) в исходное уравнение (1.2), после чего придем к уравнению [c.320]

Эта формула считается универсальной, а все разнообразие гидродинамических состояний потока и соответственно режимов движения жидкости сводится к изучению различных закономерностей изменения функции Л = / (и), отвечающих ламинарному, переходному или турбулентному режимам. Так, в соответствии с формулами Шевелева [266] для металлических водопроводов, работающих в квадратичной области при и > 1,2 м/с, [c.28]

Приведенные длины участков гидродинамической и тепловой стабилизации ламинарного потока в трубах [c.164]

Характерный профиль скоростей в ламинарном потоке устанавливается не сразу после входа жидкости в трубу распределение скоростей асимптотически приближается к параболическому, описываемому уравнением (2.19). Трансформация скоростного профиля от начального до характерного происходит на участке гидродинамической стабилизации. Очевидно, что его длина / аб. в силу асимптотического характера преобразования профиля скоростей, зависит от приемлемой погрешности, с которой скоростной профиль можно считать сформировавшимся. [c.150]

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

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

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

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

Для ламинарного потока, установившегося в гидродинамическом но не установившегося в тепловом или диффузионном отношении (т.е. со сформировавшимся профилем скоростей, но не сформировавшимся профилем температур или концентраций), аналитическое решение бьшо найдено Левеком. Согласно закону Левека 118], зависимость между критериями подобия имеет вид [c.39]

Как видно из формулы (V, 81), для ламинарного потока с определенными гидродинамическими характеристиками локальное значение р может быть найдено расчетом [c.249]

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

Помимо решений, которые могут быть получены на основе понятия пограничного слоя, возможно также получение некоторых теоретических результатов в случае чисто вязкого режима взаимодействия ламинарного потока с поверхностью твердого тела, когда в уравнении Навье — Стокса (1.1) можно пренебречь инерционными членами. При этом оказывается возможным получить решение изотермической гидродинамической задачи стационарного обтекания, например тела сферической формы. Анализ такого решения показывает, что скорость вязкого потока, обтекающего сферу (или иное тело правильной геометрической формы), плавно возрастает по мере увеличения расстояния от поверхности тела, в непосредственной близости от которой при отсутствии инерционных членов в уравнении движения не существует какой-либо зоны значительных поперечных градиентов продольной скорости, т. е. гидродинамического пограничного слоя как такового не существует. Из этого обстоятельства следует, что при чисто вязком режиме обтекания не может быть аналогии между распределениями продольной скорости и концентрации целевого компонента, что рассматривалось выше в пределах пограничного слоя при Рг == 1. [c.35]

Правая часть уравнения представляет силы, действующие на частицу при ее ускоренном движении в вязком потоке. Первое слагаемое — это сила гидродинамического воздействия потока на сферическую частицу, записанное здесь для любого режима обтекания в частном случае ламинарного режима (Re < 0,5) может быть принято X, — 64/Re. Второе слагаемое включает силу тяжести и гравитационную подъемную силу. Третье слагаемое правой части называется увлекающей силой потока, которая действовала бы на массу жидкости m.f в объеме частицы [c.65]

При ламинарном режиме течения длина участка гидродинамической стабилизации потока Л может быть определена по приближенному соотношению [c.19]

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

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

Все многообразие процессов и явлений, наблюдаемых при трении твердых тел, заключено между трением ювенильных поверхностей и гидродинамическим трением. Под трением ювенильных (идеально чистых) поверхностей понимают трение поверхностей при полном отсутствии между ними третьей фазы, способной выполнять функцию смазочной среды. Термин гидродинамическое трение определяет процессы, происходящие в присутствии смазочной среды, поведение которой подчиняется законам гидродинамики ламинарного потока жидкости, в первую очередь уравнению Ньютона. Этот термин определяет процессы трения, характеризуемые вязкостью как важнейщим физико-химическим свойством смазочной среды. Между двумя указанными предельными состояниями фрикционной системы, т. е. между сухим и жидкостным трением, существует гранич1н0е трение , наблюдаемое в том случае, когда тонкий слой смазочной среды, разделяющий трущиеся поверхности, находится в границах их влияния на смазочное вещество. [c.223]

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

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

Проводя аналогию между процессами теплопередачи и диффузии, приходится отметить, что в теплопередаче гидродинамическое подобие потоков полностью характеризуется критерием Рейнольдса только при вынужденном движении с хорошо развитой турбулентностью ири отсутствип такого движ ения, а также в потоках ламинарных и переходных режимов перенос тепла за счет естеств( Нпой конвенции характеризуется критерием Грасгофа. Аналогичный по смыслу критерий введен и для диффузионных процессов [c.34]

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

Совместное влияние трех факторов (дуги 5—7)—механического перемещивания МехП), возмущений, вносимых с входными потоками (ВхП), н геометрических особенностей формы рабочего объема аппарата — приводит к формированию определенной топологии потоков в масщтабе аппарата (ФТПА). Топология потоков в масштабе аппарата, или гидродинамическая структура потоков в аппарате, определяется характером и расположением в пространстве его рабочего объема макрогидродинамических неоднородностей потоков застойных зон, байпасов, зон ламинарного и турбулентного течения, циркуляционных токов и т. д. [c.109]

Тот факт, что гидродинамические характеристики жидкостей весьма отличаются от таковых для газос, сохранения ламинарного потока в жидкостной термодиффузионной колонне удается добиться, уменьшив ншрину зазора до 0,2—0, i мм. При этом, как уже говорилось выше, такой размер должен быть строго соблюден по всей высоте колонны, стеики которой полируются . [c.413]

СЛОИ у стенок трубы обладают большей вязкостью, чем в основном ядре тотока, поэтому скоростное поле описывается кривой с. Температурное поле и теплообмен находятся в известной зависимости от изменения скоростного поля. Таким образо М, ко1эффициент теплообмена зависит. как от направления теплового потока, тш и от его величины. Расчет теплообмена в вязких жидкостях был Выполнен К. Ямагата [Л. 88]. Во-вторых, расчетные и опытные данные трудно сравнивать потому, что часто при низких скоростях, характерных для ламинарного потока, вихревые токи свободной конвекции изменяют ламинарный характер движения в результате получается сочетание свободной и вынужденной конвекции [Л. 89]. Этот вопрос будет рассматриваться в разделе 11-5. В-третьих, для масел участок полной гидродинамической и тепловой стабилизации настолько велик, [c.247]

Для частиц, размер которых превышает 0,1 мкм, (р,2)(игь > (Р12)ьголл-Приведенные в зтом разделе выражения для частот столкновения в процессах броуновской, сдвиговой и турбулентной коагуляции получены без учета гидродинамического молекулярного и электростатического взаимодействий частиц. Учет этих взаимодействий значительно осложняет задачу. В частности, в коэффициентах броуновской и турбулентной диффузии необходимо учитывать гидродинамическое сопротивление частицы с учетом искажения поля скоростей, вызванного присутствием соседних частиц, а в уравнении диффузии учитывать конвективный поток за счет сил молекулярного взаимодействия частиц. В случае градиентной коагуляции в ламинарном потоке необходимо рассматривать траектории относительного движения частиц с учетом гидродинамических и молекулярных сил взаимодействия. [c.220]

Рассмотрим движение пузырьков в ламинарном потоке. Их взаимодействие обусловлено, с одной стороны, разностью скоростей движения относительно жидкости за счет различных размеров, а с другой — молекулярными силами взаимодействия. За счет различных размеров происходит их сближение на относительно больших по сравнению с радиусами пузырьков расстояниях. На малых расстояниях возникают силы сопротивления, которые препятствуют сближению. На этих же расстояниях начинает действовать сила притяжения Ван-дер-Ваальса, которая обеспечивает эффективный захват пузырьков. Заметим, что если происходит сближение пузырьков с полностью заторможенной поверхностью, то сила гидродинамического сопротивления при малых зазорах 5 между поверхностями пузьгрьков сингулярна Ff, 8 , поэтому столкновение пузырьков невозможно без учета силы Ван-дер-Ваальса. При сближении пузырьков со свободной поверхностью В отличие от первого случая эта особенность [c.605]

Решение указанной задачи по теплообмену Лацко [364] для стабилизованного гидродинамически турбулентного потока, а также исследование Маргулиса [365], опыты Рубинштейна [101], Аладьева [162] и других показывают, что величина так называемого начального участка трубы, где происходит практическая стабилизация критерия Ни, для турбулентного потока в трубе гораздо меньше, чем для ламинарного. Поэтому предельным значением критерия Ни можно [c.294]

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

Изучалась диффузия аммиака в ламинарном потоке газовоздушной смеси к внутренней поверхности горизонтальной трубы, на которой был нанесен слой фосфорной кислоты. В опытах использовалась трубка с внутренним диаметром (1 = 30 мм, которая имела очень плавйый вход и. участок гидродинамической стабилизации длиной 67(1, т. е. поток газовоздушной смеси был полностью стабилизирован в гидродинамическом отношении. Длина исследуемого участка была принята равной 40д. Для э.тнх условий В = Ре, а коэффициент сопротивления Я рассчитывается по формуле (16). , [c.93]

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

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