All issues

Разработан метод преобразования поверхностного представления расчетной области в формат воксельных данных. Предложен алгоритм генерации расчетной шестигранной сетки на основе объемного формата данных.

The method of conversion of surface representation to the voxel data is developed. Voxel data based algorithm of hexahedral mesh generation is proposed.

В настоящее время для численного моделирования начально-краевых задач для систем гиперболических уравнений в частных производных (например, уравнения газовой динамики, МГД, деформируемого твердого тела и т. д.) применяются различные нелинейные численные схемы пространственной аппроксимации. Это связано с необходимостью повышения порядка аппроксимации и расчета разрывных решений, часто возникающих в таких системах. Необходимость в нелинейных схемах связана с ограничением, следующим из теоремы С. К. Годунова о невозможности построения линейной схемы порядка больше первого для монотонной аппроксимации уравнений такого типа. Одними из наиболее точных нелинейных схем являются схемы типа ENO (существенно не осциллирующие схемы и их модификации), в том числе схемы WENO (взвешенные, существенно не осциллирующие схемы). Последние получили наибольшее распространение, поскольку при одинаковой ширине шаблона имеют более высокий порядок аппроксимации чем ENO-схемы. Плюсом ENO- и WENO-схем является сохранение высокого порядка аппроксимации на немонотонных участках решения. Исследование данных схем затруднительно в связи с тем, что сами схемы нелинейны и применяются для аппроксимации нелинейных уравнений. В частности, условие линейной устойчивости ранее было получено только для схемы WENO5 (пятого порядка аппроксимации на гладких решениях) и является приближенным. В настоящей работе рассматриваются вопросы построения и устойчивости схем WENO5, WENO7, WENO9, WENO11 и WENO13 для конечно-объемной схемы для уравнения Хопфа. В первой части статьи рассмотрены методы WENO в общем случае и приведены явные выражения для коэффициентов полиномов и весов линейных комбинаций, необходимых для построения схем. Доказывается ряд утверждений, позволяющих сделать выводы о порядках аппроксимации в зависимости от локального вида решения. Проводится анализ устойчивости на основе принципа замороженных коэффициентов. Рассматриваются случаи гладкого и разрывного поведения решения в области линеаризации при замороженных коэффициентах на гранях конечного объема и анализируется спектр схем для этих случаев. Доказываются условия линейной устойчивости для различных методов Рунге–Кутты при применении со схемами WENO. В результате приводятся рекомендации по выбору максимально возможного параметра устойчивости, которое наименьшим образом влияет на нелинейные свойства схем. Следуя полученным ограничениям, делается вывод о сходимости схем.

Currently, different nonlinear numerical schemes of the spatial approximation are used in numerical simulation of boundary value problems for hyperbolic systems of partial differential equations (e. g. gas dynamics equations, MHD, deformable rigid body, etc.). This is due to the need to improve the order of accuracy and perform simulation of discontinuous solutions that are often occurring in such systems. The need for non-linear schemes is followed from the barrier theorem of S. K. Godunov that states the impossibility of constructing a linear scheme for monotone approximation of such equations with approximation order two or greater. One of the most accurate non-linear type schemes are ENO (essentially non oscillating) and their modifications, including WENO (weighted, essentially non oscillating) scemes. The last received the most widespread, since the same stencil width has a higher order of approximation than the ENO scheme. The benefit of ENO and WENO schemes is the ability to maintain a high-order approximation to the areas of non-monotonic solutions. The main difficulty of the analysis of such schemes comes from the fact that they themselves are nonlinear and are used to approximate the nonlinear equations. In particular, the linear stability condition was obtained earlier only for WENO5 scheme (fifth-order approximation on smooth solutions) and it is a numerical one. In this paper we consider the problem of construction and stability for WENO5, WENO7, WENO9, WENO11, and WENO13 finite volume schemes for the Hopf equation. In the first part of this article we discuss WENO methods in general, and give the explicit expressions for the coefficients of the polynomial weights and linear combinations required to build these schemes. We prove a series of assertions that can make conclusions about the order of approximation depending on the type of local solutions. Stability analysis is carried out on the basis of the principle of frozen coefficients. The cases of a smooth and discontinuous behavior of solutions in the field of linearization with frozen coefficients on the faces of the final volume and spectra of the schemes are analyzed for these cases. We prove the linear stability conditions for a variety of Runge-Kutta methods applied to WENO schemes. As a result, our research provides guidance on choosing the best possible stability parameter, which has the smallest effect on the nonlinear properties of the schemes. The convergence of the schemes is followed from the analysis.

