Home – Home – Jornals – Numerical Analysis and Applications 2017 number 4

In this paper we propose and study the flux predictor-corrector scheme in the three-dimensional case. This scheme is depleted of drawbacks of that constructed on the basis of the Douglas-Gunn prototype-scheme. The scheme proposed demonstrates the second order of accuracy.

In Part 1 of this paper, we consider the web-pages ranking problem also known as the problem of finding the PageRank vector or Google problem. We discuss the connection of this problem with the ergodic theorem and describe different numerical methods to solve this problem together with their theoretical background, such as Markov Chain Monte Carlo and equilibrium in a macrosystem.

We consider a strongly -hard Euclidean problem of finding a subsequence in a finite sequence. The criterion of the solution is a minimum sum of squared distances from the elements of a sought subsequence to its geometric center (centroid). It is assumed that a sought subsequence contains a given number of elements. In addition, a sought subsequence should satisfy the following condition: the difference between the indices of each previous and subsequent points is bounded with given lower and upper constants. We present an approximation algorithm of solving the problem and prove that it is a fully polynomial-time approximation scheme (FPTAS) when the space dimension is bounded by a constant.

Currently, there are international microbarograph networks, with high resolution recording the wave pressure variations on the Earth's surface. This increases the interest in the problems of wave propagation in the atmosphere from variations in the atmospheric pressure. A complete system of nonlinear hydrodynamic equations for an atmospheric gas with lower boundary conditions in the form of wavelike variations on the Earth's surface is considered. Since the wave amplitudes near the Earth's surface are small, linearized equations are used in the analysis of the problem correctness. With the help of the wave energy functional method, it is shown that in the non-dissipative case, the solution of the boundary value problem is uniquely determined by the variable pressure field on the Earth's surface. The corresponding dissipative problem is correct if, in addition to the pressure field, suitable conditions on the velocity and temperature on the Earth's surface are given. In the case of an isothermal atmosphere, the problem admits analytical solutions that are harmonic in the variables x and t . A good agreement between numerical solutions and analytical ones is shown. The study has shown that in the boundary value problem, the temperature and density can rapidly vary near the lower boundary. An example of the solution of a three-dimensional problem with variable pressure on the Earth's surface, taken from experimental observations, is given. The developed algorithms and computer programs can be used to simulate the atmospheric waves from pressure variations on the Earth's surface.

An algorithm of the search for the best (in a sense) cubature formulas on a sphere that are invariant with respect to the transformations of the icosahedral group of rotations with inversion is described. This algorithm is applied to finding the parameters of all the best cubature formulas of this symmetry type up to the 79th order of accuracy. The parameters of the new cubature formulas of the 21st, 25th and 29th orders of accuracy to 16 significant digits are given.

In this paper we investigate the two-velocity stationary hydrodynamics system with a single pressure and inhomogeneous divergent and boundary conditions for the two velocities. This system is overdetermined. By replacing the unknown functions, the problem is reduced to a homogeneous one. The solution of the resulting system is reduced to the consecutive solutions of the two boundary value problems: the Stokes problem for a single velocity and pressure, and overdetermined system for the other velocity. We present the generalized statements of these problems and their discrete approximation using the finite element method. To solve the overdetermined problem we apply a version of the regularization methods.

We say that the square matrices A and B are of the same order quasi-commute if AB = σ BA for some scalar σ. Classical relations of commutation and anti-commutation are particular cases of this definition. We give a complete description of pairs of the quasi-commuting Toeplitz and Hankel matrices for σ ≠ ± 1.

It is proved that the uniform Lebesgue constant (the norm of a linear operator from C to C) of local cubic splines with equally-spaced knots, which preserve cubic polynomials, is equal to 11/9.