All issues

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

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

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

The article contains the description of developed mathematical models of fractures which can be used for numerical solution of exploration seismology problems with use of grid-characteristic method on unstructured triangular and tetrahedral meshes. The base of developed models is the concept of infinitely thin fracture. This fracture is represented by contact boundary. Such approach significantly reduces the consumption of computer resources by the absence of the mesh definition inside of fracture necessity. By the other side it lets state the fracture discretely in integration domain, therefore one can observe qualitative new effects which are not available to observe by use of effective models of fractures, actively used in computational seismic.

The main target in the development of models have been getting the most accurate result. Developed models thet can receive the response close to the actual response of the existing fracture in geological environment. We considered fluid-filled fractures, glued and partially glued fractures, and also fractures with dynamical friction force. Fracture behavior determinated by the nature of condition on the border.

Empty fracture was represented as free boundary condition. This condition give us opportunity for total reflection of wave fronts from fracture. Fluid-filling provided the condition for sliding on the border. Under this condition, there was a passage of longitudinal and total reflection of converted waves. For the real fractures, which has unequal distance between the borders has been proposed the model of partially glued fracture. At different points of the fracture's boundary were sat different conditions. Almost the same effect is achieved by using a fracture model of dynamic friction condition. But its disadvantage is the inabillity to specify the proportion of fracture's glued area due to the friction factor can take values from zero to infinity. The model of partially glued fracture is devoid of this disadvantage.

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

In the paper the statistical relationships between the size and production characteristics of phytoplankton and zooplankton of the Vistula and Curonian lagoons, the Baltic Sea, were investigated. Research phytoplankton and zooplankton within the Russian part of the area of the Vistula and the Curonian lagoon was carried out on the monthly basis (from April to November) within the framework of long-term monitoring program on evaluating of ecological status of the lagoons. The size structure of plankton is the basis for understanding of the development of production processes, mechanisms of formation of the plankton species diversity and functioning of the lagoon ecosystems. As results of the work it was found that the maximum rate of photosynthesis and the integral value of the primary production with a change in cell volume of phytoplankton are changed according to a power law. The result shows that the smaller the size of algal cells in phytoplankton communities the more actively occur metabolism and the more effective they assimilate the solar energy. It is shown that the formation of plankton species diversity in ecosystems of lagoons is closely linked with the size structure of plankton communities and with features of development of the production processes. It is proposed the structure of a spatially homogenous mathematical model of the plankton food chain for the lagoon ecosystems taking into account the size spectrum and the characteristics of phytoplankton and zooplankton. The model parameters are the sizedependent indicators allometrically linked with average volumes of cells and organisms in different ranges of their sizes. In the model the algorithm for changes over time the coefficients of food preferences in the diet of zooplankton was proposed. Developed the size-dependent mathematical model of aquatic ecosystems allows to consider the impact of turbulent exchange on the size structure and temporal dynamics of the plankton food chain of the Vistula and Curonian lagoons. The model can be used to study the different regimes of dynamic behavior of plankton systems depending on the changes in the values of its parameters and external influences, as well as to quantify the redistribution of matter flows in ecosystems of the lagoons.

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

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

We study the influence of fishery on a structured fish population under random changes of habitat conditions. The population parameters correspond to dominant pelagic fish species of Far-Eastern seas of the northwestern part of the Pacific Ocean (pollack, herring, sardine). Similar species inhabit various parts of the Word Ocean. The species body size distribution was chosen as a main population feature. This characteristic is easy to measure and adequately defines main specimen qualities such as age, maturity and other morphological and physiological peculiarities. Environmental fluctuations have a great influence on the individuals in early stages of development and have little influence on the vital activity of mature individuals. The fishery revenue was chosen as an optimality criterion. The main control characteristic is fishing effort. We have chosen quadratic dependence of fishing revenue on the fishing effort according to accepted economic ideas stating that the expenses grow with the production volume. The model study shows that the population structure ensures the increased population stability. The growth and drop out of the individuals’ due to natural mortality smoothens the oscillations of population density arising from the strong influence of the fluctuations of environment on young individuals. The smoothing part is played by diffusion component of the growth processes. The fishery in its turn smooths the fluctuations (including random fluctuations) of the environment and has a substantial impact upon the abundance of fry and the subsequent population dynamics. The optimal time-dependent fishing effort strategy was compared to stationary fishing effort strategy. It is shown that in the case of quickly changing habitat conditions and stochastic dynamics of population replenishment there exists a stationary fishing effort having approximately the same efficiency as an optimal time-dependent fishing effort. This means that a constant or weakly varying fishing effort can be very efficient strategy in terms of revenue.

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

