Все выпуски

Разработан метод и комплекс программ для численного моделирования процесса формирования поляронных состояний в конденсированных средах. Проведено численное исследование этого процесса для водной среды при воздействии лазерного облучения в ультрафиолетовом диапазоне. Показано, что в рамках предложенного подхода удается численно воспроизвести экспериментальные данные по формированию гидратированных электронов. Представлена схема численного решения системы нелинейных дифференциальных уравнений в частных производных, описывающих динамическую модельпо лярона. Программная реализация выполнена с использованием технологии параллельного программирования 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.

Построена статическая трехуровневая теоретико-игровая модель системы управления судовыми балластными водами. Используются методы иерархического управления при одновременном учете условий поддержания системы в заданном состоянии. Проводится сравнение результатов исследования модели с точки зрения игр Гермейера $\Gamma_1$ и $\Gamma_2$. Приведены примеры численных расчетов в ряде характерных случаев.

The static three-level game-theoretic model of three-level control system of the ship’s water ballast is built. The methods of hierarchical control in view of requirements of keeping the system in the given state are used. A comparison of the results of study of the model in terms of $\Gamma_1$ and $\Gamma_2$ Germeyer’s games is conducted. Numerical calculations for some typical cases are given.

В различных областях науки, техники, защиты окружающей среды, в строительстве актуальными являются вопросы изучения процессов взаимодействия пористых материалов с веществами, находящимися в различных агрегатных состояниях. Особенно актуальными с точки зрения экологии и защиты окружающей среды являются исследования процессов взаимодействия пористых материалов с водой в жидкой и газообразной фазе. Поскольку в одном моле воды содержится 6.022140857 · 10²³ молекул H₂O, для описания свойств, например, водяного пара в поре в основном используются макроскопические подходы, рассматривающие водяной пар как сплошную среду в рамках аэродинамики. В данной работе построена и использовалась для моделирования макроскопическая двумерная диффузионная модель [Bitsadze, Kalinichenko, 1980] поведения водяного пара внутри изолированной поры. Наряду с макроскопической моделью в работе предложена микроскопическая модель поведения водяного пара внутри изолированной поры, построенная в рамках молекулярно-динамического подхода [Gould et al., 2005]. В данной модели на основе классической механики Ньютона описывается движение каждой молекулы воды, взаимодействующей как с другими молекулами воды, так и со стенками поры. Рассматривается эволюция системы «водяной пар – пора» с течением времени. В зависимости от внешних по отношению к поре условий система эволюционирует к различным состояниям равновесия, которые характеризуются различными значениями макроскопических характеристик, таких как температура, плотность, давление. Сравнение результатов молекулярно-динамического моделирования с результатами вычислений на основе макроскопической диффузионной модели и экспериментальными данными позволяет сделать вывод о необходимости сочетания макроскопического и микроскопического подхода для адекватного и более точного описания процессов взаимодействия водяного пара с пористыми материалами.

In various areas of science, technology, environment protection, construction, it is very important to study processes of porous materials interaction with different substances in different aggregation states. From the point of view of ecology and environmental protection it is particularly actual to investigate processes of porous materials interaction with water in liquid and gaseous phases. Since one mole of water contains 6.022140857 · 10²³ molecules of H₂O, macroscopic approaches considering the water vapor as continuum media in the framework of classical aerodynamics are mainly used to describe properties, for example properties of water vapor in the pore. In this paper we construct and use for simulation the macroscopic two-dimensional diffusion model [Bitsadze, Kalinichenko, 1980] describing the behavior of water vapor inside the isolated pore. Together with the macroscopic model it is proposed microscopic model of the behavior of water vapor inside the isolated pores. This microscopic model is built within the molecular dynamics approach [Gould et al., 2005]. In the microscopic model a description of each water molecule motion is based on Newton classical mechanics considering interactions with other molecules and pore walls. Time evolution of “water vapor – pore” system is explored. Depending on the external to the pore conditions the system evolves to various states of equilibrium, characterized by different values of the macroscopic characteristics such as temperature, density, pressure. Comparisons of results of molecular dynamic simulations with the results of calculations based on the macroscopic diffusion model and experimental data allow to conclude that the combination of macroscopic and microscopic approach could produce more adequate and more accurate description of processes of water vapor interaction with porous materials.

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

Для осуществления полноволнового моделирования в данной работе был применен сеточно-характеристический метод на структурированных сетках и решалась гиперболическая система уравнений, описывающая распространение упругих волн в материале рассматриваемой пластины конечной толщины на конкретном примере отношения толщины к ширине 1:10.

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

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

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

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

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

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