В работе представлена новая версия программного комплекса FlowVision, предназначенного для автоматизации инженерных расчетов в области вычислительной гидродинамики: FlowVision 3.09.05. Программный комплекс (ПК) FlowVision используется для решения различных прикладных задач в различных областях промышленности. Его популярность основана на том, что он позволяет решать сложные нетрадиционные задачи, находящиеся на стыке различных дисциплин, с одной стороны, и, с другой стороны, на парадигме полной автоматизации таких трудоемких для инженера процессов, как построение расчетной сетки. FlowVision — это программный комплекс, полностью отчуждаемый от разработчиков. Он имеет развитый графический интерфейс, систему задания расчетного проекта и систему визуализации течений различными методами — от построения контуров (для скалярных переменных) и векторов (для векторных переменных) на плоскостях и поверхностях до объемной визуализации расчетных данных. Кроме этого, ПК FlowVision предоставляет пользователю возможность вычислять интегральные характеристики на поверхностях и в ограниченных объемах.

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

FlowVision позволяет решать задачи движения жидкостей со скоростями, соответствующими несжимаемому или гиперзвуковому режимам за счет использования все-скоростного метода расщепления по физическим переменным для решения уравнений Навье–Стокса. FlowVision позволяет решать междисциплинарные задачи с использованием различных средств моделирования, например: моделировать многофазные течения методом VOF, обтекание подвижных тел с помощью эйлерова подхода при неподвижной расчетной сетке, моделировать вращающиеся машины с использованием метода скользящей сетки, решать задачи взаимодействия жидкости и конструкций методом двухстороннего сопряжения FlowVision с конечно-элементными кодами. В данной работе показаны примеры решения задач-вызовов: a) посадка космического корабля на воду при торможении ракетными двигателями, где есть граница раздела «воздух–вода», подвижные тела и взаимодействие сверхзвуковой струи газа с границей раздела «вода–воздух»; б) моделирование работы человеческого сердца с искусственными и живыми клапанами, спроектированными на базе томографических исследований, с использованием двухстороннего сопряжения «жидкостной» расчетной области с конечно-элементной моделью мышц сердца.

The work submits new release of the FlowVision software designed for automation of engineering calculations in computational fluid dynamics: FlowVision 3.09.05. The FlowVision software is used for solving different industrial problems. Its popularity is based on the capability to solve complex non-tradition problems involving different physical processes. The paradigm of complete automation of labor-intensive and time-taking processes like grid generation makes FlowVision attractive for many engineers. FlowVision is completely developer-independent software. It includes an advanced graphical interface, the system for specifying a computational project as well as the system for flow visualization on planes, on curvilinear surfaces and in volume by means of different methods: plots, color contours, iso-lines, iso-surfaces, vector fields. Besides that, FlowVision provides tools for calculation of integral characteristics on surfaces and in volumetric regions.

The software is based on the finite-volume approach to approximation of the partial differential equations describing fluid motion and accompanying physical processes. It provides explicit and implicit methods for time integration of these equations. The software includes automated generator of unstructured grid with capability of its local dynamic adaptation. The solver involves two-level parallelism which allows calculations on computers with distributed and shared memory (coexisting in the same hardware). FlowVision incorporates a wide spectrum of physical models: different turbulence models, models for mass transfer accounting for chemical reactions and radioactive decay, several combustion models, a dispersed phase model, an electro-hydrodynamic model, an original VOF model for tracking moving interfaces. It should be noted that turbulence can be simulated within URANS, LES, and ILES approaches. FlowVision simulates fluid motion with velocities corresponding to all possible flow regimes: from incompressible to hypersonic. This is achieved by using an original all-speed velocity-pressure split algorithm for integration of the Navier-Stokes equations.

