Обобщенная неустойчивость Рэлея-Тейлора для суспензий, седиментационная конвекция



Extended Rayleigh-Taylor instability for suspension, stdimentatation convection

The process of suspension sedimentation with heavy particles, including the interaction between each other and with the basic liquid, is studied. The case of small Prandtl number is considered only, without inertial terms, for slow flow in high viscous liquid. Similar to Rayleigh-Benar instability theory for convective heat transfer or diffusive mass transfer in viscous liquid it is possible to create the instability analysis of the convective mass transfer in viscous liquid with settling particles, including the transition of Stokes sinking of heavy layer to convection of whole mixture, when the heavy layer deforms like jets. In limit case of zero settling velocity (i.e. frozen particles) this area of liquid can be considered as the heavy liquid. In this case sedimentation convection turns into flow of Rayleigh-Taylor theory. The analytical analysis of instability can be created for small concentration and mean description of the particle distribution. The critical modes and critical number of sedimentation convection onset are calculated. They give simple criterion. In case of slow flow (independently of viscosity and density jump) the layer of heavy particles sinks by Stokes if radius of particles is grater than ~0.15C1/2, where C is concentration. For smaller size of particles the layer deforms and Rayleigh-Taylor flows appear. The results of analytical analysis are tested by numerical solution of the system of Stokes equation and equation of concentration transfer in the same approximation of mean particle distribution both as well as for linear and finite-amplitude flows. Later the numerical experiments for sinking layer of solid interactive cylinders were done by direct solution of Stokes equation for discrete particles, without averaging-out. The particles are described by active markers with higher density and anomaly high viscosity.


Том 13, выпуск 2, 2012 год



Изучается процесс осаждения суспензий с тяжелыми частицами, взаимодействующими между собой и с основной вязкой жидкостью. Рассматривается случай малого числа Прандтля, без инерциальных членов, что соответствует медленным течениям в высоковязкой жидкости. Аналогично теории неустойчивости Рэлея-Бенара для конвективного переноса тепла и теории диффузионной конвекции можно провести анализ неустойчивости переноса массы в вязкой жидкости с оседающими частицами, включая переход от кондуктивного переноса массы в виде стокосового осаждения слоя как целого к седиментационной конвекции всей смеси, когда слой тяжелых частиц деформируется в струи. В пределе при стремлении к нулю скорости осаждения частиц (т.е. вмороженных в жидкость) область жидкости, содержащую тяжелые частицы, можно рассматривать как более плотную жидкость. В этом пределе седиментационная конвекция переходит в классические течения неустойчивости Рэлея-Тейлора. В приближении малой концентрации и усредненного описания распределения частиц удается провести аналитический анализ неустойчивости. Найдены критические моды течений и критические числа порога (начала) седиментационной конвекции, из которых выводится простой критерий. В случае медленных течений независимо от вязкости и скачка плотности слой тяжелых частиц с концентрацией C при радиусах более ~0.15C1/2 от размера области падает как целое по Стоксу, при меньших размерах частиц слой деформируется, и возникают струйные течения, в пределе переходящие в течения рэлей-тейлоровской неустойчивости. Результаты линейного аналитического анализа проверены численными решениями системы уравнений Стокса и уравнения переноса концентрации в том же приближении с усреднением распределением частиц. При этом были рассчитаны как малые, так и конечно-амплитудные течения. Далее были проведены численные эксперименты для оседания взаимодействующих частиц (твердых цилиндров) в вязкой жидкости путем численного решения непосредственно уравнения Стокса для конечных концентраций и без усреднения, с заданиями дискретных частиц в виде маркируемых перемещающихся областей повышенной плотности и аномально высокой вязкости.


Том 13, выпуск 2, 2012 год



1. Трубицын В.П., Харыбин Е.В. Конвективная неустойчивость режима седиментации в мантии // Физика Земли, 1987, № 7, С.21-30.
2. Трубицын В.П., Харыбин Е.В. Термоседиментационная конвективная неустойчивость двухкомпонентной вязкой жидкости // Физика Земли, 1991, № 2, С.3-17.
3. Chandraskhar S. Hydrodynamic and hydromagnetic stability// Oxford.Claredon Press, 1961, P.652.
4. Joseph D.D. Viscous potential flow// J. Fluid. Mech. 2003, V.479, P.191-197.
5. Трубицын В.П., Харыбин Е.В. Гидродинамическая модель дифференциации вещества в Земле // Физика Земли, 1988, №4, С.83-86.
6. Трубицын В.П., Харыбин Е.В. Конвекция в магматических камерах, вызванная инверсией распределения по глубине осаждающихся кристаллов // Физика Земли, 1997, № 5, С.47-52.
7. Трубицын В.П., Харыбин Е.В. Конвекция в вязкой жидкости с оседающими частицами// Фи-зика Земли, 2005, №12, С.3-11
8. Hill W.D., Rothhfus R.R., Kun Li. Boundary enhanced sedimentation due to settling convection// Intern. J. Multiphase Flow, 1977, V.3, N6, P. 561-583
9. McCaffrey S.J., Elliott L., Ingham D.B., Enhanced Sedimentation in Inclined Fracture Channels// Topics in Engineering. 1997, V. 32, P. 280-291
10. Невский Ю.А., Осипцов А.Н. Медленная гравитационная конвекция дисперсных систем в областях с наклонными границами//Изв. РАН, Механика жидкости и газа, 2011, №2, с. 65-81
11. Невский Ю.А. Моделирование гравитационной конвекции в дисперсных системах // Автореф. дис. канд. физ.-мат. наук. Москов. гос. ун-т им. М. В. Ломоносова. Москва, 2010, с. 1-22
12. Rudman M., Two-phase natural convection for crystal settling in magma chambers// Phys. Earth Planet. Inter., 1992, V.72, P.153-172.
13. Ландау,Л.Д., Лифшиц Е.М., Гидродинамика. <<Наука>> Москва 1986.
14. Слезкин Н.А. Динамика вязкой несжимаемой вязкости, ГИТТЛ, Москва, 1955, С. 514.
15. Richardson, J.F., Zaki, W.N. Sedimentation and fluidization// I. Trans. Inst. Chem. Eng. , 1954, V.32, P.35–53.
16. Moresi L.N., Solomatov V. S. Numerical investigation of 2D convection with extremely large vis-cosity variations// Phys. Fluids, 1995, V.7, P. 2154-2162.
17. Moresi L., Zhong S. J., Gurnis M. The accuracy of finite element solutions of Stokes' flow with strongly varying viscosity // Phys. Earth Planet. Inter., 1996, V. 97, P. 83-94.
18. Трубицын В.П., Евсеев А.Н., Баранов А.А., Трубицын А.П. Структура конвекции при
различной ширине зон фазовых переходов // Физика Земли. 2008. №8 . С. 3-14.