Synchronization of piecewise continuous systems of fractional order

This paper proves analytically that synchronization of a class of piecewise continuous fractional-order systems can be achieved. Since there are no dedicated numerical methods to integrate differential equations with discontinuous right-hand sides for fractional-order models, Filippov’s regularization (Filippov, Differential Equations with Discontinuous Right-Hand Sides, 1988) is applied, and Cellina’s Theorem (Aubin and Cellina, Differential Inclusions Set-valued Maps and Viability Theory, 1984; Aubin and Frankowska, Set-valued Analysis, 1990) is used. It is proved that the corresponding initial value problem can be converted to a continuous problem of fractional-order systems, to which numerical methods can be applied. In this way, the synchronization problem is transformed into a standard problem for continuous fractional-order systems. Three examples are presented: the Sprott’s system, Chen’s system, and Shimizu–Morioka’s system.     

Stability and Freezing of Nonlinear Waves in First Order Hyperbolic PDEs

It is a well-known problem to derive nonlinear stability of a traveling wave from the spectral stability of a linearization. In this paper we prove such a result for a large class of hyperbolic systems. To cope with the unknown asymptotic phase, the problem is reformulated as a partial differential algebraic equation for which asymptotic stability becomes usual Lyapunov stability. The stability proof is then based on linear estimates from (Rottmann-Matthes, J Dyn Diff Equat 23:365–393, 2011) and a careful analysis of the nonlinear terms. Moreover, we show that the freezing method (Beyn and Thümmler, SIAM J Appl Dyn Syst 3:85–116, 2004; Rowley et al. Nonlinearity 16:1257–1275, 2003) is well-suited for the long time simulation and numerical approximation of the asymptotic behavior. The theory is illustrated by numerical examples, including a hyperbolic version of the Hodgkin–Huxley equations.     

New results on stability and stabilization of a class of nonlinear fractional-order systems

The asymptotic stability and stabilization problem of a class of fractional-order nonlinear systems with Caputo derivative are discussed in this paper. By using of Mittag–Leffler function, Laplace transform, and the generalized Gronwall inequality, a new sufficient condition ensuring local asymptotic stability and stabilization of a class of fractional-order nonlinear systems with fractional-order α:1<α<2 is proposed. Then a sufficient condition for the global asymptotic stability and stabilization of such system is presented firstly. Finally, two numerical examples are provided to show the validity and feasibility of the proposed method.     

OGY method for a class of discontinuous dynamical systems

In this paper, we prove that the OGY method to control unstable periodic orbits (UPOs) of continuous-time systems can be applied to a class of systems discontinuous with respect the state variable, by using a generalized derivative. Because the discontinuous problem may have not classical solutions, the initial value problem is transformed into a set-valued problem via Filippov regularization. The existence of the ingredients necessary to apply OGY method (UPO, Poincaré map and stable and unstable directions) is proved and the numerically implementation is explained. Another possible way analyzed in this paper is the continuous approximation of the underlying initial value problem, via Cellina??s theorem for differential inclusions. Thus, the problem is approximated by a continuous initial value problem, and the OGY method can be applied as usual.     

Stabilization using fractional-order PI and PID controllers

This paper presents a solution to the problem of stabilizing a given fractional dynamic system using fractional-order PIλ and PI^λD^μ controllers. It is based on plotting the global stability region in the (k _p, k _i)-plane for the PI^λ controller and in the (k _p , k _i , k _d)-space for the PIλDμ controller. Analytical expressions are derived for the purpose of describing the stability domain boundaries which are described by real root boundary, infinite root boundary and complex root boundary. Thus, the complete set of stabilizing parameters of the fractional-order controller is obtained. The algorithm has a simple and reliable result which is illustrated by several examples, and hence is practically useful in the analysis and design of fractional-order control systems.     

Stability and stabilization of a class of fractional-order nonlinear systems for $$\varvec{0<}\,{\varvec{\alpha }} \,\varvec{< 2}$$

