Все выпуски

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

For the first time, numerical modeling has determined the mechanism of the initiated-by-cavitation rise of the liquid level in tubes and capillaries, known as the laser-induced optocapillary effect, as well as its analogues — the acoustocapillary and plasmocapillary effects. It is shown that the key condition for the occurrence of the liquid rise is the asymmetric collapse of a single relatively large cavitation bubble inside a vertically oriented tube or capillary. The proximity of boundaries (the tube wall, the fiber optic tip, and others) disrupts the spherical symmetry of the bubble during its collapse, leading to the appearance of a liquid flow that rolls up into a long-lived toroidal vortex ring. Due to viscous entrainment of the surrounding medium, the vortex generates a directed liquid flow upward and also ensures the suction of a new portion of liquid through the open lower end of the tube. The simulation results show that the characteristic lifetime of the toroidal vortex significantly exceeds the duration of the growth and collapse stages of the cavitation bubble that generated it. It is demonstrated that the rise of the liquid level in the tube does not begin at the moment of bubble expansion, but after its complete disappearance, and continues over a relatively long period due to the inertia of the vortex motion. This result is in complete agreement with experimental data, confirming the validity of the proposed mechanism.

The study investigated the practically significant configuration of the laser-induced optocapillary effect using an optical fiber. This configuration opens broad prospects for technical and medical applications, particularly in laser surgery. The investigated mechanism can be used to create cavitation pumps — effective tools for cleaning technical surfaces and wound surfaces, where the process of removing damaged tissue and foreign bodies due to thermal exposure will be accompanied by the removal of debris through the tube, significantly increasing the efficiency and safety of the procedure.

The obtained results represent the first consistent explanation for a class of cavitation-induced capillary phenomena and create a foundation for their controlled application in biomedical and microfluidic technologies.

Рассматриваются вопросы моделирования течения в дальних областях осесимметричного безымпульсного турбулентного следа, возникающего при обтекании тела. Основными особенностями исследуемого течения являются его практически бессдвиговый характер и близость к автомодельному режиму, что подтверждается имеющимися экспериментальными данными. Для описания течения в дальнем следе привлекается полуэмпирическая $e-\varepsilon$-модель турбулентности, основанная на алгебраической модели рейнольдсовых напряжений Роди. С помощью теоретико-группового анализа получена автомодельная редукция привлекаемой модели турбулентности к системе обыкновенных дифференциальных уравнений. Проведено построение приближенного автомодельного решения второго рода соответствующей краевой задачи для редуцированной системы уравнений с помощью метода, основанного на сшивании усечений асимптотических разложений ее решения, найденных в окрестностях граничных точек. Привлекаемый подход также позволил осуществить выделение приближенного значения показателя автомодельности и других неизвестных параметров задачи в процессе выполнения процедуры сшивания асимптотических разложений. Сопоставление сконструированного в виде рядов Пюизё приближенного решения поставленной краевой задачи, включающего 16 аппроксимирующих слагаемых, с численным решением, найденным с помощью метода стрельбы, показало, что максимальная относительная погрешность между решениями не превосходит 2%. Установленное в данной работе приближенное значение показателя автомодельности задачи хорошо согласуется с точным значением, полученным ранее другими авторами в результате аналитического решения соответствующей нелинейной задачи на собственные значения. Проведено сопоставление приближенного решения с доступными экспериментальными данными по осесимметричному безымпульсному турбулентному следу, полученными в ИГиЛ СО РАН, имеет место удовлетворительное соответствие. Построенные в работе приближенные формулы имеют практическое значение и могут быть использованы для быстрой оценки значений гидродинамических величин в дальних областях безымпульсного осесимметричного турбулентного следа. Развиваемый метод построения приближенных решений краевых задач свободной турбулентности может найти применение и в других прикладных областях, которые не связаны с теорией турбулентности.

