Все выпуски

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

Рассматривается нелинейное эволюционное операторное уравнение второго рода $\varphi=\mathcal{F}\bigl[f[u]\varphi\bigr]$, $\varphi\in W[0;T]\subset L_q\bigl([0;T];X\bigr)$, в произвольном банаховом пространстве $X$, с эволюционными (вольтерровыми) операторами $\mathcal{F}\colon L_p\bigl([0;\tau];Y\bigr)\to W[0;T]$, $f[u]\colon W[0;T]\to L_p\bigl([0;T];Y\bigr)$ общего вида, $Y$ - произвольное банахово пространство, $u\in\mathcal{D}$ - управляющий параметр. Для указанного уравнения доказываются теорема единственности решения, а также теорема о достаточных условиях тотально (по множеству допустимых управлений) глобальной разрешимости при варьировании управления. При некоторых естественных предположениях, связанных с поточечными по времени $t$ оценками, заключение об однозначной тотально глобальной разрешимости делается, исходя из факта глобальной разрешимости системы сравнения, в качестве которой выступает система функционально-интегральных неравенств (можно заменить ее системой уравнений аналогичного типа, а в некоторых случаях - системой обыкновенных дифференциальных уравнений) относительно функций одного переменного $t\in[0;T]$ со значениями в пространстве $\mathbb{R}$. В качестве примера устанавливаются условия однозначной тотально глобальной разрешимости управляемой нелинейной нестационарной системы уравнений Навье-Стокса.

Пусть $U$ — множество допустимых управлений, $T>0$ и задана шкала банаховых пространств $W[0;\tau]$, $\tau\in(0;T]$, такая, что множество сужений функций из $W=W[0;T]$ на $[0;\tau]$ совпадает с $W[0;\tau]$, $F[\cdot;u]\colon W\to W$ — управляемый вольтерров оператор, $u\in U$. Для операторного уравнения $x=F[x;u]$, $x\in W$, вводится система сравнения в форме функционально-интегрального уравнения в пространстве $\mathbf{C}[0;T]$. Устанавливается, что при естественных предположениях относительно оператора $F$ для сохранения (относительно малых вариаций правой части) глобальной разрешимости операторного уравнения достаточно сохранения глобальной разрешимости указанной системы сравнения. Сам по себе этот факт аналогичен некоторым результатам, установленным автором ранее. Центральный результат статьи составляет ряд признаков устойчивой глобальной разрешимости функционально-интегральных уравнений, упомянутых выше, без предположения типа локальной липшицевости правой части. В качестве содержательного примера, представляющего самостоятельный интерес, рассматривается нелинейная нестационарная система Навье–Стокса в пространстве $\mathbb{R}^3$.

Численно исследуются газодинамические процессы, протекающие в начальный момент работы сверхзвукового сопла с высокой степенью геометрического расширения. Основное внимание уделяется изучению механизмов потери течением осевой симметрии за счет неустойчивости образующихся в сверхзвуковой части сопла зон отрывного течения. Модель нестационарного течения вязкого теплопроводного сжимаемого газа по соплу основана на системе уравнений сохранения в форме Навье-Стокса. Турбулентность исследуемого течения моделируется методом отсоединенных вихрей DES и его модификацией DDES с привлечением полуэмпирической модели Спаларта-Аллмараса. Выполнено сравнение распределения давления на стенке сопла, проекции годографа вектора тяги, мгновенных и осредненных картин течения с экспериментальными данными и численными результатами других авторов. Показано, что применение вихреразрешающего моделирования DES и DDES позволяет адекватно описать основные особенности течения и воспроизвести феномен возникновения боковой составляющей тяги сверхзвукового сопла при приемлемом уровне вычислительных затрат.

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

В статье рассматривается модельная задача несжимаемого течения жидкости и переноса тепла в коротком плоском канале с обратным уступом. Цель работы состоит в исследовании влияния граничного условия для потока тепла (температуры) на выходе из канала на характеристики теплопереноса внутри канала. Система уравнений Навье-Стокса и баланса тепла решаются численно с использованием равномерной сетки разрешением $6001\times301$ узлов. Для разностной аппроксимации пространственных производных используется метод контрольного объема второго порядка. Достоверность получаемых решений подтверждена для широкого диапазона числа Рейнольдса $(100 \leqslant \text{Re} \leqslant 1000)$ и числа Прандтля $\text{Pr} = 0.71$ путем сравнения с экспериментальными и теоретическими результатами, найденными в литературе. Анализируются картины течения, поля изотерм перегрева потока и поведение локального числа Нуссельта вдоль нагретой нижней стенки канала в зависимости от выбора выходного граничного условия для потока тепла (температуры). Показано, что этот выбор может оказать существенное влияние на характер прогрева течения внутри всего канала. По результатам исследования выбор сделан в пользу нелинейного граничного условия.

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

В работе рассматриваются результаты решения задачи стационарного течения вязкой несжимаемой жидкости в плоском канале с обратным уступом и прогреваемой нижней стенкой в широком диапазоне числа Рейнольдса $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$.

Представлена методика моделирования турбулентного течения вязкого газа, основанная на схеме высокого порядка аппроксимации WENO (взвешенная существенно неосциллирующая схема). Данная схема характеризуется значительной устойчивостью при выполнении расчетов, так как WENO позволяет устранять нефизичные осцилляции численного решения, которые могут возникнуть в ходе вычислений. Приведена система определяющих уравнений, описывающая поток вязкого газа, основанная на системе уравнений Навье-Стокса. Разработаны и реализованы алгоритмы 3-го и 5-го порядков точности. Приведено описание численных методик использованных в расчетах потока газа. Моделирование турбулентности производилось с применением метода крупных вихрей. Предложенные алгоритмы были использованы для исследования течения вязкого газа в канале с обратным уступом. Число Рейнольдса потока газа в канале составляло Re=15000. Проведено сравнение результатов численного моделирования с экспериментальными данными.

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