FlowVision enables solving multi-physic problems with use of different modeling tools. For instance, one can simulate multi-phase flows with use of the VOF method, flows past bodies moving across a stationary grid (within Euler approach), flows in rotary machines with use of the technology of sliding grid. Besides that, the software solves fluid-structure interaction problems using the technology of two-way coupling of FlowVision with finite-element codes. Two examples of solving challenging problems in the FlowVision software are demonstrated in the given article. The first one is splashdown of a spacecraft after deceleration by means of jet engines. This problem is characterized by presence of moving bodies and contact surface between the air and the water in the computational domain. The supersonic jets interact with the air-water interphase. The second problem is simulation of the work of a human heart with artificial and natural valves designed on the basis of tomographic investigations with use of a finite-element model of the heart. This problem is characterized by two-way coupling between the “liquid” computational domain and the finite-element model of the hart muscles.

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

The paper presents results of experimental, computational and analytical study of the effect of nonequilibrium chemical activation of air-fuel mixture on effectiveness of Diesel process. The generation of a high-voltage multi-streamer discharge in combustion chamber at the compression phase is considered as the method of the activation. The description of electrical discharge system, results of measurement and visualization are presented. The plasma-chemical kinetics of nonequilibrium ignition is analyzed to establish a passway for a proper reduction of chemical kinetics scheme. The results of numerical simulation of gas dynamic processes at presence of plasma-assisted combustion in a geometrical configuration close to the experimental one are described.

В работе представлены результаты численного моделирования влияния протяженного искрового разряда на динамику перемешивания инжектируемой газовой струи со сверхзвуковым воздушным потоком. Расчеты проводились в программном комплексе FlowVision. Подача топлива осуществляется при помощи инжектора, расположенного на стенке канала, а разряд организован вблизи стенки ниже по потоку относительно инжектора. Моделирование электрического искрового разряда выполнено при помощи объемного источника тепла. С целью описания принципиального вида плазменного актуатора для ускорения перемешивания в сверхзвуковом потоке (число Маха М = 2) в ходе исследования выполнено варьирование энерговклада в разряд в диапазоне 100–500 мДж на один импульс, а также определено влияние формы и местоположения разряда относительно топливного инжектора. Проведено исследование режимов инжекции топлива в сверхзвуковой воздушный поток и найден оптимальный режим истечения струи газа для исследования влияния искрового разряда на смешение. Разработан метод анализа картины возмущений границы раздела «топливо–окислитель», вызванных работой импульсного искрового разряда. Подготовлена программа в среде LabView для получения количественной характеристики для дальнейшего сравнения полученных результатов с экспериментальными данными.

Результаты моделирования позволяют сделать вывод, что протяженный искровой разряд, расположенный ниже по потоку относительно инжектора и расположенный вдоль потока, обеспечивает максимальное увеличение границы раздела между струей топлива и основным потоком. Типичная частота повторения импульсов разряда в импульсно-периодическом режиме должна составлять более 6 кГц при длине разряда ~10 мм, чтобы обеспечить постоянное влияние на смешение в потоке со скоростью 500 м/с.

The paper presents the results of numerical simulation of the effect of a long spark discharge on the mixing dynamics of an injected gas jet with supersonic air flow. The calculations were performed using the CFD software package FlowVision. The fuel was supplied using an injector located on the channel wall, and the discharge was organized near the wall downstream of the injector. Simulation of electrical spark discharge was performed using a volumetric heat source. In order to describe the principal specifications of a plasma actuator to accelerate mixing in a supersonic flow (Mach number M = 2), the research involved varying the energy impact to the discharge in the range of 100–500 mJ per pulse, determining the influence of the shape and location of the discharge. A study of the fuel injection modes in a supersonic air flow has been carried out and an optimal gas jet outflow regime has been found to study the effect of a spark discharge. A method has been developed for analyzing the disturbance pattern of the fuel-oxidant interface caused by the operation of a pulsed spark discharge. A program was prepared in the LabView software environment for obtaining quantitative characteristics for further comparison with the results obtained in the experiment.

