Все выпуски

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

The paper presents the model of anisotropic growth of dendritic crystallization of chemically pure and binary liquid (solution or melt) based on forced convection of the liquid phase. The dependencies of the growth rate and the radius of the top of a dendrite from under-cooling fluid in cases of a chemically pure material and alloys are presented. A comparative analysis of the influence of forced convection on the dendrite growth kinetics is carried out. Evaluation of growth rate and morphology of dendrite by high-speed crystal growth model was done. The contribution of convective flow and the anisotropic properties of the liquid-crystal boundary were taking into account. The model is also used hyperbolic diffusion equation to describe the non-equilibrium impurity capture by crystal surface, which occurs under the rapid crystals growth.

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

The formation process of a spherical shock impulse in gas and its interaction with a protective aqueous foam barrier, accompanied by formation of vortex flows, are numerically investigated. The problem is solved in a two-dimensional axisymmetric formulation using a two-phase model of a gas-liquid mixture, which includes the laws of conservation of mass, momentum and energy of the mixture and an equation for the dynamics of volume content of phases.The numerical implementation is carried out on the basis of the OpenFOAM package using the standard compressibleMultiphaseInterFoam solver, modified in accordance with the conditions of the problem and model representations. The discretization of the system of equations in the chosen solver is carried out by the method of finite volumes using the computational Pimple algorithm. A significant decrease in the intensity of the shock wave in its interaction with the aqueous foam barrier is shown and the causes leading to vortex formation in the gas region are revealed. The reliability of the results obtained is estimated by comparison with solutions of a similar problem by other numerical methods.

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

Numerical simulation flow of viscoelastic fluid with free surface, which is realized in entrance and output flow  in extrusion die was performed. The flow of liquid is described by equations of conservation of mass, momentum and thermal energy with rheological constitutive equation of Giesekesus.  On basis of finite element method the  stable numerical scheme was developed  to solve  this problem. Different numerical experiments was performed to define the configuration of outflow jet  in various regimes and construction of die. The distribution of flow velocity fields, pressure and temperature  are investigated on dependence of heating the walls. The ratio of extrusion in dependence of parameters the rheological model are investigated.

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

The problem of choice of boundary conditions is discussed for the case of numerical integration of the shallow water equations on a substantially irregular relief. While modeling unsteady surface water flows there is a dynamic boundary that partitions liquid and dry bottom. The situation is complicated by the emergence of sub- and supercritical flow regimes for the problems of seasonal floodplain flooding, flash floods, tsunami landfalls. Analysis of the use of various methods of setting conditions for the physical quantities of liquid at the settlement of the boundary shows the advantages of using the waterfall type conditions in the presence of strong heterogeneities of landforms. When there is a waterfall on the border of computational domain and heterogeneity of the relief in the vicinity of the boundary, a portion may occur which is formed by the region of critical flow with the formation of a hydraulic jump, which greatly weakens the effect of the waterfall on the flow pattern upstream.

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

The results of numerical simulation of the structure of a two-phase flow of a gas-liquid bubble mixture in a vertical ascending flow in a pipe are presented. The mathematical model is based on the use of the two-fluid Eulerian approach taking into account the inverse influence of bubbles on averaged characteristics and turbulence of the carrying phase. The turbulent kinetic energy of a liquid is calculated using equations for the transfer of Reynolds stresses. To describe the dynamics of bubble size distribution, the equations of particle number conservation for individual groups of bubbles with different constant diameters for each fraction are used taking into account the processes of breakup and coalescence. The influence of changes in the degree of dispersion of the gas phase, volume flow gas content and the velocity of the dispersed phase on the local structure and surface friction in the two-phase flow is numerically investigated. Comparison of simulation results with experimental data has shown that the developed approach allows an adequate description of turbulent gas-liquid flows in a wide range of changes in gas content and initial bubble sizes.

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

The results of the modeling of hydrodynamic processes in the operating space of 3 different types of fixed flow capillary viscometers are represented in the article. The results were obtained from computational solution of the Navier-Stokes equation for laminar flow with the use of finite-element method. The influence of capillary tube and cylinder bottom shape on the metrological performance of viscometer was established.

Рассматриваются методы моделирования взаимодействия потока несжимаемой жидкости и преграды в рамках эйлерова (метод объема жидкости в ячейке, Volume of Fluid - VOF) и лагранжева (метод гидродинамики сглаженных частиц, Smoothed Particle Hydrodynamics - SPH) описаний. На примере решения задач о движении потока жидкости, вызванного распадом начального уровня жидкости (задача о разрушении плотины), оцениваются преимущества и недостатки применения метода SPH для моделирования гидродинамических нагрузок на преграду, развитой свободной поверхности и каплеобразования. Определяется влияние способа конкретной численной реализации граничных условий Дирихле на твердых стенках на величину давления и характер ее изменения во времени. Численные результаты, полученные с использованием методов VOF и SPH, сопоставляются с известными экспериментальными данными.

The paper considers the methods of modeling of the incompressible fluid flow interaction with barriers in Euler formulation (volume of fluid - VOF) and Lagrangian (smoothed particle hydrodynamics - SPH) description. By the example of solving the problems of motion of the fluid flow caused by the collapse of the initial liquid level (dam break problem), the authors estimate advantages and disadvantages of using the SPH method for the simulation of hydrodynamic loads, free-form surface and formation of drops. The influence of the specific numerical implementation of the Dirichlet boundary conditions on solid walls on both the pressure magnitude and its time behavior is determined. Numerical results obtained by the methods of VOF and SPH are compared with known experimental data.