Existence of quasibounded solutions for the higher order dynamic equations on measure chains

By using the exponential dichotomy and Schauder’s fixed point theorem, some new criteria are established for the existence of quasibounded solutions of the inhomogeneous system x^Δ = A(t)x + g(t, x) + h(t), which generalize the previous results in [15], [19].     

Entropy estimation of the Hénon attractor

The topological entropy of the Hénon attractor is estimated using a function that describes the stable and unstable manifolds of the Hénon map. This function provides an accurate estimate of the length of curves in the attractor. The estimation method presented here can be applied to cases in which the invariant set is not hyperbolic. From the result of the length calculation, we have estimated the topological entropy h as h ～ 0.49703 for the original parameters a = 1.4 and b = 0.3 adopted by Hénon.     

Inverse problems for parabolic equations 2

Let u_t − u_xx = h(t) in 0 ⩽ x ⩽ π, t ⩾ 0. Assume that u(0, t) = v(t), u(π, t) = 0, and u(x, 0) = g(t). The problem is: what extra data determine the three unknown functions {h, v, g} uniquely? This question is answered and an analytical method for recovery of the above three functions is proposed.     

Discrepancy principle for DSM II

Let Ay = f, A is a linear operator in a Hilbert space H, y ⊥ N(A) ≔ {u : Au = 0}, R(A) ≔ {h : h = Au, u ∈ D(A)} is not closed, ∥f_δ − f∥ ⩽ δ. Given f_δ, one wants to construct u_δ such that lim_δ→0∥u_δ − y∥ = 0. Two versions of discrepancy principles for the DSM (dynamical systems method) for finding the stopping time and calculating the stable solution u_δ to the original equation Ay = f are formulated and mathematically justified.     

On tropical matrices of small factor rank

The notion of the factor rank of tropical matrices is considered. We construct a linear-time algorithm that either finds a full-rank 3 × 3 submatrix of a given matrix A or concludes that the factor rank of A is less than 3. We show that there exist matrices of factor rank 4 whose 4 × 4 submatrices are all rank deficient.     

On the positive definite solutions of nonlinear matrix equation X + A⋆ X−δA = Q

In this paper we consider the positive definite solutions of nonlinear matrix equation X + A^⋆X^−δA = Q, where δ ∈ (0, 1], which appears for the first time in [S.M. El-Sayed, A.C.M. Ran, On an iteration methods for solving a class of nonlinear matrix equations, SIAM J. Matrix Anal. Appl. 23 (2001) 632–645]. The necessary and sufficient conditions for the existence of a solution are derived. An iterative algorithm for obtaining the positive definite solutions of the equation is discussed. The error estimations are found.     

Conditional contractivity of Runge–Kutta methods for nonlinear differential equations with many variable delays

This paper deals with conditional contractivity properties of Runge–Kutta (RK) methods with variable step-size applied to nonlinear differential equations with many variable delays (MDDEs). The concepts of CRN_m(ω, H)- and BN_f(μ, ?)-stability are introduced. It is shown that the numerical solution produced by a BN_f(μ, ?)-stable Runge–Kutta method with an appropriate interpolation is contractive. In particular, these results are also novel for nonlinear differential equations with many constant delays or single variable delay. To obtain BN_f(μ, ?)-stable methods, (k, l)-algebraically stable Runge–Kutta methods are also investigated.     

Integrals,Quantum Galois Extensions,and the Affineness Criterion for Quantum Yetter–Drinfel'd Modules

In this paper we shall generalize the notion of an integral on a Hopf algebra introduced by Sweedler, by defining the more general concept of an integral of a threetuple (H, A, C), where H is a Hopf algebra coacting on an algebra A and acting on a coalgebra C. We prove that there exists a total integral γ: C → Hom(C, A) of (H, A, C) if and only if any representation of (H, A, C) is injective in a functorial way, as a corepresentation of C. In particular, the quantum integrals associated to Yetter–Drinfel'd modules are defined. Let now A be an H-bicomodule algebra, ^H$\text{Y}$$\text{D}$_A the category of quantum Yetter–Drinfel'd modules, and B = {a ∈ A|∑S^− 1(a_〈1〉)a_{〈 − 1〉} ⊗ a_〈0〉 = 1_H ⊗ a}, the subalgebra of coinvariants of the Verma structure A ∈ ^H$\text{Y}$$\text{D}$_A. We shall prove the following affineness criterion: if there exists γ: H → Hom(H, A) a total quantum integral and the canonical map β: A ⊗ _B A → H ⊗ A, β(a ⊗ _B b) = ∑S^− 1(b_〈1〉)b_{〈 − 1〉} ⊗ ab_〈0〉 is surjective (i.e., A/B is a quantum homogeneous space), then the induction functor – ⊗ _B A: $\text{M}$_B → ^H$\text{Y}$$\text{D}$_A is an equivalence of categories. The affineness criteria proven by Cline, Parshall, and Scott, and independently by Oberst (for affine algebraic groups schemes) and Schneider (in the noncommutative case), are recovered as special cases.     

