All issues

В рамках проблемы предотвращения астероидно-кометной угрозы выполнен физический и теоретический анализ процессов воздействия различных факторов надповерхностного ядерного взрыва достаточно высокой энергии на астероид во внеатмосферных условиях космического пространства. Показано, что в соответствии с энергией и проницаемой способностью плазмы продуктов взрыва, рентгеновского и гамма-нейтронного излучения на поверхности астероида, обращенной к взрыву, образуется слоистая структура с разной плотностью энергии, зависящей от угловых координат. Для каждого слоя выяснен временной характер трансформации энергии внутри него и определены роли различных фото- и столкновительных процессов. Воздействие высокоскоростного потока плазмы носит эрозионный характер, при этом импульс плазмы передается астероиду. Показано, что в тонком слое поглощения рентгеновского излучения вещество астероида разогревается до высоких температур, и в результате его расширения формируется импульс отдачи, который не является определяющим из-за малой массы расширяющейся высокотемпературной плазмы. Расчеты показали, что основной импульс, полученный астероидом, связан с уносом разогретого слоя вещества, образованного нейтронным потоком (7.5 · 10¹⁴ г · см/с). Показано, что астероид с радиусом ~100 м приобретает при этом скорость ≈ 100 см/с. Расчеты выполнены с учетом затрат энергии взрыва на разрушение аморфной структуры вещества астероида (~1 эВ/атом = 3.8 · 10¹⁰ эрг/г) и на ионизацию в области высокотемпературного слоя. На основе аналогичного анализа получено приближенное выражение для оценки среднего размера осколков при возможном разрушении астероида ударными волнами, образующимися внутри него под действием импульсов давления. Выполнен физический эксперимент в лабораторных условиях, имитирующий фрагментацию каменного астероида и подтвердивший справедливость полученной зависимости от выбранных значений определенных параметров. В результате численных исследований воздействия взрыва, произведенных на различном расстоянии от поверхности астероида, показано, что учет реальной геометрии отколочного слоя дает оптимальную высоту для формирования максимального импульса астероида примерно в 1.5 раза большую, чем аналогичные оценки по упрощенной модели. Предложена двухэтапная концепция воздействия ядерных взрывов на астероид с использованием радиолокационных средств наведения. Проанализировано возможное влияние возникающих ионизационных помех на радиолокационное слежение за разлетом крупных осколков астероида в условиях пространственно-временной эволюции всех элементов исследуемой динамической системы.

As part of the paper, a physical and theoretical analysis of the impact processes of various factors of a highaltitude and high-energy nuclear explosion on the asteroid in extra-atmospheric conditions of open space is done. It is shown that, in accordance with the energy and permeability of the plasma of explosion products, X-ray and gamma-neutron radiation, a layered structure with a different energy density depending on angular coordinates is formed on the surface of the asteroid. The temporal patterns of the energy transformation for each layer is clarified and the roles of various photo- and collision processes are determined. The effect of a high-speed plasma flow is erosive in nature, and the plasma pulse is transmitted to the asteroid. The paper presents that in a thin layer of x-ray absorption, the asteroid substance is heated to high temperatures and as a result of its expansion, a recoil impulse is formed, which is not decisive due to the small mass of the expanding high-temperature plasma. Calculations shows that the main impulse received by an asteroid is associated with the entrainment of a heated layer of a substance formed by a neutron flux (7.5 E 10¹⁴ g E cm/s). It is shown that an asteroid with a radius of ~100 m acquires a velocity of . 100 cm/s. The calculations were performed taking into account the explosion energy spent on the destruction of the amorphous structure of the asteroid material (~1 eV/atom = 3.8 E 10¹⁰ erg/g) and ionization in the region of the high-temperature layer. Based on a similar analysis, an approximation is obtained for estimating the average size of fragments in the event of the possible destruction of the asteroid by shock waves generated inside it under the influence of pressure impulses. A physical experiment was conducted in laboratory conditions, simulating the fragmentation of a stone asteroid and confirming the validity of the obtained dependence on the selected values of certain parameters. As a result of numerical studies of the effects of the explosion, carried out at different distances from the surface of the asteroid, it is shown that taking into account the real geometry of the spallation layer gives the optimal height for the formation of the maximum asteroid momentum by a factor of 1.5 greater than similar estimates according to the simplified model. A two-stage concept of the impact of nuclear explosions on an asteroid using radar guidance tools is proposed. The paper analyzes the possible impact of the emerging ionization interference on the radar tracking of the movement of large fragments of the asteroid in the space-time evolution of all elements of the studied dynamic system.

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

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

