Bundle bispectrality for matrix differential equations

We consider the fundamental solutions of a wide class of first order systems with polynomial dependence on the spectral parameter and rational matrix potentials. Such matrix potentials are rational solutions of a large class of integrable nonlinear equations, which play an important role in different mathematical physics problems. The concept of bispectrality, which was originally introduced by Grünbaum, is extended in a natural way for the systems under consideration and their bispectrality is derived via the representation of the fundamental solutions. This bispectrality is preserved under the flows of the corresponding integrable nonlinear equations. For the case of Dirac type (canonical) systems the complete characterization of the bispectral potentials under consideration is obtained in terms of the system's spectral function.     

Mechanical oscillators described by a system of differential-algebraic equations

The classical framework for studying the equations governing the motion of lumped parameter systems presumes one can provide expressions for the forces in terms of kinematical quantities for the individual constituents. This is not possible for a very large class of problems where one can only provide implicit relations between the forces and the kinematical quantities. In certain special cases, one can provide non-invertible expressions for a kinematical quantity in terms of the force, which then reduces the problem to a system of differential-algebraic equations.     

On the asymptotic solutions of the equations of motion of mechanical systems

The existence of solutions of the equations of motion of non-natural mechanical systems of a certain form which tend to a position of equilibrium when the time increases without limit is proved by the methods described in /1, 2/. The corresponding instability theorem is proved.     

非线性振动系统周期解的数值分析

用直接数值积分法求非线性振动系统的周期解,求解时对初始条件进行迭代,使它与终点条件相一致.积分时间区间(即周期)或运动方程中的某些参数,也可在迭代过程中随同变化,积分方法是变步长的. 用这种“打靶”法求周期解,所需计算工作量相对较少.其中误差主要来源于数值积分,故不难估计并控制它足够小.这种方法可处理各种类型的振动问题,如单自由度和多自由度系统的自由无阻尼振动、强迫振动、自激振动和参数振动等等;也能求得不稳定解和那些对参数变动十分敏感的解.解的稳定性根据相关的周期系数微分方程来研究.求共振曲线或其他振动特性曲线时,利用插值方法并自动调节步长来定出迭代始值. 为了阐明这种方法的通用性,计算了若干例子.非线性的描述可用解析函数或任何其他形式,例如分段线性函数.文中还就所得周期解指出了非线性振动的一些值得注意的性质.部分计算结果与已有的近似解或实验结果作了比较.     

The dynamics of systems with large gyroscopic forces and the realization of constraints

Lagrangian systems with a large multiplier N on the gyroscopic terms are considered. Simplified equations of motion of general form with holonomic constraints are obtained in the first approximation with respect to the small parameter ɛ = 1/N. The structure of the solutions of the precessional equations is examined.     

Investigation of the problem of the flow of a compressible viscous fluid past a bounded body in the hölder classes of functions

It is shown that the stationary exterior problem for the equations of motion of a viscous compressible fluid is uniquely solvable in weighted Hölder spaces if the exterior forces and the value of the velocity at infinity are sufficiently small. As a weight function, a power function (1+|x|)^?m, m>0, is taken. The proof, which is carried out by the method of decomposition, relies on the estimates of singular integrals and solutions of the transport equations in weighted spaces. Bibliography: 5 titles.     

A computational method for solving optimal control of a system of parallel beams using Legendre wavelets

The optimal control of transverse vibration of two Euler–Bernoulli beams coupled in parallel by discrete springs is considered. An index of performance is formulated which consists of a modified energy functional of two coupled structures at a specified time and penalty functions involving the point control forces. The minimization of the performance index over these forces is subject to the equation of motion governing the structural vibrations, the imposed initial condition as well as the boundary conditions. By use of the modal space technique, the optimal control of distributed parameter systems is simplified into the optimal control of a linear time-invariant lumped-parameter systems. A computationally attractive method based on Legendre wavelets in time domain for solving the optimal control of the lumped parameter systems for any finite interval is proposed. Legendre wavelet integral operational matrix and the properties of a Kronecker product are used to find the approximated optimal trajectory and optimal law of the linear systems with respect to a quadratic cost function by only solving a linear system of algebraic equations. This method provides a straightforward and convenient approach for digital computation. A numerical example is provided to demonstrate the applicability and effectiveness of the proposed method.     

On approximate methods of analysing certain singularly-perturbed systems

