Search

The mixed problem for the telegraph equation well-known in electrical engineering and electronics, provided that the line is free from distortions, is reduced to a similar problem for one-dimensional inhomogeneous wave equation. An effective way to solve this problem is based on the use of special functions – polylogarithms, which are complex power series with power coefficients, converging in the unit circle. The exact solution of the problem is expressed in integral form in terms of the imaginary part of the first-order polylogarithm on the unit circle, and the approximate one – in the form of a finite sum in terms of the real part of the dilogarithm and the imaginary part of the third-order polylogarithm. All the indicated parts of the polylogarithms are periodic functions that have polynomial expressions of the corresponding degrees on an interval of length in the period, which makes it possible to obtain a solution to the problem in elementary functions. In the paper, we study a mixed problem for the telegrapher’s equation which is well-known in applications. This problem of linear substitution of the desired function witha time-exponential coefficient is reduced to a similar problem for the Klein – Gordon equation. The solution of the latter can be found by dividing the variables in the form of a series of trigonometric functions of a line point with time-dependent coefficients. Such a solution is of little use for practical application, since it requires the calculation of a large number of coefficients-integrals and it is difficult to estimate the error of calculations. In the present paper, we propose another way to solve this problem, based on the use of special He-functions, which are complex power series of a certain type that converge in the unit circle. The exact solution of the problem is presented in integral form in terms of second-order He-functions on the unit circle. The approximate solution is expressed in the final form in terms of third-order He-functions. The paper also proposes a simple and effective estimate of the error of the approximate solution of the problem. It is linear in relation to the line splitting step with a time-exponential coefficient. An example of solving the problem for the Klein – Gordon equation in the way that has been developed is given, and the graphs of exact and approximate solutions are constructed.

WENO schemes (weighted, essentially non oscillating) are currently having a wide range of applications as approximate high order schemes for discontinuous solutions of partial differential equations. These schemes are used for direct numerical simulation (DNS) and large eddy simmulation in the gas dynamic problems, problems for DNS in MHD and even neutron kinetics. This work is dedicated to clarify some characteristics of WENO schemes and numerical simulation of specific tasks. Results of the simulations can be used to clarify the field of application of these schemes. The first part of the work contained proofs of the approximation properties, stability and convergence of WENO5, WENO7, WENO9, WENO11 and WENO13 schemes. In the second part of the work the modified wave number analysis is conducted that allows to conclude the dispersion and dissipative properties of schemes. Further, a numerical simulation of a number of specific problems for hyperbolic equations is conducted, namely for advection equations (one-dimensional and two-dimensional), Hopf equation, Burgers equation (with low dissipation) and equations of non viscous gas dynamics (onedimensional and two-dimensional). For each problem that is implying a smooth solution, the practical calculation of the order of approximation via Runge method is performed. The influence of a time step on nonlinear properties of the schemes is analyzed experimentally in all problems and cross checked with the first part of the paper. In particular, the advection equations of a discontinuous function and Hopf equations show that the failure of the recommendations from the first part of the paper leads first to an increase in total variation of the solution and then the approximation is decreased by the non-linear dissipative mechanics of the schemes. Dissipation of randomly distributed initial conditions in a periodic domain for one-dimensional Burgers equation is conducted and a comparison with the spectral method is performed. It is concluded that the WENO7-WENO13 schemes are suitable for direct numerical simulation of turbulence. At the end we demonstrate the possibility of the schemes to be used in solution of initial-boundary value problems for equations of non viscous gas dynamics: Rayleigh-Taylor instability and the reflection of the shock wave from a wedge with the formation a complex configuration of shock waves and discontinuities.

The total probability formula for continuous random variables is the integral of product of two probability density functions that defines the unconditional probability density function from the conditional one. The need for calculation of such integrals arises in many applications, for instant, in statistical decision theory. The statistical decision theory attracts attention due to the ability to formulate the problems in a strict mathematical form. One of the technical problems solved by the statistical decision theory is the problem of dual control that requires calculation of integrals connected with the multivariate probability distributions. The necessary integrals are not available in the literature. One theorem on the total probability formula for vector Gaussian distributions was published by the authors earlier. In this paper we repeat this theorem and prove a new theorem that uses more familiar form of the initial data and has more familiar form of the result. The new form of the theorem allows us to obtain the unconditional mathematical expectation and the unconditional variance-covariance matrix very simply. We also confirm the new theorem by direct calculation for the case of the simple linear regression.

