Stability Analysis for Digital Redesign of Discrete-Time Switched Systems Using H∞ Linear Matrix Inequality

In this paper, stability is studied for a class of switched singular systems containing both stable and unstable modes. By introducing a time-varying piecewise Lyapunov function (TVPLF) and a mode-dependent average dwell time (ADT) switching rule, the computable sufficient conditions for system stability are derived. The time-varying piecewise Lyapunov functions are piecewise continuously differentiable on every mode (but may not be differentiable at the interpolating points of the dwell time). This Lyapunov function method is particularly advantageous in overcoming the limitations of traditional multiple Lyapunov function (MLF) methods, which may not have a feasible solution when dealing with switched systems containing only unstable modes. As such, the TVPLF offers greater flexibility in application. Compared with the conventional ADT switching rule, the mode-dependent ADT switching rule not only enables each mode to have its own ADT but also allows for its own switching strategy. Specifically, the stable mode adopts a slow switching strategy while the unstable mode adopts a fast one, thereby reducing the conservatism of the ADT switching rule. Furthermore, based on the stability analysis, the time-varying controllers are proposed to stabilize the switched singular system, which can be expressed as the sequential linear combination of a series of linear state feedback on each mode. The proposed controllers are continuous for each mode, which are different from the controllers designed through the traditional MLF and MDLF methods, where the controllers designed by traditional MLF are the time-invariant linear state feedback in each mode while the controllers designed by the MDLF are piecewise continuous for each mode.

While shear flows in laboratory and geo- and astrophysical fluids abound, are often unstable, and have been studied both experimentally and numerically, a reliable understanding of the nonlinear saturation of the instability remains elusive. The common assumption of instability-saturation is that the unstable-mode energy, in its entirety, is cascaded to small scales where it is dissipated. Though the cascade exists, I will show, using novel eigenmode-coupled energy transfer analyses, that in both two- and three-dimensional driven shear-flow turbulence almost all the unstable-mode energy is transferred to large-scale linearly stable modes at the instability scale itself. These stable modes then return the turbulent energy to the background gradient. This finding is consequential also to magnetized turbulence where removing the stable modes of Kelvin-Helmholtz-instability in numerical simulations enhances, by an order of magnitude, both the small-scale cascade and dissipation, including that of magnetic fluctuations. Detailed examination shows that the flow composed of the unstable modes, in stable modes’ absentia, rapidly distorts and folds the magnetic field lines [1]. In the shear layer, up-gradient momentum transport by the stable modes cancels 70%-90% of down-gradient momentum transport by the unstable modes [2]. A model based on only a few unstable and stable modes accurately reproduces the transport rates from full-scale turbulence simulations. With this new instability-saturation paradigm, inroads to old and important problems of shear-driven dynamos and reconnection-driven outflows may be possible.

Visco-resistive magnetohydrodynamic turbulence, driven by a two-dimensional unstable shear layer that is maintained by an imposed body force, is examined by decomposing it into dissipationless linear eigenmodes of the initial profiles. The down-gradient momentum flux, as expected, originates from the large-scale instability. However, continual up-gradient momentum transport by large-scale linearly stable but nonlinearly excited eigenmodes is identified and found to nearly cancel the down-gradient transport by unstable modes. The stable modes effectuate this by depleting the large-scale turbulent fluctuations via energy transfer to the mean flow. This establishes a physical mechanism underlying the long-known observation that coherent vortices formed from nonlinear saturation of the instability reduce turbulent transport and fluctuations, as such vortices are composed of both the stable and unstable modes, which are nearly equal in their amplitudes. The impact of magnetic fields on the nonlinearly excited stable modes is then quantified. Even when imposing a strong magnetic field that almost completely suppresses the instability, the up-gradient transport by the stable modes is at least two-thirds of the down-gradient transport by the unstable modes, whereas for weaker fields, this fraction reaches up to 98%. These effects are persistent with variations in magnetic Prandtl number and forcing strength. Finally, continuum modes are shown to be energetically less important, but essential for capturing the magnetic fluctuations and Maxwell stress. A simple analytical scaling law is derived for their saturated turbulent amplitudes. It predicts the falloff rate as the inverse of the Fourier wavenumber, a property which is confirmed in numerical simulations.

In this work, we study finite-time stability of hybrid systems with unstable modes. We present sufficient conditions in terms of multiple Lyapunov functions for the origin of a class of hybrid systems to be finite-time stable. More specifically, we show that even if the value of the Lyapunov function increases during continuous flow, i.e., if the unstable modes in the system are active for some time, finite-time stability can be guaranteed if the finite-time convergent mode is active for a sufficient amount of cumulative time. This is the first work on finite-time stability of hybrid systems using multiple Lyapunov functions. Prior work uses a common Lyapunov function approach, and requires the Lyapunov function to be decreasing during the continuous flows and non-increasing at the discrete jumps, thereby, restricting the hybrid system to only have stable modes, or to only evolve along the stable modes. In contrast, we allow Lyapunov functions to increasebothduring the continuous flows and the discrete jumps. As thus, the derived stability results are less conservative compared to the earlier results in the related literature, and in effect allow the hybrid system to have unstable modes.