In this paper we consider the problem of modeling of flow in the far regions of an axisymmetric momentumless turbulent wake that occurs when a body is flowed. The main features of the studied flow are its almost shear-free nature and its proximity to the self-similar regime, which is confirmed by the available experimental data. A semi-empirical $e-\varepsilon$ model of turbulence based on the Rodi algebraic model of Reynolds stresses is invoked to describe a flow in the far wake. Using group-theoretic analysis, similarity reduction of the invoked model of turbulence to a system of ordinary differential equations is obtained. An approximate self-similar solution of the second kind of the corresponding boundary value problem for the reduced system of equations was constructed using a method based on matching the truncations of asymptotic expansions of its solution found in neighborhoods of the boundary points. The involved approach also allowed for the extraction of an approximate value of the self-similarity exponent and other unknown parameters of the problem during the procedure of matching of the asymptotic expansion. Comparison of the approximate solution of the given boundary value problem constructed in the form of Puiseux series, which includes 16 approximating terms, with the numerical solution found using the shooting method, showed that the maximum relative error between the solutions does not exceed 2%. The approximate value of the self-similarity exponent of the problem established in this work is in a good agreement with the exact value obtained earlier by other authors as a result of the analytical solution of the corresponding nonlinear eigenvalue problem. The approximate solution was compared with the available experimental data on the axisymmetric momentumless turbulent wake obtained at the Lavrentyev Institute of Hydrodynamics of the Siberian Branch of the Russian Academy of Sciences, and a satisfactory agreement was found. The approximate formulas developed in this work have practical significance and can be used for quick estimation of the values of hydrodynamic quantities in the far regions of the axisymmetric momentumless turbulent wake. The developed method for constructing approximate solutions to the boundary value problems of free turbulence can also be used in other applied fields that are not related to turbulence theory.

Предложена модель расчета межслоевой магнитостатической энергии и оценка ее вклада в общую магнитную энергию многослойных пленок. Данная модель была применена для расчета трехслойной структуры Co/Cu/Co. Сравнение полученных теоретических результатов с экспериментальными данными позволило оценить влияние магнитостатического взаимодействия на величину и характер поведения поля насыщения подобных структур.

The computation model of the interlayer magnetostatic energy and its contribution into the total magnetic energy of multilayer films has been proposed. The model has been used for the energy computation of three-layer Co/Cu/Co structure. The comparison of obtained theoretical results with experimental data has been allowed made an estimation of magnetostatic interaction influence on a value and a character of saturation field behavior for such structures.

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

The results of the experimental and theoretical researches of kinetics of erythrosine penetration into collagenase molecules have represented in this paper. The case with introduction of the compound (fluorescent probe) which has dimers to enzyme solution as an injection has been considered. It was shown that monomers and dimers can penetrate into enzyme molecules with formation complexes monomer — enzyme, dimer- enzyme. Moreover, transformation of probe fluorescence spectra is at each time moment. Spectrum maximum shift, and its form change. At a time, the immobilized dye dimers greatly impact to formation of end fluorescence spectrum. Well correlation between experimental and theoretical results confirms reality of the obtained data.

В работе рассматриваются различные модели гидравлического гистерезиса, возникающего при инвазивной ртутной порометрии. Для моделирования гидравлического гистерезиса используется изотропная перколяция узлов на трехмерных квадратных решетках с $(1,\,\pi)$-окрестностью. Феноменологически исследуется взаимосвязь данных инвазивной порометрии с параметрами перколяционной модели. Реализация перколяционной модели основана на библиотеках SPSL и SECP, выпущенных под лицензией GNU GPL-3 с использованием свободного языка программирования R.

In this paper we consider various models of hydraulic hysteresis in invasive mercury porosimetry. For simulating the hydraulic hysteresis is used isotropic site percolation on three-dimensional square lattices with $(1,\,\pi)$-neighborhood. The relationship between the percolation model parameters and invasive porosimetry data is studied phenomenologically. The implementation of the percolation model is based on libraries SPSL and SECP, released under license GNU GPL-3 using the free programming language R.

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