Ultrasound examination of material properties is a precision method for determining their elastic and strength properties in connection with the small wavelength formed in the material after impact of a laser beam. In this paper, the wave processes arising during these measurements are considered in detail. It is shown that full-wave numerical modeling allows us to study in detail the types of waves, topological characteristics of their profile, speed of arrival of waves at various points, identification the types of waves whose measurements are most optimal for examining a sample made of a specific material of a particular shape, and to develop measurement procedures.

To carry out full-wave modeling, a grid-characteristic method on structured grids was used in this work and a hyperbolic system of equations that describes the propagation of elastic waves in the material of the thin plate under consideration on a specific example of a ratio of thickness to width of 1:10 was solved.

To simulate an elastic front that arose in the plate due to a laser beam, a model of the corresponding initial conditions was proposed. A comparison of the wave effects that arise during its use in the case of a point source and with the data of physical experiments on the propagation of laser ultrasound in metal plates was made.

A study was made on the basis of which the characteristic topological features of the wave processes under consideration were identified and revealed. The main types of elastic waves arising due to a laser beam are investigated, the possibility of their use for studying the properties of materials is analyzed. A method based on the analysis of multiple waves is proposed. The proposed method for studying the properties of a plate with the help of multiple waves on synthetic data was tested, and it showed good results.

It should be noted that most of the studies of multiple waves are aimed at developing methods for their suppression. Multiple waves are not used to process the results of ultrasound studies due to the complexity of their detection in the recorded data of a physical experiment.

Due to the use of full wave modeling and analysis of spatial dynamic wave processes, multiple waves are considered in detail in this work and it is proposed to divide materials into three classes, which allows using multiple waves to obtain information about the material of the plate.

The main results of the work are the developed problem statements for the numerical simulation of the study of plates of a finite thickness by laser ultrasound; the revealed features of the wave phenomena arising in plates of a finite thickness; the developed method for studying the properties of the plate on the basis of multiple waves; the developed classification of materials.

The results of the studies presented in this paper may be of interest not only for developments in the field of ultrasonic non-destructive testing, but also in the field of seismic exploration of the earth's interior, since the proposed approach can be extended to more complex cases of heterogeneous media and applied in geophysics.

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

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

We present a physico-mathematical statement of coupled geometrical and gas dynamics problem of intrachamber processes simulation and calculation of main internal ballistics characteristics of solid rocket motors in axisymmetric approximation. Method and numerical algorithm of solving the problem are described in this paper. We track the propellant burning surface using the level set method. This method allows us to implicitly represent the surface on a fixed Cartesian grid as zero-level of some function. Two-dimensional gas-dynamics equations describe a flow of combustion products in a solid rocket motor. Due to inconsistency of domain boundaries and nodes of computational grid, presence of ghost points lying outside the computational domain is taken into account. For setting the values of flow parameters in ghost points, we use the inverse Lax – Wendroff procedure. We discretize spatial derivatives of level set and gas-dynamics equations with standard WENO schemes of fifth and third-order respectively and time derivatives using total variation diminishing Runge –Kutta methods. We parallelize the presented numerical algorithm using CUDA technology and further optimize it with regard to peculiarities of graphics processors architecture.

Created software package is used for calculating internal ballistics characteristics of nozzleless solid rocket motor during main firing phase. On the base of obtained numerical results, we discuss efficiency of parallelization using CUDA technology and applying considered optimizations. It has been shown that implemented parallelization technique leads to a significant acceleration in comparison with central processes. Distributions of key parameters of combustion products flow in different periods of time have been presented in this paper. We make a comparison of obtained results between quasione-dimensional approach and developed numerical technique.

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

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

The quasistatic deformation problem for shape memory alloys is reviewed within the phenomenological mechanics of solids without microphysics analysis. The phenomenological approach is based on comparison of two material deformation diagrams. The first diagram corresponds to the active proportional loading when the alloy behaves as an ideal elastoplastic material; the residual strain is observed after unloading. The second diagram is relevant to the case when the deformed sample is heated to a certain temperature for each alloy. The initial shape is restored: the reverse distortion matches deformations on the first diagram, except for the sign. Because the first step of distortion can be described with the variational principle, for which the existence of the generalized solutions is proved under arbitrary loading, it becomes clear how to explain the reverse distortion within the slightly modified theory of plasticity. The simply connected surface of loading needs to be replaced with the doubly connected one, and the variational principle needs to be updated with two laws of thermodynamics and the principle of orthogonality for thermodynamic forces and streams. In this case it is not difficult to prove the existence of solutions either. The successful application of the theory of plasticity under the constant temperature causes the need to obtain a similar result for a more general case of variable external forces and temperatures. The paper studies the ideal elastoplastic von Mises model at linear strain rates. Taking into account hardening and arbitrary loading surface does not cause any additional difficulties.