Studies on the valveless scheme to produce high-frequency detonations with different purge methods

We examine the geometry of the finite-dimensional attractor associated with fluid flows described by Navier–Stokes equations and relate its nonlinear dimensionality to energy exchanges between dynamical components (modes) of the flow. Specifically, we use a stochastic framework based on the dynamically orthogonal equations to perform efficient order-reduction and describe the stochastic attractor in the reduced-order phase space in terms of the associated probability measure. We introduce the notion of local fractal dimensionality to describe the geometry of the attractor and we establish a connection with the number of positive finite-time Lyapunov exponents. Subsequently, we illustrate in specific fluid flows that the low dimensionality of the stochastic attractor is caused by the synergistic activity of linearly unstable and stable modes as well as the action of the quadratic terms. In particular, we illustrate the connection of the low-dimensionality of the attractor with the circulation of energy: (i) from the mean flow to the unstable modes (due to their linearly unstable character), (ii) from the unstable modes to the stable ones (due to a nonlinear energy transfer mechanism) and (iii) from the stable modes back to the mean (due to the linearly stable character of these modes).

Decreasing skin friction in boundary layers attached to aircraft wings can have an impact in both fuel consumption and pollutant production, which are becoming crucial to reduce operation costs and meet environmental regulations, respectively. Skin friction in turbulent boundary layers is about ten times that of laminar boundary layers. Thus, an obvious method to reduce friction drag is to delay transition to turbulence, which is a fairly involved process in real aircraft wings [J98]. Transition sis promoted either by Tollmien—Schlichting (TS) and Klebanov (K) modes [K94], with the former playing an essential role. Various methods (e.g., suction [SG00,ZLB04], wave cancellation [WAA01,LG06]) have been proposed to reduce TS modes in laminar boundary layers. Mode interaction methods have been successfully used in fluid systems to control related instabilities, such as the Rayleigh—Taylor instability [LMV01]. Here, we present some recent results on using these methods to control TS modes in a compressible, 2D boundary layer over a flat plate at zero incidence. A given unstable TS mode can be stabilized by coupling its spatial evolution with that of a second selected stable TS mode, in such a way that the stable mode takes energy from the unstable one and gives a stable coupled evolution of both modes. The coupling device is a wavetrain in the boundary layer, with appropriate wavenumber and frequency, which can be created by an array of oscillators on the wall, and promotes both (i) parametric coupling between the stable and unstable TS modes and (ii) a mean flow that is also stabilizing. Three differences with wave cancelation methods are relevant. Namely, (a) nonlinear terms play an essential role in the process; (b) the unstable TS mode is stabilized (its growth rate is decreased), not just canceled; and (c) stabilization does not depend on the phase of the incoming wave, which implies that active control is not necessary.

Limit-cycle oscillations are studied for ion temperature gradient turbulence, which, in the absence of large diamagnetic (mean) shear flows, saturates through energy transfer from unstable modes to large-scale stable modes via zonal-flow intermediary modes. Oscillations of zonal flow and turbulence levels are strongly constrained by the reactive, largely non-dissipative character of the zonal flows. Since existing predator–prey models for observed oscillations in experiments do not include energy transfer through zonal flows to stable modes, low-order fluid models with this physics are constructed and investigated. A simple three-wave truncation produces low-amplitude zonal flows that slowly oscillate around a zero mean, with turbulence oscillations between coupled wavenumbers that exceed linear frequencies by orders of magnitude. This inconsistency with experimental observations is caused by the weak non-linear drive of zonal flows in three-wave systems and the lack of multiple-wavenumber turbulent interactions. A more comprehensive model that preserves multiple wavenumber interactions within the context of conservative zonal-flow-mediated energy transfer to stable modes accurately reflects observed dynamics when the phase between stable and unstable modes is occasionally randomized.

The longitudinal dynamic flight stability of a hovering bumblebee was studied using the method of computational fluid dynamics to compute the aerodynamic derivatives and the techniques of eigenvalue and eigenvector analysis for solving the equations of motion. For the longitudinal disturbed motion, three natural modes were identified: one unstable oscillatory mode, one stable fast subsidence mode and one stable slow subsidence mode. The unstable oscillatory mode consists of pitching and horizontal moving oscillations with negligible vertical motion. The period of the oscillations is 0.32 s (approx. 50 times the wingbeat period of the bumblebee). The oscillations double in amplitude in 0.1 s; coupling of nose-up pitching with forward horizontal motion (and nose-down pitching with backward horizontal motion) in this mode causes the instability. The stable fast subsidence mode consists of monotonic pitching and horizontal motions, which decay to half of the starting values in 0.024 s. The stable slow subsidence mode is mainly a monotonic descending (or ascending) motion, which decays to half of its starting value in 0.37 s. Due to the unstable oscillatory mode, the hovering flight of the bumblebee is dynamically unstable. However, the instability might not be a great problem to a bumblebee that tries to stay hovering: the time for the initial disturbances to double (0.1 s) is more than 15 times the wingbeat period (6.4 ms), and the bumblebee has plenty of time to adjust its wing motion before the disturbances grow large.