The paper discusses the possibility of modeling the strength anisotropy of layered elastic medium using a scalar damage parameter. Thermodynamically consistent constitutive equations are formulated. Using SIMULIA / Abaqus we numerically simulated the stretching and compression of the samples. The results of calculation using the proposed model are compared with the known experimental data from the literature and the predictions of traditional models.

В основу системы интегрированной интерпретации геофизических данных при изучении глубинного строения Земли положена система ГИС ИНТЕГРО, являющаяся геоинформационной системой функционирования разнообразных вычислительных и аналитических приложений при решении различных геологических задач. ГИС ИНТЕГРО включает в себя многообразные интерфейсы, позволяющие изменять форму представления данных (растр, вектор, регулярная и нерегулярная сеть наблюдений), блок преобразования картографических проекций, а также прикладные блоки, включающие блок интегрированного анализа данных и решения прогнозно-диагностических задач.

Методический подход базируется на интеграции и комплексном анализе геофизических данных по региональным профилям, геофизических потенциальных полей и дополнительной геологической информации на изучаемую территорию.

Аналитическое обеспечение включает пакеты трансформаций, фильтрации, статистической обработки полей, расчета характеристик, выделения линеаментов, решения прямых и обратных задач, интегрирования геоинформации.

Технология и программно-аналитическое обеспечение апробировались при решении задач тектонического районирования в масштабах 1:200000, 1:1000000 в Якутии, Казахстане, Ростовской области, изучения глубинного строения по региональным профилям 1:ЕВ, 1-СБ, 2-СБ, 3-СБ и 2-ДВ, прогноза нефтегазоносности в районах Восточной Сибири, Бразилии.

GIS INTEGRO is the geo-information software system forming the basis for the integrated interpretation of geophysical data in researching a deep structure of Earth. GIS INTEGRO combines a variety of computational and analytical applications for the solution of geological and geophysical problems. It includes various interfaces that allow you to change the form of representation of data (raster, vector, regular and irregular network of observations), the conversion unit of map projections, application blocks, including block integrated data analysis and decision prognostic and diagnostic tasks.

The methodological approach is based on integration and integrated analysis of geophysical data on regional profiles, geophysical potential fields and additional geological information on the study area. Analytical support includes packages transformations, filtering, statistical processing, calculation, finding of lineaments, solving direct and inverse tasks, integration of geographic information.

Technology and software and analytical support was tested in solving problems tectonic zoning in scale 1:200000, 1:1000000 in Yakutia, Kazakhstan, Rostov region, studying the deep structure of regional profiles 1:S, 1-SC, 2-SAT, 3-SAT and 2-DV, oil and gas forecast in the regions of Eastern Siberia, Brazil.

The article describes two possible approaches of parallel calculations for data processing 2D or 3D nets in the field of geophysical research. As an example presented realization in the environment of GRID of the application software ZondGeoStat (statistical sensing), which create 3D net model on the basis of data 2d net. The experience has demonstrated the high efficiency of the use of environment of GRID during realization of calculations in field of geophysical researches.

В работе на основе предложенной ранее русловой модели решена одномерная задача устойчивости песчаного дна напорного канала. Особенностью исследуемой задачи является используемое оригинальное уравнение русловых деформаций, учитывающее влияние физико-механических и гранулометрических характеристик донного материала и неровности донной поверхности при русловом анализе. Еще одной особенностью рассматриваемой задачи является учет влияния не только придонного касательного, но и нормального напряжения при изучении русловой неустойчивости. Из решения задачи устойчивости песчаного дна для напорного канала получена аналитическая зависимость, определяющая длину волны для быстрорастущих донных возмущений. Выполнен анализ полученной аналитической зависимости, показано, что она обобщает ряд известных эмпирических формул: Коулмана, Шуляка и Бэгнольда. Структура полученной аналитической зависимости указывает на существование двух гидродинамических режимов, характеризуемых числом Фруда, при которых рост донных возмущений может сильно или слабо зависеть от числа Фруда. Учитывая природную стохастичность процесса движения донных волн и наличие области определения решения со слабой зависимостью от чисел Фруда, можно сделать вывод о том, что экспериментальное наблюдение за процессом развития движения донных волн в данной области должно приводить к получению данных, имеющих существенную дисперсию, что и происходит в действительности.