The paper provides the mathematical and numerical models of the interrelated thermo- and hydrodynamic processes in the operational mode of development the unified oil-producing complex during the hydrogel flooding of the non-uniform oil reservoir exploited with a system of arbitrarily located injecting wells and producing wells equipped with submersible multistage electrical centrifugal pumps. A special feature of our approach is the modeling of the special ground-based equipment operation (control stations of submersible pumps, drossel devices on the head of producing wells), designed to regulate the operation modes of both the whole complex and its individual elements.

The complete differential model includes equations governing non-stationary two-phase five-component filtration in the reservoir, quasi-stationary heat and mass transfer in the wells and working channels of pumps. Special non-linear boundary conditions and dependencies simulate, respectively, the influence of the drossel diameter on the flow rate and pressure at the wellhead of each producing well and the frequency electric current on the performance characteristics of the submersible pump unit. Oil field development is also regulated by the change in bottom-hole pressure of each injection well, concentration of the gel-forming components pumping into the reservoir, their total volume and duration of injection. The problem is solved numerically using conservative difference schemes constructed on the base of the finite difference method, and developed iterative algorithms oriented on the parallel computing technologies. Numerical model is implemented in a software package which can be considered as the «Intellectual System of Wells» for the virtual control the oil field development.

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

The article focuses on the problem of mechanisms’ co-design for robotic systems to perform adaptive physical interaction with an unstructured environment, including physical human robot interaction. The co-design means simultaneous optimization of mechanics and control system, ensuring optimal behavior and performance of the system. Mechanics optimization refers to the search for optimal structure, geometric parameters, mass distribution among the links and their compliance; control refers to the search for motion trajectories for mechanism’s joints. The paper presents a generalized method of structural-parametric synthesis of underactuated mechanisms with closed kinematics for robotic systems for various purposes, e. g., it was previously used for the co-design of fingers’ mechanisms for anthropomorphic gripper and legs’ mechanisms for galloping robots. The method implements the concept of morphological computation of control laws due to the features of mechanical design, minimizing the control effort from the algorithmic component of the control system, which reduces the requirements for the level of technical equipment and reduces energy consumption. In this paper, the proposed method is used to optimize the structure and geometric parameters of the passive mechanism of the back support module of an industrial exosuit. Human movements are diverse and non-deterministic when compared with the movements of autonomous robots, which complicates the design of wearable robotic devices. To reduce injuries, fatigue and increase the productivity of workers, the synthesized industrial exosuit should not only compensate for loads, but also not interfere with the natural human motions. To test the developed exosuit, kinematic datasets from motion capture of an entire human body during industrial operations were used. The proposed method of structural-parametric synthesis was used to improve the ergonomics of a wearable robotic device. Verification of the synthesized mechanism was carried out using simulation: the passive module of the back is attached to two geometric primitives that move the chest and pelvis of the exosuit operator in accordance with the motion capture data. The ergonomics of the back module is quantified by the distance between the joints connecting the upper and bottom parts of the exosuit; minimizing deviation from the average value corresponds to a lesser limitation of the operator’s movement, i. e. greater ergonomics. The article provides a detailed description of the method of structural-parametric synthesis, an example of synthesis of an exosuit module and the results of simulation.

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