The simulation results allow us to conclude that the long spark discharge located along the flow downstream of the injector provides the maximum increase in the interface between the jet of fuel and the main flow. A typical repetition frequency of discharge pulses in a pulse-periodic mode should be more than 6 kHz with a discharge length of ~10 mm to ensure a continuous effect on the mixing at a flow velocity of 500 m/s.

Методом молекулярной динамики изучен эффект нелинейной супратрансмиссии в кристалле стехиометрии А₃В, на примере Pt₃Al, заключающийся в передаче энергии на частотах вне фононного спектра кристалла. Исследование механизмов транспорта энергии от поверхности материала вглубь является важной задачей как с теоретической точки зрения, так и с точки зрения перспектив практического применения при модификации приповерхностных слоев обработкой интенсивными внешними воздействиями различного характера. Модель представляла собой объемный гранецентрированный кубический кристалл, атомы которого взаимодействовали посредством многочастичного потенциала, полученного методом погруженного атома, что обеспечивает большую реалистичность модели по сравнению с применением парных потенциалов. Рассмотрены разные формы осцилляции области внешнего воздействия. Показана возможность транспорта энергии от поверхности кристалла вглубь посредством возбуждения квазибризеров вблизи области воздействия и последующего их разрушения в кристалле и рассеяния запасенной на них энергии. Отметим, что под квазибризерами понимаются высокоамплитудные нелинейные колебания атомов легкого компонента сплава на частотах вне фононного спектра кристалла. При этом установлено, что не при любой форме осцилляции области воздействия наблюдался данный эффект. Наиболее интенсивно квазибризеры возникали вблизи области воздействия при синусоидальной форме колебаний. Полученные результаты свидетельствуют, что вклад квазибризеров в передачу энергии по кристаллу возрастает при увеличении амплитуды воздействия. Рассмотрен диапазон амплитуд от 0.05 до 0.5 Å. Частота воздействия варьировалась от 0.2 до 15 ТГц, что обеспечивало охват всего спектра малоамплитудных колебаний для данной модели кристалла. Установлена минимальная величина амплитуды внешнего воздействия, при которой наблюдался данный эффект, которая составила 0.2 Å. При амплитудах более 0.5 Å происходит быстрое разрушение ячейки для частот, близких к оптической ветви фононного спектра. Результаты проведенного исследования могут быть полезны при лазерной обработке материалов и обработке поверхности низкоэнергетической плазмой, а также в радиационном материаловедении.

The effect of the nonlinear supratransmission in crystal of A₃B stoichiometry is studied by molecular dynamics on the example of Pt₃Al alloy. This effect is the transfer of energy at frequencies outside the phonon spectrum of the crystal. Research of the mechanisms of energy transport from the material surface to the interior is the important task, both from the theoretical point of view and from the prospects for practical application in the modification of near-surface layers by treatment with intense external influence of various types. The model was a three-dimensional face-centered cubic crystal whose atoms interact by means of the multiparticle potential obtained by the embedded atom method, which provides greater realism of the model in comparison with the use of pair potentials. Various forms of oscillation of the external influence region are considered. The possibility of energy transport from the crystal surface to the interior is shown by excitation of quasi-breathers near the region of influence and their subsequent destruction in the crystal and scattering of the energy stored on them. The quasibreathers are high-amplitude nonlinear atoms' oscillations of the alloy lightweight component at frequencies outside the phonon spectrum of the crystal. This effect was observed not with every oscillation's form of the region of influence. Quasi-breathers appeared most intensely near the region of influence with sinusoidal form oscillations. The results obtained indicate that the contribution of quasi-breathers to the energy transfer through the crystal increases with increasing amplitude of the influence. The range of amplitudes from 0.05 to 0.5 Å is considered. The frequency of the influence varied from 0.2 to 15 THz, which ensured the coverage of the entire spectrum of lowamplitude oscillations for this crystal's model. The minimum magnitude of the external effect amplitude at which this effect was observed was found to be 0.15 Å. At amplitudes greater than 0.5 Å, the cell rapidly decays for frequencies close to the optical branch of the phonon spectrum. The results of the study can be useful for laser processing of materials, surface treatment by low-energy plasma, and also in radiation materials science.

