Все выпуски

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

Данная работа посвящена систематическому исследованию равновесных конфигураций кластеров золота и серебра размером от 2 до 9 атомов. Для найденных самых низких по энергии конфигураций проводится расчет электронных спектров возбуждения. Все расчеты проводятся в рамках теории функционала плотности (ТФП) с использованием гибридного функционала B3LYP и псевдопотенциала LANL2DZ. На основании анализа полученных электронных спектров возбуждения делается попытка подтвердить предположение, что небольшие кластеры серебра и золота, способные к люминесценции, могут образовываться только в присутствии стабилизирующей матрицы. Авторы показывают, что без стабилизации в растворе образуются «плоские» и «сферические» конфигурации, спектры которых не соответствуют экспериментальным спектрам возбуждения люминесценции для комплексов «ДНК-кластер».

Present work is dedicated to systematic study of equilibrium configurations of gold and silver clusters in size from 2 to 9 atoms. To the lowest energy configurations is performed calculation of electronic excitation spectra. All calculations are carried out within the framework of density functional theory (DFT) using the hybrid B3LYP functional and ECP LANL2DZ. Based on the analysis of the received electronic excitation spectra we try to confirm the assumption that small gold and silver clusters, capable of luminescence can be produced only in the presence of a stabilizing template. The authors show that without stabilization in the solution formed “flat” and “spherical” configuration whose spectra do not correspond to the experimental luminescence excitation spectra for DNA cluster complexes.

Предлагается нелинейная математическая модель динамики численности занятого населения разных возрастных групп с учетом их взаимодействий, которые рассматриваются по аналогии с популяционными взаимодействиями (конкуренция, дискриминация, помощь, угнетение и т. п.). Под взаимодействиями понимаются такие обобщенные социально-экономические механизмы, которые вызывают взаимосвязанные изменения численности занятых различных возрастных групп. Рассматриваются три возрастные группы занятого населения: молодые специалисты (15–29 лет), с опытом работы (30–49 лет), работники предпенсионного и пенсионного возраста (50 и старше). На основе статистических данных выполнена оценка параметров предложенной модели для южных регионов Дальневосточного федерального округа (ДФО). Анализ модели и модельных сценариев позволяет заключить, что наблюдаемые колебания численности разновозрастных работников на фоне стабильной общей численности занятого населения могут быть следствием сложных взаимодействий этих групп между собой. Вычислительные эксперименты, проведенные при полученных значениях параметров, позволили рассчитать темпы снижения численности и старения занятого населения, а также определить характер взаимодействий между возрастными группами занятых, прямо не отраженный в статистических данных. Установлено, что в целом по ДФО занятые 50 лет и старше находятся с работающей молодежью до 29 лет в отношениях дискриминации, занятые до 29 лет и 30–49 лет — в отношениях партнерства. Наиболее развитые регионы (Приморский край и Хабаровский край) демонстрируют «равномерную» конкуренцию среди разных возрастных групп занятого населения. Для Приморского края удалось выявить эффект перемешивания сценариев динамики, что характерно для систем, находящихся в состоянии структурной перестройки. Этот эффект выражается в том, что при значительном уменьшении миграционного притока занятых 30–49 лет будут формироваться длинные циклы занятости. Кроме того, изменение миграции сопровождается сменой типа взаимодействия — с дискриминации старшего поколения средним на дискриминацию среднего возраста старшим. Для менее развитых регионов Дальнего Востока (Амурская, Магаданская и Еврейская автономная области) характерны более низкие значения миграционного сальдо почти всех возрастов, а также дискриминация со стороны занятой молодежи до 29 лет других возрастных групп и дискриминация занятыми 30–49 лет старшего поколения.

The article focuses on a nonlinear mathematical model that describes the interaction of the different age groups of the employed population. The interactions are treated by analogy with population relationship (competition, discrimination, assistance, oppression, etc). Under interaction of peoples we mean the generalized social and economic mechanisms that cause related changes in the number of employees of different age groups. Three age groups of the employed population are considered. It is young specialists (15–29 years), workers with experience (30–49 years), the employees of pre-retirement and retirement age (50 and older). The estimation of model’s parameters for the southern regions of the Far Eastern Federal District (FEFD) is executed by statistical data. Analysis of model scenarios allows us to conclude the observed number fluctuations of the different ages employees on the background of a stable total employed population may be a consequence of complex interactions between these groups of peoples. Computational experiments with the obtained values of the parameters allowed us to calculate the rate of decline and the aging of the working population and to determine the nature of the interaction between the age groups of employees that are not directly as reflected in the statistics. It was found that in FEFD the employed of 50 years and older are discriminated against by the young workers under 29, employed up to 29 and 30–49 years are in a partnership. It is shown in most developed regions (Primorsky and Khabarovsk Krai) there is “uniform” competition among different age groups of the employed population. For Primorsky Krai we were able to identify the mixing effect dynamics. It is a typical situation for systems in a state of structural adjustment. This effect is reflected in the fact the long cycles of employed population form with a significant decrease in migration inflows of employees 30–49 years. Besides, the change of migration is accompanied by a change of interaction type — from employment discrimination by the oldest of middle generation to discrimination by the middle of older generation. In less developed regions (Amur, Magadan and Jewish Autonomous Regions) there are lower values of migration balance of almost all age groups and discrimination by young workers up 29 years of other age groups and employment discrimination 30–49 years of the older generation.

