


Том 58, № 7 (2018)
- Жылы: 2018
- Мақалалар: 15
- URL: https://journal-vniispk.ru/0965-5425/issue/view/11179
Article
Integro-Differential Polynomial and Trigonometrical Splines and Quadrature Formulas
Аннотация
This work is one of many that are devoted to the further investigation of local interpolating polynomial splines of the fifth order approximation. Here, new polynomial and trigonometrical basic splines are presented. The main features of these splines are the following: the approximation is constructed separately for each grid interval (or elementary rectangular), the approximation constructed as the sum of products of the basic splines and the values of function in nodes and/or the values of its derivatives and/or the values of integrals of this function over subintervals. Basic splines are determined by using a solving system of equations which are provided by the set of functions. It is known that when integrals of the function over the intervals is equal to the integrals of the approximation of the function over the intervals then the approximation has some physical parallel. The splines which are constructed here satisfy the property of the fifth order approximation. Here, the one-dimensional polynomial and trigonometrical basic splines of the fifth order approximation are constructed when the values of the function are known in each point of interpolation. For the construction of the spline, we use the discrete analogues of the first derivative and quadrature with the appropriate order of approximation. We compare the properties of these splines with splines which are constructed when the values of the first derivative of the function are known in each point of interpolation and the values of integral over each grid interval are given. The one-dimensional case can be extended to multiple dimensions through the use of tensor product spline constructs. Numerical examples are represented.



The Bauer-Type Factorization of Matrix Polynomials Revisited and Extended
Аннотация
For a Laurent polynomial \(a(\lambda )\), which is Hermitian and positive definite on the unit circle, the Bauer method provides the spectral factorization \(a(\lambda ) = p(\lambda )p{\kern 1pt} {\text{*}}({{\lambda }^{{ - 1}}})\), where \(p(\lambda )\) is a polynomial having all its roots outside the unit circle. Namely, as the size of the banded Hermitian positive definite Toeplitz matrix associated with the Laurent polynomial increases, the coefficients at the bottom of its Cholesky lower triangular factor tend to the coefficients of \(p(\lambda )\). We study extensions of the Bauer method to the non-Hermitian matrix case. In the Hermitian case, we give new convergence bounds with computable coefficients.



A New Proof of the Kuhn–Tucker and Farkas Theorems
Аннотация
For the minimization problem for a differentiable function on a set defined by inequality constraints, a simple proof of the Kuhn–Tucker theorem in the Fritz John form is presented. Only an elementary property of the projection of a point onto a convex closed set is used. The approach proposed by the authors is applied to prove Farkas’ theorem. All results are extended to the case of Banach spaces.



Method for Constructing Optimal Dark Coverings
Аннотация
The problem of constructing metric ε-nets and corresponding coverings by balls for compact sets with a probability measure is considered. In the case of sets having metrically significant parts with a small measure (dark sets), methods for constructing ε-nets are combined with the deep holes method in a unified approach. According to this approach, a constructed metric net is supplemented with its deep hole (the most distant element of the set) until the required accuracy is achieved. An existing implementation of the method for a metric set with a given probability measure is based on a pure global search for deep holes. To construct dark coverings, the method is implemented on the basis of a random multistart. For the resulting nets, the logarithm of the number of their elements is shown to be close to ε-entropy, which means that they are optimal. Techniques for estimating the reliability and completeness of constructed (ε, δ)-coverings in the sense of C.E. Shannon are described. The methods under consideration can be used to construct coverings of implicitly given sets with a measure defined on the preimage and to recover compact supports of multidimensional random variables with an unknown distribution law.



Affine Controlled Systems and t-Systems of Pfaffian Equations
Аннотация



Development and Application of the Fourier Method for the Numerical Solution of Ito Stochastic Differential Equations
Аннотация
This paper is devoted to the development and application of the Fourier method to the numerical solution of Ito stochastic differential equations. Fourier series are widely used in various fields of applied mathematics and physics. However, the method of Fourier series as applied to the numerical solution of stochastic differential equations, which are proper mathematical models of various dynamic systems affected by random disturbances, has not been adequately studied. This paper partially fills this gap.



Dynamic Reconstruction of Disturbances in a Quasilinear Stochastic Differential Equation
Аннотация
The problem of reconstructing unknown inputs in a first-order quasilinear stochastic differential equation is studied by applying dynamic inversion theory. The disturbances in the deterministic and stochastic terms of the equation are simultaneously reconstructed using discrete information on some realizations of the stochastic process. The problem is reduced to an inverse one for ordinary differential equations satisfied by the expectation and variance of the original process. A finite-step software implementable solution algorithm is proposed, and its accuracy with respect to the number of measured realizations is estimated. An illustrative example is given.






Estimation of Two Error Components in the Numerical Solution to the Problem of Nonisothermal Flow of Polymer Fluid between Two Coaxial Cylinders
Аннотация
An algorithm for solving a stationary nonlinear problem of a nonisothermal flow of an incompressible viscoelastic polymer fluid between two coaxial cylinders is developed on the basis of Chebyshev approximations and the collocation method. In test calculations, the absence of saturation of the algorithm is shown. A posteriori estimates of two error components in the numerical solution—the error of approximation method and the round-off error—are obtained. The behavior of these components as a function of the number of nodes in the spatial grid of the algorithm and the radius of the inner cylinder is analyzed. The calculations show exponential convergence, stability to rounding errors, and high time efficiency of the algorithm developed.