This paper investigates the stability and stabilization problem of fractional-order nonlinear systems for \(0<\alpha <2\). Based on the fractional-order Lyapunov stability theorem, S-procedure and Mittag–Leffler function, the stability conditions that ensure local stability and stabilization of a class of fractional-order nonlinear systems under the Caputo derivative with \(0<\alpha <2\) are proposed. Finally, typical instances, including the fractional-order nonlinear Chen system and the fractional-order nonlinear Lorenz system, are implemented to demonstrate the feasibility and validity of the proposed method.     

A practical synchronization approach for fractional-order chaotic systems

A practical synchronization approach is proposed for a class of fractional-order chaotic systems to realize perfect \(\delta \)-synchronization, and the nonlinear functions in the fractional-order chaotic systems are all polynomials. The \(\delta \)-synchronization scheme in this paper means that the origin in synchronization error system is stable. The reliability of \(\delta \)-synchronization has been confirmed on a class of fractional-order chaotic systems with detailed theoretical proof and discussion. Furthermore, the \(\delta \)-synchronization scheme for the fractional-order Lorenz chaotic system and the fractional-order Chua circuit is presented to demonstrate the effectiveness of the proposed method.     

Decay Structure for Symmetric Hyperbolic Systems with Non-Symmetric Relaxation and its Application

This paper is concerned with the decay structure for linear symmetric hyperbolic systems with relaxation. When the relaxation matrix is symmetric, the dissipative structure of the systems is completely characterized by the Kawashima–Shizuta stability condition formulated in Umeda et al. (Jpn J Appl Math 1:435–457, 1984) and Shizuta and Kawashima (Hokkaido Math J 14:249–275, 1985) and we obtain the asymptotic stability result together with the explicit time-decay rate under that stability condition. However, some physical models which satisfy the stability condition have non-symmetric relaxation term (for example, the Timoshenko system and the Euler–Maxwell system). Moreover, it had been already known that the dissipative structure of such systems is weaker than the standard type and is of the regularity-loss type (see Duan in J Hyperbolic Differ Equ 8:375–413, 2011; Ide et al. in Math Models Meth Appl Sci 18:647–667, 2008; Ide and Kawashima in Math Models Meth Appl Sci 18:1001–1025, 2008; Ueda et al. in SIAM J Math Anal 2012; Ueda and Kawashima in Methods Appl Anal 2012). Therefore our purpose in this paper is to formulate a new structural condition which includes the Kawashima–Shizuta condition, and to analyze the weak dissipative structure for general systems with non-symmetric relaxation.     

Application of the method of direct separation of motions to the parametric stabilization of an elastic wire

The paper considers the application of the method of direct separation of motions to the investigation of distributed systems. An approach is proposed which allows one to apply the method directly to the initial equation of motion and to satisfy all boundary conditions, arising for both slow and fast components of motion. The methodology is demonstrated by means of a classical problem concerning the so-called Indian magic rope trick (Blekhman et al. in Selected topics in vibrational mechanics, vol. 11, pp. 139–149, [2004]; Champneys and Fraser in Proc. R. Soc. Lond. A 456:553–570, [2000]; in SIAM J. Appl. Math. 65(1):267–298, [2004]; Fraser and Champneys in Proc. R. Soc. Lond. A 458:1353–1373, [2002]; Galan et al. in J. Sound Vib. 280:359–377, [2005]), in which a wire with an unstable upper vertical position is stabilized due to vertical vibration of its bottom support point. The wire is modeled as a heavy Bernoulli–Euler beam with a vertically vibrating lower end. As a result of the treatment, an explicit formula is obtained for the vibrational correction to the critical flexural stiffness of the nonexcited system.     

Stability analysis of Caputo fractional-order nonlinear systems revisited

In this paper stability analysis of fractional-order nonlinear systems is studied. An extension of Lyapunov direct method for fractional-order systems using Bihari’s and Bellman–Gronwall’s inequality and a proof of comparison theorem for fractional-order systems are proposed.     

A numerical simulation of external heat transfer around turbine blades

