Все выпуски

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

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

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

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.

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

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

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

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

We study an appearance of gravitational convection in a porous medium saturated by the double-diffusive fluid. The rectangle heated from below is considered with anisotropy of media properties. We analyze Darcy – Boussinesq equations for a binary fluid with Soret effect.

Resulting system for the stream function, the deviation of temperature and concentration is cosymmetric under some additional conditions for the parameters of the problem. It means that the quiescent state (mechanical equilibrium) loses its stability and a continuous family of stationary regimes branches off. We derive explicit formulas for the critical values of the Rayleigh numbers both for temperature and concentration under these conditions of the cosymmetry. It allows to analyze monotonic instability of mechanical equilibrium, the results of corresponding computations are presented.

A finite-difference discretization of a second-order accuracy is developed with preserving of the cosymmetry of the underlying system. The derived numerical scheme is applied to analyze the stability of mechanical equilibrium.

The appearance of stationary and nonstationary convective regimes is studied. The neutral stability curves for the mechanical equilibrium are presented. The map for the plane of the Rayleigh numbers (temperature and concentration) are displayed. The impact of the parameters of thermal diffusion on the Rayleigh concentration number is established, at which the oscillating instability precedes the monotonic instability. In the general situation, when the conditions of cosymmetry are not satisfied, the derived formulas of the critical Rayleigh numbers can be used to estimate the thresholds for the convection onset.

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

Численные исследования нестационарных режимов смешанной конвекции жидкости с переменной вязкостью в горизонтальном канале с тепловыделяющим источником были проведены при следующих значениях безразмерных параметров: $\mathrm{Pr} = 7.0$, $\varepsilon = 0.8$, $\mathrm{Gr} = 10^5$, $C = 0-1$, $10^{-5} < \mathrm{Da} < 10^{-1}$, $50 < \mathrm{Re} < 500$, $\delta = l/H = 0.6-3$. Все распределения изолиний функции тока и температуры, а также зависимости среднего числа Нуссельта и средней температуры были получены в стационарном режиме, когда наблюдается установление картины течения и теплопереноса. В результате анализа установлено, что введение пористой вставки позволяет интенсифицировать теплосъем с поверхности источника энергии. Увеличение размеров пористой ставки, а также использование рабочих сред с разными теплофизическими характеристиками приводят к снижению температуры в источнике энергии.

Numerical study of unsteady mixed convection in an open partially porous horizontal channel with a heatgenerating source was performed. The outer surfaces of horizontal walls of finite thickness were adiabatic. In the channel there was a Newtonian heat-conducting fluid with a temperature-dependent viscosity. The discrete heatconducting and heat-generating source is located inside the bottom wall. The temperature of the fluid phase was equal to the temperature of the porous medium, and calculations were performed using the local thermal equilibrium model. The porous insertion is isotropic, homogeneous and permeable to fluid. The Darcy–Brinkman model was used to simulate the transport process within the porous medium. Governing equations formulated in dimensionless variables “stream function – vorticity – temperature” using the Boussinesq approximation were solved numerically by the finite difference method. The vorticity dispersion equation and energy equation were solved using locally one-dimensional Samarskii scheme. The diffusive terms were approximated by central differences, while the convective terms were approximated using monotonic Samarskii scheme. The difference equations were solved by the Thomas algorithm. The approximated Poisson equation for the stream function was solved separately by successive over-relaxation method. Optimal value of the relaxation parameter was found on the basis of computational experiments. The developed computational code was tested using a set of uniform grids and verified by comparing the results obtained of other authors.

