Derivatives of (modified) Fredholm determinants and stability of standing and traveling waves

Continuing a line of investigation initiated in [F. Gesztesy, Y. Latushkin, K.A. Makarov, Evans functions, Jost functions, and Fredholm determinants, Arch. Rat. Mech. Anal. 186 (2007) 361–421] exploring the connections between Jost and Evans functions and (modified) Fredholm determinants of Birman–Schwinger type integral operators, we here examine the stability index, or sign of the first nonvanishing derivative at frequency zero of the characteristic determinant, an object that has found considerable use in the study by Evans function techniques of stability of standing and traveling wave solutions of partial differential equations (PDE) in one dimension. This leads us to the derivation of general perturbation expansions for analytically-varying modified Fredholm determinants of abstract operators. Our main conclusion, similarly in the analysis of the determinant itself, is that the derivative of the characteristic Fredholm determinant may be efficiently computed from first principles for integral operators with semi-separable integral kernels, which include in particular the general one-dimensional case, and for sums thereof, which appears to offer applications in the multi-dimensional case.A second main result is to show that the multi-dimensional characteristic Fredholm determinant is the renormalized limit of a sequence of Evans functions defined in [G.J. Lord, D. Peterhof, B. Sandstede, A. Scheel, Numerical computation of solitary waves in infinite cylindrical domains, SIAM J. Numer. Anal. 37 (2000) 1420–1454] on successive Galerkin subspaces, giving a natural extension of the one-dimensional results of [F. Gesztesy, Y. Latushkin, K.A. Makarov, Evans functions, Jost functions, and Fredholm determinants, Arch. Rat. Mech. Anal. 186 (2007) 361–421] and answering a question of [J. Niesen, Evans function calculations for a two-dimensional system, presented talk, SIAM Conference on Applications of Dynamical Systems, Snowbird, UT, USA, May 2007] whether this sequence might possibly converge (in general, no, but with renormalization, yes). Convergence is useful in practice for numerical error control and acceleration.     

Mathematical modeling of time fractional reaction–diffusion systems

We study a fractional reaction–diffusion system with two types of variables: activator and inhibitor. The interactions between components are modeled by cubical nonlinearity. Linearization of the system around the homogeneous state provides information about the stability of the solutions which is quite different from linear stability analysis of the regular system with integer derivatives. It is shown that by combining the fractional derivatives index with the ratio of characteristic times, it is possible to find the marginal value of the index where the oscillatory instability arises. The increase of the value of fractional derivative index leads to the time periodic solutions. The domains of existing periodic solutions for different parameters of the problem are obtained. A computer simulation of the corresponding nonlinear fractional ordinary differential equations is presented. For the fractional reaction–diffusion systems it is established that there exists a set of stable spatio-temporal structures of the one-dimensional system under the Neumann and periodic boundary conditions. The characteristic features of these solutions consist of the transformation of the steady-state dissipative structures to homogeneous oscillations or space temporary structures at a certain value of fractional index and the ratio of characteristic times of system.     

Case of a generating family of quasi-periodic solutions in the theory of small parameter : PMM vol. 37, n6, 1973, pp. 990–998

We consider the problem of the existence and the stability in-the-small of periodic solutions of systems of ordinary differential equations with a small parameter μ, which in the generating approximation (μ = 0) admit of a family of quasi-periodic solutions (we are concerned only with the solutions belonging to the indicated family when μ = 0). The case to be investigated is in a specific sense a more general case of the unisolated generating solution in the small parameter theory and, therefore, includes everything previously treated by Malkin [1], Blekhman [2], and others. The main difficulty in the investigation is the presence of a multiple zero root in the characteristic determinant of the problem's generating system, to which both simple as well as quadratic elementary divisors [3] correspond. This fact predestines the presence of three groups of stability criteria for the solution being examined. The method for constructing these criteria, proposed here, assumes, in contrast to a previous one [1], the preliminary determination of not only the generating approximation but also the first one to the desired periodic solution. Particular aspects of the general “mixed” problem treated here were studied earlier in [4, 5].     

On existence and nonexistence of global solutions of Cauchy–Goursat problem for nonlinear wave equations

We consider the Cauchy–Goursat initial characteristic problem for nonlinear wave equations with power nonlinearity. Depending on the power of nonlinearity and the parameter in an equation we investigate the problem on existence and nonexistence of global solutions of the Cauchy–Goursat problem. The question on local solvability of the problem is also considered.     

On the theory of layered anisotropic cylindrical shells stiffened with ribs

The deformation, stability and vibration equations for anisotropic cylindrical shells stiffened with individual longitudinal and circumferential ribs are derived without introducing the hypothesis of nondeformable normals. The more general assumption adopted for layered materials (for example, glass-reinforced plastics) of a linear variation of the displacements over the thickness of the shell and the height of the ribs is used; in this case for the points of contact of the shell and the ribs after deformation the common normals form broken lines. The solution of the problem of the stability of a cylindrical shell stiffened with circumferential ribs is examined. For a shell with different, arbitrarily located ribs the problem is reduced to a homogeneous algebraic system of equations equal in number to three times the number of ribs.Moscow. Translated from Mekhanika Polimerov, No. 4, pp. 647–654, July–August, 1974.     

