Все выпуски

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

This paper covers deployable systems assembled from trapezium plates. When the plate package is unwrapped, a net shell with six loop cells is formed. It is proved that additional degrees of freedom appear in case of certain correlation between the sizes of the six loop faces. When thin plates were used, the continuum approximation of the deployed net could be interpreted as a shell with a wide variety of local curvatures. Kinematics of the continuum model is analyzed by the method of Cartan moving hedron. Mechanical behavior of continuum nets is studied when cylindrical hinges between the plates are completed of shape memory plastic materials. The paper researches into shell transformations from one stable form to the other. Various practical applications of the continuum nets are demonstrated.

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

In the current paper, we study a Galerkin–Petrov method for a parabolic equations of higher order in domain with a moving boundary. Asymptotic estimates for the convergence rate of approximate solutions are obtained.

Для многомерного нелокального уравнения Фишера–Колмогорова–Петровского–Пискунова в классе траекторно-сосредоточенных функций построены квазиклассические асимптотики с точностью $O(D^{N/2})$, $N\geqslant3$. С помощью операторов симметрии получен счетный набор асимптотических решений исходного уравнения с точностью $O(D^{3/2})$. В явном виде построены асимптотические решения двумерного уравнения Фишера–Колмогорова–Петровского–Пискунова.

Semiclassical asymptotic solutions with accuracy $O(D^{N/2})$, $N\geqslant3$ are constructed for the multidimensional Fisher–Kolmogorov–Petrovskii–Piskunov equation in the class of trajectory-concentrated functions. Using the symmetry operators a countable set of asymptotic solutions with accuracy $O(D^{3/2})$ is obtained. Asymptotic solutions of two-dimensional Fisher–Kolmogorov–Petrovskii–Piskunov equation are found in explicit
form.

Настоящая работа посвящена применению теории недавно разработанных полулокальных сглаживающих сплайнов, или $S$-сплайнов высоких степеней, к решению задач теории упругости. $S$-сплайн — кусочно-полиномиальная функция, коэффициенты полиномов которой определяются из двух условий: первая часть коэффициентов определяется условиями гладкой склейки, остальные определяются методом наименьших квадратов. Мы рассмотрим, каким образом могут быть применены сплайны 7-ой степени класса $C^4$ при решении бигармонического уравнения на круге.

This article is dedicated to the use of higher degree $S$-splines for solving equations of the elasticity theory. As an example we consider the solution to the equation of bending of a plate on a circle. $S$-spline is a piecewise-polynomial function. Its coefficients are defined by two conditions. The first part of the coefficients are defined by the smoothness of the spline. The rest are determined using the least-squares method. We consider class $C^4$ 7th degree $S$-splines.

Рассматривается модель стека длинных джозефсоновских переходов (ДДП), состоящего из чередующихся сверхпроводящих слоев и слоев диэлектрика, с учетом индуктивной и емкостной связи между слоями. Модель описывается системой нелинейных дифференциальных уравнений в частных производных относительно разности фаз и напряжения между соседними сверхпроводящими слоями в стеке ДДП, с соответствующими начальными и граничными условиями. Численное решение этой системы уравнений основано на использовании стандартных трехточечных конечно-разностных формул для дискретной аппроксимации по пространственной координате и применении четырехшагового метода Рунге–Кутты для решения полученной задачи Коши. Разработанный параллельный алгоритм реализован на основе технологии MPI (Message Passing Interface). В работе дана математическая постановка задачи в рамках рассматриваемой модели, описаны вычислительная схема и методика расчета вольт-амперных характеристик системы ДДП, представлены два варианта параллельной реализации. Продемонстрировано влияние индуктивной и емкостной связи между ДДП на структуру вольт-амперной характеристики в рамках рассматриваемой модели. Представлены результаты методических расчетов с различными параметрами длины и количества джозефсоновских переходов в стеке ДДП в зависимости от количества задействованных параллельных вычислительных узлов. Расчеты выполнены на многопроцессорных кластерах HybriLIT и ЦИВК Многофункционального информационно-вычислительного комплекса Лаборатории информационных технологий Объединенного института ядерных исследований (Дубна). На основе полученных численных результатов обсуждается эффективность рассмотренных вариантов распределения вычислений для численного моделирования системы ДДП в параллельном режиме. Показано, что один из предложенных подходов приводит к ускорению вычислений до 9 раз по сравнению с расчетами в однопроцессорном режиме.