The article presented the computational algorithm developed to study chemical processes in the internal flows of a multicomponent gas under the influence of laser radiation. The mathematical model is the gas dynamics’ equations with chemical reactions at low Mach numbers. It takes into account dissipative terms that describe the dynamics of a viscous heat-conducting medium with diffusion, chemical reactions and energy supply by laser radiation. This mathematical model is characterized by the presence of several very different time and spatial scales. The computational algorithm is based on a splitting scheme by physical processes. Each time integration step is divided into the following blocks: solving the equations of chemical kinetics, solving the equation for the radiation intensity, solving the convection-diffusion equations, calculating the dynamic component of pressure and calculating the correction of the velocity vector. The solution of a stiff system of chemical kinetics equations is carried out using a specialized explicit second-order accuracy scheme or a plug-in RADAU5 module. Numerical Rusanov flows and a WENO scheme of an increased order of approximation are used to find convective terms in the equations. The code based on the obtained algorithm has been developed using MPI parallel computing technology. The developed code is used to calculate the pyrolysis of ethane with radical reactions. The superequilibrium concentrations’ formation of radicals in the reactor volume is studied in detail. Numerical simulation of the reaction gas flow in a flat tube with laser radiation supply is carried out, which is in demand for the interpretation of experimental results. It is shown that laser radiation significantly increases the conversion of ethane and yields of target products at short lengths closer to the entrance to the reaction zone. Reducing the effective length of the reaction zone allows us to offer new solutions in the design of ethane conversion reactors into valuable hydrocarbons. The developed algorithm and program will find their application in the creation of new technologies of laser thermochemistry.

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

The work is devoted to the problem of creating a model with stationary parameters using historical data under conditions of unknown disturbances. The case is considered when a representative sample of object states can be formed using historical data accumulated only over a significant period of time. It is assumed that unknown disturbances can act in a wide frequency range and may have low-frequency and trend components. In such a situation, including data from different time periods in the sample can lead to inconsistencies and greatly reduce the accuracy of the model. The paper provides an overview of approaches and methods for data harmonization. In this case, the main attention is paid to data sampling. An assessment is made of the applicability of various data sampling options as a tool for reducing the level of uncertainty. We propose a method for identifying a self-leveling object model using data accumulated over a significant period of time under conditions of unknown disturbances with a wide frequency range. The method is focused on creating a model with stationary parameters that does not require periodic reconfiguration to new conditions. The method is based on the combined use of sampling and presentation of data from individual periods of time in the form of increments relative to the initial point in time for the period. This makes it possible to reduce the number of parameters that characterize unknown disturbances with a minimum of assumptions that limit the application of the method. As a result, the dimensionality of the search problem is reduced and the computational costs associated with setting up the model are minimized. It is possible to configure both linear and, in some cases, nonlinear models. The method was used to develop a model of closed cooling of steel on a unit for continuous hot-dip galvanizing of steel strip. The model can be used for predictive control of thermal processes and for selecting strip speed. It is shown that the method makes it possible to develop a model of thermal processes from a closed cooling section under conditions of unknown disturbances, including low-frequency components.

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

This paper describes the design of an aquatic robot moving on the surface of a fluid and driven by two internal moving masses. The body of the aquatic robot in cross section has the shape of a symmetrical airfoil with a sharp edge. In this prototype, two internal masses move in circles and are rotated by a single DC motor and a gear mechanism that transmits torque from the motor to each mass. Angular velocities of moving masses are used as a control action, and the developed kinematic scheme for transmitting rotation from the motor to the moving masses allows the rotation of two masses with equal angular velocities in magnitude, but with a different direction of rotation. It is also possible to install additional tail fins of various shapes and sizes on the body of this robot. Also in the work for this object, the equations of motion are presented, written in the form of Kirchhoff equations for the motion of a solid body in an ideal fluid, which are supplemented by terms of viscous resistance. A mathematical description of the additional forces acting on the flexible tail fin is presented. Experimental studies on the influence of various tail fins on the speed of motion in the fluid were carried out with the developed prototype of the robot. In this work, tail fins of the same shape and size were installed on the robot, while having different stiffness. The experiments were carried out in a pool with water, over which a camera was installed, on which video recordings of all the experiments were obtained. Next processing of the video recordings made it possible to obtain the object’s movements coordinates, as well as its linear and angular velocities. The paper shows the difference in the velocities developed by the robot when moving without a tail fin, as well as with tail fins having different stiffness. The comparison of the velocities developed by the robot, obtained in experimental studies, with the results of mathematical modeling of the system is given.

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