We present results of a study aimed at identification of a controlled object’s channels based on postprocessing of measurements with development of a model of a multiple-input controlled object and subsequent active modelling experiment. The controlled object model is developed using approximation of its behavior by a neural network model using trends obtained during a passive experiment in the mode of normal operation. Recurrent neural network containing feedback elements allows to simulate behavior of dynamic objects; input and feedback time delays allow to simulate behavior of inertial objects with pure delay. The model was taught using examples of the object’s operation with a control system and is presented by a dynamic neural network and a model of a regulator with a known regulation function. The neural network model simulates the system’s behavior and is used to conduct active computing experiments. Neural network model allows to obtain the controlled object’s response to an exploratory stimulus, including a periodic one. The obtained complex frequency response is used to evaluate parameters of the object’s transfer system using the least squares method. We present an example of identification of a channel of the simulated control system. The simulated object has two input ports and one output port and varying transport delays in transfer channels. One of the input ports serves as a controlling stimulus, the second is a controlled perturbation. The controlled output value changes as a result of control stimulus produced by the regulator operating according to the proportional-integral regulation law based on deviation of the controlled value from the task. The obtained parameters of the object’s channels’ transfer functions are close to the parameters of the input simulated object. The obtained normalized error of the reaction for a single step-wise stimulus of the control system model developed based on identification of the simulated control system doesn’t exceed 0.08. The considered objects pertain to the class of technological processes with continuous production. Such objects are characteristic of chemical, metallurgic, mine-mill, pulp and paper, and other industries.

В статье представлены результаты численного моделирования и экспериментальные данные по высокоскоростному внедрению ударника в преграду. В расчетах использовалась дискретно-элементная модель, основанная на представлении ударника и преграды совокупностью плотно упакованных взаимосвязанных частиц. Данный класс моделей находит все более широкое применение в задачах высокоскоростного взаимодействия тел. В предыдущих работах авторов рассмотрены вопросы применения дискретно-элементной модели к задаче внедрения металлических шаров в массивные преграды. На основе сравнительного анализа данных вычислительных и физических экспериментов было показано, что для широкого класса задач высокоскоростного внедрения достаточно высокая точность дискретно-элементного моделирования может быть достигнута с использованием двухпараметрического потенциала Леннарда–Джонса. При этом была идентифицирована зависимость энергии межэлементной связи от динамической твердости материалов. Использование построенной таким образом дискретно-элементной модели позволило достаточно точно описать наблюдаемые в экспериментах процессы внедрения ударника в массивную преграду в диапазоне скоростей взаимодействия 500–2500 м/c.

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

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

The paper presents the results of numerical simulation and experimental data on the high-speed penetration of the impactor into the obstacle. In the calculations, a discrete-element model has been used, based on the representation of the impactor and the target by a set of close packed interconnected particles. This class of models finds an increasingly wide application in the problems of high-speed interaction of bodies. In the previous works of the authors, the questions of application of the discrete-element model to the problem of the penetration of spherical impactors into massive targets were considered. On the basis of a comparative analysis of the data of computational and physical experiments, it was found out that for a wide class of high-speed penetration problems, a high accuracy of discrete-element modeling can be achieved using the two-parameter Lennard–Jones potential. The binding energy was identified as a function of the dynamic hardness of materials. It was shown that the use of this approach makes it possible to describe accurately the penetration process in the range of impactor velocities 500–2500 m/c.

In this paper, we compare the results of discrete-element modeling with experimental data on penetration of high-strength targets of different thickness by steel impactors. The use of computational parallelization technologies on graphic processors in combination with 3D visualization and animation of the results makes it possible to obtain detailed spatio-temporal patterns of the penetration process and compare them with experimental data.

A comparative analysis of the experimental and calculated data has shown a sufficiently high accuracy of discrete-element modeling for a wide range of target thicknesses: for thin targets pierced with preservation of the integrity of the deformed impactor, for targets of medium thickness, pierced with practically complete fragmentation of the impactor at the exit from the target, and for thick impenetrable targets.

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

Оценка надежности программного обеспечения (ПО) — важная составляющая процесса разработки современного программного обеспечения. Многие исследования направлены на улучшение моделей для измерения и прогнозирования надежности программных продуктов. Однако мало внимания уделяется подходам к сопоставлению существующих систем с точки зрения надежности ПО. Несмотря на огромное значение для практики (и для управления разработкой ПО), полной и проверенной методологии сравнения не существует. В этой статье мы предлагаем методологию сравнения надежности программного обеспечения, в которой широко применяются модели роста надежности программного обеспечения. Методология была оценена на примере трех мобильных операционных систем с открытым исходным кодом: Sailfish, Tizen, CyanogenMod.