We consider a model of stacked long Josephson junctions (LJJ), which consists of alternating superconducting and dielectric layers. The model takes into account the inductive and capacitive coupling between the neighbor junctions. The model is described by a system of nonlinear partial differential equations with respect to the phase differences and the voltage of LJJ, with appropriate initial and boundary conditions. The numerical solution of this system of equations is based on the use of standard three-point finite-difference formulae for discrete approximations in the space coordinate, and the applying the four-step Runge-Kutta method for solving the Cauchy problem obtained. Designed parallel algorithm is implemented by means of the MPI technology (Message Passing Interface). In the paper, the mathematical formulation of the problem is given, numerical scheme and a method of calculation of the current-voltage characteristics of the LJJ system are described. Two variants of parallel implementation are presented. The influence of inductive and capacitive coupling between junctions on the structure of the current-voltage characteristics is demonstrated. The results of methodical calculations with various parameters of length and number of Josephson junctions in the LJJ stack depending on the number of parallel computing nodes, are presented. The calculations have been performed on multiprocessor clusters HybriLIT and CICC of Multi-Functional Information and Computing Complex (Laboratory of Information Technologies, Joint Institute for Nuclear Research, Dubna). The numerical results are discussed from the viewpoint of the effectiveness of presented approaches of the LJJ system numerical simulation in parallel. It has been shown that one of parallel algorithms provides the 9 times speedup of calculations.

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

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

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

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

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

We consider the approaches to the construction of methods for solving four-dimensional programming problems for calculating directions for multiple minimizations of smooth functions on a set of a given set of linear equalities. The approach consists of two stages.

At the first stage, the problem of quadratic programming is transformed by a numerically stable direct multiplicative algorithm into an equivalent problem of designing the origin of coordinates on a linear manifold, which defines a new mathematical formulation of the dual quadratic problem. For this, a numerically stable direct multiplicative method for solving systems of linear equations is proposed, taking into account the sparsity of matrices presented in packaged form. The advantage of this approach is to calculate the modified Cholesky factors to construct a substantially positive definite matrix of the system of equations and its solution in the framework of one procedure. And also in the possibility of minimizing the filling of the main rows of multipliers without losing the accuracy of the results, and no changes are made in the position of the next processed row of the matrix, which allows the use of static data storage formats.

At the second stage, the necessary and sufficient optimality conditions in the form of Kuhn–Tucker determine the calculation of the direction of descent — the solution of the dual quadratic problem is reduced to solving a system of linear equations with symmetric positive definite matrix for calculating of Lagrange's coefficients multipliers and to substituting the solution into the formula for calculating the direction of descent.

It is proved that the proposed approach to the calculation of the direction of descent by numerically stable direct multiplicative methods at one iteration requires a cubic law less computation than one iteration compared to the well-known dual method of Gill and Murray. Besides, the proposed method allows the organization of the computational process from any starting point that the user chooses as the initial approximation of the solution.

Variants of the problem of designing the origin of coordinates on a linear manifold, a convex polyhedron and a vertex of a convex polyhedron are presented. Also the relationship and implementation of methods for solving these problems are described.

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

Non-collapsing soliton-like wave functions are shown to exist in semiclassical approximation for the Bose-Einstein condensate model based on the Gross–Pitaevskii equation with attractive nonlinearity and external field of magnetic trap of special form.

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

The a priori analysis of approximation of magnetohydrodynamic equations on irregular quadrangular analysis was performed. The values of coefficients wich determine the misalignment norm for difference analogs of operators gradient and divergence were calculated. Was studied the influence of properties of grid cells on misalignment. For the numerical confirmation of obtained estimations were cited the examples of calculations with specifying identical initial data on different grids.

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

The modified two-sweep iteration procedure was proposed, built according to the least-squares method scheme, which determines progressive approximations to the periodic solution of Mathieu’s equation and his own function, considerably superior according to the accuracy earlier well-known results.