A mathematical model of predator – prey microecosystem with lower critical population number of prey is considered. The predator – prey system is assumed to be under harvesting. Harvesting intensity variations generate changes in two model parameters which are considered as controllable. Bifurcation diagram in control-lable parameters plane is constructed and corresponding phase portraits are represented.

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

The paper presents results that confirm the ability to control the movement of the kink of the modified sine-Gordon equation with variable external force parameters. Three types of external influences have been considered: permanent action, periodic action with a constant frequency and a frequency-modulated periodic exposure. The dependences of the position and velocity of the kink on time for various values of the parameters of external influence were obtained using the method of McLaughlin and Scott. It is shown that by changing the settings, one can adjust the velocity and direction of movement of the kink.

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

A method and a complex of computer programs are developed for the numerical simulation of the polaron states excitation process in condensed media. A numerical study of the polaron states formation in water under the action of the ultraviolet range laser irradiation is carried out. Our approach allows to reproduce the experimental data of the hydrated electrons formation. A numerical scheme is presented for the solution of the respective system of nonlinear partial differential equations. Parallel implementation is based on the MPI technique. The numerical results are given in comparison with the experimental data and theoretical estimations.

Волны возбуждения являются прообразом самоорганизующихся динамических структур в неравновесных системах. Они характеризуются своей собственной внутренней динамикой, приводящей к формированию бегущих волн различных типов и форм. Яркие примеры — это вращающиеся спирали и скрученные свитки. Интересная и сложная задача — найти способы управления их поведением, применяя внешние сигналы, влияющие на распространяющиеся волны. В качестве такого воздействия мы используем внешние электрические поля, наложенные на возбудимую реакцию Белоусова–Жаботинского (БЖ). Существенные эффекты влияния полей на волны включают изменение скорости волны, обращение направления распространения, взаимное уничтожение вращающихся в противоположных направлениях спиральных волн и переориентацию нитей скрученных свитков. Эти эффекты могут быть объяснены в численных экспериментах, при этом существенную роль играет отрицательно заряженный ингибиторбромид. Эффекты электрического поля также были исследованы в биологических возбудимых средах, таких как социальные амебы Dictyostelium discoideum. Совсем недавно мы начали исследовать влияние электрического поля на реакцию БЖ, протекающую в водно-масляной микроэмульсии. Удалось наблюдать дрейф сложных структур, а также изменение вязкости и электрической проводимости. Мы обсуждаем предположение, что эта система может выступать в качестве модели для дальнодействующего взаимодействия между нейронами.

Excitation waves are a prototype of self-organized dynamic patterns in non-equilibrium systems. They develop their own intrinsic dynamics resulting in travelling waves of various forms and shapes. Prominent examples are rotating spirals and scroll waves. It is an interesting and challenging task to find ways to control their behavior by applying external signals, upon which these propagating waves react. We apply external electric fields to such waves in the excitable Belousov–Zhabotinsky (BZ) reaction. Remarkable effects include the change of wave speed, reversal of propagation direction, annihilation of counter-rotating spiral waves and reorientation of scroll wave filaments. These effects can be explained in numerical simulations, where the negatively charged inhibitor bromide plays an essential role. Electric field effects have also been investigated in biological excitable media such as the social amoebae Dictyostelium discoideum. Quite recently we have started to investigate electric field effect in the BZ reaction dissolved in an Aerosol OT water-in-oil microemulsion. A drift of complex patterns can be observed, and also the viscosity and electric conductivity change. We discuss the assumption that this system can act as a model for long range communication between neurons.