Stability of incompressible current-vortex sheets

We revisit the study in [Y. Trakhinin, On the existence of incompressible current-vortex sheets: study of a linearized free boundary value problem, Math. Methods Appl. Sci. 28 (2005) 917–945] where an energy a priori estimate for the linearized free boundary value problem for planar current-vortex sheets in ideal incompressible magnetohydrodynamics was proved for a part of the whole stability domain found a long time ago in [S.I. Syrovatskij, The stability of tangential discontinuities in a magnetohydrodynamic medium, Zh. Eksper. Teor. Fiz. 24 (1953) 622–629 (in Russian); W.I. Axford, Note on a problem of magnetohydrodynamic stability, Canad. J. Phys. 40 (1962) 654–655]. In this paper we derive an a priori estimate in the whole stability domain. The crucial point in deriving this estimate is the construction of a symbolic symmetrizer for a nonstandard elliptic problem for the small perturbation of total pressure. This symmetrizer is an analogue of Kreiss' type symmetrizers. As in hyperbolic theory, the failure of the uniform Lopatinski condition, i.e., the fact that current-vortex sheets are only weakly (neutrally) stable yields loss of derivatives in the energy estimate. The result of this paper is a necessary step to prove the local-in-time existence of stable nonplanar incompressible current-vortex sheets by a suitable Nash–Moser type iteration scheme.     

Thermodynamic stability of crystalline polymers in the transition region

The type of phase transition can be determined from the temperature dependence of the coefficient (–P/av)_s and the stability determinant D. We have calculated these quantities for polyethylene from data on the speed of ultrasound, density, and the specific heat C_p.Krupskaya Moscow Regional Pedagogical Institute. Translated from Mekhanika Polimerov, No. 4, pp. 724–726, July–August, 1971.     

Asymptotic behavior of the Toeplitz determinant

One considers the asymptotic behavior of the Toeplitz determinant for a nonnegative function. One proves that under certain conditions on the function the determinant admits the asymptotic representation where Translated from Zapiski Nauchnykh Seminarov Leningradskogo Otdeleniya Matematicheskogo Instituta im. V. A. Steklova AN SSSR, Vol. 97, pp. 22–31, 1980.In conclusion I express my gratitude to V. N. Solev for the formulation of the problem and for his help in the preparation of the paper.     

Experimental-theoretical studies of the stability of orthotropic shells with a filler in axial compression

Conclusions 1. An analysis has been made of the solution to the problem of the stability of multilayer cylindrical shells having a filler and simple calculation formulas have been obtained for determining the critical forces.2. The stability of fiberglass-plastic shells with rubber-like fillers has been studied experimentally.3. Comparative experimental-theoretical studies of critical forces have been made, and the stability coefficients have been ascertained for the shell class under consideration.Translated from Mekhanika Polimerov, No. 3, pp. 485–489, May–June, 1978.     

Perturbation analysis of the orthogonal procrustes problem

Given two arbitrary real matricesA andB of the same size, the orthogonal Procrustes problem is to find an orthogonal matrixM such that the Frobenius norm MA – B is minimized. This paper treats the common case when the orthogonal matrixM is required to have a positive determinant. The stability of the problem is studied and supremum results for the perturbation bounds are derived.     

Polynomial Fibonacci–Hessenberg matrices

A Fibonacci–Hessenberg matrix with Fibonacci polynomial determinant is referred to as a polynomial Fibonacci–Hessenberg matrix. Several classes of polynomial Fibonacci–Hessenberg matrices are introduced. The notion of two-dimensional Fibonacci polynomial array is introduced and three classes of polynomial Fibonacci–Hessenberg matrices satisfying this property are given.     

Computing Stability of Differential Equations with Bounded Distributed Delays

This paper deals with the stability analysis of scalar delay integro-differential equations (DIDEs). We propose a numerical scheme for computing the stability determining characteristic roots of DIDEs which involves a linear multistep method as time integration scheme and a quadrature method based on Lagrange interpolation and a Gauss–Legendre quadrature rule. We investigate to which extent the proposed scheme preserves the stability properties of the original equation. We derive and prove a sufficient condition for (asymptotic) stability of a DIDE (with a constant kernel) which we call RHP-stability. Conditions are obtained under which the proposed scheme preserves RHP-stability. We compare the obtained results with corresponding ones using Newton–Cotes formulas. Results of numerical experiments on computing the stability of DIDEs with constant and nonconstant kernel functions are presented.     

Stability of a sandwich plate composed of rubbery material