Nonlinear normal modes of a fixed-fixed buckled beam about its first post-buckling configuration are investigated. The cases of three-to-one and one-to-one internal resonances are analyzed. Approximate solutions for the nonlinear normal modes are computed by applying the method of multiple scales directly to the governing integral-partial-differential equation and associated boundary conditions. Curves displaying variation of the amplitude of one of the modes with the internal-resonance-detuning parameter are generated. It is shown that, for a three-to-one internal resonance between the first and third modes, the beam may possess one stable uncoupled mode (high-frequency mode) and either (a) one stable coupled mode, (b) three stable coupled modes, or (c) two stable and one unstable coupled modes. For the same resonance, the beam possesses one degenerate mode (with a multiplicity of two) and two stable and one unstable coupled modes. On the other hand, for a one-to-one internal resonance between the first and second modes, the beam possesses (a) two stable uncoupled modes and two stable and two unstable coupled modes; (b) one stable and one unstable uncoupled modes and two stable and two unstable coupled modes; and (c) two stable uncoupled and two unstable coupled modes (with a multiplicity of two). For a one-to-one internal resonance between the third and fourth modes, the beam possesses (a) two stable uncoupled modes and four stable coupled modes; (b) one stable and one unstable uncoupled modes and four stable coupled modes; (c) two unstable uncoupled modes and four stable coupled modes; and (d) two stable uncoupled modes and two stable coupled modes (each with a multiplicity of two).

The modes in the electrocardiogram (ECG) signal can be divided into stable modes and unstable modes. The unstable modes are of great significance for signal analysis and classification. When the traditional dynamic mode decomposition (DMD) method is directly applied to these signals, it is difficult to effectively extract these key modes due to the rank mismatch problem of the data matrix. In order to better capture and analyze unstable modes, this study introduced the Hankel matrix to expand ECG data and used it as the input of DMD to propose the Hankel dynamic mode decomposition (HDMD) method. Although HDMD has been applied in other physiological signal processing, this study is the first to successfully apply it to multi-lead ECG signal analysis. By optimizing the delay parameters of the Hankel matrix and retaining the modes that account for 90% of the singular value decomposition energy, we significantly improve the effectiveness of feature extraction. Each type of abnormal ECG signal data is classified on the PTB Diagnostic ECG database. The experimental results of the proposed model are compared with the direct use of DMD for modal extraction. The highest classification accuracy of multi-modal ECG signal data extracted by HDMD is more than 10% higher than that using DMD. At the same time, the mean squared error (MSE) of the reconstruction using HDMD is 0.282 lower than that of DMD, indicating a significant improvement in reconstruction accuracy. Further illustrating the proposed method, the HDMD method can better capture and analyze the unstable modes in the ECG signal, thereby significantly improving the accuracy and robustness of signal classification. Future work will focus on further optimizing the parameter selection of the HDMD model and exploring its application potential in real-time heart disease monitoring and early warning systems.

The overall pattern of normal modes of parallel flow of inviscid stratified fluid is examined. For a given flow and wavenumber the modes are divided into five classes, some of which may be empty: (i) a finite class of non-singular unstable modes; (ii) a conjugate finite class of non-singular damped stable modes; (iii) a finite class of singular stable modes, each of these having a branch point and being the limit of unstable modes; (iv) a discrete class of modified internal gravity waves, these being non-singular stable modes (if the density decreases with height everywhere); (v) a continuous class of singular stable modes. The modified internal gravity waves are described asymptotically for large values of the Richardson number. These asymptotic results are related to and extended by numerical calculations for a sinusoidal basic velocity profile and a Bickley jet. The wave speeds for small values of the Richardson number are found to depend only upon the local behaviour of the mean flow near an overall simple maximum or minimum of the velocity profile. Finally some difficulties in the use of the Howard formula for perturbation at a curve of marginal stability are elucidated.

Using weakly non-linear theory of oscillation, we estimate the amplitudes of low frequency modes in a slowly pulsating B (SPB) star, taking account of the effects of rotation on the modes. Applying the formulation by Schenk et al (2002), we compute non-linear coupling coefficient between the low frequency modes and estimate the equilibrium amplitudes of the modes excited in the star, assuming the amplitudes of the unstable modes are saturated as a result of non-linear coupling with stable modes, that is, as a result of parametric instability expected between one unstable mode and two stable modes. We use the traditional approximation to calculate adiabatic and non-adiabatic oscillations in a rotating star. We find $r$-modes in a rapidly rotating star play a significant role in the amplitude determination through non-linear coupling. We also find that for low $m$ modes, the fractional amplitudes of the radiative luminosity caused by the low frequency modes are of order $10^{-4}$ to $10^{-3}$ at the surface.

Kink mode m = 1 in magnetic tube with discontinuous magnetic field

This technical note focuses on a class of switched nonlinear systems with unstable modes and switching uncertainties. Both prescribed switching instants and switching sequence of the nominal switching signal would change. Two novel indexes named time changing ratios and mode changing ratios are proposed, based on which several conditions that fully utilize the trade-off among stable and unstable modes are established to achieve the robust stability of the switched system. An example of multi-agent systems illustrates the efficiency of the new results.