In this paper on the basis of the riverbed model proposed earlier the one-dimensional stability problem of closed flow channel with sandy bed is solved. The feature of the investigated problem is used original equation of riverbed deformations, which takes into account the influence of mechanical and granulometric bed material characteristics and the bed slope when riverbed analyzing. Another feature of the discussed problem is the consideration together with shear stress influence normal stress influence when investigating the riverbed instability. The analytical dependence determined the wave length of fast-growing bed perturbations is obtained from the solution of the sandy bed stability problem for closed flow channel. The analysis of the obtained analytical dependence is performed. It is shown that the obtained dependence generalizes the row of well-known empirical formulas: Coleman, Shulyak and Bagnold. The structure of the obtained analytical dependence denotes the existence of two hydrodynamic regimes characterized by the Froude number, at which the bed perturbations growth can strongly or weakly depend on the Froude number. Considering a natural stochasticity of the waves movement process and the presence of a definition domain of the solution with a weak dependence on the Froude numbers it can be concluded that the experimental observation of the of the bed waves movement development should lead to the data acquisition with a significant dispersion and it occurs in reality.

Исследуются динамические теоретико-игровые модели инновационного развития корпорации. Предлагаемые модели основаны на согласовании частных и общественных интересов агентов. Предполагается, что структура интересов каждого агента включает как частную (личные интересы), так и общественную (интересы компании в целом, в первую очередь отражающие необходимость ее инновационного развития) составляющие. Агенты могут делить персональные ресурсы между этими направлениями. Динамика системы описывается не дифференциальным, а разностным уравнением. При исследовании предложенной модели инновационного развития используются имитация и метод перебора областей допустимых управлений субъектов с некоторым шагом. Основной вклад работы — сравнительный анализ эффективности методов иерархического управления для информационных регламентов Штакельберга/Гермейера при принуждении/побуждении (четыре регламента) с помощью индексов системной согласованности. Предлагаемая модель носит универсальный характер и может быть использована для научно обоснованной поддержки ПИР компаний всех отраслей экономики. Специфика конкретной компании учитывается в ходе идентификации модели (определения конкретных классов ис- пользуемых в модели функций и числовых значений параметров), которая представляет собой отдельную сложную задачу и предполагает анализ системы официальной отчетности компании и применение экспертных оценок ее специалистов. Приняты следующие предположения относительно информационного регламента иерархической игры: все игроки используют программные стратегии; ведущий выбирает и сообщает ведомым экономические управления либо административные управления, которые могут быть только функциями времени (игры Штакельберга) либо зависеть также от управлений ведомых (игры Гермейера); при известных стратегиях ведущего ведомые одновременно и независимо выбирают свои стратегии, что приводит к равновесию Нэша в игре ведомых. За конечное число итераций предложенный алгоритм имитационного моделирования позволяет построить приближенное решение модели или сделать вывод, что равновесия не существует. Достоверность и эффективность предложенного алгоритма следуют из свойств методов сценариев и прямого упорядоченного перебора с постоянным шагом. Получен ряд содержательных выводов относительно сравнительной эффективности методов иерархического управления инновациями.

