Стационарная и колебательная конвекция бидисперсной коллоидной суспензии

В рамках бидисперсной модели с помощью метода конечных разностей проведено численное моделирование и изучены нелинейная эволюция и свойства двумерных режимов конвекции коллоидной суспензии, возникающих в горизонтальном слое при подогреве снизу. Слой имеет твердые теплопроводные и непроницаемые для суспензии горизонтальные границы. На боковых границах расчетной области используются периодические граничные условия, позволяющие исследовать не только стационарную конвекцию и стоячие волны, но и бегущие волны. Рассматривается случай, когда в состоянии механического равновесия термодиффузионный транспорт и гравитационная стратификация наночастиц направлены противоположно, что является источником оседания тяжелой примеси и возникновения колебательной неустойчивости. При достижении числом Релея некоторого критического значения в слое образуется стоячая волна, которая оказывается неустойчивой. Ее разрушение, сопровождающееся перестройкой полей концентрации наночастиц, приводит к длительному переходному режиму бегущих волн, характеристики которого изучены и проанализированы. Построена бифуркационная диаграмма (зависимость максимальной функции тока от числа Релея), изображающая нелинейные режимы течения жидкости в зависимости от интенсивности нагрева. Показано, что при высокой интенсивности нагрева конвективное перемешивание размывает гравитационную седиментацию наночастиц и приводит к сглаживанию концентрационных неоднородностей. В результате нелинейной эволюции колебательных возмущений устанавливается режим стационарной конвекции, для которого свойственна зеркальная симметрия решений. Устойчивые режимы бегущих волн обнаружены в узкой подкритической области значений числа Релея. Получены поля функции тока, температуры и концентраций примеси малых и больших наночастиц.