Оптимизация судовых обводов для снижения сопротивления движению является актуальной задачей гидродинамики судна. Однако вопросы проектирования и совершенствования обводов в современной практике все еще слабо обобщены и формализованы. Они решаются с помощью комбинации научных знаний, инженерного опыта и критериев из области искусства. Практическое использование экспериментального и численного моделирования в задачах оптимизации формы корпуса обычно сводится к методу проб и ошибок. В статье представлен новый метод оптимизации обводов, предназначенный для детального совершенствования формы корпуса, концепция которого использует теоретические закономерности формирования волновой системы судна. Метод предусматривает систематическое варьирование продольного распределения полноты корпуса при фиксации или контроле ее вертикального распределения. Как известно, вертикальное распределение водоизмещения не имеет оптимума по волновому сопротивлению, которое является основным активным компонентом, особенно в отношении формы носовой части. Варьирование продольного распределения водоизмещения предусмотрено путем задания конечных приращений водоизмещения на строевой по шпангоутам, которые затем переносятся на теоретический чертеж с помощью специальных методов трансформации шпангоутов и реализуются в 3D-моделях корпуса. Для оценки влияния модификаций геометрии на сопротивление используется численное моделирование буксировки полученных моделей. Дальнейшие оптимизационные процедуры базируются на выдвинутой гипотезе о независимости влияния различных участков корпуса, выделенных по длине, на буксировочное сопротивление. В результате применения метода к форме корпуса хорошо известного судна KCS, рекомендованного конференцией «Гетеборг-2000» в качестве эталонного объекта для тестирования численных методов, получены оптимальное продольное распределение полноты и соответствующие обводы корпуса, которые позволили снизить его сопротивление на 8.9 %. Оптимизация выполнена на базе результатов по шести моделям с вариациями формы, которые обусловили колебания полного сопротивления корпуса разного знака, величиной 1.3–6.5 %. Визуализация волновых систем показала, что при снижении сопротивления происходит заметное ослабление поперечных волн и усиление расходящихся.

Optimization of hull lines for the minimum resistance to movement is a problem of current interest in ship hydrodynamics. In practice, lines design is still to some extent an art. The usual approaches to decrease the ship resistance are based on the model experiment and/or CFD simulation, following the trial and error method. The paper presents a new method of in-detail hull form design based on the wave-based optimization approach. The method provides systematic variation of the hull geometrical form, which corresponds to alteration of longitudinal distribution of the hull volume, while its vertical volume distribution is fixed or highly controlled. It’s well known from the theoretical studies that the vertical distribution can't be optimized by condition of minimum wave resistance, thus it can be neglected for the optimization procedures. The method efficiency was investigated by application to the foreship of KCS, the well-known test object from the workshop Gothenburg-2000. The variations of the longitudinal distribution of the volume were set on the sectional area curve as finite volume increments and then transferred to the lines plan with the help of special frame transformation methods. The CFD towing simulations were carried out for the initial hull form and the six modified variants. According to the simulation results, examined modifications caused the resistance increments in the range 1.3–6.5 %. Optimization process was underpinned with the respective data analysis based on the new hypothesis, according to which, the resistance increments caused by separate longitudinal segments of hull form meet the principle of superposition. The achieved results, which are presented as the optimum distribution of volume present in the optimized designed hull form, which shows the interesting characteristics that its resistance has decrease by 8.9 % in respect to initial KCS hull form. Visualization of the wave patterns showed an attenuation of the transversal wave components, and the intensification of the diverging wave components.