Exponential stability of delayed fuzzy cellular neural networks with diffusion

The exponential stability of delayed fuzzy cellular neural networks (FCNN) with diffusion is investigated. Exponential stability, significant for applications of neural networks, is obtained under conditions that are easily verified by a new approach. Earlier results on the exponential stability of FCNN with time-dependent delay, a special case of the model studied in this paper, are improved without using the time-varying term condition: dτ(t)/dt < μ.     

A semi-Markov process with an inverse Gaussian distribution as sojourn time

Let X_n denote the state of a device after n repairs. We assume that the time between two repairs is the time τ taken by a Wiener process {W(t), t ? 0}, starting from w₀ and with drift μ < 0, to reach c ∈ [0, w₀). After the nth repair, the process takes on either the value X_n?1 + 1 or X_n?1 + 2. The probability that X_n = X_n?1 + j, for j = 1, 2, depends on whether τ ? t₀ (a fixed constant) or τ > t₀. The device is considered to be worn out when X_n ? k, where k ∈ {1, 2, …}. This model is based on the ones proposed by Rishel (1991) [1] and Tseng and Peng (2007) [2]. We obtain an explicit expression for the mean lifetime of the device. Numerical methods are used to illustrate the analytical findings.     

On “P” property and the column-W property

P-matrices play an important role in the well-posedness of a linear complementarity problem (LCP). Similarly, the well-posedness of a horizontal linear complementarity problem (HLCP) is closely related to the column-W property of a matrix k-tuple.In this paper we first consider the problem of generating P-matrices from a given pair of matrices. Given a matrix pair (D, F) where D is a square matrix of order m and matrix F has m rows, “what are the conditions under which there exists a matrix G such that (D + FG) is a P-matrix?”. We obtain necessary and sufficient conditions for the special case when the column rank of F is m ? 1. A decision algorithm of complexity O(m²) to check whether the given pair of matrices (D, F) is P-matrisable is obtained. We also obtain a necessary and an independent sufficient condition for the general case when rank(F) is less than m ? 1.We then generalise the P-matrix generating problem to the generation of matrix k-tuples satisfying the column-W property from a given matrix (k + 1)-tuple. That is, given a matrix (k + 1)-tuple (D₁, … ,D_k, F), where D_js are square matrices of order m and F is a matrix having m rows, we determine the conditions under which the matrix k-tuple (D₁ + FG₁, … ,D_k + FG_k) satisfies the column-W property. As in the case of P-matrices we obtain necessary and sufficient conditions for the case when rank(F) = m ? 1. Using these conditions a decision algorithm of complexity O(km²) to check whether the given matrix (k + 1)-tuple is column-W matrisable is obtained. Then for the case when rank(F) is less than m ? 1, we obtain a necessary and an independent sufficient condition.For a special sub-class of P-matrices we give a polynomial time decision algorithm for P-matrisability. Finally, we obtain a geometric characterisation of column-W property by generalising the well known separation theorem for P-matrices.     

The stochastic soliton-like solutions of stochastic KdV equations

By means of a generalized method and symbolic computation, we consider a stochastic KdV equation U_t + f(t)U ♢ U_x + g(t)U_xxx = W(t) ♢ R^♢(t, U, U_x, U_xxx). We construct new and more general formal solutions. At the same time, we recover all the solutions found by Xie [Phys. Lett. A 310 (2003) 161]. The solutions obtained include the nontravelling wave and coefficient function’s stochastic soliton-like solutions, singular stochastic soliton-like solutions, stochastic triangular functions solutions.     

