Все выпуски

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

В работе рассмотрена интегрируемая гамильтонова система на алгебре Ли $so(4)$ с дополнительным интегралом четвертой степени - интегрируемый случай Адлера-ван Мёрбеке. Рассмотрены классические работы, посвященные, с одной стороны, динамике твердого тела, содержащего полости, полностью заполненные идеальной жидкостью, совершающей однородное вихревое движение, а с другой стороны, изучению геодезических потоков левоинвариантных метрик на группах Ли. Приведены уравнения движения, функция Гамильтона, скобки Ли-Пуассона, функции Казимира и фазовое пространство рассматриваемого случая. В предыдущих работах начато исследование фазовой топологии интегрируемого случая Адлера-ван Мёрбеке: приводятся в явном виде спектральная кривая, дискриминантное множество, бифуркационная диаграмма отображения момента, предъявлены характеристические показатели для определения типа критических точек ранга 0 и 1 отображения момента. В данной работе излагается алгоритм построения торов Лиувилля. Рассмотрены примеры перестроек лиувиллиевых торов при пересечении бифуркационных кривых для перестроек одного тора в два и двух торов в два.

Представлены результаты математического моделирования течения степенной жидкости в смесительном аппарате лопастного типа. Численное решение задачи основано на методе контрольного объема и процедуре SIMPLE, реализованных на неравномерной ортогональной сетке. Исследована картина распределения кинематических характеристик в объеме и изучена структура потока с выделением характерных циркуляционных зон для псевдопластичной, ньютоновской и дилатантной реологий. Для оценки процесса смешения во времени решена задача о перераспределении ансамбля маркерных частиц, в начальный момент расположенных в выделенной области и движущихся со скоростью среды. Проведены параметрические исследования двух технологических критериев: интеграла от диссипативной функции, показывающего энергетические затраты на организацию течения; оригинальной характеристики, рассчитываемой по текущему распределению маркерных частиц, позволяющей оценить однородность перемешивания.

Представлены результаты численных исследований собственных колебаний усеченных прямых конических оболочек вращения, полностью заполненных идеальной сжимаемой жидкостью. Толщина оболочек непостоянна вдоль образующей и изменяется по различным законам. Поведение упругой конструкции и жидкой среды описывается в рамках классической теории оболочек, основанной на гипотезах Кирхгофа–Лява, и уравнений Эйлера. Уравнения движения оболочки совместно с соответствующими геометрическими и физическими соотношениями сводятся к системе обыкновенных дифференциальных уравнений относительно новых неизвестных. Акустическое волновое уравнение, записанное относительно гидродинамического давления, преобразуется к системе дифференциальных уравнений с помощью метода обобщенных дифференциальных квадратур. Решение сформулированной краевой задачи осуществляется методом ортогональной прогонки Годунова и сводится к вычислению собственных частот колебаний. Для этой цели используется сочетание пошаговой процедуры с последующим уточнением найденных значений в полученном диапазоне методом Мюллера. Достоверность получаемых результатов подтверждена сравнением с известными численными решениями. Для оболочек с различными углами конусности и комбинациями граничных условий (свободное опирание, жесткое и консольное закрепления) исследованы зависимости низших частот колебаний, полученных при степенном (линейном и квадратичном, имеющих симметричную и несимметричную формы) и гармоническом (с положительной и отрицательной кривизной) изменении толщины. Оценено влияние граничных условий на возможность существования конфигураций (угол конусности, закон изменения толщины, отношение максимальной и минимальной толщины профиля), обеспечивавших повышение фундаментальной частоты по сравнению с оболочками постоянной толщины при ограничениях на вес конструкции.

Работа посвящена моделированию ползущего движения вязкоупругой жидкости со свободной поверхностью, реализующейся при входе полимерной жидкости в формующий канал и выходе из него.  Движение жидкости описывается  уравнениями сохранения массы, импульса и энергии, дополненное определяющим реологическим уравнением состояния среды Гиезекуса. На основе метода конечных элементов разработан устойчивый численный алгоритм решения задачи. Проведены численные исследования по определению формы выходной струи для различных режимов течения и формы насадки. Исследована картина распределения скоростей жидкости, давления, напряжений и температуры при увеличении степени нагрева стенки формующего канала. Получены численные результаты зависимости эффекта разбухания полимерной жидкости от параметров реологической модели и температурных факторов.

В работе представлены результаты расчетного исследования локальной структуры восходящего газожидкостного потока в вертикальной трубе. Математическая модель основана на использовании двухжидкостного эйлерова подхода с учетом обратного влияния пузырьков на осредненные характеристики и турбулентность несущей фазы. Турбулентная кинетическая энергия жидкости рассчитывается с применением двухпараметрической изотропной модели турбулентности $k - \varepsilon$, модифицированной для двухфазных сред. Для описания динамики распределения пузырьков по размерам используются уравнения сохранения количества частиц для отдельных групп пузырьков с различными диаметрами для каждой фракции с учетом процессов дробления и коалесценции. Численно исследовано влияние изменения степени дисперсности газовой фазы, объемного расходного газосодержания, скорости дисперсной фазы на локальную структуру и поверхностное трение в двухфазном потоке. Сравнение результатов моделирования с экспериментальными данными показало, что разработанный подход позволяет адекватно описывать турбулентные газожидкостные течения в широком диапазоне изменения газосодержания и начальных размеров пузырьков.

В работе рассматриваются результаты решения задачи стационарного течения вязкой несжимаемой жидкости в плоском канале с обратным уступом и прогреваемой нижней стенкой в широком диапазоне числа Рейнольдса $100\leqslant \text{Re}\leqslant 1000$ и параметра расширения потока $1.11 \leqslant ER \leqslant 10$. Исследование выполнено путем численного интегрирования системы двумерных уравнений Навье-Стокса в переменных «скорость-давление» на равномерных сетках с шагом 1/300. Достоверность полученных результатов подтверждается их сравнением с литературными данными. Приводятся подробные картины течения и перегрева жидкости в зависимости от двух основных параметров задачи: $\text{Re}$ и $ER$. Показывается, что с одновременным ростом параметров $\text{Re}$ и $ER$ существенно усложняется структура течения - увеличиваются количество вихрей и их размеры вплоть до образования вихря за уступом с двумя центрами вращения. Также показывается, что характерная высота зоны прогрева течения слабо зависит от $\text{Re}$ и $ER$ и в конечном счете ближе к выходу из канала составляет приблизительно половину его высоты. Для всех центров вихрей определяются их основные характеристики: координаты, экстремумы функции тока, завихренности. Анализируется сложное немонотонное поведение профилей коэффициентов трения, сопротивления и теплоотдачи (числа Нуссельта) по длине канала. Показывается, что эти коэффициенты в одинаковой степени сильно зависят как от числа Рейнольдса, так и от параметра расширения канала, достигая своих максимальных значений при максимальных значениях $\text{Re}$ и $ER$.

В работе предложен общий топологический подход к исследованию устойчивости периодических решений интегрируемых динамических систем с двумя степенями свободы. Развиваемые методы проиллюстрированы на примерах нескольких интегрируемых задач, связанных с классическими уравнениями Эйлера—Пуассона, движением твердого тела в жидкости, а также динамикой газообразных расширяющихся эллипсоидов. Данные топологические методы позволяют также отыскивать невырожденные периодические решения интегрируемых систем, что является особенно актуальным в тех случаях, когда общее решение, например, при помощи разделения переменных неизвестно.