Numerical analysis of unsteady mixed convection of variable viscosity fluid in the horizontal channel with a heat-generating source was performed for the following parameters: $\mathrm{Pr} = 7.0$, $\varepsilon = 0.8$, $\mathrm{Gr} = 10^5$, $C = 0-1$, $10^{-5} < \mathrm{Da} < 10^{-1}$, $50 < \mathrm{Re} < 500$, $\delta = l/H = 0.6-3$. Distributions of the isolines of the stream function, temperature and the dependences of the average Nusselt number and the average temperature inside the heater were obtained in a steady-state regime, when the stationary picture of the flow and heat transfer is observed. As a result we showed that an addition of a porous insertion leads to an intensification of heat removal from the surface of the energy source. The increase in the porous insertion sizes and the use of working fluid with different thermal characteristics, lead to a decrease in temperature inside the source.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Статья посвящена моделированию гидродинамических процессов мелководных водоемов на примере Азовского моря. В статье приведена математическая модель гидродинамики мелководного водоема, позволяющая вычислить трехмерные поля вектора скорости движения водной среды. Применение регуляризаторов по Б.Н. Четверушкину в уравнении неразрывности привело к изменению способа расчета поля давления, базирующегося на решении волнового уравнения. Построена дискретная конечно-разностная схема для расчета давления в области, линейные размеры которой по вертикали существенно меньше размеров по горизонтальным координатным направлениям, что является характерным для геометрии мелководных водоемов. Описаны метод и алгоритм решения сеточных уравнений с предобуславливателем трехдиагонального вида. Предложенный метод применен для решения сеточных уравнений, возникающих при расчете давления для трехмерной задачи гидродинамики Азовского моря. Показано, что предложенный метод сходится быстрее модифицированного попеременно-треугольного метода. Представлена параллельная реализация предложенного метода решения сеточных уравнений и проведены теоретические и практические оценки ускорения алгоритма с учетом времени латентности вычислительной системы. Приведены результаты вычислительных экспериментов для решения задач гидродинамики Азовского моря с использованием гибридной технологии MPI + OpenMP. Разработанные модели и алгоритмы применялись для реконструкции произошедшей в 2001 году в Азовском море экологической катастрофы и решения задачи движения водной среды в устьевых районах. Численные эксперименты проводились на гибридном вычислительном кластере К-60 ИПМ им. М.В. Келдыша РАН.

The article is devoted to modeling the shallow water hydrodynamic processes using the example of the Azov Sea. The article presents a mathematical model of the hydrodynamics of a shallow water body, which allows one to calculate three-dimensional fields of the velocity vector of movement of the aquatic environment. Application of regularizers according to B.N.Chetverushkin in the continuity equation led to a change in the method of calculating the pressure field, based on solving the wave equation. A discrete finite-difference scheme has been constructed for calculating pressure in an area whose linear vertical dimensions are significantly smaller than those in horizontal coordinate directions, which is typical for the geometry of shallow water bodies. The method and algorithm for solving grid equations with a tridiagonal preconditioner are described. The proposed method is used to solve grid equations that arise when calculating pressure for the three-dimensional problem of hydrodynamics of the Azov Sea. It is shown that the proposed method converges faster than the modified alternating triangular method. A parallel implementation of the proposed method for solving grid equations is presented and theoretical and practical estimates of the acceleration of the algorithm are carried out taking into account the latency time of the computing system. The results of computational experiments for solving problems of hydrodynamics of the Sea of Azov using the hybrid MPI + OpenMP technology are presented. The developed models and algorithms were used to reconstruct the environmental disaster that occurred in the Sea of Azov in 2001 and to solve the problem of the movement of the aquatic environment in estuary areas. Numerical experiments were carried out on the K-60 hybrid computing cluster of the Keldysh Institute of Applied Mathematics of Russian Academy of Sciences.

Представлена программа NINE (Newtonian Iteration for Nonlinear Equation) численного решения граничных задач для нелинейных дифференциальных уравнений второго порядка на основе непрерывного аналога метода Ньютона (НАМН) с использованием нумеровской конечно-разностной аппроксимации четвертого порядка относительно шага дискретизации по пространственной переменной. Обсуждаются алгоритмы вычисления ньютоновского итерационного параметра. Выполнены методические расчеты, демонстрирующие влияние выбора итерационного параметра на сходимость итерационного процесса. Представлены результаты проведенного с помощью программы NINE численного исследования положительных частицеподобных решений уравнения скалярного поля.

The computer code NINE (Newtonian Iteration for Nonlinear Equation) for numerical solution of the boundary problems for nonlinear differential equations on the basis of continuous analogue of the Newton method (CANM) is presented. Numerov’s finite-difference appproximation is applied to provide the fourth accuracy order with respect to the discretization stepsize. Algorithms of calculating the Newtonian iterative parameter are discussed. A convergence of iteration process in dependence on choice of the iteration parameter has been studied. Results of numerical investigation of the particle-like solutions of the scalar field equation are given.

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

This paper presents a method for obtaining fractal signals using fractal splines similar to signals generated by fractal functions. The hypothesis about the identity of discrete fractal functions and linear fractal splines is justified. There are considered the features of planning matrix calculation of cumulative fractal spline, examples of generated curves are shown.