The paper studies linear differential operators in derivatives with respect to one variable. Such operators include, in particular, operators defined on infinite prolongations of evolutionary systems of differential equations with one spatial variable. In this case, differential operators in total derivatives with respect to the spatial variable are considered. In parallel, linear differential operators with one independent variable are investigated. The known algorithms for reducing the matrix to a stepwise or diagonal form are generalized to the operator matrices of both types. These generalizations are useful at points, where the functions, into which the matrix components are divided when applying the algorithm, are nonzero.In addition, the integral operator is defined as a multi-valued operator that is the right inverse of the total derivative. Linear operators that involve both the total derivatives and the integral operator are called integro-differential. An invertible operator in the integro-differential sense is an operator for which there exists a two-sided inverse integro-differential operator. A description of scalar differential operators that are invertible in this sense is obtained. An algorithm for checking the invertibility in the integro-differential sense of a differential operator and for constructing the inverse integro-differential operator is formulated.The results of the work can be used to solve linear equations for matrix differential operators arising in the theory of evolutionary systems with one spatial variable. Such operator equations arise when describing systems that are integrable by the inverse scattering method, when calculating recursion operators, higher symmetries, conservation laws and symplectic operators, and also when solving some other problems. The proposed method for solving operator equations is based on reducing the matrices defining the operator equation to a stepwise or diagonal form and solving the resulting scalar operator equations.

The creation of methods for the study of nonlinear phase systems has a long history, since the 60s of the last century (V.I. Tikhonov, V. Lindsay, M.V. Kapranov, B.I. Shakhtarin, etc.). By now, rigorous and approximate analysis methods of such systems have been developed. However, most methods are limited to the analysis of low order systems. Only in recent years attempts have been made to create methods, which allow to carry out the analysis of high-order phase systems. The material of this article deals with these methods. The article considers the construction of solutions of phase systems on the example of phase-locked frequency of arbitrary dimension with piecewise linear approximation of the nonlinear function. This approximation allowed to use an explicit form of solutions in the linearity and to obtain analytical conditions for the existence of various types of system behavior. The analytical conditions for the existence of solutions leading to the emergence of complex limit sets of the trajectories of phase systems and their bifurcations are obtained. These are homoclinic trajectories in the case of the saddle-focus equilibrium state, which play a decisive role in the occurrence of chaos. It is also shown that it is possible to obtain analytical conditions for the bifurcation of the birth and the existence of multi-pass rotational cycles in a piecewise linear phase system, on the basis of which a criterion for the transition to chaos through bifurcations cascade of doubling the stable cycle period can be obtained; in accordance with the Sharkovsky theorem it ends with the bifurcation of the cycle birth of the period three and the occurrence of developed chaos. It should be noted that the research methods of piecewise linear systems described in the paper were applied by the authors not only to phase systems, but, for example, to the Chua system, allowing various chaotic behavior.

In this paper we introduce a new class of foliations on Rie-mannian 3-manifolds, called B-foliations, generalizing the class of foliations of non-negative curvature. The leaves of B-foliations have bounded total absolute curvature in the induced Riemannian metric. We describe several topological and geometric properties of B-foliations and the structure of closed oriented 3-dimensional manifolds admitting B-foliations with non-positive curvature of leaves.

Let $M$ be a compact two-dimensional manifold and, $f \in C^{\infty}(M,\mathbb{R})$ be a Morse function, and $\Gamma_f$ be its Kronrod-Reeb graph. Denote by $\mathcal{O}_{f}=\{f \circ h \mid h \in \mathcal{D}\}$ the orbit of $f$ with respect to the natural right action of the group of diffeomorphisms $\mathcal{D}$ on $C^{\infty}(M,\mathbb{R})$, and by $\mathcal{S}(f)=\{h\in\mathcal{D} \mid f \circ h = f\}$ the corresponding stabilizer of this function. It is easy to show that each $h\in\mathcal{S}(f)$ induces a homeomorphism of $\Gamma_f$. Let also $\mathcal{D}_{\mathrm{id}}(M)$ be the identity path component of $\mathcal{D}(M)$, $\mathcal{S}'(f)= \mathcal{S}(f) \cap \mathcal{D}_{\mathrm{id}}(M)$ be group of diffeomorphisms of $M$ preserving $f$ and isotopic to identity map, and $G_f$ be the group of homeomorphisms of the graph $\Gamma_f$ induced by diffeomorphisms belonging to $\mathcal{S}'(f)$. This group is one of the key ingredients for calculating the homotopy type of the orbit $\mathcal{O}_{f}$. Recently the authors described the structure of groups $G_f$ for Morse functions on all orientable surfaces distinct from $2$-torus $T^2$ and $2$-sphere $S^2$. The present paper is devoted to the case $M=S^{2}$. In this situation $\Gamma_f$ is always a tree, and therefore all elements of the group $G_f$ have a common fixed subtree $\mathrm{Fix}(G_f)$, which may even consist of a unique vertex. Our main result calculates the groups $G_f$ for all Morse functions $f:S^{2}\to\mathbb{R}$ whose fixed subtree $\mathrm{Fix}(G_f)$ consists of more than one point., Comment: 5 pages