В статье проведен анализ исторического процесса с использованием методов синергетики (науки о нелинейных развивающихся системах в природе и обществе), развитых в работах Д. С. Чернавского применительно к экономическим и социальным системам. Показано, что социальная самоорганизация в зависимости от условий приводит к формированию как обществ с сильной внутренней конкуренцией (Y-структуры), так и обществ кооперативного типа (Х-структуры). Y-структуры характерны для стран Запада, Х-структуры характерны для стран Востока. Показано, что в XIX и XX веках имело место ускоренное формирование и усиление Y-структур. Однако в настоящее время мировая система вошла в период серьезных структурных перемен в экономической, политической, идеологической сферах: доминирование Y-структур заканчивается. Рассмотрены возможные пути дальнейшего развития мировой системы, связанные с изменением режимов самоорганизации и ограничением внутренней конкуренции. Этот переход будет длительным и сложным. В этих условиях объективно будет возрастать ценность цивилизационного опыта России, на основе которого в ней была сформирована социальная система комбинированного типа. Показано, что в конечном итоге неизбежен переход от нынешнего доминирования Y-структур к абсолютно новой глобальной системе, устойчивость которой будет основана на новой идеологии, новой духовности (то есть новой «условной информации», по Д. С. Чернавскому), делающей разворот от принципов конкуренции к принципам сотрудничества.

In the article is carried out the analysis of historical process with the use of methods of synergetics (science about the nonlinear developing systems in nature and the society), developed in the works of D. S. Chernavskii in connection with to economic and social systems. It is shown that social self-organizing depending on conditions leads to the formation of both the societies with the strong internal competition (Y-structures) and cooperative type societies (X-structures). Y-structures are characteristic for the countries of the West, X-structure are characteristic for the countries of the East. It is shown that in XIX and in XX centuries occurred accelerated shaping and strengthening of Y-structures. However, at present world system entered into the period of serious structural changes in the economic, political, ideological spheres: the domination of Y-structures concludes. Are examined the possible ways of further development of the world system, connected with change in the regimes of self-organizing and limitation of internal competition. This passage will be prolonged and complex. Under these conditions it will objectively grow the value of the civilizational experience of Russia, on basis of which was formed combined type social system. It is shown that ultimately inevitable the passage from the present do-mination of Y-structures to the absolutely new global system, whose stability will be based on the new ideology, the new spirituality (i.e., new “conditional information” according D. S. Chernavskii), which makes a turn from the principles of competition to the principles of collaboration.

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

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

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

Выполнен анализ диссипативных свойств метода с использованием известных ограничителей вязкости и потоков, а также их линейной комбинации. Разрешающая способность схемы и качество численных решений продемонстрированы на примерах двумерных тестов с обтеканием ступеньки потоком газа с числами Маха 3, 10 и 20, двойным маховским отражением сильной ударной волны и с импульсным сжатием газа. Изучено влияние схемной вязкости метода на численное воспроизведение неустойчивости на контактных поверхностях газов. Установлено, что уменьшение уровня диссипативных свойств схемы в задаче с импульсным сжатием газа приводит к разрушению симметричного решения и формированию хаотической неустойчивости на контактной поверхности.

Численные решения сопоставлены с результатами других авторов, полученных по схемам повышенного порядка аппроксимации: КАБАРЕ, HLLC (Harten Lax van Leer Contact), CFLFh (CFLF hybrid scheme), JT (centered scheme with limiter by Jiang and Tadmor), PPM (Piecewise Parabolic Method), WENO5 (weighted essentially non-oscillatory scheme), RKGD (Runge–Kutta Discontinuous Galerkin), с гибридной взвешенной нелинейной интерполяцией CCSSR-HW4 и CCSSR-HW6. К достоинствам гибридного метода крупных частиц относятся расширенные возможности решения задач гиперболического и смешанного типов, хорошее соотношение диссипативных и дисперсионных свойств, сочетание алгоритмической простоты и высокой разрешающей способности в задачах со сложной ударно-волновой структурой, развитием неустойчивости и вихреобразованием на контактных границах.

We study the computational properties of a parametric class of finite-volume schemes with customizable dissipative properties with splitting by physical processes into Lagrangian, Eulerian, and the final stages (the hybrid large-particle method). The method has a second-order approximation in space and time on smooth solutions. The regularization of a numerical solution at the Lagrangian stage is performed by nonlinear correction of artificial viscosity. Regardless of the grid resolution, the artificial viscosity value tends to zero outside the zone of discontinuities and extremes in the solution. At Eulerian and final stages, primitive variables (density, velocity, and total energy) are first reconstructed by an additive combination of upwind and central approximations weighted by a flux limiter. Then numerical divergent fluxes are formed from them. In this case, discrete analogs of conservation laws are performed.

The analysis of dissipative properties of the method using known viscosity and flow limiters, as well as their linear combination, is performed. The resolution of the scheme and the quality of numerical solutions are demonstrated by examples of two-dimensional benchmarks: a gas flow around the step with Mach numbers 3, 10 and 20, the double Mach reflection of a strong shock wave, and the implosion problem. The influence of the scheme viscosity of the method on the numerical reproduction of a gases interface instability is studied. It is found that a decrease of the dissipation level in the implosion problem leads to the symmetric solution destruction and formation of a chaotic instability on the contact surface.

