Vol 25, No 4 (2023)
- Year: 2023
- Published: 24.12.2025
- Articles: 10
- URL: https://journal-vniispk.ru/2079-6900/issue/view/24372
Full Issue
Mathematics
On Logarithmic Hölder Condition and Local Extrema of Power Takagi Functions
Abstract
This paper studies one class of real functions, which we call Takagi power functions. Such functions have one positive real parameter; they are continuous, but nowhere differentiable, and are given on a real line using functional series. These series are similar to the series defining the continuous, nowhere differentiable Takagi function described in 1903. For each parameter value, we derive a functional equation for functions related to Takagi power functions. Then, using this equation, we obtain an accurate two-sides estimate for the functions under study. Next, we prove that for parameter values not exceeding 1, Takagi power functions satisfy the Hölder logarithmic condition, and find the smallest value of the constant in this condition. As a result, we get the usual Hölder condition, which follows from the logarithmic Hölder condition. Moreover, for parameter values ranging from 0 to 1, we investigate the behavior of Takagi power functions in the neighborhood of their global maximum points. Then we show that the functions under study reach a strict local minimum on the real axis at binary-rational points, and only at them. Finally, we describe the set of points at which our functions reach a strict local maximum. The benefit of our research lies in the development of methods applicable to continuous functions that cannot be differentiated anywhere. This can significantly expand the set of functions being studied.
223-241
Estimating the Lebesgue constant for the Chebyshev distribution of nodes
Abstract
In this paper an approach to estimation of the Lebesgue constant for the Lagrange interpolation process with nodes in the zeros of Chebyshev polynomials of the first kind is done. Two-sided estimation of this constant is carried out by using the logarithmic derivative of the Euler gamma function and of the Riemann zeta function. The choice of interpolation nodes is due to the fact that with a fixed number of Chebyshev nodes, the Lebesgue constant tends to its minimum value, thus reducing the error of algebraic interpolation and providing less sensitivity to rounding errors. The expressions for the upper and the lower bounds of this constant are represented as finite sums of an asymptotic alternating series. Based on the expressions obtained, these boundaries are calculated depending on the number of nodes of the interpolation process. The error of each of the boundaries’ value is estimated based on the first discarded term in the corresponding asymptotic series. The results of the calculations are presented in tables showing deviations of the Lebesgue constant from its lower and upper estimated bounds. Dependence of the values’ errors on the number of Chebyshev nodes is depicted in these tables as well. It is numerically shown that with an increase in the number of these nodes, the estimation boundaries rapidly get close to each other. The presented results can be used in the theory of interpolation to estimate the norm of the operator matching a function to its interpolation polynomial and to estimate a deviation of the constructed perturbed polynomial from the unperturbed one.
242-254
Numerical Study of the Rate of Convergence of Chernoff Approximations to Solutions of the Heat Equation
Abstract
The article is devoted to construction of examples illustrating (using computer calculation) the convergence rate of Chernoff approximations to the solution of the Cauchy problem for the heat equation. Two Chernoff functions (of the first and second order of Chernoff tangency to the double differentiation operator) and several initial conditions of different smoothness are considered. As an illustration for the initial condition equal to the absolute value of the sine function to the power of five over two, a graph of the exact solution of the Cauchy problem and graphs of the tenth Chernoff approximations given by two different Chernoff functions are plotted. It is visually determined that the approximations are close to the solution. For each of the two Chernoff functions, for several initial conditions of different smoothness and for the approximation numbers up to 11, the error corresponding to each approximation is numerically found. This error is understood as the supremum of the absolute value of the difference between the exact solution and its approximating function. As it turned out, in all the cases studied, the error dependence on the approximation number nearly power-law form. This follows from the fact that the dependence of the error logarithm on the logarithm of the approximation number is close to linear. By finding the equation of the approximating line using linear regression, we find the exponent in the power dependence of the error on the approximation number and call it the order of convergence. These orders for all studied initial conditions are compiled in a table. The empirical dependence of the convergence order on the smoothness class of the initial condition is found on the considered family of initial conditions.
255-272
Superstructures over Cartesian products of orientation-preserving rough circle transformations
Abstract
One of the constructions for obtaining flows on a manifold is building a superstructure over a cascade. In this case, the flow is non-singular, that is, it has no fixed points. C. Smale showed that superstructures over conjugate diffeomorphisms are topologically equivalent. The converse statement is not generally true, but under certain assumptions the conjugacy of diffeomorphisms is tantamount to equivalence of superstructures. Thus, J. Ikegami showed that the criterion works in the case when a diffeomorphism is given on a manifold whose fundamental group does not admit an epimorphism into the group Z. He also constructed examples of non-conjugate diffeomorphisms of a circle whose superstructures are equivalent. In the work of I. V. Golikova and O. V. Pochinka superstructures over diffeomorphisms of circles are examined. It is also proven in this paper that the complete invariant of the equivalence of superstructures over orientation-preserving diffeomorphisms is the equality of periods for periodic points generating their diffeomorphisms. For the other side, it is known from the result of A.G. Mayer that the coincidence of rotation numbers is also necessary for conjugacy of orientation-preserving diffeomorphisms. At the same time, superstructures over orientation-changing diffeomorphisms of circles are equivalent if and only if the corresponding diffeomorphisms of circles are topologically conjugate. Work of S. Kh. Zinina and P. I. Pochinka proved that superstructures over orientation-changing Cartesian products of diffeomorphisms of circles are equivalent if and only if the corresponding diffeomorphisms of tori are topologically conjugate. In this paper a classification result is obtained for superstructures over Cartesian products of orientation-preserving diffeomorphisms of circles.
273-283
On the Similarity over the Ring of Integers of Certain Nilpotent Matrices of Maximal Rank
Abstract
This paper is devoted to the problem of matrix similarity recognition over the ring of integers for some families of matrices. Namely, nilpotent upper triangular matrices of maximal rank are considered such that only first and second superdiagonals of these matrices are non-zero. Several necessary conditions are obtained for similarity of such matrices to matrices of the form superdiag(a1,a2,...,an-1) with a single non-zero superdiagonal, that is a generalization of the Jordan cell Jn(0)=superdiag(1,1,...,1). These conditions are formulated in simple terms of divisibility and greatest common divisors of matrix elements. The result is obtained by reducing the problem of similarity recognition to the problem of solving in integers a system of linear equations and applying the known necessary similarity conditions for arbitrary matrices. Under some additional conditions on the elements a1,a2,...,an-1 of the first superdiagonal of matrix A, it is proven that the matrix A is similar to matrix superdiag(a1,a2,...,an-1) regardless of the values of the elements of the second superdiagonal. Moreover, for the considered matrices of the third and the fourth orders, easily verifiable necessary and sufficient similarity conditions are obtained describing their similarity to a matrix of the form superdiag(a1,a2,...,an-1).
284-298
On the stability of a nonlinear nonautonomous scalar equation with variable delay
Abstract
The stability problem of a scalar functional differential equation is a classical one. It has been most fully studied for linear equations. Modern research on modeling biological, infectious and other processes leads to the need to determine the qualitative properties of the solutions for more general equations. In this paper we study the stability and the global limit behavior of solutions to a nonlinear one-dimensional (scalar) equation with variable delay with unbounded and bounded right-hand sides. In particular, our research is reduced to a problem on the stability of a non-stationary solution of a nonlinear scalar Lotka-Volterra-type equation, on the stabilization and control of a non-stationary process described by such an equation. The problem posed is considered depending on the delay behavior: is it a bounded differentiable function or a continuous and bounded one. The study is based on the application of the Lyapunov-Krasovsky functionals method as well as the corresponding theorems on the stability of non-autonomous functional differential equations of retarded type with finite delay. Sufficient conditions are derived for uniform asymptotic stability of the zero solution, including global stability, for every continuous initial function. Using the theorem proven by one of the co-authors on the limiting behavior of solutions to a non-autonomous functional differential equation based on the Lyapunov functional with a semidefinite derivative, the properties of the solutions’ attraction to the set of equilibrium states of the equation under study are obtained. In addition, illustrative examples are provided.
299-312
Applied mathematics and mechanics
Nonlinear mathematical model of pressure measurement systems in gas-liquid media
Abstract
The primary element of the instrumentation for measuring the pressure of a gas-liquid medium is a sensor that supplies data on the pressure of the working medium. It determines the proper functioning of machines, mechanisms, and systems. Increasing the service life, reducing development time, and reducing the cost of sensors is one of the important tasks. Mathematical modeling of pressure measurement systems’ functioning plays an important role at the design stage of such systems. This article examines a nonlinear one-dimensional model of a mechanical system “pipeline – pressure sensor” designed to measure and control the pressure of the working gas-liquid medium in the combustion chambers of engines. In such a system, the sensor is connected to the engine via a pipeline and is located at some distance from it to reduce the impact of vibration accelerations and high temperatures. The purpose of the work is to study the dynamics and stability of joint oscillations of the elastic sensitive element in the pressure sensor and of the working medium in the pipeline for a given law of pressure change in the combustion chamber. The study is provided under the assumption that the working medium is ideal and compressible. To describe the movement of the working medium (gas or liquid), a nonlinear model of fluid and gas mechanics is used. Mathematical description of the process of interest includes an initial boundary value problem, whose formulation contains a nonlinear partial differential equation. To solve it, numerical-analytical method of solution based on the Galerkin method is proposed, which makes it possible to reduce the study of the problem to solving a system of ordinary differential equations. A numerical experiment is carried out and examples of calculating the dynamics of the sensor's sensitive element are presented. The proposed mathematical model makes it possible to determine the law of change in the deviation of the sensor's sensitive element depending on the law of change in pressure in the combustion chamber. The research results are intended for use at the design stage of pressure measurement systems.
313-325
Nonholonomic mechanical systems on a plane with a variable slope
Abstract
This paper considers such nonholonomic mechanical systems as Chaplygin skate, inhomogeneous Chaplygin sleigh and Chaplygin sphere moving in the gravity field along an oscillating plane with a slope varying with the periodic law. By explicit integration of the equations of motion, analytical expressions for the velocities and trajectories of the contact point for Chaplygin skate and Chaplygin sleigh are obtained. Numerical parameters of the periodic law for the inclination angle change are found, such that the velocity of Chaplygin skate will be unbounded, that is, an acceleration will take place. In the case of inhomogeneous Chaplygin sleigh, on the contrary, numerical parameters of the periodic law of the inclination angle change are found, for which the sleigh velocity is bounded and there is no drift of the sleigh. For similar numerical parameters and initial conditions, when the sleigh moves along a horizontal or inclined plane with the constant slope, the velocity and trajectory of the contact point are unbounded, that is, there is a drift of the sleigh. A similar problem is solved for the Chaplygin sphere; its trajectories are constructed on the basis of numerical integration. The results are illustrated graphically. The control of the slope of the plane, depending on the angular momentum of the sphere, is proposed for discussion. Regardless of the initial conditions, such control can almost always prevent the drift of the sphere in one of the directions.
326-341
Mathematical modeling and computer science
Investigation of different influence functions in peridynamics
Abstract
Peridynamics is a non–local numerical method for solving fracture problems based on integral equations. It is assumed that particles in a continuum are endowed with volume and interact with each other at a finite distance, as in molecular dynamics. The influence function in peridynamic models is used to limit the force acting on a particle and to adjust the bond strength depending on the distance between the particles. It satisfies certain continuity conditions and describes the behavior of non-local interaction. The article investigates various types of influence function in peridynamic models on the example of three-dimensional problems of elasticity and fracture. In the course of the work done, the bond-based and state-based fracture models used in the Sandia Laboratory are described, 6 types of influence functions for the bond-based model and 2 types of functions for the state-based model are presented, and the corresponding formulas for calculating the stiffness of the bond are obtained. For testing, we used the problem of propagation of a spherically symmetric elastic wave, which has an analytical solution, and a qualitative problem of destruction of a brittle disk under the action of a spherical impactor. Graphs of radial displacement are given, raster images of simulation results are shown.
342-360
Numerical model of vapor-air-diesel autothermal reformer
Abstract
The use of fuel cell power plants is a promising area in the generation of electricity. The path to their widespread use is hindered by their high cost and the availability of the fuel used. To solve this problem, effective energy conversion systems operating on diesel fuel are being developed. The main goal is to create a device (a fuel processor), which would convert diesel fuel into a hydrogen-containing gas. The device consists of several components: a nozzle for injecting liquid fuel having the form of drops into superheated steam, a mixing and vaporization zone for diesel fuel, an air supply area, and a reaction zone including a catalyst. The selection of temperature for the vaporization process should be made in such a way that, on the one hand, liquid droplets do not come into contact with the catalyst surface, and, on the other hand, gas-phase reactions are not initiated in the mixing zone. Developing such a device requires not only conducting laboratory experiments and studying the process catalyst, but also optimizing the basic physical parameters of the device. These parameters are its linear dimensions, operational temperature, reactant flow rates, and many others. Carrying out such a study is impossible without using methods of mathematical modeling. This significantly reduces the time and cost of work. This paper presents a digital model of an air-hydrocarbon mixture generator in an axisymmetric formulation. The dynamics of subsonic multiphase flow of water vapor carrying drops of liquid diesel fuel, the process of diesel fuel evaporation and mixing with water vapor and air are studied. The mathematical model is implemented in the ANSYS Fluent package (academic license of SSCC SBRAS). A series of calculations for various mixture feed temperatures are performed to optimize the main parameters. For the established optimal temperature, modeling of the mixture mixing process with air is carried out.
361-374