Dynamic game theoretic models of the corporative innovative development are investigated. The proposed models are based on concordance of private and public interests of agents. It is supposed that the structure of interests of each agent includes both private (personal interests) and public (interests of the whole company connected with its innovative development first) components. The agents allocate their personal resources between these two directions. The system dynamics is described by a difference (not differential) equation. The proposed model of innovative development is studied by simulation and the method of enumeration of the domains of feasible controls with a constant step. The main contribution of the paper consists in comparative analysis of efficiency of the methods of hierarchical control (compulsion or impulsion) for information structures of Stackelberg or Germeier (four structures) by means of the indices of system compatibility. The proposed model is a universal one and can be used for a scientifically grounded support of the programs of innovative development of any economic firm. The features of a specific company are considered in the process of model identification (a determination of the specific classes of model functions and numerical values of its parameters) which forms a separate complex problem and requires an analysis of the statistical data and expert estimations. The following assumptions about information rules of the hierarchical game are accepted: all players use open-loop strategies; the leader chooses and reports to the followers some values of administrative (compulsion) or economic (impulsion) control variables which can be only functions of time (Stackelberg games) or depend also on the followers’ controls (Germeier games); given the leader’s strategies all followers simultaneously and independently choose their strategies that gives a Nash equilibrium in the followers’ game. For a finite number of iterations the proposed algorithm of simulation modeling allows to build an approximate solution of the model or to conclude that it doesn’t exist. A reliability and efficiency of the proposed algorithm follow from the properties of the scenario method and the method of a direct ordered enumeration with a constant step. Some comprehensive conclusions about the comparative efficiency of methods of hierarchical control of innovations are received.

Модернизация методики аэродинамического эксперимента на современном уровне подразумевает создание математических моделей аэродинамических труб (электронных АДТ), предназначенных для вычислительного сопровождения экспериментальных исследований. Применение электронных АДТ в перспективе способно обеспечить получение достоверных аэродинамических характеристик летательных аппаратов по результатам исследования их моделей в аэродинамических трубах, согласования результатов, полученных на разных экспериментальных установках, сравнения расчетов моделей в безграничном потоке с учетом влияния подвесных устройств и границ потока в рабочей части экспериментальной установки.

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

В данной работе представлены результаты расчетов аэродинамических характеристик тематической модели летательного аппарата схемы «летающее крыло» в безграничном потоке при разных углах атаки, полученные в рамках первого этапа работы по созданию математической модели рабочей части аэродинамической трубы Т-102 ЦАГИ. Расчеты выполнялись с использованием двухпараметрической k–ε модели турбулентности со специальными пристеночными функциями, приспособленными для расчета отрывных течений. В рамках данной работы исследовались основные продольные аэродинамические характеристики, было выполнено сравнение с результатами экспериментальных исследований в аэродинамической трубе Т-102 ЦАГИ с учетом погрешностей.

Modern approach to modernization of the experimental techniques involves design of mathematical models of the wind-tunnel, which are also referred to as Electronic of Digital Wind-Tunnels. They are meant to supplement experimental data with computational analysis. Using Electronic Wind-Tunnels is supposed to provide accurate information on aerodynamic performance of an aircraft basing on a set of experimental data, to obtain agreement between data from different test facilities and perform comparison between computational results for flight conditions and data with the presence of support system and test section.

Completing this task requires some preliminary research, which involves extensive wind-tunnel testing as well as RANS-based computational research with the use of supercomputer technologies. At different stages of computational investigation one may have to model not only the aircraft itself but also the wind-tunnel test section and the model support system. Modelling such complex geometries will inevitably result in quite complex vertical and separated flows one will have to simulate. Another problem is that boundary layer transition is often present in wind-tunnel testing due to quite small model scales and therefore low Reynolds numbers.

In the current article the first stage of the Electronic Wind-Tunnel design program is covered. This stage involves computational investigation of aerodynamic characteristics of the generic flying-wing UAV model previously tested in TsAGI T-102 wind-tunnel. Since this stage is preliminary the model was simulated without taking test-section and support system geometry into account. The boundary layer was considered to be fully turbulent.

For the current research FlowVision computational code was used because of its automatic grid generation feature and stability of the solver when simulating complex flows. A two-equation k–ε turbulence model was used with special wall functions designed to properly capture flow separation. Computed lift force and drag force coefficients for different angles-of-attack were compared to the experimental data.