The extended variational principle of the Reissner type is defined. Together with the laws of thermal plasticity it enables to prove the existence of the generalized solutions for three-dimensional bodies made of shape memory materials. The main issue to resolve is a challenge to choose a functional space for the rates and deformations of the continuum points. The space of bounded deformation, which is the main instrument of the mathematical theory of plasticity, serves this purpose in the paper. The proving process shows that the choice of the functional spaces used in the paper is not the only one. The study of other possible problem settings for the extended variational principle and search for regularity of generalized solutions seem an interesting challenge for future research.

Анализируются варианты учета неоднородности среды при компьютерном моделировании динамики хищника и жертвы на основе системы уравнений реакции–диффузии–адвекции. Локальное взаимодействие видов (члены реакции) описывается логистическим законом роста для жертвы и соотношениями Беддингтона – ДеАнгелиса, частными случаями которых являются функциональный отклик Холлинга второго рода и модель Ардити – Гинзбурга. Рассматривается одномерная по пространству задача для неоднородного ресурса (емкости среды) и трех видов таксиса (жертвы на ресурс и от хищника, хищника к жертве). Используется аналитический подход для исследования устойчивости стационарных решений в случае локального взаимодействия (бездиффузионный подход) и вычисления на основе метода прямых для учета диффузионных и адвективных процессов. Сравнение критических значений параметра смертности хищников показало, что при постоянных коэффициентах в соотношениях Беддингтона – ДеАнгелиса получаются переменные по пространственной координате критические величины, а для модели Ардити – Гинзбурга данный эффект не наблюдается. Предложена модификация членов реакции, позволяющая учесть неоднородность ресурса. Представлены численные результаты по динамике видов для больших и малых миграционных коэффициентов, демонстрирующие снижение влияния вида локальных членов на формирующиеся пространственно-временные распределения популяций. Проанализированы бифуркационные переходы при изменении параметров диффузии–адвекции и членов реакции.

We analyze variants of considering the inhomogeneity of the environment in computer modeling of the dynamics of a predator and prey based on a system of reaction-diffusion–advection equations. The local interaction of species (reaction terms) is described by the logistic law for the prey and the Beddington –DeAngelis functional response, special cases of which are the Holling type II functional response and the Arditi – Ginzburg model. We consider a one-dimensional problem in space for a heterogeneous resource (carrying capacity) and three types of taxis (the prey to resource and from the predator, the predator to the prey). An analytical approach is used to study the stability of stationary solutions in the case of local interaction (diffusionless approach). We employ the method of lines to study diffusion and advective processes. A comparison of the critical values of the mortality parameter of predators is given. Analysis showed that at constant coefficients in the Beddington –DeAngelis model, critical values are variable along the spatial coordinate, while we do not observe this effect for the Arditi –Ginzburg model. We propose a modification of the reaction terms, which makes it possible to take into account the heterogeneity of the resource. Numerical results on the dynamics of species for large and small migration coefficients are presented, demonstrating a decrease in the influence of the species of local members on the emerging spatio-temporal distributions of populations. Bifurcation transitions are analyzed when changing the parameters of diffusion–advection and reaction terms.

При распознавании рабочих на изображениях строительной площадки, получаемых с камер наблюдения, типичной является ситуация, при которой объекты детекции имеют сильно различающийся пространственный масштаб относительно друг друга и других объектов. Повышение точности детекции мелких объектов может быть обеспечено путем использования Feature-Fused модификации детектора SSD (Single Shot Detector). Вместе с применением на инференсе нарезки изображения с перекрытием такая модель хорошо справляется с детекцией мелких объектов. Однако при практическом использовании данного подхода требуется ручная настройка параметров нарезки. При этом снижается точность детекции объектов на сценах, отличающихся от сцен, использованных при обучении, а также крупных объектов. В данной работе предложен алгоритм автоматического выбора оптимальных параметров нарезки изображения в зависимости от соотношений характерных геометрических размеров объектов на изображении. Нами разработан двухпроходной вариант детектора Feature-Fused SSD для автоматического определения параметров нарезки изображения. На первом проходе применяется усеченная версия детектора, позволяющая определять характерные размеры объектов интереса. На втором проходе осуществляется финальная детекция объектов с параметрами нарезки, выбранными после первого прохода. Был собран датасет с изображениями рабочих на строительной площадке. Датасет включает крупные, мелкие и разноплановые изображения рабочих. Для сравнения результатов детекции для однопроходного алгоритма без разбиения входного изображения, однопроходного алгоритма с равномерным разбиением и двухпроходного алгоритма с подбором оптимального разбиения рассматривались тесты по детекции отдельно крупных объектов, очень мелких объектов, с высокой плотностью объектов как на переднем, так и на заднем плане, только на заднем плане. В диапазоне рассмотренных нами случаев наш подход превосходит подходы, взятые в сравнение, позволяет хорошо бороться с проблемой двойных детекций и демонстрирует качество 0,82–0,91 по метрике mAP (mean Average Precision).