The stability of a sandwich plate is investigated within the framework of the plane problem of the stability of nonlinear elastic incompressible bodies for finite subcritical strains and an arbitrary form of the potential. Numerical examples are considered for rubbery bodies with Treloar (neo-Hookean body) and Mooney potentials.Institute of Mechanics, Academy of Sciences of the Ukrainian SSR, Kiev. Kiev Institute of the National Economy. Translated from Mekhanika Polimerov, No. 2, pp. 334–337, March–April, 1971.     

The number of roots of a polynomial outside a circle

We obtain a new criterion in terms of determinant inequalities that all the roots of a real polynomial should lie inside the unit circle, i.e., a criterion for the stability of periodic motions. In comparison with the Shur-Kon criterion, the number of determinants is halved.The results of this paper were published without proof in [2].Translated from Matematicheskie Zametki, Vol. 13, No. 1, pp. 3–12, January, 1973.     

Surface instability of high-elastic materials

The problem of the surface instability of a lower half-plane in compression is investigated for an arbitrary form of the potential within the framework of the plane problem of the deformational stability of incompressible elastic bodies at finite subcritical strains [1, 4]. Numerical examples are considered for high-elastic bodies with potentials of the Mooney and Treloar types.Institute of Mechanics, Academy of Sciences of the Ukrainian SSR, Kiev. Translated from Mekhanika Polimerov, No. 6, pp. 1107–1110, November–December, 1970.     

Fuzzy control of thyristor-controlled series compensator in power system transients

Power system transient stability is one of the most challenging technical areas in electric power industry. Thyristor-controlled series compensation (TCSC) is expected to improve transient stability and damp power oscillations. TCSC control in power system transients is a nonlinear control problem. This paper presents a T–S-model-based fuzzy control scheme and a systematic design method for the TCSC fuzzy controller. The nonlinear power system containing TCSC is modelled as a fuzzy “blending” of a set of locally linearized models. A linear optimal control is designed for each local linear model. Different control requirements at different stages during power system transients can be considered in deriving the linear control rules. The resulting fuzzy controller is then a fuzzy “blending” of these linear controllers. Quadratic stability of the overall nonlinear controlled system can be checked and ensured using H^∞ control theory. Digital simulation with NETOMAC software has verified that the fuzzy control scheme can improve power system transient stability and damp power swings very quickly.     

Distribution of integral points on the determinant surface

One obtains an asymptotic formula with remainder term for the number of second-order integral matrices with an increasing determinant, belonging to a given region of the discriminant surface and to a given residue class. The results are more accurate than in A. M. Istamov's paper (this issue, pp. 14–17) and are obtained in a somewhat different manner. The presentation is more detailed.Translated from Zapiski Nauchnykh Seminarov Leningradskogo Otdeleniya Matematicheskogo Instituta im. V. A. Steklova AN SSSR, Vol. 93, pp. 30–40, 1980.     

Dynamic stability of glass-reinforced plastic shells

The problem of the dynamic stability of circular-cylindrical glass-reinforced plastic shells subjected to external transverse pressure is examined in the nonlinear formulation. After the Lagrange equations have been constructed, the problem reduces to the integration of a system of ordinary differential equations with aperiodic coefficients. The integration has been carried out numerically on a computer for various loading rates and shell parameters. Analogous problems for isotropic metal shells were examined in [1–4]. A review of the subject may be found in [5].Mekhanika Polimerov, Vol. 4, No. 1, pp. 109–115, 1968     

Determinant of the Schrödinger operator

For the example of the nonrelativistic Schrödinger operator, methods are formulated for calculating the determinant of an elliptic operator on the basis of scattering theory. It is shown that such a determinant is identical to the Jost determinant at zero energy. In the centrally symmetric case, it reduces to ordinary Jost functions and ultimately to the values of the zero-energy wave functions at the origin. The relationship between the determinant of the Schrödinger operator and the characteristics of the scattering resonances and the number of bound states in a field of opposite sign is noted. This makes it possible to find the first terms in the gradient expansion of the determinant as a functional of the potential. The problem of the correlation free energy of a classical plasma serves as a physical illustration.P. N. Lebedev Physics Institute, Russian Academy of Sciences. Translated from Teoreticheskaya i Matematicheskaya Fizika, Vol. 96, No. 3, pp. 373–384, September, 1993.     

Mathematical modeling and control of population systems: Applications in biological pest control

The aim of this paper is to apply methods from optimal control theory, and from the theory of dynamic systems to the mathematical modeling of biological pest control. The linear feedback control problem for nonlinear systems has been formulated in order to obtain the optimal pest control strategy only through the introduction of natural enemies. Asymptotic stability of the closed-loop nonlinear Kolmogorov system is guaranteed by means of a Lyapunov function which can clearly be seen to be the solution of the Hamilton–Jacobi–Bellman equation, thus guaranteeing both stability and optimality. Numerical simulations for three possible scenarios of biological pest control based on the Lotka–Volterra models are provided to show the effectiveness of this method.