Побочным продуктом исследования является сравнение трех мобильных операционных систем с открытым исходным кодом. Целью данного исследования является определение того, какая ОС является более надежной. Для этого были определены 3 вопроса и 8 метрик. С учетом сравнения этих метрик оказалось, что Sailfish в большинстве случаев является самой эффективной операционной системой. Напротив, Tizen показывает лучшее в 3 случаях из 8, но оказывается хуже других систем только в одном случае из 8.

Evaluation of software reliability is an important part of the process of developing modern software. Many studies are aimed at improving models for measuring and predicting the reliability of software products. However, little attention is paid to approaches to comparing existing systems in terms of software reliability. Despite the enormous importance for practice (and for managing software development), a complete and proven comparison methodology does not exist. In this article, we propose a software reliability comparison methodology in which software reliability growth models are widely used. The proposed methodology has the following features: it provides certain level of flexibility and abstraction while keeping objectivity, i.e. providing measurable comparison criteria. Also, given the comparison methodology with a set of SRGMs and evaluation criteria it becomes much easier to disseminate information about reliability of wide range of software systems. The methodology was evaluated on the example of three mobile operating systems with open source: Sailfish, Tizen, CyanogenMod.

A byproduct of our study is a comparison of the three analyzed Open Source mobile operating systems. The goal of this research is to determine which OS is stronger in terms of reliability. To this end we have performed a GQM analysis and we have identified 3 questions and 8 metrics. Considering the comparison of metrics, it appears that Sailfish is in most case the best performing OS. However, it is also the OS that performs the worst in most cases. On the contrary, Tizen scores the best in 3 cases out of 8, but the worst only in one case out of 8.

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

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

Приведены результаты решения типовых задач разрушения твердых деформируемых тел для медленных (квазистатических) и быстрых (динамических) процессов деформации. На основании опыта рас- четов даны рекомендации по моделированию процессов разрушения и обеспечению достоверности численных решений.

The paper reviews the existing approaches to calculating the destruction of solids. The main attention is paid to algorithms using a unified approach to the calculation of deformation both for nondestructive and for the destroyed states of the material. The thermodynamic derivation of the unified rheological relationships taking into account the elastic, viscous and plastic properties of materials and describing the loss of the deformation resistance ability with the accumulation of microdamages is presented. It is shown that the mathematical model under consideration provides a continuous dependence of the solution on input parameters (parameters of the material medium, initial and boundary conditions, discretization parameters) with softening of the material.

Explicit and implicit non-matrix algorithms for calculating the evolution of deformation and fracture development are presented. Non-explicit schemes are implemented using iterations of the conjugate gradient method, with the calculation of each iteration exactly coinciding with the calculation of the time step for two-layer explicit schemes. So, the solution algorithms are very simple.

The results of solving typical problems of destruction of solid deformable bodies for slow (quasistatic) and fast (dynamic) deformation processes are presented. Based on the experience of calculations, recommendations are given for modeling the processes of destruction and ensuring the reliability of numerical solutions.

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

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

Большой объем вычислений требует эффективного распараллеливания алгоритма, что реализовано на основе методики Message Passing Interface (MPI). Передача информации от одного узла MPI к другому осуществляется как разновидность граничного условия — таким образом, каждый MPI узел может хранить только ту часть сетки, которая имеет отношение конкретно к нему.

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

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

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

The purpose of this work is to develop a universal computer program (solver) which solves kinetic Boltzmann equation for simulations of rarefied gas flows in complexly shaped devices. The structure of the solver is described in details. Its efficiency is demonstrated on an example of calculations of a modern many tubes Knudsen pump. The kinetic Boltzmann equation is solved by finite-difference method on discrete grid in spatial and velocity spaces. The differential advection operator is approximated by finite difference method. The calculation of the collision integral is based on the conservative projection method.

In the developed computational program the unstructured spatial mesh is generated using GMSH and may include prisms, tetrahedrons, hexahedrons and pyramids. The mesh is denser in areas of flow with large gradients of gas parameters. A three-dimensional velocity grid consists of cubic cells of equal volume.

A huge amount of calculations requires effective parallelization of the algorithm which is implemented in the program with the use of Message Passing Interface (MPI) technology. An information transfer from one node to another is implemented as a kind of boundary condition. As a result, every MPI node contains the information about only its part of the grid.

The main result of the work is presented in the graph of pressure difference in 2 reservoirs connected by a multitube Knudsen pump from Knudsen number. This characteristic of the Knudsen pump obtained by numerical methods shows the quality of the pump. Distributions of pressure, temperature and gas concentration in a steady state inside the pump and the reservoirs are presented as well.

The correctness of the solver is checked using two special test solutions of more simple boundary problems — test with temperature distribution between 2 planes with different temperatures and test with conservation of total gas mass.

The correctness of the obtained data for multitube Knudsen pump is checked using denser spatial and velocity grids, using more collisions in collision integral per time step.