The study of non-stationary rarefied gas flows is, currently, attracting a great deal of attention. Such an interest arises from creating the pulsed jets used for deposition of thin films and special coatings on the solid surfaces. However, the problems of non-stationary rarefied gas flows are still understudied because of their large computational complexity. The paper considers the computational aspects of investigating non-stationary movement of gas reflected from a wall and flowing through a suddenly formed gap. The study objective is to analyse the possible numerical kinetic approaches to solve such problems and identify the difficulties in their solving. When modeling the gas flows in strong rarefaction one should consider the Boltzmann kinetic equation, but its numerical implementation is rather time-consuming. In order to use more simple approaches based, for example, on approximation kinetic equations (Ellipsoidal-Statistical model, Shakhov model), it is important to estimate the difference between the solutions of the model equations and of the Boltzmann equation. For this purpose, two auxiliary problems are considered, namely reflection of the gas flow from the wall and outflow of the free jet into the rarefied background gas.A numerical solution of these problems shows a weak dependence of the solution on the type of the collision operator in the rarefied region, but at the same time a strong dependence of a behavior of the macro-parameters on the velocity grid step. The detailed velocity grid is necessary to avoid a non-monotonous behavior of the macro-parameters caused by so-called ray effect. To reduce computational costs of the detailed velocity grid solution, a hybrid method based on the synthesis of model equations and the Boltzmann equation is proposed. Such an approach can be promising since it reduces the domain in which the Boltzmann collision integral should be used.The article presents the results obtained using two different software packages, namely a Unified Flow Solver (UFS) [13] and a Nesvetay 3D software complex [14-15]. Note that the UFS uses the discrete ordinate method for velocity space on a uniform grid and a hierarchical adaptive mesh refinement in physical space. The possibility to calculate both the Boltzmann equation and the model equations is realized. The Nesvetay 3D software complex was created to solve the Shakhov model equation (S-model) for calculations based on non-structured non-uniform grids, both in velocity space and in physical one.

A numerically stable direct multiplicative method for solving systems of linear equations that takes into account the sparseness of matrices presented in a packed form is considered. The advantage of the method is the calculation of the Cholesky factors for a positive definite matrix of the system of equations and its solution within the framework of one procedure. And also in the possibility of minimizing the filling of the main rows of multipliers without losing the accuracy of the results, and no changes are made to the position of the next processed row of the matrix, which allows using static data storage formats. The solution of the system of linear equations by a direct multiplicative algorithm is, like the solution with LU-decomposition, just another scheme for implementing the Gaussian elimination method. The calculation of the Cholesky factors for a positive definite matrix of the system and its solution underlies the construction of a new mathematical formulation of the unconditional problem of quadratic programming and a new form of specifying necessary and sufficient conditions for optimality that are quite simple and are used in this paper to construct a new mathematical formulation for the problem of quadratic programming on a polyhedral set of constraints, which is the problem of finding the minimum distance between the origin ordinate and polyhedral boundary by means of a set of constraints and linear algebra dimensional geometry. To determine the distance, it is proposed to apply the known exact method based on solving systems of linear equations whose dimension is not higher than the number of variables of the objective function. The distances are determined by the construction of perpendiculars to the faces of a polyhedron of different dimensions. To reduce the number of faces examined, the proposed method involves a special order of sorting the faces. Only the faces containing the vertex closest to the point of the unconditional extremum and visible from this point are subject to investigation. In the case of the presence of several nearest equidistant vertices, we investigate a face containing all these vertices and faces of smaller dimension that have at least two common nearest vertices with the first face.