The maximal and minimal ranks of a quaternion matrix expression with applications

In this paper, we establish the formulas of the extermal ranks of the quaternion matrix expression f(X₁, X₂) = C₇ ? A₄X₁B₄ ? A₅X₂B₅ where X₁, X₂ are variant quaternion matrices subject to quaternion matrix equations A₁X₁ = C₁, A₂X₁ = C₂, A₃X₁ = C₃, X₂B₁ = C₄, X₂B₂ = C₅, X₂B₃ = C₆. As applications, we give a new necessary and sufficient condition for the existence of solutions to some systems of quaternion matrix equations. Some results can be viewed as special cases of the results of this paper.     

Magnetohydrodynamic stationary and oscillatory convective stability in a mushy layer during binary alloy solidification

For the case of solidification of a bottom cooled binary alloy, the magnetohydrodynamic stationary and oscillatory convective stability in the mushy layer is investigated analytically using normal mode linear stability analysis. In the limit of large Stefan number (S_t), a near–eutectic approximation with large far field temperature is considered in the present research. To ascertain the instability in the mushy layer, the strength of the superimposed magnetic field is so chosen that it corresponds to a given mush Hartmann number (Ha_m) of the problem. The results are presented for various values of mush Hartmann numbers in the range, 0 ≤ Ha_m ≤ 50. The critical Rayleigh number for stationary convection shows a linear relationship with increasing Ha_m. The magnetohydrodynamic effect imparts a stabilizing influence during stationary convection. In comparison to that of the stationary convective mode, the oscillatory mode appears to be critically susceptible at higher values of β (β = S_t/℘² ϒ², ℘ is the compositional ratio, ϒ = 1 + S_t/℘), and vice versa for lower β values. Analogous to the behavior for stationary convection, the magnetic field also offers a stabilizing effect in oscillatory convection and thus influences global stability of the mushy layer. Increasing magnetic strength shows reduction in the wavenumber and in the number of rolls formed in the mushy layer.     

Time-delayed control to suppress the nonlinear vibrations of a horizontally suspended Jeffcott-rotor system

Time-delay is an unavoidable phenomenon in active control systems. Measuring of the system states, processing of the measured signals, executing the control laws, conditioning and enforcing the control actions are the main reasons of time-delayed systems. This paper studies the vibration control of a horizontally suspended Jeffcott-rotor system having cubic and quadratic nonlinearities via time-delayed position-velocity controller. The intervals of the time-delays (τ₁ and τ₂) at which the system response is stable has been studied. The τ₁ − τ₂ plane is constructed to illustrate the area at which the system solutions are stable. The influences of the controller gains on the stable-solutions area in τ₁ − τ₂ plane are explored. The analysis revealed that the time-delay increases the vibration amplitudes and can destabilize the system solution in the case of negative position feedback control, while at positive position feedback control it improves the vibration suppression performance. The time-delays mechanism in stabilizing and destabilizing the dynamical systems is explained. Then, we proposed a simple and concrete method to determine the optimal value for time-delays that can improve the vibrations suppression efficiency. The acquired analytical results are confirmed numerically and the optimal working conditions of the system are concluded. Finally, a comparison with the papers that published previously is included.     

Unsteady magnetohydrodynamic Hartmann–Couette flow and heat transfer in a Darcian channel with Hall current,ionslip, viscous and Joule heating effects: Network numerical solutions

