

Home – Home – Jornals – Siberian Journal of Numerical Mathematics 2019 number 4
2019 year, number 4
A.B. Bakushinsky^{1}, A.S. Leonov^{2}
^{1}Institute for Systems Analysis, Moscow, Russia ^{2}National Nuclear Research University В«MEPHI», Moscow, Russia
Keywords: трехмерное волновое уравнение, обратная коэффициентная задача, регуляризующий алгоритм, быстрое преобразование Фурье, threedimensional wave equation, wave field, inverse coefficient problem, regularizing algorithm, fast Fourier transform
Abstract >>
A threedimensional coefficient inverse problem for the wave equation (with losses) in a cylindrical domain is under consideration. The data for its solution are special time integrals of the wave field measured in a cylindrical layer. We present and substantiate an efficient algorithm for solving such a threedimensional problem based on the fast Fourier transform. The algorithm proposed makes possible to obtain a solution on grids of 512x512x512 size in a time of about 1.4 hours on a typical PC without parallelizing the calculations. The results of the numerical experiments for solving the corresponding model inverse problems are presented.

V.G. Borisov^{1}, Y.N. Zakharov^{1}, Y.I. Shokin^{2}, E.A. Ovcharenko^{3}, K.Y. Klyshnikov^{3}, I.N. Sizova^{3}, A.V. Batranin^{4}, Y.A. Kudryavtseva^{3}, P.S. Onishchenko^{2,3}
^{1}Institute of Computational Technologies SB RAS, Kemerovo, Russia ^{2}Institute of Computational Technologies SB RAS, Novosibirsk, Russia ^{3}Research Institute for Complex Issues of Cardiovascular Diseases under the Siberian Branch of the Russian Academy of Medical Sciences, Kemerovo, Russia ^{4}Tomsk Polytechnic University, Tomsk, Russia
Keywords: компьютерное моделирование, течение крови, биопротезы, пристеночное напряжение сдвига, computer modeling, blood flow, bioprostheses, wall shear stress
Abstract >>
The threedimensional unsteadystate periodic flow of blood in xenogenic vascular bioprostheses is simulated using computational fluid dynamics methods. The geometry of the computational domain is based on microtomographic scanning of bioprostheses. To set a variable pressure gradient causing a nonstationary flow in the prostheses, personalspecific data of the Dopplerechography of the blood flow of a particular patient are used. A comparative analysis of the velocity fields in the flow areas corresponding to three real samples of bioprostheses with multiple stenoses is carried out. In the zones of stenosis and outside of them, the distribution of the nearwall shear stress, which influences the risk factors for thrombosis in the prostheses, is analyzed. An algorithm for predicting the hemodynamic effects arising in vascular bioprostheses, based on the numerical modeling of a blood flow in them, is proposed.

E.A. Vorontsova^{1,2}, A.V. Gasnikov^{3,4,5}, A.C. Ivanova^{3}, E.A. Nurminsky^{1}
^{1}Far Eastern Federal University, Vladivostok, Russia ^{2}Universite de GrenobleAlpes, SaintMartind'Heres, France ^{3}Moscow Institute of Physics and Technology, Dolgoprudny, Russia ^{4}Institute for Information Transmission Problems RAS, Moscow, Russia ^{5}Adyghe State University, Maikop, Russia
Keywords: вальрасов механизм, децентрализация цен, прямодвойственный метод, субградиентный метод, условие Слейтера, Walrasian equilibrium, decentralized pricing, primaldual method, subgradient method, Slater condition
Abstract >>
We consider the resource allocation problem and its numerical solution. The following is demonstrated: 1) the Walrasian priceadjustment mechanism for determining the equilibrium; 2) the decentralized role of the prices; 3) Slater's method for price restrictions (dual Lagrange multipliers); 4) a new mechanism for determining equilibrium prices, in which prices are fully controlled not by Center (Government), but by economic agents  nodes (factories). In the economic literature, only the convergence of the methods considered is proved. In contrast, this paper provides an accurate analysis of the convergence rate of the described procedures for determining the equilibrium. The analysis is based on the primaldual nature of the algorithms proposed. More precisely, in this paper, we propose the economic interpretation of the following numerical primaldual methods of the convex optimization: dichotomy and subgradient projection method.

M.I. Ivanov^{1}, I.A. Kremer^{1,2}, M.V. Urev^{1,2}
^{1}Institute of Computational Mathematics and Mathematical Geophysics SB RAS, Novosibirsk, Russia ^{2}Novosibirsk State University, Novosibirsk, Russia
Keywords: вырожденная задача Неймана, условия согласования, ортогонализация правой части, конечные элементы, degenerate Neumann problem, matching conditions, orthogonalization of the righthand side, finite elements
Abstract >>
This paper deals with the solution of the degenerate Neumann problem for the diffusion equation by the finite element method. First, an extended generalized formulation of the Neumann problem in the Sobolev space H^{1}(Ω) is derived and investigated. Then a discrete analogue of this problem is formulated using standard finite element approximations of the space H^{1}(Ω). An iterative method for solving the corresponding SLAE is proposed. Some examples of solving the model problems are used to discuss the numerical peculiarities of the algorithm proposed.