Numerical solutions are compared with the results of other authors obtained using higher-order approximation schemes: CABARET, HLLC (Harten Lax van Leer Contact), CFLFh (CFLF hybrid scheme), JT (centered scheme with limiter by Jiang and Tadmor), PPM (Piecewise Parabolic Method), WENO5 (weighted essentially non-oscillatory scheme), RKGD (Runge –Kutta Discontinuous Galerkin), hybrid weighted nonlinear schemes CCSSR-HW4 and CCSSR-HW6. The advantages of the hybrid large-particle method include extended possibilities for solving hyperbolic and mixed types of problems, a good ratio of dissipative and dispersive properties, a combination of algorithmic simplicity and high resolution in problems with complex shock-wave structure, both instability and vortex formation at interfaces.

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

One of the numerical methods for solving problems of scattering of electromagnetic waves by systems formed by parallel oriented cylindrical elements — two-dimensional photonic crystals — is considered. The method is based on the classical method of separation of variables for solving the wave equation. Тhe essence of the method is to represent the field as the sum of the primary field and the unknown secondary scattered on the elements of the medium field. The mathematical expression for the latter is written in the form of infinite series in elementary wave functions with unknown coefficients. In particular, the field scattered by N elements is sought as the sum of N diffraction series, in which one of the series is composed of the wave functions of one body, and the wave functions in the remaining series are expressed in terms of the eigenfunctions of the first body using addition theorems. From satisfying the boundary conditions on the surface of each element we obtain systems of linear algebraic equations with an infinite number of unknowns — the required expansion coefficients, which are solved by standard methods. A feature of the method is the use of analytical expressions describing diffraction by a single element of the system. In contrast to most numerical methods, this approach allows one to obtain information on the amplitude-phase or spectral characteristics of the field only at local points of the structure. The absence of the need to determine the field parameters in the entire area of space occupied by the considered multi-element system determines the high efficiency of this method. The paper compares the results of calculating the transmission spectra of two-dimensional photonic crystals by the considered method with experimental data and numerical results obtained using other approaches. Their good agreement is demonstrated.

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

A method for studying the dispersion characteristics of photonic crystals — media with a dielectric constant that varies periodically in space — is considered. The method is based on the representation of the wave functions and permittivity of a periodic medium in the form of Fourier series and their subsequent substitution into the wave equation, which leads to the formulation of the dispersion equation. Using the latter, for each value of the wave vector it is possible determined a set of eigen frequencies. Each of eigen frequency forms a separate dispersion curve as a continuous function of the wave number. The Fourier expansion coefficients of the permittivity, which depend on the vectors of the reciprocal lattice of the photonic crystal, are determined on the basis of data on the geometric characteristics of the elements that form the crystal, their electrophysical properties and the density of the crystal. The solution of the dispersion equation found makes it possible to obtain complete information about the number of modes propagating in a periodic structure at different frequencies, and about the possibility of forming band gaps, i.e. frequency ranges within which wave propagation through a photonic crystal is impossible. The focus of this work is on the application of this method to the analysis of the dispersion properties of metallic photonic crystals. The difficulties that arise in this case due to the presence of intrinsic dispersion properties of the metals that form the elements of the crystal are overcome by an analytical description of their permittivity based on the model of free electrons. As a result, a dispersion equation is formulated, the numerical solution of which is easily algorithmized. That makes possible to determine the dispersion characteristics of metallic photonic crystals with arbitrary parameters. Obtained by this method the results of calculation of dispersion diagrams, which characterize two-dimensional metal photonic crystals, are compared with experimental data and numerical results obtained using the method of self-consistent equations. Their good agreement is demonstrated.

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

A mathematical model of two-phase capillary-nonequilibrium isothermal flows of incompressible phases in a double porosity medium is constructed. A double porosity medium is considered, which is a composition of two porous media with contrasting capillary properties (absolute permeability, capillary pressure). One of the constituent media has high permeability and is conductive, the second is characterized by low permeability and forms an disconnected system of matrix blocks. A feature of the model is to take into account the influence of capillary nonequilibrium on mass transfer between subsystems of double porosity, while the nonequilibrium properties of two-phase flow in the constituent media are described in a linear approximation within the Hassanizadeh model. Homogenization by the method of formal asymptotic expansions leads to a system of partial differential equations, the coefficients of which depend on internal variables determined from the solution of cell problems. Numerical solution of cell problems for a system of partial differential equations is computationally expensive. Therefore, a thermodynamically consistent kinetic equation is formulated for the internal parameter characterizing the phase distribution between the subsystems of double porosity. Dynamic relative phase permeability and capillary pressure in the processes of drainage and impregnation are constructed. It is shown that the capillary nonequilibrium of flows in the constituent subsystems has a strong influence on them. Thus, the analysis and modeling of this factor is important in transfer problems in systems with double porosity.