Reflection of a Rarefaction Wave from the Center of Symmetry: Theoretical Analysis of the Flow Features and Calculation by the Method of Characteristics
Аннотация
For an unsteady spherical rarefaction wave in an ideal (inviscid and non-heat-conducting) gas, the features of the flow near the reflection of the first characteristic from the center of symmetry are investigated. Computations performed by the method of characteristics on nearly uniform grids usually used in such problems reveal sawtooth irregularities in the parameter distributions near the reflection point, whereas similar irregularities in the cylindrical and plane cases are absent. The amplitudes of the irregularities and the sizes of the domains where they are observed remain nearly unchanged when the number of points of the characteristic grid is increased by many times. Away from the reflection point in both time and space, the numerical solution completely “forgets” about the irregularities. This finding explains why these irregularities were ignored or overlooked earlier, but the nature of this phenomenon remains an open question. The present study has established that the spherical rarefaction flow near the reflection point differs fundamentally in structure from its plane and cylindrical counterparts. In the spherical case, the rarefaction flow near the reflection point was found to be nearly conical (entirely conical in the linear approximation). Allowance for this feature in the method of characteristics led to continuous regular distributions of the parameters. The performed analysis and computations revealed that a spherical rarefaction wave is strengthened (cumulates) theoretically unlimitedly in a small neighborhood of the reflection point (center of symmetry) of the first characteristic. Moreover, the claim that a gradient catastrophe occurs in this neighborhood was found to be untenable.



Solution of a Boundary Value Problem for Velocity-Linearized Navier–Stokes Equations in the Case of a Heated Spherical Solid Particle Settling in Fluid
Аннотация
Assuming that the fluid viscosity is an exponential-power function of temperature, a boundary value problem for the Navier–Stokes equations linearized with respect to velocity is solved and the uniqueness of the solution is proved. The problem of a nonuniformly heated spherical solid particle settling in fluid is considered as an application.



On the Dispersion Curves of Anisotropic Waveguides
Аннотация
Behavior of the dispersion curves of waveguides with anisotropic filling is considered. The existence of backward and complex waves at certain values of the anisotropy coefficients is established. The region of localization of the singular points of the dispersion curves in which the generation of backward waves takes place is found.



Asymptotics of the Deflection of a Cruciform Junction of Two Narrow Kirchhoff Plates
Аннотация
Two two-dimensional plates with bending described by Sophie Germain’s equation with the biharmonic operator are joined in the form of a cross with clamped ends, but with free lateral sides outside the cross core. Asymptotics of the deflection of the junction with respect to the relative width of the plates regarded as a small parameter is constructed and justified. The variational-asymptotic model includes a system of two ordinary differential equations of the fourth and second orders with Dirichlet conditions at the endpoints of the one-dimensional cross and the Kirchhoff transmission conditions at its center. They are derived by analyzing the boundary layer near the crossing of the plates and mean that the deflection and the angles of rotation at the central point are continuous and that the total bending force and the principal torques vanish. Possible generalizations of the asymptotic analysis are discussed.



Dynamics and Stability of Air Bubbles in a Porous Medium
Аннотация
A numerical method is developed for reliably computing the evolution of the boundary of a multiply connected water-saturated domain with air bubbles in the case when the pressure inside them depends on the bubble volume. It is assumed that the distance between the gas bubbles is comparable with their size. Gas bubbles can be near an extended phase transition boundary separating a porous medium flow and a domain saturated with a mixture of air and water vapor. The numerical method is verified by comparing the numerical solution of a test problem with its analytical solution. Caused by finite-amplitude perturbations of the phase interface, the deformation of an air bubble in an extended horizontal water-saturated porous layer with a constant pressure gradient is studied. It is shown that the instability of the bubble boundary with respect to finite perturbations leads to the splitting of the bubble. An analysis of the numerical solution shows that, although all circular bubbles move at the same velocity irrespective of their size, nevertheless, due to instability, a portion of the bubble boundary where the air displaces the fluid moves faster than an opposite portion where the fluid displaces the air. As a result, nearby bubbles are capable of merging before splitting.



On the Number of Roots of Boolean Polynomials
Аннотация
This paper continues the series of studies of Boolean polynomials (Zhegalkin polynomials or algebraic normal forms (ANFs)). The investigation of properties of Boolean polynomials is an important topic of discrete mathematics and combinatorial analysis. Theoretical results in this field have wide practical applications. For example, a number of popular public-key cryptosystems are based on Reed–Muller codes, and the representation of these codes, as well as encryption and decryption algorithms are based on Boolean polynomials. The spectral properties of such codes are determined by the number of zeros of Boolean polynomials and are analyzed using the randomization lemma. Since the problem of determining the number of zeros Zg of the Boolean polynomial g(x) is NP-hard, the algorithms taking into account the “combinatorial structure of the polynomial” (even though they are search algorithms) are of practical interest. In this paper, such an algorithm based on the properties of the monomial matrix is proposed. For various types of polynomials, exact formulas for the number of roots and the expectation of the number of roots are obtained. A subclass of Boolean polynomials consisting of polynomials with separating variables is considered. A result that can be considered as a generalization of the randomization lemma is proved. The theoretical results provide a basis for estimating the applicability of polynomials for solving various practical problems.