External heat transfer prediction is performed in two-dimensional turbine blade cascades using the Reynolds-averaged Navier–Stokes equations. For this purpose, six different turbulence models including the algebraic Baldwin–Lomax (AIAA paper 78-257, 1978), three low-Re k−ɛ models (Chien in AIAA J 20:33–38, 1982; Launder and Sharma in Lett Heat Mass Transf 1(2):131–138, 1974; Biswas and Fukuyama in J Turbomach 116:765–773, 1994), and two k−ω models (Wilcox in AIAA J 32(2):247–255, 1994) are taken into account. The computer code developed employs a finite volume method to solve governing equations based on an explicit time marching approach with capability to simulate subsonic, transonic and supersonic flows. The Roe method is used to decompose the inviscid fluxes and the gradient theorem to decompose viscous fluxes. The performance of different turbulence models in prediction of heat transfer is examined. To do so, the effect of Reynolds and Mach numbers along with the turbulent intensity are taken into account, and the numerical results obtained are compared with the experimental data available.     

Bursting generation mechanism in a three-dimensional autonomous system,chaos control,and synchronization in its fractional-order form

The bifurcation mechanism of bursting oscillations in a three-dimensional autonomous slow-fast Kingni et al. system (Nonlinear Dyn. 73, 1111–1123, 2013) and its fractional-order form are investigated in this paper. The stability analysis of the system is carried out assuming that the slow subsystem evolves on quasi-static state. It is reveaved that the bursting oscillations found in the system result from the system switching between the unstable and the stable states of the only equilibrium point of the fast subsystem. We refer this class of bursting to “source/bursting.” The coexistence of symmetrical bursting limit cycles and chaotic bursting attractors is observed. In addition, the fractional-order chaotic slow-fast system is studied. The lowest order of the commensurate form of this system to exhibit chaotic behavior is found to be 2.199. By tuning the commensurate fractional-order, the chaotic slow-fast system displays Chen- and Lorenz-like chaotic attractors, respectively. The stability analysis of the controlled fractional-order-form of the system to its equilibria is undertaken using Routh–Hurwitz conditions for fractional-order systems. Moreover, the synchronization of chaotic bursting oscillations in two identical fractional-order systems is numerically studied using the unidirectional linear error feedback coupling scheme. It is shown that the system can achieve synchronization for appropriate coupling strength. Furthermore, the effect of fractional derivatives orders on chaos control and synchronization is analyzed.     

Explicit Rate–Time Solutions for Modeling Matrix–Fracture Flow of Single Phase Gas in Dual-Porosity Media

One of the widely used methods for modeling matrix–fracture fluid exchange in naturally fractured reservoirs is dual porosity approach. In this type of modeling, matrix blocks are regarded as sources/sinks in the fracture network medium. The rate of fluid transfer from matrix blocks into fracture medium may be modeled using shape factor concept (Warren and Root, SPEJ 3:245–255, 1963); or the rate–time solution is directly derived for the specific matrix geometry (de Swaan, SPEJ 16:117–122, 1976). Numerous works have been conducted to study matrix–fracture fluid exchange for slightly compressible fluids (e.g. oil). However, little attention has been taken to systems containing gas (compressible fluid). The objective of this work is to develop explicit rate–time solutions for matrix–fracture fluid transfer in systems containing single phase gas. For this purpose, the governing equation describing flow of gas from matrix block into fracture system is linearized using pseudopressure and pseudotime functions. Then, the governing equation is solved under specific boundary conditions to obtain an implicit relation between rate and time. Since rate calculations using such an implicit relation need iterations, which may be computationally inconvenient, an explicit rate–time relation is developed with the aid of material balance equation and several specific assumptions. Also, expressions are derived for average pseudopressure in matrix block. Furthermore, simplified solutions (originated from the complex general solutions) are introduced applicable in infinite and finite acting flow periods in matrix. Based on the derived solutions, expressions are developed for shape factor. An important observation is that the shape factor for gas systems is the same as that of oil bearing matrix blocks. Subsequently, a multiplier is introduced which relates rate to matrix pressure instead of matrix pseudopressure. Finally, the introduced equations are verified using a numerical simulator.     

On the motion of a gyrostat similar to Lagrange’s gyroscope under the influence of a gyrostatic moment vector