A certainclass of sigularly-perturbed systems which have a variety of m-dimensional stationary positions is considered. When a small parameter disappears, the system also has an m-dimensional manifold of stationary positions and, therefore, the corresponding characteristic equation has m zero roots. The conditions under which the solution of a stability problem reduces to the same problem for a degenerate system are defined. As an application in practice gyroscopic stabilizing systems (the critical case corresponds to such systems) with elastic elements of high stiffness are discussed. The conditions under which the solution of the problem of the stability of steady motion follows from the solution of this problem for an ideal system (with absolutely rigid elements) are obtained. The problem of the closeness of the corresponding solutions of the complete and a simplified system of differential equations over an infinite time interval is discussed.     

Stabilization of the motions of non-autonomous mechanical systems

The problem of stabilizing the motions of mechanical systems that can be described by non-autonomous systems of ordinary differential equations is considered. The sufficient conditions for stabilizing of the motions of mechanical systems with assigned forces due to forces of another structure are obtained by constructing a vector Lyapunov function and a reference system. Examples of the solution of the problems of stabilizing the rotational motion of an axisymmetric satellite in an elliptic orbit, a non-tumbling gyro horizon, etc. are considered ©2009     

Time periodic solutions to the Navier–Stokes equations in the rotational framework

We consider the Navier–Stokes equations in the rotational framework with the time periodic external force. We give sufficient conditions on the size of the external forces for the existence of time periodic solutions in terms of the Coriolis parameter. It follows from our conditions that the unique existence of time periodic solutions is guaranteed for large external forces provided the speed of rotation is sufficiently fast.     

Existence and zero‐electron‐mass limit of strong solutions to the stationary compressible Navier‐Stokes‐Poisson equation with large external force

In this study, we consider a viscous compressible model of plasma and semiconductors, which is expressed as a compressible Navier‐Stokes‐Poisson equation. We prove that there exists a strong solution to the boundary value problem of the steady compressible Navier‐Stokes‐Poisson equation with large external forces in bounded domain, provided that the ratio of the electron/ions mass is appropriately small. Moreover, the zero‐electron‐mass limit of the strong solutions is rigorously verified. The main idea in the proof is to split the original equation into 4 parts, a system of stationary incompressible Navier‐Stokes equations with large forces, a system of stationary compressible Navier‐Stokes equations with small forces, coupled with 2 Poisson equations. Based on the known results about linear incompressible Navier‐Stokes equation, linear compressible Navier‐Stokes, linear transport, and Poisson equations, we try to establish uniform in the ratio of the electron/ions mass a priori estimates. Further, using Schauder fixed point theorem, we can show the existence of a strong solution to the boundary value problem of the steady compressible Navier‐Stokes‐Poisson equation with large external forces. At the same time, from the uniform a priori estimates, we present the zero‐electron‐mass limit of the strong solutions, which converge to the solutions of the corresponding incompressible Navier‐Stokes‐Poisson equations.     

On the Bifurcation and Stability of Rigidly Rotating Inviscid Liquid Bridges

Summary. We consider a mathematical model that describes the motion of an ideal fluid of finite volume that forms a bridge between two fixed parallel plates. Most importantly, this model includes capillarity effects at the plates and surface tension at the free surface of the liquid bridge. We point out that the liquid can stick to the plates due to the inner pressure even in the absence of adhesion forces. We use both the Hamiltonian structure and the symmetry group of this model to perform a bifurcation and stability analysis for relative equilibrium solutions. Starting from rigidly rotating, circularly cylindrical fluid bridges, which exist for arbitrary values of the angular velocity and vanishing adhesion forces, we find various symmetry-breaking bifurcations and prove corresponding stability results. Either the angular velocity or the angular momentum can be used as a bifurcation parameter. This analysis reduces to find critical points and corresponding definiteness properties of a potential function involving the respective bifurcation parameter. Received June 21, 1996; revision received October 2, 1997, and accepted for publication October 9, 1997     

An elliptic averaging method for harmonically excited oscillators with a purely non-linear non-negative real-power restoring force

Harmonically excited oscillators with a purely non-linear non-negative real-power restoring force are considered in this paper. The solution for motion is assumed in the form of a Jacobi cn elliptic function, the frequency and the parameter of which are obtained from the exact period of motion calculated from the energy conservation law. The parameter of the elliptic function is found to be negative for under-linear restoring forces and positive for over-linear restoring forces. By taking into account the time variation of the parameter of the elliptic function, a new elliptic averaging method is developed. The use of the equations derived is illustrated on the examples of oscillators with van der Pol damping and linear viscous damping with various integer and non-integer powers of the restoring force. New insights into dynamics of these oscillators are engendered. Numerical confirmations of analytical results are provided.     

