Все выпуски

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

We consider a problem with data on the characteristics for a loaded system of hyperbolic equations of the second order on a rectangular domain. The questions of the existence and uniqueness of the classical solution of the considered problem, as well as the continuity dependence of the solution on the initial data, are investigated. We propose a new approach to solving the problem with data on the characteristics for the loaded system of hyperbolic equations second order based on the introduction new functions. By introducing new unknown functions the problem is reduced to an equivalent family of Cauchy problems for a loaded system of differential with a parameters and integral relations. An algorithm for finding an approximate solution to the equivalent problem is proposed and its convergence is proved. Conditions for the unique solvability of the problem with data on the characteristics for the loaded system of hyperbolic equations of the second order are established in the terms of coefficient's system.

Проведено математическое моделирование конвективно-кондуктивно-радиационного теплообмена в кубической полости, заполненной прозрачной для излучения средой. Анализируемый объект представлял собой замкнутый объем с теплопроводными стенками конечной толщины, имеющими диффузно-серые внутренние поверхности. Внешние поверхности двух вертикальных стенок являлись изотермическими, а остальные внешние грани области решения - адиабатическими. Краевая задача сформулирована в безразмерных переменных «векторный потенциал-вектор завихренности-температура» в приближении Буссинеска и с учетом диатермичности сплошной среды. Анализ радиационного теплообмена проведен с использованием метода сальдо в варианте Поляка. Сформулированная нестационарная краевая задача реализована численно методом конечных разностей в широком диапазоне изменения числа Рэлея, коэффициента теплопроводности материала ограждающих твердых стенок и коэффициента излучения. Получены корреляционные соотношения для средних конвективного и радиационного чисел Нуссельта на характерной внутренней границе раздела сред. Проведено сравнение полученных результатов с данными двумерной модели. Установлено, что при рассмотрении трехмерной задачи можно оценить формирование интенсивных поперечных перетоков среды со стороны двух вертикальных поверхностей, которые отсутствуют в двумерной постановке. Показано, что решение задач конвективно-радиационного теплопереноса в сопряженной постановке приводит к существенным изменениям в распределениях локальных и интегральных характеристик по сравнению с несопряженной моделью, что в первую очередь связано с более корректным описанием механизма теплового излучения в диатермичных средах за счет учета теплопроводности ограждающих твердых стенок.

Mathematical simulation of convective-conductive-radiative heat transfer in a cubical cavity filled with diathermanous medium has been carried out. The domain of interest is a closed volume having heat-conducting solid walls of finite thickness with diffuse grey inner surfaces. The outer surfaces of two vertical walls are isothermal while the other walls are adiabatic. The boundary-value problem has been formulated in dimensionless variables such as “vector potential-vorticity vector-temperature’’ in the Boussinesq approximation and taking into account the diathermancy of the continuous medium. An analysis of surface thermal radiation has been conducted on the basis of the net-radiation method in the form of Poljak. The formulated transient boundary-value problem has been solved by finite difference method in a wide range of the Rayleigh number, thermal conductivity ratio and surface emissivity. Correlations for the average convective and radiative Nusselt numbers at the characteristic internal solid-fluid interface have been obtained. The comparison between the obtained three-dimensional results and the two-dimensional data has been conducted. It has been found, that on the basis of a three-dimensional model it is possible to analyze the formation of intensive transverse flows from two vertical surfaces which are absent in a two-dimensional model. It has been also shown, that the solution of convective-radiative heat transfer problems in the conjugate statement leads to essential changes in distributions of local and integral parameters in comparison with the non-conjugate model, which first of all is related to a more correct description of the thermal radiation in diathermanous media due to taking into account the thermal conduction of the solid walls.

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

The paper deals with the results of solving the problem of steady-state flow of a viscous incompressible fluid in a plane channel with a backward-facing step and a heated bottom wall for the Reynolds number in the range $100\leqslant \text{Re}\leqslant1000$ and the expansion ratio of a plane channel in the range $1.11 \leqslant ER \leqslant 10$. The study was carried out by numerical integration of the 2-D Navier-Stokes equations in velocity-pressure formulation on uniform grids with a step which equals to 1/300. Correction of the results is confirmed by comparing them with the literature data. Detailed flow patterns and fields of stream overheating depending on two basic parameters of the problem $\text{Re}$ and $ER$ are demonstrated. It is shown that with the increase of parameters $\text{Re}$ and $ER$ the structure of flow becomes much more complicated, that is, there is an increase of the number of vortices and their sizes up to the formation of a vortex behind the backward-facing step with two centers of rotation. It is also stated that the typical height of the heating zone of the flow depends weakly on $\text{Re}$ and $ER$ and eventually, near the exit of the channel, equals approximately half of the channel height. For all centers of vortices their main characteristics are defined: location, extremums of stream function, vorticity. Complex nonmonotonic behaviors of the coefficients of friction, hydrodynamic resistance and heat transfer (Nusselt number) along the channel are analyzed. It is shown that these coefficients strongly depend both on Reynolds number and on expansion ratio, reaching the maximum values at the maximum values of $\text{Re}$ and $ER$.