In this paper, the problem of the motion of a gyrostat fixed at one point under the action of a gyrostatic moment vector whose components are ℓ _i (i=1,2,3) about the axes of rotation, similar to a Lagrange gyroscope is investigated. We assume that the center of mass G of this gyrostat is displaced by a small quantity relative to the axis of symmetry, and that quantity is used to obtain the small parameter ε (Elfimov in PMM, 42(2):251–258, [1978]). The equations of motion will be studied under certain initial conditions of motion. The Poincaré small parameter method (Malkin in USAEC, Technical Information Service, ABC. Tr-3766, [1959]; Nayfeh in Perturbation methods, Wiley-Interscience, New York, [1973]) is applied to obtain the periodic solutions of motion. The periodic solutions for the case of irrational frequencies ratio are given. The periodic solutions are analyzed geometrically using Euler’s angles to describe the orientation of the body at any instant t of time. These solutions are performed by our computer programs to get their graphical representations.     

Development and characterization of a variable turbulence generation system

Experimental turbulent combustion studies require systems that can simulate the turbulence intensities [u′/U ₀ ~ 20–30% (Koutmos and McGuirk in Exp Fluids 7(5):344–354, 1989)] and operating conditions of real systems. Furthermore, it is important to have systems where turbulence intensity can be varied independently of mean flow velocity, as quantities such as turbulent flame speed and turbulent flame brush thickness exhibit complex and not yet fully understood dependencies upon both U ₀ and u′. Finally, high pressure operation in a highly pre-heated environment requires systems that can be sealed, withstand high gas temperatures, and have remotely variable turbulence intensity that does not require system shut down and disassembly. This paper describes the development and characterization of a variable turbulence generation system for turbulent combustion studies. The system is capable of a wide range of turbulence intensities (10–30%) and turbulent Reynolds numbers (140–2,200) over a range of flow velocities. An important aspect of this system is the ability to vary the turbulence intensity remotely, without changing the mean flow velocity. This system is similar to the turbulence generators described by Videto and Santavicca (Combust Sci Technol 76(1):159–164, 1991) and Coppola and Gomez (Exp Therm Fluid Sci 33(7):1037–1048, 2009), where variable blockage ratio slots are located upstream of a contoured nozzle. Vortical structures from the slots impinge on the walls of the contoured nozzle to produce fine-scale turbulence. The flow field was characterized for two nozzle diameters using three-component Laser Doppler velocimetry (LDV) and hotwire anemometry for mean flow velocities from 4 to 50 m/s. This paper describes the key design features of the system, as well as the variation of mean and RMS velocity, integral length scales, and spectra with nozzle diameter, flow velocity, and turbulence generator blockage ratio.     

Application of Rheological Models in Prediction of Turbulent Slurry Flow

The paper deals with fully developed steady turbulent flow of slurry in a circular straight and smooth pipe. The Kaolin slurry consists of very fine solid particles, so the solid particles concentration, and density, and viscosity are assumed to be constant across the pipe. The mathematical model is based on the time averaged momentum equation. The problem of closure was solved by the Launder and Sharma k-ε turbulence model (Launder and Sharma, Lett Heat Mass Transf 1:131–138, 1974) but with a different turbulence damping function. The turbulence damping function, used in the mathematical model in the present paper, is that proposed by Bartosik (1997). The mathematical model uses the apparent viscosity concept and the apparent viscosity was calculated using two- and three-parameter rheological models, namely Bingham and Herschel–Bulkley. The main aim of the paper is to compare measurements and predictions of the frictional head loss and velocity distribution, taking into account two- and three-parameter rheological models, namely Bingham and Herschel–Bulkley, if the Kaolin slurry possesses low, moderate, and high yield stress. Predictions compared with measurements show an observable advantage of the Herschel–Bulkley rheological model over the Bingham model particularly if the bulk velocity decreases.     

A Fractional Model of Continuum Mechanics