Darboux theorem to more general context of Frechet manifolds we face an obstacle: in general vector fields do not have local flows. Recently, Fr\'{e}chet geometry has been developed in terms of projective limit of Banach manifolds. In this framework under an appropriate Lipchitz condition The Darboux theorem asserts that a symplectic manifold $(M^{2n},\omega)$ is locally symplectomorphic to $(R^{2n}, \omega_0)$, where $\omega_0$ is the standard symplectic form on $R^{2n}$. This theorem was proved by Moser in 1965, the idea of proof, known as the Moser’s trick, works in many situations. The Moser tricks is to construct an appropriate isotopy $ \ff_t $ generated by a time-dependent vector field $ X_t $ on $M$ such that $ \ff_1^{*} \omega = \omega_0$. Nevertheless, it was showed by Marsden that Darboux theorem is not valid for weak symplectic Banach manifolds. However, in 1999 Bambusi showed that if we associate to each point of a Banach manifold a suitable Banach space (classifying space) via a given symplectic form then the Moser trick can be applied to obtain the theorem if the classifying space does not depend on the point of the manifold and a suitable smoothness condition holds. If we want to try to generalize the local flows exist and with some restrictive conditions the Darboux theorem was proved by Kumar. In this paper we consider the category of so-called bounded Fr\'{e}chet manifolds and prove that in this category vector fields have local flows and following the idea of Bambusi we associate to each point of a manifold a Fr\'{e}chet space independent of the choice of the point and with the assumption of bounded smoothness on vector fields we prove the Darboux theorem.

In 1994, in [13], N. Papaghiuc introduced the notion of semi-slant submanifold in a Hermitian manifold which is a generalization of CR- and slant-submanifolds. In particular, he considered this submanifold in Kaehlerian manifolds, [13]. Then, in 2007, V. A. Khan and M. A. Khan considered this submanifold in a nearly Kaehler manifold and obtained interesting results, [11]. Recently, we considered semi-slant submanifolds in a locally conformal Kaehler manifold and gave a necessary and sufficient conditions for two distributions (holomorphic and slant) to be integrable. Moreover, we considered these submanifolds in a locally conformal Kaehler space form, [4]. In this paper, we define 2-kind warped product semi-slant submanifolds in a locally conformal Kaehler manifold and consider some properties of these submanifolds.

This paper considers spectral-difference methods of a high-order of accuracy for solving the one-way wave equation using the Laguerre integral transform with respect to time as the base. In order to provide a high spatial accuracy and calculation stability, the Richardson method can be employed. However such an approach requires high computer costs, therefore we consider alternative algorithms based on the Adams multistep schemes. To reach the stability first for the 1D one-way equation and then for the 2D case, the stabilizing procedures using the spline interpolation were developed. This made possible to efficiently implement a predictor-corrector type method in terms of which a boundary value problem for high-order elliptic equations is substituted for a sequence of inversions of second order elliptic operators thus decreasing computer costs. To assess the accuracy and stability of difference approximations for the 1D one-way wave equation, the analytical solution based on the double Laguerre transform was obtained. This solution can be efficiently calculated if for summation one makes use of fast algorithms of computing a discrete linear convolution. For the 2D one-way wave equation the stability and accuracy of the procedures proposed have been studied on implementing the migration algorithm within a problem of seismic prospecting., In Russian. arXiv admin note: substantial text overlap with arXiv:1605.03292

The paper introduces and investigates partial orders that are finer than the minus partial order and are based on inverses along an element and other specific outer inverses. It turns out that in this way a number of classical partial orders can be equivalently defined.

International audience; http://mech.math.msu.su/~fpm/ps/k14/k143/k14303.pdfIn this paper, a computerized model for morphological analysis of languages with word-formation based on affixation processes is proposed. The main idea consists in defining structural patterns of words and corresponding lists of suffixes. First, a detaileddescription of a stemming algorithm, its modification, and the technique of determining grammatical characteristics of word-forms are given. The next part of this work focuses on the application of the proposed algorithms for the French language. Finally, some results of execution of these algorithms are provided.; http://mech.math.msu.su/~fpm/ps/k14/k143/k14303.pdfВ работе предлагается метод структурных схем в качестве модели морфологического анализа словоформ естественного языка с развитым аффиксальным словообразованием и словоизменением. Дано описание алгоритма выделения псевдоосновы, его модификация, а также алгоритм восстановления грамматических характеристик словоформ. Описано применение предложенного метода для анализа словоформ французского языка. Представлены результаты работы предложенных алгоритмов.