A theoretical study of unsteady magnetohydrodynamic viscous Hartmann–Couette laminar flow and heat transfer in a Darcian porous medium intercalated between parallel plates, under a constant pressure gradient is presented. Viscous dissipation, Joule heating, Hall current and ionslip current effects are included as is lateral mass flux at both plates. The dimensionless conservation equations for the primary (x¹-direction), secondary (z¹-direction) momentum and also energy conservation equation are derived and solved using a computational technique known as Network Simulation Methodology (NSM). Velocity distributions (u¹, w¹) and temperature distribution (T¹) at the channel centre (y¹ = 0) over time (t¹) are studied graphically for the effects of Darcy number (Da), Hartmann number (Ha), transpiration (N_t), Hall current parameter (Be), ionslip parameter (Bi), pressure gradient parameter (dP/dx¹) with Prandtl number prescribed at 7.0 (electrically conducting water), Eckert number held constant at 0.25 (heat convection from the plates to the fluid) and Reynolds number (Re) fixed at 5.0 (for Re < 10, Darcian model is generally valid). Increasing Darcy number causes an increase in temperature, T¹; values are however significantly reduced for the higher Hartmann number case (Ha = 10). For the case of low transpiration (i.e. N_t = 1 which corresponds to weak suction at the upper plate and weak injection at the lower plate), both primary velocity (u¹) and secondary velocity (w¹) are increased with a rise in Darcy number (owing to a simultaneous decrease in Darcian porous drag); temperature T¹ is also increased considerably with increasing Da. However, for stronger transpiration (N_t = 10), magnitudes of u¹, w¹ and T¹ are significantly reduced and also significant overshoots are detected prior to the establishment of steady state flow. With increasing Hall current parameter, Be, (for the purely fluid regime i.e. Da → ∞), primary velocity is considerably increased, whereas secondary velocity is reduced; temperatures are decreased in the early stages of flow but effectively increased in the steady state with increasing Be. With strong Darcian drag present (Da = 0.01 i.e. very low permeability), magnitudes of u¹, w¹ and T¹ are considerably reduced and temperatures are found to be reduced for all t¹, with increasing Hall current effect (Be). Increasing ionslip current parameter (Bi) increases primary velocity (u¹), decreases secondary velocity (w¹) and also temperature (T¹) for all time (t¹), in the infinite permeability case (Da → ∞). For weakly Darcian flow, ionslip parameter (Bi) has a much reduced effect on the velocity distributions. Temperature, T¹ is strongly increased with a rise in pressure gradient parameter, dP¹/dx¹, as is primary velocity (u¹); however, secondary velocity (w¹) is reduced. The present study has applications in hybrid magnetohydrodynamic (MHD) energy generators, materials processing, geophysical hydromagnetics, etc.     

Dynamics of landslide model with time delay and periodic parameter perturbations

In present paper, we analyze the dynamics of a single-block model on an inclined slope with Dieterich–Ruina friction law under the variation of two new introduced parameters: time delay T_d and initial shear stress μ. It is assumed that this phenomenological model qualitatively simulates the motion along the infinite creeping slope. The introduction of time delay is proposed to mimic the memory effect of the sliding surface and it is generally considered as a function of history of sliding. On the other hand, periodic perturbation of initial shear stress emulates external triggering effect of long-distant earthquakes or some non-natural vibration source. The effects of variation of a single observed parameter, T_d or μ, as well as their co-action, are estimated for three different sliding regimes: β < 1, β = 1 and β > 1, where β stands for the ratio of long-term to short-term stress changes. The results of standard local bifurcation analysis indicate the onset of complex dynamics for very low values of time delay. On the other side, numerical approach confirms an additional complexity that was not observed by local analysis, due to the possible effect of global bifurcations. The most complex dynamics is detected for β < 1, with a complete Ruelle–Takens–Newhouse route to chaos under the variation of T_d, or the co-action of both parameters T_d and μ. These results correspond well with the previous experimental observations on clay and siltstone with low clay fraction. In the same regime, the perturbation of only a single parameter, μ, renders the oscillatory motion of the block. Within the velocity-independent regime, β = 1, the inclusion and variation of T_d generates a transition to equilibrium state, whereas the small oscillations of μ induce oscillatory motion with decreasing amplitude. The co-action of both parameters, in the same regime, causes the decrease of block’s velocity. As for β > 1, highly-frequent, limit-amplitude oscillations of initial stress give rise to oscillatory motion. Also for β > 1, in case of perturbing only the initial shear stress, with smaller amplitude, velocity of the block changes exponentially fast. If the time delay is introduced, besides the stress perturbation, within the same regime, the co-action of T_d (T_d < 0.1) and small oscillations of μ induce the onset of deterministic chaos.     

Solution of the Burgers equation using an implicit linearizing transformation

The initial-boundary-value problem on the semi-infinite interval and on a finite interval for the Burgers equation u_t = u_xx + 2u_xu is solved using a stream function ? and a linearizing transformation w = e^?. The transformation reduces the equation to a heat equation with appropriate initial and homogeneous time-dependent linear boundary conditions. One advantage of this method is that we never need to find an explicit expression for ? in our computations.