Численное моделирование процесса перемешивания взаиморастворимых жидкостей в роторно-пульсационном аппарате

Rotor-pulse apparatus are successfully used for the creation of homogeneous mixtures from high-viscosity components. The fields of concentration of components are defined as a result of the numerical simulation of the liquid flow and the process of mass transfer in the working space of an apparatus....

Full description

Saved in:
Bibliographic Details
Date:2007
Main Authors: Басок, Б.И., Давыденко, Б.В., Ободович, А.Н., Пироженко, И.А.
Format: Article
Language:Ukrainian
Published: Видавничий дім "Академперіодика" НАН України 2007
Subjects:
Online Access:https://nasplib.isofts.kiev.ua/handle/123456789/1724
Tags: Add Tag
No Tags, Be the first to tag this record!
Journal Title:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Cite this:Численное моделирование процесса перемешивания взаиморастворимых жидкостей в роторно-пульсационном аппарате / Б.И. Басок, Б.В. Давыденко, А.Н. Ободович, И.А. Пироженко // Доп. НАН України. — 2007. — N 7. — С. 79–86. — Бібліогр.: 5 назв. — рос.

Institution

Digital Library of Periodicals of National Academy of Sciences of Ukraine
_version_ 1860268555152392192
author Басок, Б.И.
Давыденко, Б.В.
Ободович, А.Н.
Пироженко, И.А.
author_facet Басок, Б.И.
Давыденко, Б.В.
Ободович, А.Н.
Пироженко, И.А.
citation_txt Численное моделирование процесса перемешивания взаиморастворимых жидкостей в роторно-пульсационном аппарате / Б.И. Басок, Б.В. Давыденко, А.Н. Ободович, И.А. Пироженко // Доп. НАН України. — 2007. — N 7. — С. 79–86. — Бібліогр.: 5 назв. — рос.
collection DSpace DC
description Rotor-pulse apparatus are successfully used for the creation of homogeneous mixtures from high-viscosity components. The fields of concentration of components are defined as a result of the numerical simulation of the liquid flow and the process of mass transfer in the working space of an apparatus. The time dependence of the average concentrations in the input and output sections of the working space is investigated. The estimation of the efficiency of a rotor-pulse apparatus as a mixing device is made on the basis of the proposed criterion.
first_indexed 2025-12-07T19:03:55Z
format Article
fulltext оповiдi НАЦIОНАЛЬНОЇ АКАДЕМIЇ НАУК УКРАЇНИ 7 • 2007 ТЕПЛОФIЗИКА УДК 66.021.2:663.913.7 © 2007 Член-корреспондент НАН Украины Б.И. Басок, Б.В. Давыденко, А.Н. Ободович, И.А. Пироженко Численное моделирование процесса перемешивания взаиморастворимых жидкостей в роторно-пульсационном аппарате Rotor-pulse apparatus are successfully used for the creation of homogeneous mixtures from high- viscosity components. The fields of concentration of components are defined as a result of the numerical simulation of the liquid flow and the process of mass transfer in the working space of an apparatus. The time dependence of the average concentrations in the input and output sections of the working space is investigated. The estimation of the efficiency of a rotor-pulse apparatus as a mixing device is made on the basis of the proposed criterion. Роторно-пульсационные аппараты (РПА) как устройства для создания однородных смесей различных жидкостей широко применяются в пищевой, фармацевтической и других отрас- лях промышленности [1]. Примером успешного использования РПА как смесителя может быть технология гидролиза сахарозы для получения глюкозы и фруктозы. Как показали результаты исследований, представленные в [2], включение в технологическую схему гид- ролиза сахарозы аппарата роторного типа приводит к значительному ускорению данного процесса. Одна из причин указанного явления связана с созданием условий для интенсив- ного перемешивания в роторном аппарате компонентов, принимающих участие в рассмат- риваемом физико-химическом процессе. Компонентами в данном случае являются сахарный сироп и концентрированный водный раствор кислоты. По отношению к массе сахарного сиропа количество кислоты составляет от одного до нескольких процентов. В связи с тем, что компонентами, как правило, являются высоковязкие среды, созда- ние однородной смеси с целью обеспечения условий для интенсивного протекания реакции гидролиза оказывается задачей достаточно сложной. При использовании для этой цели традиционных перемешивающих устройств, время, затрачиваемое на процесс гидролиза, составляет около двух часов. Обычно в качестве такого устройства применяется аппарат с пропеллерной мешалкой [3]. В отличие от упомянутых устройств РПА обеспечивают достаточно быстрое и каче- ственное перемешивание компонентов. Этому способствуют высокочастотные пульсации ISSN 1025-6415 Доповiдi Нацiональної академiї наук України, 2007, №7 79 Рис. 1. Расчетная область роторно-пульсационного аппарата скорости, ускорения и давления в потоке жидкости, протекающей через рабочие органы аппарата. Рабочими органами РПА являются коаксиальные цилиндрические тела, разде- ленные зазорами и имеющие радиальные прорези. Часть цилиндрических тел находится в покое, а другая часть вращается с высокой угловой скоростью. Неподвижные и вращаю- щиеся цилиндрические тела (статоры и роторы) располагаются последовательно. Движение жидкости через рабочие органы обеспечивается как за счет перепада давления между вхо- дным и выходным сечениями рабочего пространства аппарата, так и за счет центробежной силы, вызванной вращением роторов (рис. 1). Для выяснения вопроса о количественных характеристиках работы РПА как перемеши- вающего устройства рассмотрим модельную задачу массопереноса одного из компонентов смеси в потоке другого компонента при их совместном движении через рабочее пространст- во данного аппарата. Рассматривается случай аппарата, рабочая зона которого включает два статора и один ротор, расположенный между статорами. Течение вязких жидкостей предполагается двумерным, ламинарным и рассматривается в горизонтальном сечении ап- парата, перпендикулярном общей оси цилиндрических рабочих элементов. Задача о течении вязких жидкостей, сопровождаемом массообменном, описывается сис- темой дифференциальных уравнений, имеющей в полярных координатах вид: уравнение неразрывности ∂ρ ∂τ + 1 r ∂(ρrv) ∂r + ∂(ρω) ∂θ = 0; (1) уравнения сохранения импульса ∂(ρv) ∂τ + 1 r ∂(ρrv2) ∂r + ∂(ρrvω) ∂θ − ρω2r = − ∂p ∂r + 2 r ∂ ∂r ( rµ ∂v ∂r ) + 2µ r ( ∂ω ∂θ + v r ) + + ∂ ∂θ [ µ ( 1 r2 ∂v ∂θ + ∂ω ∂r )] ; (2) ∂(ρω) ∂τ + 1 r2 ∂(ρr2ωv) ∂r + ∂(ρω2) ∂θ + ρvω r = − 1 r2 ∂p ∂θ + 2 r2 ∂ ∂θ [ µ ( ∂ω ∂θ + v r )] + + 1 r3 ∂ ∂r [ rµ ( r2 ∂ω ∂r + ∂v ∂θ )] ; (3) 80 ISSN 1025-6415 Reports of the National Academy of Sciences of Ukraine, 2007, №7 уравнение массопереноса ∂(ρC1) ∂τ + 1 r ∂(rvρC1) ∂r + ∂(ωρC1) ∂θ = 1 r ∂ ∂r ( rρD ∂C1 ∂r ) + 1 r2 ∂ ∂θ ( ρD ∂C1 ∂θ ) , (4) где τ — время; r — радиус; θ — полярный угол; v — радиальная скорость; ω — угловая скорость; p — давление; µ — динамический коэффициент вязкости; ρ — плотность смеси; C1 — локальная концентрация одного из компонентов раствора в смеси; D — коэффициент взаимной диффузии компонентов. Величина C1 задается, как отношение массы одного из компонентов m1 к массе смеси m1+m2, сосредоточенной в элементарном объеме. Очевидно, что концентрация второго компонента C2 = 1−C1, а плотность смеси может быть получена из выражения ρ = ρ1ρ2 ρ2C1 + ρ1(1 − C1) , где ρ1 и ρ2 — плотности отдельно взятых компонентов раствора. Таким образом, ρ = ρ(C1). От локальных значений концентрации C1 зависят также коэффициенты µ = µ(C1) и D = = D(C1). Их зависимости от температуры в данной задаче не учитываются. Расчетная область включает в себя участки рабочего пространства, занятые статора- ми 1, 3, ротором 2, зазорами между ними δ1 = δ2 = δ, а также небольшие участки пространс- тва, прилегающие к внутренней поверхности внутреннего статора и к внешней поверхности внешнего статора (см. рис. 1). Поскольку рабочие элементы представляют собой тонкие ци- линдрические тела с периодически повторяющимися прорезями, есть смысл рассматривать только один периодически повторяющийся сегмент с углом раскрытия ∆θ, включающий одну прорезь и две половины прилегающих к ней стенок, отделяющих данную прорезь от соседних. Граничными условиями для динамической задачи (1)–(3) будут: перепад давления ∆p = = pвх−pвых между входным и выходным сечениями расчетной области, которые характери- зуются минимальной (r0) и максимальной (r1) радиальными координатами (в этих сечени- ях тангенциальная составляющая скорости принимается равной нулю); угловая скорость ω0 вращения ротора; условия периодичности по углу ∆θ для всех неизвестных величин, вхо- дящих в исходную систему уравнений. Расчетная схема и результаты численного решения системы уравнений (1)–(3) представ- лены в [4]. Для решения уравнения массопереноса (4) сформулируем граничные условия для концентрации первого компонента C1. Будем при этом предполагать, что в получае- мой смеси массовая доля второго компонента значительно преобладает над массовой долей первого. Исходя из этого, примем в качестве условия на входе в рабочую зону следующие выражения: для r = r0 при 0 < θ < θ∗: C1 = 0 (C2 = 1); при θ∗ < θ < θ∗ + δθ∗: C1 = 1 (C2 = 0); при θ∗ + δθ∗ < θ < ∆θ: C1 = 0 (C2 = 1), (5) где θ∗ — в определенной степени произвольное значение угловой координаты (0 < θ∗ < ∆θ), относящееся к расчетной области, а δθ∗ ≪ ∆θ — участок входного сечения в рабочую зону, на котором, как предполагается, в основной поток преобладающего второго компонента вводится первый компонент раствора. ISSN 1025-6415 Доповiдi Нацiональної академiї наук України, 2007, №7 81 Как видно из условия (5), постановка данной задачи предполагает, что концентрация первого компонента скачкообразно изменяется от 0 до 1 на достаточно малом участке вход- ного сечения рабочей зоны. Строго говоря, условие (5) лишь приближенно моделирует ре- альную картину процесса, протекающего в РПА. В действительности, загрузка различных компонентов в аппарат производится через специальную загрузочную емкость или через штуцеры, расположенные на его корпусе. Естественно, что двухкомпонентный поток на подходе к рабочей зоне аппарата оказывается до некоторой степени уже перемешанным. Условие же в форме (5) будет в этом смысле более “жестким”, чем реальное. Его использо- вание позволит лучше оценить вклад собственно активных элементов аппарата в процесс перемешивания, поскольку оно является наименее “благоприятным” для начала процесса перемешивания в рабочей зоне. На выходе из рабочей зоны условия формулируются в виде для r = r1 при 0 < θ < ∆θ: ∂C1 ∂r = 0. (6) В качестве примера рассматривается аппарат со следующими геометрическими разме- рами [5]: внутренний статор: r1вн = 55,0 мм; r1нар = 58,9 мм; ротор: r2вн = 59,1 мм; r2нар = 62,9 мм; внешний статор: r3вн = 63,1 мм; r3нар = 67,0 мм. При этом r0 = 54,0 мм, а r1 = 68,0 мм. Каждый цилиндрический элемент имеет 60 прорезей. Угловой размер прорези соответствует угловому размеру стенки, разделяющей соседние прорези. Скорость враще- ния ротора — 48 оборотов в секунду. Перепад давления ∆p = 50 кПа. Решение модельной задачи о перемешивании двух высоковязких жидкостей выполня- ется при условии, что их плотности одинаковые: ρ1 = ρ2 = 1135 кг/м3, а коэффициенты переноса — постоянные: µ = 5 Па · с; D = 0,63 · 10−6 м2/с. Численное решение задачи (1)–(4) с указанными граничными условиями выполнено на разностной сетке с количеством шагов n = 48 по угловой координате, и m = 155 — по радиальной. Радиальные шаги разностной сетки сгущались около поверхностей рабочих элементов, а шаги по угловой координате были постоянными: dθ = ∆θ/n. Шаг по времени dτ выбирался так, чтобы за интервал dτ ротор поворачивался на угол dθ, т. е. dτ = dθ/ω0. Аппараты роторного типа относятся к устройствам периодического действия. Геомет- рическая периодичность конструкций рабочих элементов аппарата, о которой говорилось выше, предопределяет периодичность изменения во времени динамических характеристик течения жидкости через рабочую зону. Другими словами, при установившемся режиме ра- боты аппарата картина течения жидкости повторяется после очередного поворота ротора на периодический угол ∆θ. Указанным периодом времени будет интервал ∆τ = ∆θ/ω0. Учитывая вышесказанное, будем в дальнейшем рассматривать результаты расчета хара- ктеристик работы аппарата на указанном интервале времени 0< τ < ∆τ . Результаты расчетов полей скорости и концентрации, выполненных для случая θ∗ = = ∆θ/2; δθ∗ = dθ, представлены на рис. 2. Рассмотрены два характерных момента вре- мени: момент совпадения каналов ротора и статоров (τ = 0; ∆τ) и момент их полного взаимного перекрытия (τ = ∆τ/2). Векторы скорости на рис. 2 представлены так, что на участке прорези ротора они рассматриваются в подвижной системе координат, связанной с вращающимся ротором, а на остальной части пространства, включающей прорези двух статоров и зазоры между рабочими элементами, — в неподвижной системе. Как видно из рис. 2, интенсивное перемешивание компонентов раствора начинается во второй половине участка прорези внутреннего статора, в которой располагается циркуляционная зона по- 82 ISSN 1025-6415 Reports of the National Academy of Sciences of Ukraine, 2007, №7 Рис. 2. Поля скорости и концентрации в рабочем пространстве РПА при двух положениях ротора: а — полное открытие прорезей (τ = 0); б — полное взаимное перекрытие (τ = ∆τ/2) тока, способствующая процессу перемешивания. Далее по тракту рабочего пространства аппарата концентрация C1 изменяется слабо и находится в пределах C1 = 0,0541. . .0,0545. Изменение радиальной составляющей скорости v во входном сечении расчетной зоны (r = r0) на интервале 0 < θ < ∆θ для указанных выше характерных моментов време- ни представлены на рис. 3, откуда видно, что форма профиля скорости претерпевает за указанный период времени существенные изменения. Это связано с изменением площади проходного сечения для потока жидкости вследствие периодического взаимного перекрытия прорезей в роторе и в статорах. Кроме того, профили радиальной скорости v(τ, r0, θ) нерав- номерны по углу θ. При θ = ∆θ/2 указанная функция имеет максимум, а в окрестностях ISSN 1025-6415 Доповiдi Нацiональної академiї наук України, 2007, №7 83 Рис. 3. Изменение радиальной скорости v вдоль входного сечения (r = r0) рабочего пространства РПА при двух положениях ротора: 1 — полное открытие прорезей; 2 — полное взаимное перекрытие точек θ = 0 и θ = ∆θ она отрицательна. В данном случае область входного сечения, на котором радиальная скорость v максимальна, совпадает с участком подачи в раствор пер- вого (меньшего по количеству) компонента. Этим объясняется тот факт, что концентрация данного компонента в растворе C1(τ, r1, θ) на выходе из рабочей зоны аппарата (r = r1) оказывается больше величины δθ∗/∆θ = 1/n = 0,0208, характеризующей ту часть площади входного сечения рабочей зоны, на которой производится ввод этого компонента в раствор. Для оценки качественных показателей работы роторного аппарата как перемешиваю- щего устройства определим вначале величину среднемассовой концентрации компонента в смеси по формуле C1(τ, r) = ∆θ ∫ 0 v(τ, r, θ)C1(τ, r, θ)rρdθ ∆θ ∫ 0 v(τ, r, θ)rρdθ , (7) выражающей отношение массового расхода первого компонента к массовому расходу всей смеси. Поскольку по условию рассматриваемой задачи плотность смеси ρ — величина по- стоянная, значение массового расхода смеси ∆G = ∆θ ∫ 0 v(τ, r, θ)rρdθ, стоящее в знаменателе выражения (7), оказывается не зависящим от радиуса r, что следует из закона сохранения массы (1). Вместе с тем, как видно из рис. 3, в течение периода ∆τ значение расхода ∆G изменяется во времени. Числитель же выражения (7) изменяется как вдоль радиуса, так и во времени. Связано это с тем, что масса первого компонента вводится в рабочую зону изменяющимися во времени порциями. Поскольку для продвижения указанной массы по всему тракту рабочей зоны требуется определенное время, расходы первого компонента во входном и выходном сечениях рабочей зоны оказываются в отличие от суммарного массо- вого расхода ∆G во времени не совпадающими. Поэтому величина C1(τ, r) зависит как от времени, так и от радиальной координаты. Характер ее изменения за период ∆τ во вход- ном и выходном сечениях рабочей зоны представлены на рис. 4, откуда видно, что средняя 84 ISSN 1025-6415 Reports of the National Academy of Sciences of Ukraine, 2007, №7 Рис. 4. Изменение во времени средней концентрации C1 в рабочем пространстве РПА: 1 — во входном сечении; 2 — в выходном сечении; 3 — C1 по сечению концентрация C1(τ, r) существенно изменяется во времени во входном сечении (кривая 1 ) и слабо изменяется во времени в выходном сечении (кривая 2 ). Оценивая работу РПА как перемешивающего устройства, целесообразно исключить вли- яние времени τ на его динамические и масоообменные характеристики. Для этого уместно кроме величины C1(τ, r) рассмотреть также значение концентрации, осредненной как по сечению, так и за период времени ∆τ C1 = ∆τ ∫ 0 ∆θ ∫ 0 v(τ, r, θ)C1(τ, r, θ)rρdθdτ ∆τ ∫ 0 ∆θ ∫ 0 v(τ, r, θ)rρdθdτ . Следует отметить, что величина C1 оказывается не зависящей от радиуса r. Связано это с тем, что при установившемся периодическом режиме работы аппарата вся масса пер- вого компонента, вошедшая с радиальным потоком жидкости в рабочую зону, из этой зоны за период времени ∆τ удаляется (если бы этого не происходило, процесс не был бы пери- одическим). Поэтому осредненные по сечению и во времени расходы первого компонента будут одинаковыми как во входном, так и в выходном сечениях рабочей зоны. Как показали результаты вычислений, в рассматриваемом случае C1 = 0,054. Степени неоднородности σC1 (r) смеси двух жидкостей во входом и выходном сечениях рабочей зоны РПА можно оценить по осредненному за период ∆τ среднеквадратичному отклонению локальных концентраций C1(τ, r, θ) в рассматриваемых сечениях от средней концентрации C1. Эта величина будет зависеть как от качества работы перемешивающего устройства, так и от соотношения масс взаимнорастворяемых компонентов. Поэтому, оце- нивая работоспособность РПА как смесителя, целесообразно данную величину разделить на среднюю концентрацию C1. В итоге получим выражение σC1 (r) = ∆τ ∫ 0 { 1 ∆θ ∆θ ∫ 0 [C1(τ, r, θ) − C1] 2dθ }0,5 dτ C1∆τ . ISSN 1025-6415 Доповiдi Нацiональної академiї наук України, 2007, №7 85 Рассчитанная по приведенной формуле степень неоднородности смеси на входе в рабочую зону (r = r0) составляет σC1 (r0) ≈ 2,71, а на выходе из нее — σC1 (r1) ≈ 0,032. Как следует из результатов решения динамической задачи (1)–(3), в рабочей зоне аппа- рата среднемассовая скорость движения жидкости в радиальном направлении составляет приблизительно 0,128 м/с. При этом длина рабочей зоны вдоль радиуса равна ∆r = r1−r0 = = 0,014 м. Следовательно, обрабатываемая в аппарате смесь находится под воздействием рабочих органов примерно 0,11 с. За это время степень неоднородности раствора в рабо- чей зоне уменьшается в 84,7 раза. Данный показатель работы роторного аппарата можно считать достаточно высоким. Следует также отметить, что указанное снижение степени не- однородности смеси достигается за единичный проход жидкости через рабочую зону. Если по техническим условиям требуемая степень неоднородности смеси окажется не достигну- той за один проход, смесь можно неоднократно пропустить через рабочие органы аппарата до достижения им необходимых количественных показателей. 1. Долинский А.А., Павленко А.М., Басок Б.И. Теплофизические процессы в эмульсиях: получение, использование, утилизация. – Киев: Наук. думка, 2005. – 264 с. 2. Ободович А.Н., Хибина М.А., Боряк Л.А. и др. Интенсификация производства глюкозо-фруктозных сиропов за счет механохимической деструкции // Промышл. теплотехника. – 2006. – 28, № 3. – С. 44–49. 3. Мальцев П.М. Технология бродильных производств. – Москва: Пищевая промышленность, 1980. – 560 с. 4. Басок Б.И., Давыденко Б.В., Кравченко Ю.С., Пироженко И.А. Исследование микроструктуры потока жидкости в роторно-пульсационном аппарате // Доп. НАН України. – 2003. – № 11. – С. 71–76. 5. Басок Б.И., Гартвиг А.П., Коба А.Р., Горячев О.А. Оборудование для получения и обработки высо- ковязких дисперсных сред // Промышл. теплотехника. – 1996. – 18, № 1. – С. 50–56. Поступило в редакцию 05.03.2007Институт технической теплофизики НАН Украины, Киев 86 ISSN 1025-6415 Reports of the National Academy of Sciences of Ukraine, 2007, №7
id nasplib_isofts_kiev_ua-123456789-1724
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
issn 1025-6415
language Ukrainian
last_indexed 2025-12-07T19:03:55Z
publishDate 2007
publisher Видавничий дім "Академперіодика" НАН України
record_format dspace
spelling Басок, Б.И.
Давыденко, Б.В.
Ободович, А.Н.
Пироженко, И.А.
2008-09-02T16:56:24Z
2008-09-02T16:56:24Z
2007
Численное моделирование процесса перемешивания взаиморастворимых жидкостей в роторно-пульсационном аппарате / Б.И. Басок, Б.В. Давыденко, А.Н. Ободович, И.А. Пироженко // Доп. НАН України. — 2007. — N 7. — С. 79–86. — Бібліогр.: 5 назв. — рос.
1025-6415
https://nasplib.isofts.kiev.ua/handle/123456789/1724
66.021.2:663.913.7
Rotor-pulse apparatus are successfully used for the creation of homogeneous mixtures from high-viscosity components. The fields of concentration of components are defined as a result of the numerical simulation of the liquid flow and the process of mass transfer in the working space of an apparatus. The time dependence of the average concentrations in the input and output sections of the working space is investigated. The estimation of the efficiency of a rotor-pulse apparatus as a mixing device is made on the basis of the proposed criterion.
uk
Видавничий дім "Академперіодика" НАН України
Теплофізика
Численное моделирование процесса перемешивания взаиморастворимых жидкостей в роторно-пульсационном аппарате
Article
published earlier
spellingShingle Численное моделирование процесса перемешивания взаиморастворимых жидкостей в роторно-пульсационном аппарате
Басок, Б.И.
Давыденко, Б.В.
Ободович, А.Н.
Пироженко, И.А.
Теплофізика
title Численное моделирование процесса перемешивания взаиморастворимых жидкостей в роторно-пульсационном аппарате
title_full Численное моделирование процесса перемешивания взаиморастворимых жидкостей в роторно-пульсационном аппарате
title_fullStr Численное моделирование процесса перемешивания взаиморастворимых жидкостей в роторно-пульсационном аппарате
title_full_unstemmed Численное моделирование процесса перемешивания взаиморастворимых жидкостей в роторно-пульсационном аппарате
title_short Численное моделирование процесса перемешивания взаиморастворимых жидкостей в роторно-пульсационном аппарате
title_sort численное моделирование процесса перемешивания взаиморастворимых жидкостей в роторно-пульсационном аппарате
topic Теплофізика
topic_facet Теплофізика
url https://nasplib.isofts.kiev.ua/handle/123456789/1724
work_keys_str_mv AT basokbi čislennoemodelirovanieprocessaperemešivaniâvzaimorastvorimyhžidkosteivrotornopulʹsacionnomapparate
AT davydenkobv čislennoemodelirovanieprocessaperemešivaniâvzaimorastvorimyhžidkosteivrotornopulʹsacionnomapparate
AT obodovičan čislennoemodelirovanieprocessaperemešivaniâvzaimorastvorimyhžidkosteivrotornopulʹsacionnomapparate
AT piroženkoia čislennoemodelirovanieprocessaperemešivaniâvzaimorastvorimyhžidkosteivrotornopulʹsacionnomapparate