All issues

Теоретические и экспериментальные исследования взаимодействия водяного пара с пористыми материалами проводятся как на макро-, так и на микроуровне. На макроуровне исследуется влияние структуры расположения индивидуальных пор на процессы взаимодействия водяного пара с пористым материалом как сплошной средой. На микроуровне исследуется зависимость характеристик взаимодействия водяного пара с пористой средой от геометрии и размеров индивидуальной поры.

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

Theoretical and experimental investigations of water vapor interaction with porous materials are carried out both at the macro level and at the micro level. At the macro level, the influence of the arrangement structure of individual pores on the processes of water vapor interaction with porous material as a continuous medium is studied. At the micro level, it is very interesting to investigate the dependence of the characteristics of the water vapor interaction with porous media on the geometry and dimensions of the individual pore.

In this paper, a study was carried out by means of mathematical modelling of the processes of water vapor interaction with suffering pore of the cylindrical type. The calculations were performed using a model of a hybrid type combining a molecular-dynamic and a macro-diffusion approach for describing water vapor interaction with an individual pore. The processes of evolution to the state of thermodynamic equilibrium of macroscopic characteristics of the system such as temperature, density, and pressure, depending on external conditions with respect to pore, were explored. The dependence of the evolution parameters on the distribution of the diffusion coefficient in the pore, obtained as a result of molecular dynamics modelling, is examined. The relevance of these studies is due to the fact that all methods and programs used for the modelling of the moisture and heat conductivity are based on the use of transport equations in a porous material as a continuous medium with known values of the transport coefficients, which are usually obtained experimentally.

Проведен сравнительный анализ двух моделей пористой среды (Дарси и Бринкмана) на примере математического моделирования нестационарных режимов термогравитационной конвекции в пористой вертикальной цилиндрической полости с теплопроводной оболочкой конечной толщины в условиях конвективного охлаждения со стороны окружающей среды. Краевая задача математической физики, сформулированная в безразмерных переменных «функция тока — завихренность — температура», реализована численно неявным методом конечных разностей. Представлены результаты тестовых расчетов и влияния сеточных параметров, отражающие правомерность применения предлагаемого численного подхода. Установлены особенности класса сопряженных задач при использовании рассматриваемых моделей пористой среды.

Comparative analysis of two models of porous medium (Dacry and Brinkman) on an example of mathematical simulation of transient natural convection in a porous vertical cylindrical cavity with heat-conducting shell of finite thickness in conditions of convective cooling from an environment has been carried out. The boundary-value problem of mathematical physics formulated in dimensionless variables such as stream function, vorticity and temperature has been solved by implicit finite difference method. The presented verification results validate used numerical approach and also confirm that the solution is not dependent on the mesh size. Features of the conjugate heat transfer problems with considered models of porous medium have been determined.

Проведено компьютерное моделирование динамики нагрева воздушной фурмы доменной печи с помощью вычислительной среды конечно-элементного анализа DEFORM-2D. Исследовано влияние теплоизолирующей вставки, установленной в дутьевой канал с воздушным зазором и без зазора, а также газотермического покрытия на температурное поле воздушной фурмы доменной печи. Результаты моделирования показали значительное влияние теплоизолирующей вставки в дутьевой канал и воздушного зазора, отделяющего ее от внутреннего стакана, на температурное поле фурмы. При наличии вставки наблюдается градиент температуры по ее толщине до 540–555 °С, причем максимального значения температура вставки достигает на поверхности со стороны дутьевого канала. В то же время температура внутреннего стакана снижается на 35–40 °С по сравнению с фурмой без вставки. При наличии вставки с воздушным зазором градиент температуры вставки по ее толщине снижается до 160–250 °С по сравнению с вариантом без воздушного зазора, причем максимальное значение температуры поверхности вставки со стороны дутьевого канала также увеличивается. Температура внутреннего стакана также снижается еще на 15–20 °С по сравнению с вариантом без воздушного зазора. Однако наблюдается резкий градиент температуры воздушного зазора по его толщине до 760 °С из-за низкой теплопроводности воздуха. При наличии газотермического покрытия максимальная температура нагрева торца рыльной части снизилась до 326 °С, а максимальный градиент температуры по его толщине также снизился до 67 °С по сравнению с вариантом без покрытия. С помощью программного комплекса DEFORM-2D создана модель, имитирующая прогар фурмы вследствие контака с жидким чугуном. Показано, что через 40 с контакта с чугуном температура на поверхности рыльной части со стороны воды достигает 1050 °С, а через 100 с — 1060 °С, что практически равносильно прогару.

Study of work of heating equipment is an actual issue because it allows determining optimal regimes to reach highest efficiency. At that it is very helpful to use computer simulation to predict how different heating modes influence the effectiveness of the heating process and wear of heating equipment. Computer simulation provides results whose accuracy is proven by many studies and requires costs and time less than real experiments. In terms of present research, computer simulation of heating of air tuyere of blast furnace was realized with the help of FEM software. Background studies revealed possibility to simulate it as a flat, axisymmetric problem and DEFORM-2D software was used for simulation. Geometry, necessary for simulation, was designed with the help of SolidWorks, saved in .dxf format. Then it was exported to DEFORM-2D pre-processor and positioned. Preliminary and boundary conditions were set up. Several modes of operating regimes were under analysis. In order to demonstrate influence of eah of the modes and for better visualization point tracking option of the DEFORM-2D post-processor was applied. Influence of thermal insulation box plugged into blow channel, with and without air gap, and thermal coating on air tuyere’s temperature field was investigated. Simulation data demonstrated significant effect of thermal insulation box on air tuyere’s temperature field. Designed model allowed to simulate tuyere’s burnout as a result of interaction with liquid iron. Conducted researches have demonstrated DEFORM-2D effectiveness while using it for simulation of heat transfer and heating processes. DEFORM-2D is about to be used in further studies dedicated to more complex process connected with temperature field of blast furnace’s air tuyere.

В статье рассматривается устойчивость эффективной численной схемы, предложенной С.В. Патанкаром. Численная схема Патанкара нашла широкое применение в решении разнообразных прикладных задач, поэтому вопросы, связанные с математическим обоснованием этой схемы, являются достаточно актуальными.

In this paper we consider the issues of Patankar's numerical scheme stability. The Patankar’s numerical scheme is applied in the most number of the applications. So, the issues of Patankar's numerical scheme stability are very important question for the applications.