Although there has been renewed interest in the use of fractional models in many application areas, in reality fractional analysis has a long and distinguished history and can be traced back to the likes of Leibniz (Letter to L’Hospital, 1695), Liouville (J. éc. Polytech. 13:71, 1832), and Riemann (Gesammelte Werke, p. 62, 1876). Recent publications (Podlubny in Math. Sci. Eng. 198, 1999; Sabatier et al. in Advances in fractional calculus: theoretical developments and applications in physics and engineering, Springer, Berlin, 2007; Das in Functional fractional calculus for system identification and controls, Springer, Berlin, 2007) demonstrate that fractional derivative models have found widespread applications in science and engineering. Late fundamental considerations have led to the introduction of fractional calculus in continuum mechanics in an attempt to develop non-local constitutive relations (Lazopoulos in Mech. Res. Commun. 33:753–757, 2006). Attempts have also been made to model microscopic forces using fractional derivatives (Vazquez in Nonlinear waves: classical and quantum aspects, pp. 129–133, 2004). Our approach in this paper differs from previous theoretical work, in that we develop a general framework directly from the classical continuum mechanics, by defining the laws of motion and the stresses using fractional derivatives. The timeliness and relevance of this work is justified by the surge in interest in applications of fractional order models to biological, physical and economic systems. The aim of the present paper is to lay the foundations for a new non-local model of continuum mechanics based on fractional order derivatives which we will refer to as the fractional model of continuum mechanics. Following the theoretical development, we apply this framework to two one-dimensional model problems: the deformation of an infinite bar subjected to a self-equilibrated load distribution, and the propagation of longitudinal waves in a thin finite bar.     

Synchronization between integer-order chaotic systems and a class of fractional-order chaotic system based on fuzzy sliding mode control

In this paper, we focus on the synchronization between integer-order chaotic systems and a class of fractional-order chaotic system using the stability theory of fractional-order systems. A new fuzzy sliding mode method is proposed to accomplish this end for different initial conditions and number of dimensions. Furthermore, three examples are presented to illustrate the effectiveness of the proposed scheme, which are the synchronization between a fractional-order Lü chaotic system and an integer-order Liu chaotic system, the synchronization between a fractional-order hyperchaotic system based on Chen??s system and an integer-order hyperchaotic system based upon the Lorenz system, and the synchronization between a fractional-order hyperchaotic system based on Chen??s system, and an integer-order Liu chaotic system. Finally, numerical results are presented and are in agreement with theoretical analysis.     

On the Amplitude Equations for Weakly Nonlinear Surface Waves

Nonlocal generalizations of Burgers’ equation were derived in earlier work by Hunter (Contemp Math, vol 100, pp 185–202. AMS, 1989), and more recently by Benzoni-Gavage and Rosini (Comput Math Appl 57(3–4):1463–1484, 2009), as weakly nonlinear amplitude equations for hyperbolic boundary value problems admitting linear surface waves. The local-in-time well-posedness of such equations in Sobolev spaces was proved by Benzoni-Gavage (Differ Integr Equ 22(3–4):303–320, 2009) under an appropriate stability condition originally pointed out by Hunter. The same stability condition has also been shown to be necessary for well-posedness in Sobolev spaces in a previous work of the authors in collaboration with Tzvetkov (Benzoni-Gavage et al. in Adv Math 227(6):2220–2240, 2011). In this article, we show how the verification of Hunter’s stability condition follows from natural stability assumptions on the original hyperbolic boundary value problem, thus avoiding lengthy computations in each particular situation. We also show that the resulting amplitude equation has a Hamiltonian structure when the original boundary value problem has a variational origin. Our analysis encompasses previous equations derived for nonlinear Rayleigh waves in elasticity.     

Degenerate Regularization of Forward–Backward Parabolic Equations: The Regularized Problem

We study a quasilinear parabolic equation of forward–backward type in one space dimension, under assumptions on the nonlinearity which hold for a number of important mathematical models (for example, the one-dimensional Perona–Malik equation), using a degenerate pseudoparabolic regularization proposed in Barenblatt et al. (SIAM J Math Anal 24:1414–1439, 1993), which takes time delay effects into account. We prove existence and uniqueness of positive solutions of the regularized problem in a space of Radon measures. We also study qualitative properties of such solutions, in particular concerning their decomposition into an absolutely continuous part and a singular part with respect to the Lebesgue measure. In this respect, the existence of a family of viscous entropy inequalities plays an important role.