When recognizing workers on images of a construction site obtained from surveillance cameras, a situation is typical in which the objects of detection have a very different spatial scale relative to each other and other objects. An increase in the accuracy of detection of small objects can be achieved by using the Feature-Fused modification of the SSD detector. Together with the use of overlapping image slicing on the inference, this model copes well with the detection of small objects. However, the practical use of this approach requires manual adjustment of the slicing parameters. This reduces the accuracy of object detection on scenes that differ from the scenes used in training, as well as large objects. In this paper, we propose an algorithm for automatic selection of image slicing parameters depending on the ratio of the characteristic geometric dimensions of objects in the image. We have developed a two-pass version of the Feature-Fused SSD detector for automatic determination of optimal image slicing parameters. On the first pass, a fast truncated version of the detector is used, which makes it possible to determine the characteristic sizes of objects of interest. On the second pass, the final detection of objects with slicing parameters selected after the first pass is performed. A dataset was collected with images of workers on a construction site. The dataset includes large, small and diverse images of workers. To compare the detection results for a one-pass algorithm without splitting the input image, a one-pass algorithm with uniform splitting, and a two-pass algorithm with the selection of the optimal splitting, we considered tests for the detection of separately large objects, very small objects, with a high density of objects both in the foreground and in the background, only in the background. In the range of cases we have considered, our approach is superior to the approaches taken in comparison, allows us to deal well with the problem of double detections and demonstrates a quality of 0.82–0.91 according to the mAP (mean Average Precision) metric.

Дискретизация задач по методу гидродинамики сглаженных частиц (SPH) предполагает присутствие в решении нескольких констант — параметров дискретизации. Среди них особо следует отметить модельную скорость звука $c_0$, которая связывает мгновенную плотность в SPH-частице с возникающим давлением через замыкающее уравнение состояния.

В работе изложен подход к точному определению необходимого значения модельной скорости звука, имеющий в своей основе анализ изменения плотностей в SPH-частицах при их относительном смещении. Примером движения сплошной среды принята задача о плоском сдвиговом течении; объектом анализа является функция относительного уплотнения $\varepsilon_\rho$ в SPH-частице, определяемая формой ядра сглаживания. Идеальный плоскопараллельный относительный сдвиг частиц в области сглаживания определяет периодическое изменение их плотностей. Исследование функций $\varepsilon_\rho$, получаемых от использования различных ядер сглаживания в аппроксимации плотности с учетом такого сдвига, позволило установить пульсационный характер возникновения давлений в частицах. Кроме того, определен случай расположения соседей в области сглаживания, обеспечивающий максимум уплотнения в частице.

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

The problem discrete statement by the smoothed particle hydrodynamics method (SPH) include a discretization constants parameters set. Of them particular note is the model sound speed $c_0$, which relates the SPH-particle instantaneous density to the resulting pressure through the equation of state.

The paper describes an approach to the exact determination of the model sound speed required value. It is on the analysis based, how SPH-particle density changes with their relative shift. An example of the continuous medium motion taken the plane shear flow problem; the analysis object is the relative compaction function $\varepsilon_\rho$ in the SPH-particle. For various smoothing kernels was research the functions of $\varepsilon_\rho$, that allowed the pulsating nature of the pressures occurrence in particles to establish. Also the neighbors uniform distribution in the smoothing domain was determined, at which shaping the maximum of compaction in the particle.

Through comparison the function $\varepsilon_\rho$ with the SPH-approximation of motion equation is defined associate the discretization parameter $c_0$ with the smoothing kernel shape and other problem parameters. As a result, an equation is formulated that the necessary and sufficient model sound speed value provides finding. For such equation the expressions of root $c_0$ are given for three different smoothing kernels, that simplified from polynomials to numerical coefficients for the plane shear flow problem parameters.

Предложена модель расчета межслоевой магнитостатической энергии и оценка ее вклада в общую магнитную энергию многослойных пленок. Данная модель была применена для расчета трехслойной структуры 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.