非完整系统有耗散和循环力时的平衡不稳定性

讨论定常非完整系统在耗散、保守、循环力作用下的不稳定平衡问题．应用方法是基于运动微分方程解的存在性,当t→-∞时,系统渐近地趋于平衡状态．假定在平衡位置附近,动能、Reyleigh耗散函数、位置力都是无限可微函数．结果将通过一个实例说明．部分结果参见Kozlov V V. On the asymptotic motions of systems with dissipation. Prikl Math Mekh, 1994, 58(4): 31-36. (in Russian); Merkin D R. Introduction to the Theory of the Stability of Motion. Moscow: Nauka, 1987. (in Russian); Thomson W, Tait P. Treatise on Natural Philosophy, Part Ⅰ. Cambridge: Cambridge University Press, 1879．     

Elasticity problems involving coupled half-planes

A general method is proposed for reducing problems concerning cracks, cuts, inclusions and interacting blocks in coupled half-planes to complex integral equations, both singular and hyper-singular. The method is based on the fact that if the Kolosov-Muskhelishvili functions are known for a whole plane, then the corresponding functions for coupled half-planes are obtained from them by simple transformations. Boundary integral equations (BIE) are presented, as well as fundamental solutions for isolated forces and periodic systems of forces, which may be used to construct new complex BIEs.     

On the relationship between solutions of delay differential equations and infinite-dimensional systems of differential equations

We study the limit properties of solutions for a class of systems of ordinary differential equations as the number of equations and a certain parameter grow unboundedly. We show that the sequence of functions formed by the last components of solutions of such systems has a repeated limit. The limit function is a solution of a delay differential equation. Estimates of the convergence rate are obtained.     

On variational principles in which the Lagrangian function involves third order derivatives

Summary Recently variational principles whose Lagrangian functions involve third order derivatives of the position vector have been considered with a view to applying them to certain aspects of elementary particle theory. It is known that definite consistency conditions arise when parameter invariance of the associated action integral is required. By invoking two additional assumptions, which have been adopted in the past, a general parameter invariant Lagrangian is deduced. The structure of the corresponding momentum vector is investigated. It is shown that a large class of the resulting equations of motion are derivable from a different approach due to Rund. Borelowski's equations of motion are also derived. Although no radiation effects are considered the parameter invariant counterpart of the Abraham vector plays an important role in the theory.     

Trajectory sensitivity analysis of hybrid systems with sliding motion

This paper concerns hybrid control systems exhibiting the sliding motion. It is assumed that the system’s motion on the switching surface is described by index-2 differential–algebraic equations (DAEs), which guarantee the accurate tracking of the sliding motion surface. For those systems the sensitivity analysis is performed with the help of solutions to system’s linearized equations. The paper states conditions under which the solutions to the linearized equations for original DAEs and the solutions to linearized equations for underlying ordinary differential equations (ODEs) exhibit similar properties. Due to the presence of sliding motion, we restrict the class of admissible control functions to piecewise differentiable functions. The presented sensitivity analysis might be useful in deriving the weak maximum principle for optimal control problems with hybrid systems exhibiting sliding motion and in establishing the global convergence of algorithms for solving those problems.     

Integration of the equations of gas dynamics for 2.5-dimensional solutions

We integrate the equations of gas dynamics in finite form for the solutions in which the thermodynamic parameters depend only on one spatial variable. The corresponding motion of gas represents the nonlinear superposition of the one-dimensional gas motion corresponding to the invariant system and the two-dimensional motion determined by noninvariant functions. These motions are called 2.5-dimensional. We reduce the invariant system to a first-order implicit ordinary differential equation. We study various solutions of the latter. We construct some continuous and discontinuous solutions to the equations of gas dynamics and give their physical interpretation.     

Integration of stiff mechanical systems by Runge-Kutta methods

The numerical integration of stiff mechanical systems is studied in which a strong potential forces the motion to remain close to a manifold. The equations of motion are written as a singular singular perturbation problem with a small stiffness parameter . Smooth solutions of such systems are characterized, in distinction to highly oscillatory general solutions. Implicit Runge-Kutta methods using step sizes larger than are shown to approximate smooth solutions, and precise error estimates are derived. As 0, Runge-Kutta solutions of the stiff system converge to Runge-Kutta solutions of the associated constrained system formulated as a differential-algebraic equation of index 3. Standard software for stiff initial-value problems does not work satisfactorily on the stiff systems considered here. The reasons for this failure are explained, and remedies are proposed.This work was supported in part by the Austrian Science Foundation, grant P8443-PHY.