E.A. Karatsuba
Federal Research Center "Computer Science and Control" of Russian Academy of Sciences, Moscow, Russia
Keywords: функции Бесселя, быстрые алгоритмы, сложность вычисления, метод БВЕ, большой аргумент, эффективное вычисление, Bessel functions, fast algorithms, computational complexity, FEE method, large argument, efficient calculation
Abstract >>
Two algorithms for an effective calculation of the Bessel function are presented: a fast algorithm with an increasing accuracy of computation and a computational algorithm for the case of a large argument of the Bessel function.

V.N. Lutay
Southern Federal University, RostovonDon, Russia
Keywords: плохо обусловленные матрицы, треугольное разложение, повышение устойчивости, отсечение младших разрядов, неполное скалярное произведение, illconditioned matrix, triangular decomposition, improving resilience, cutting off the least significant bits of partial scalar product
Abstract >>
An approach to increasing the stability of triangular decomposition of a dense positive definite matrix with a large condition number with the use of the Gauss and the Cholesky methods is considered. It is proposed to introduce additions to standard computational schemes, which consist in the use of an incomplete scalar product of two vectors, which is formed by cutting off the lower digits of the sum of the products of two numbers. Cutting off being performed in the process of factorization leads to an increase in the diagonal elements of triangular matrices to a random number and prevents the appearance of very small numbers during the decomposition according to Gauss and a negative radical expression in the Cholesky method. The number of additional operations required to obtain an accurate solution is estimated. The results of computational experiments are presented.

A.G. Megrabov^{1,2}
^{1}Institute of Computational Mathematics and Mathematical Geophysics SB RAS, Novosibirsk, Russia ^{2}Novosibirsk State Technical University, Novosibirsk, Russia
Keywords: кинематическая сейсмика, геометрическая оптика, уравнение эйконала, семейство лучей, семейство фронтов, законы сохранения, дифференциальная геометрия, геометрия векторных полей, kinematic seismic, geometric optics, eikonal equation, family of rays, family of wavefronts, conservation laws, differential geometry, geometry of vector fields
Abstract >>
In the previous studies, the author has obtained the conservation laws for the 2D eikonal equation in an inhomogeneous isotropic medium. These laws represent the divergent identities of the form div F =0. The vector field F is expressed in terms of the solution to the eikonal equation (the time field), the refractive index (the equation parameter) and their partial derivatives. Also, there were found equivalent conservation laws (divergent identities) for the families of rays and the families of wavefronts in terms of their geometric characteristics. Thus, the geometric essence (interpretation) of the abovementioned conservation laws for the 2D eikonal equation was discovered. In this paper, the 3D analogs to the results obtained are presented: differential conservation laws for the 3D eikonal equation and the conservation laws (divergent identities of the form div F =0) for the family of rays and the family of wavefronts, the vector field F is expressed in terms of classical geometric characteristics of the ray curves: their Frenet basis (unit tangent vector, a principal normal and a binormal), the first curvature and the second curvature, or in terms of the classical geometric characteristics of the wavefront surfaces, i. e. their normal, principal directions, principal curvatures, the Gaussian curvature and the mean curvature. All the results have been obtained based on the vector and geometric formulas (differential conservation laws and some formulas) obtained for the families of arbitrary smooth curves, the families of arbitrary smooth surfaces and arbitrary smooth vector fields.

S.I. Fadeev^{1,2}, V.V. Kogai^{1,2}
^{1}Sobolev Institute of Mathematics, Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russia ^{2}Novosibirsk State University, Novosibirsk, Russia
Keywords: математическая модель, микрогенератор, задача Коши, краевая задача, периодические колебания, предельный цикл, устойчивость, фазовая плоскость, продолжение решения по параметру, mathematical model, microgenerator, Cauchy problem, boundary value problem, periodic oscillations, limit cycle, phase plane, continuation of the solution with respect to the parameter
Abstract >>
In this paper, we consider a mathematical model of a new type of a microgenerator, based on generating the oscillations of a mobile electrode in a microgap due to electrostatic forces. The principle of operation of the generator is analogous to the wellknown theory of the clockescapement mechanism, with the difference that in the equation of motion the form of the righthand side corresponds to the electrostatic nature of the impulse action. The numerical analysis shows that the bounded oscillations with an increase in time tend towards a stable limit cycle in the phase plane and, thereby, the emerging oscillations are stable with respect to external perturbations. In studying periodic oscillations, depending on the parameters of a model, we use the solution of the boundary value problem for the equation with a discontinuous righthand side, transformed to a form allowing the application of the numerical continuation method. In this way, the area in the plane of the model parameters is defined, in which stable limit cycles exist.