В работе рассматривается задача стационарной смешанной конвекции и теплообмена вязкой теплопроводной жидкости в плоской квадратной каверне с подвижной верхней крышкой. Нагретая верхняя стенка каверны имеет температуру $T_{\mathrm{H}}$, холодная нижняя — $T_\mathrm{0}$ $(T_\mathrm{H} > T_\mathrm{0})$, а боковые стенки каверны теплоизолированы. Особенностью задачи является тот факт, что плотность жидкости может принимать произвольные значения в зависимости от величины перегрева крышки каверны. Математическая постановка включает в себя уравнения Навье–Стокса в переменных «скорость–давление» и баланса тепла, сформулированные с учетом несжимаемости течения жидкости и воздействия объемной силы плавучести. Разностная аппроксимация исходных дифференциальных уравнений выполнена методом контрольного объема. Численные решения задачи получены на сетке $501 \times 501$ для следующих значений параметров подобия: число Прандтля Pr = 0.70; число Рейнольдса Re = 100, 1000; число Ричардсона Ri = 0.1, 1, 10 и относительный перегрев верхней стенки $(T_\mathrm{H} − T_\mathrm{0})/T_\mathrm{0} = 0, 1, 2, 3$. Достоверность полученных результатов подтверждена их сравнением с литературными данными. Представлены подробные картины течения в виде линий тока и изотерм перегрева потока. Показано, что увеличение значения числа Ричардсона (рост влияния силы плавучести) приводит к принципиальному изменению структуры течения жидкости. Также установлено, что учет переменности плотности жидкости приводит к ослаблению влияния роста Ri на трансформацию структуры течения. Это связано с тем, что изменение плотности в замкнутом объеме всегда приводит к возникновению зон с отрицательной плавучестью. Как следствие, конкуренция положительных и отрицательных объемных сил приводит в целом к ослаблению эффекта плавучести. Также проанализировано поведение коэффициентов теплоотдачи (числа Нуссельта) и трения вдоль нижней стенки каверны в зависимости от параметров задачи. Выявлено, что влияние переменности плотности на эти коэффициенты тем больше, чем большие значения при прочих равных условиях принимает число Ричардсона.

The paper considers the problem of stationary mixed convection and heat transfer of a viscous heatconducting fluid in a plane square lid-driven cavity. The hot top cover of the cavity has any temperature $T_\mathrm{H}$ and cold bottom wall has temperature $T_\mathrm{0} (T_\mathrm{H} > T_\mathrm{0})$, whereas in contrast the side walls are insulated. The fact that the fluid density can take arbitrary values depending on the amount of overheating of the cavity cover is a feature of the problem. The mathematical formulation includes the Navier–Stokes equations in the ’velocity–pressure’ variables and the heat balance equation which take into account the incompressibility of the fluid flow and the influence of volumetric buoyancy force. The difference approximation of the original differential equations has been performed by the control volume method. Numerical solutions of the problem have been obtained on the $501 \times 501$ grid for the following values of similarity parameters: Prandtl number Pr = 0.70; Reynolds number Re = 100 and 1000; Richardson number Ri = 0.1, 1, and 10; and the relative cover overheating $(T_\mathrm{H}-T_\mathrm{0})/T_\mathrm{0} = 0, 1, 2, 3$. Detailed flow patterns in the form of streamlines and isotherms of relative overheating of the fluid flow are given in the work. It is shown that the increase in the value of the Richardson number (the increase in the influence of buoyancy force) leads to a fundamental change in the structure of the liquid stream. It is also found out that taking into account the variability of the liquid density leads to weakening of the influence of Ri growth on the transformation of the flow structure. The change in density in a closed volume is the cause of this weakening, since it always leads to the existence of zones with negative buoyancy in the presence of a volumetric force. As a consequence, the competition of positive and negative volumetric forces leads in general to weakening of the buoyancy effect. The behaviors of heat exchange coefficient (Nusselt number) and coefficient of friction along the bottom wall of the cavity depending on the parameters of the problem are also analyzed. It is revealed that the greater the values of the Richardson number are, the greater, ceteris paribus, the influence of density variation on these coefficients is.