Аннотация
Введение 3
1 Математическая постановка и метод решения 6
1.1 Математическая постановка 6
1.2 Метод контрольного объёма для структурированной сетки 8
1.2.1 Метод, основанный на теореме Грина 12
1.2.2 Метод на основе МНК 15
1.3 Структурированная сетка 23
1.3 Алгоритм SIMPLE 27
2 Результаты исследования 31
2.1 Апроксимационная сходимость 31
2.2 Случай Ньютоновской жидкости для якорных мешалок 32
2.3 Случай степенной жидкости в лопастном смесителе 34
2.4 Случай бингамовской жидкости в лопастном смесителе 37
Заключение 41
Список использованной литературы 42
Процессы перемешивания жидкостей в специализированных аппаратах, с целью получения высокооднородных по составу смесей распространены во многих отраслях промышленности , например химической, фармацевтической, металлургической и т.д. Широкое применение нашел механический метод перемешивания, в котором смешение происходи за счет вращение лопастей и мешалок. Отметим, что данный метод не единственный, однако один из самых дешевых в производстве. Решающей значение организации технического процесса в производстве, играет детальное понимание структуры течения и возможности влиять на нее. Экспериментальное исследование течений в смесителях трудоемкий и зачастую дорогостоящий процесс, так как зачастую исследуются течение неньютоновской или многофазной жидкости свойствами. Современный уровень экспериментальных исследований описан . Аналитическое решение для такого типа задач зачастую либо получить нельзя, либо с большим числом допущений, некоторые такие результаты представлены . Поэтому все большое распространение находит применение вычислительное гидродинамики, для анализа течения.
На сегодняшний день проведено много работ по математическому моделированию течений в различных смесителях. Результаты численных исследований течений ньютоновской и степенной жидкости в смесительных аппаратах представлены . Анализ процессов диспергирования твердых частиц в вязкой жидкости даны в. Существует большое количество работ, посвященных моделированию процессов смешения в турбулентном режиме Однако число исследований течений высоковязких жидкостей, к которым относятся многие полимерные материалы, ограничено. Так же сравнительно мало работ посвященных влиянию параметров смесителя на структуру потока.
Вычислительная гидродинамика позволяет решать самые разные, по сложности задачи, за счет большого арсенала разработанных методов. Так существует множество методов дискретизации уравнений, например, таких как метод контрольных объёмов, метод конечных элементов. Выбор одних из них диктуется особенностями решаемой задачи.
Наибольшее распространение в нынешнее время получил метод контрольного объёма. Так как уравнений получающиеся с помощью данного метода обладают важными свойствами, такими как консервативность, монотонность и т.д. Так же он легко распространяется на разные типы сеток, что позволяет решать задачи с разнообразной геометрией. Надо отметить, что теоретическая база МКО не так сильно развита по сравнению, например с методом конечных разностей.
Существенным в вычислительной гидродинамике является способы дискретизации области. Так как от способа дискретизации зависит не только точность отражаемой геометрии, но и точность дискретизированных уравнений. Так структурированные сетки позволяют получать более точные решения, для простых областей, однако их распространение на сложную геометрию связано с существенными математическими трудностями. Между тем не структурированные сетки, например треугольные, позволяют точнее представить сложную геометрию. Однако такие недостатки как: большая сложность получения дискретизированных уравнений, значительно более трудоемкий процесс построения самой сетки и большая ресурсоемкость. Ограничивают их применение.
...
В рамках данной работы сформулирована математическая постановка
задачи и разработан численный алгоритм решения, основанный на методе
контрольного объёма и корректирующей процедуры SIMPLE. Дискретизация
дифференциальных уравнений для лопастного смесителя производилась на
структурированной, разнесенной, с не постоянным шагом по радиусу, сетке.
Для якорного смесителя дискретизация проводилась на неструктурированной
треугольной сетке. Продемонстрированно наличие аппроксимационной
сходимости. Проведен анализ влияния параметров якорного смесителя на
течение. Исследовано влияния реологии на кинематику течения в лопастном
смесителе неньютоновской жидкости для разных безразмерных параметров
таких, как Re, Bn, n. Показан качественно процесс смешения с помощью
введения в поток специальных маркеров
Oldshue J. Y. Fluid mixing technology / J. Y. Oldshue. - New York:
McGraw-Hill Publications Co. 1983 - 574 pages.
2. Shiue A. et al. Mixing performance of a non-Newtonian fluid in a coaxial
agitated impeller reactor // J Taiwan Inst Chem Eng. Elsevier, 2023. Vol. 143. P.
104715.
3. Kazemzadeh A. et al. Investigation of hydrodynamic performances of
coaxial mixers in agitation of yield-pseudoplasitc fluids: Single and double central
impellers in combination with the anchor // Chemical Engineering Journal. Elsevier,
2016. Vol. 294. P. 417–430.
4. Климов Д.М. Вязкопластические течения: динамический хаос,
устойчивость, перемешивание / Д. М. Климов, Д. В. Георгиевский, А. Г.
Петров. – М.: Наука, 2005. – 394 с.
5. AMEUR H. Effect of some parameters on the performance of anchor
impellers for stirring shear-thinning fluids in a cylindrical vessel // J Hydrodynam B.
No longer published by Elsevier, 2016. Vol. 28, № 4. P. 669–675.
6. Patel D., Ein-Mozaffari F., Mehrvar M. Effect of rheological parameters on
non-ideal flows in the continuous-flow mixing of biopolymer solutions // Chemical
Engineering Research and Design. Elsevier, 2015. Vol. 100. P. 126–134.
7. Xu Z. et al. Numerical Study on the Synergistic Mechanism of Coaxial
Mixers in Dense Solid-Liquid Mixing Systems // Ind Eng Chem Res. American
Chemical Society, 2023. Vol. 62, № 29. P. 11744–11755.
8. Blais B. et al. Development of an unresolved CFD–DEM model for the
flow of viscous suspensions and its application to solid–liquid mixing // J Comput
Phys. Academic Press, 2016. Vol. 318. P. 201–221.
9. Charalambidou A.D., Micheletti M., Ducci A. Study of trailing vortices and
impeller jet instabilities of a flat blade impeller in small-scale reactors // AIChE
Journal. John Wiley and Sons Inc, 2023. Vol. 69, № 2. P. e17842.
10. Komoda Y., Date T. Enhancement of laminar mixing by an anchor
impeller with rotationally reciprocating motion Articles You May Be Interested In.
2022.
11. Фирсов Д.К. Метод контрольного объёма на неструктурированной
сетке в вычислительной механике : учебное пособие / Д. К. Фирсов. – Томск :
ТГУ, 2007. – 72 с.
12. Wang J.W., Liu R.X. A comparative study of finite volume methods on
unstructured meshes for simulation of 2D shallow water wave problems // Math
Comput Simul. North-Holland, 2000. Vol. 53, № 3. P. 171–184.
13. Hubbard M.E. Multidimensional Slope Limiters for MUSCL-Type Finite
Volume Schemes on Unstructured Grids // J Comput Phys. Academic Press, 1999.
Vol. 155, № 1. P. 54–74.
14. Ollivier-Gooch C.F. Quasi-ENO Schemes for Unstructured Meshes Based
on Unlimited Data-Dependent Least-Squares Reconstruction // J Comput Phys.
Academic Press, 1997. Vol. 133, № 1. P. 6–17.
15. Горбаченко В. И. Вычислительная линейная алгебра с примерами на
MATLAB / В. И. Горбаченко. – Санкт-Петербург : БХВ-Петербург, 2011. –
320 с. // БХВ-Петербург.