Stabilization of linear time-varying systems with state and input constraints using convex optimization

The stabilization problem of linear time-varying systems with both state and input constraints is considered. Sufficient conditions for the existence of the solution to this problem are derived and a gain-switched（gain-scheduled） state feedback control scheme is built to stabilize the constrained timevarying system. The design problem is transformed to a series of convex feasibility problems which can be solved efficiently. A design example is given to illustrate the effect of the proposed algorithm.

A High-Order Semi-Lagrangian Finite Difference Method for Nonlinear Vlasov and BGK Models

In this paper,we propose a new conservative high-order semi-Lagrangian finite difference(SLFD)method to solve linear advection equation and the nonlinear Vlasov and BGK models.The finite difference scheme has better computational flexibility by working with point values,especially when working with high-dimensional problems in an operator splitting setting.The reconstruction procedure in the proposed SLFD scheme is motivated from the SL finite volume scheme.In particular,we define a new sliding average function,whose cell averages agree with point values of the underlying function.By developing the SL finite volume scheme for the sliding average function,we derive the proposed SLFD scheme,which is high-order accurate,mass conservative and unconditionally stable for linear problems.The performance of the scheme is showcased by linear transport applications,as well as the nonlinear Vlasov-Poisson and BGK models.Furthermore,we apply the Fourier stability analysis to a fully discrete SLFD scheme coupled with diagonally implicit Runge-Kutta(DIRK)method when applied to a stiff two-velocity hyperbolic relaxation system.Numerical stability and asymptotic accuracy properties of DIRK methods are discussed in theoretical and computational aspects.

Influence of variable viscosity and double diffusion on the convective stability of a nanofluid flow in an inclined porous channel

The influence of variable viscosity and double diffusion on the convective stability of a nanofluid flow in an inclined porous channel is investigated.The DarcyBrinkman model is used to characterize the fluid flow dynamics in porous materials.The analytical solutions are obtained for the unidirectional and completely developed flow.Based on a normal mode analysis,the generalized eigenvalue problem under a perturbed state is solved.The eigenvalue problem is then solved by the spectral method.Finally,the critical Rayleigh number with the corresponding wavenumber is evaluated at the assigned values of the other flow-governing parameters.The results show that increasing the Darcy number,the Lewis number,the Dufour parameter,or the Soret parameter increases the stability of the system,whereas increasing the inclination angle of the channel destabilizes the flow.Besides,the flow is the most unstable when the channel is vertically oriented.

Implications for fault reactivation and seismicity induced by hydraulic fracturing

Evaluating the physical mechanisms that link hydraulic fracturing(HF) operations to induced earthquakes and the anticipated form of the resulting events is significant in informing subsurface fluid injection operations. Current understanding supports the overriding role of the effective stress magnitude in triggering earthquakes, while the impact of change rate of effective stress has not been systematically addressed. In this work, a modified critical stiffness was brought up to investigate the likelihood, impact,and mitigation of induced seismicity during and after hydraulic fracturing by developing a poroelastic model based on rate-and-state fraction law and linear stability analysis. In the new criterion, the change rate of effective stress was considered a key variable to explore the evolution of this criterion and hence the likelihood of instability slip of fault. A coupled fluid flow-deformation model was used to represent the entire hydraulic fracturing process in COMSOL Multiphysics. The possibility of triggering an earthquake throughout the entire hydraulic fracturing process, from fracturing to cessation, was investigated considering different fault locations, orientations, and positions along the fault. The competition between the effects of the magnitude and change rate of effective stress was notable at each fracturing stage. The effective stress magnitude is a significant controlling factor during fracturing events, with the change rate dominating when fracturing is suddenly started or stopped. Instability dominates when the magnitude of the effective stress increases(constant injection at each fracturing stage) and the change rate of effective stress decreases(the injection process is suddenly stopped). Fracturing with a high injection rate, a fault adjacent to the hydraulic fracturing location and the position of the junction between the reservoir and fault are important to reduce the Coulomb failure stress(CFS) and enhance the critical stiffness as the significant disturbance of stresses at these positions in the coupled process. Therefore,notable attention should be given to the injection rate during fracturing, fault position, and position along faults as important considerations to help reduce the potential for induced seismicity. Our model was verified and confirmed using the case of the Longmaxi Formation in the Sichuan Basin, China, in which the reported microseismic data were correlated with high critical stiffness values. This work supplies new thoughts of the seismic risk associated with HF engineering.

Effects of connected automated vehicle on stability and energy consumption of heterogeneous traffic flow system

With the development of intelligent and interconnected traffic system,a convergence of traffic stream is anticipated in the foreseeable future,where both connected automated vehicle(CAV)and human driven vehicle(HDV)will coexist.In order to examine the effect of CAV on the overall stability and energy consumption of such a heterogeneous traffic system,we first take into account the interrelated perception of distance and speed by CAV to establish a macroscopic dynamic model through utilizing the full velocity difference(FVD)model.Subsequently,adopting the linear stability theory,we propose the linear stability condition for the model through using the small perturbation method,and the validity of the heterogeneous model is verified by comparing with the FVD model.Through nonlinear theoretical analysis,we further derive the KdV-Burgers equation,which captures the propagation characteristics of traffic density waves.Finally,by numerical simulation experiments through utilizing a macroscopic model of heterogeneous traffic flow,the effect of CAV permeability on the stability of density wave in heterogeneous traffic flow and the energy consumption of the traffic system is investigated.Subsequent analysis reveals emergent traffic phenomena.The experimental findings demonstrate that as CAV permeability increases,the ability to dampen the propagation of fluctuations in heterogeneous traffic flow gradually intensifies when giving system perturbation,leading to enhanced stability of the traffic system.Furthermore,higher initial traffic density renders the traffic system more susceptible to congestion,resulting in local clustering effect and stop-and-go traffic phenomenon.Remarkably,the total energy consumption of the heterogeneous traffic system exhibits a gradual decline with CAV permeability increasing.Further evidence has demonstrated the positive influence of CAV on heterogeneous traffic flow.This research contributes to providing theoretical guidance for future CAV applications,aiming to enhance urban road traffic efficiency and alleviate congestion.

Linear stability of a fluid channel with a porous layer in the center

We perform a Poiseuille flow in a channel linear stability analysis of a inserted with one porous layer in the centre, and focus mainly on the effect of porous filling ratio. The spectral collocation technique is adopted to solve the coupled linear stability problem. We investigate the effect of permeability, σ, with fixed porous filling ratio ψ = 1/3 and then the effect of change in porous filling ratio. As shown in the paper, with increasing σ, almost each eigenvalue on the upper left branch has two subbranches at ψ = 1/3. The channel flow with one porous layer inserted at its middle （ψ = 1/3） is more stable than the structure of two porous layers at upper and bottom walls with the same parameters. By decreasing the filling ratio ψ, the modes on the upper left branch are almost in pairs and move in opposite directions, especially one of the two unstable modes moves back to a stable mode, while the other becomes more instable. It is concluded that there are at most two unstable modes with decreasing filling ratio ψ. By analyzing the relation between ψ and the maximum imaginary part of the streamwise phase speed, Cimax, we find that increasing Re has a destabilizing effect and the effect is more obvious for small Re, where ψ a remarkable drop in Cimax can be observed. The most unstable mode is more sensitive at small filling ratio ψ, and decreasing ψ can not always increase the linear stability. There is a maximum value of Cimax which appears at a small porous filling ratio when Re is larger than 2 000. And the value of filling ratio 0 corresponding to the maximum value of Cimax in the most unstable state is increased with in- creasing Re. There is a critical value of porous filling ratio （= 0.24） for Re = 500; the structure will become stable as ψ grows to surpass the threshold of 0.24; When porous filling ratio ψ increases from 0.4 to 0.6, there is hardly any changes in the values of Cimax. We have also observed that the critical Reynolds number is especially sensitive for small ψ where the fastest drop is observed, and there may be a wide range in which the porous filling ratio has less effect on the stability （ψ ranges from 0.2 to 0.6 at σ = 0.002）. At larger permeability, σ, the critical Reynolds number tends to converge no matter what the value of porous filling ratio is.

Linear spatial instability analysis in 3D boundary layers using plane-marching 3D-LPSE

It is widely accepted that a robust and efficient method to compute the linear spatial amplified rate ought to be developed in three-dimensional （3D） boundary layers to predict the transition with the e^N method, especially when the boundary layer varies significantly in the spanwise direction. The 3D-linear parabolized stability equation （3D- LPSE） approach, a 3D extension of the two-dimensional LPSE （2D-LPSE）, is developed with a plane-marching procedure for investigating the instability of a 3D boundary layer with a significant spanwise variation. The method is suitable for a full Mach number region, and is validated by computing the unstable modes in 2D and 3D boundary layers, in both global and local instability problems. The predictions are in better agreement with the ones of the direct numerical simulation （DNS） rather than a 2D-eigenvalue problem （EVP） procedure. These results suggest that the plane-marching 3D-LPSE approach is a robust, efficient, and accurate choice for the local and global instability analysis in 2D and 3D boundary layers for all free-stream Mach numbers.

Linear stability of plane creeping Couette flow for Burgers fluid

It is well known that plane creeping Couette flow of UCM and Oldroy-B fluids are linearly stable. However, for Burges fluid, which includes UCM and Oldroyd-B fluids as special cases, unstable modes are detected in the present work. The wave speed, critical parameters and perturbation mode are studied for neutral waves. Energy analysis shows that the sustaining of perturbation energy in Poiseuille flow and Couette flow is completely different. At low Reynolds number limit, analytical solutions are obtained for simpli- fied perturbation equations. The essential difference between Burgers fluid and Oldroyd-B fluid is revealed to be the fact that neutral mode exists only in the former.

The Linear Stability of a Solution of a Parallel Redundant Repairable System

In this paper,the solution of a parallel redundant repairable system is investigated.by using the method of functional analysis.Especially,the linear semigroups of operator theory on Banach space,we prove the existence of the strictly dominant eigenvalue,and show the linear stability of solution.

Exact Solutions of Discrete Complex Cubic Ginzburg-Landau Equation and Their Linear Stability

The discrete complex cubic Ginzburg-Landau equation is an important model to describe a number of physical systems such as Taylor and frustrated vortices in hydrodynamics and semiconductor laser arrays in optics. In this paper, the exact solutions of the discrete complex cubic Ginzburg-Landau equation are derived using homogeneous balance principle and the GI/G-expansion method, and the linear stability of exact solutions is discussed.

Exact solutions and linear stability analysis for two-dimensional Ablowitz Ladik equation

The Ablowitz-Ladik equation is a very important model in nonlinear mathematical physics. In this paper, the hyper- bolic function solitary wave solutions, the trigonometric function periodic wave solutions, and the rational wave solutions with more arbitrary parameters of two-dimensional Ablowitz-Ladik equation are derived by using the （GI/G）-expansion method, and the effects of the parameters （including the coupling constant and other parameters） on the linear stability of the exact solutions are analysed and numerically simulated.

Fundamental and dressed annular solitons in saturable nonlinearity with parity–time symmetric Bessel potential

We theoretically study the existence and stability of optical solitons in saturable nonlinearity with a two-dimensional parity–time（PT） symmetric Bessel potential.Besides the fundamental solitons,a novel type of dressed soliton,whose intensity looks like a ring dressed on an intensity hump,are presented.It is found that both the fundamental solitons and dressed solitons can exist when the propagation constant is beyond a certain critical value.The propagation stability is investigated with a linear stability analysis corroborated by a beam propagation method.All the fundamental solitons are stable,while dressed solitons are unstable for low values of saturable parameter.As the value of saturable parameter increases,the dressed solitons tend to be stable at high powers.

Optical solitons supported by finite waveguide lattices with diffusive nonlocal nonlinearity

We investigate the properties of fundamental,multi-peak,and multi-peaked twisted solitons in three types of finite waveguide lattices imprinted in photorefractive media with asymmetrical diffusion nonlinearity.Two opposite soliton selfbending signals are considered for different families of solitons.Power thresholdless fundamental and multi-peaked solitons are stable in the low power region.The existence domain of two-peaked twisted solitons can be changed by the soliton self-bending signals.When solitons tend to self-bend toward the waveguide lattice,stable two-peaked twisted solitons can be found in a larger region in the middle of their existence region.Three-peaked twisted solitons are stable in the lower(upper)cutoff region for a shallow(deep)lattice depth.Our results provide an effective guidance for revealing the soliton characteristics supported by a finite waveguide lattice with diffusive nonlocal nonlinearity.

ON THE STABILITY OF DISTORTED LAMINAR FLOW(Ⅱ)——THE LINEAR STABILITY ANALYSIS OF DISTORTED PARALLEL SHEAR FLOW

Based on the hydrodynamic stability theory of distorted laminar flow and the kind of distortion profiles on the mean velocity in parallel shear flow given in paper [1], this paper investigates the linear stability behaviour of parallel shear flow, presents unstable results of plane Couette flow and pipe Poiseuille flow to two-dimensional or axisymmetric disturbances for the first time, and obtains neutral curves of these two motions under certain definition.

Linear stability analysis of supersonic axisymmetric jets

Stabilities of supersonic jets are examined with different velocities, momentum thicknesses, and core temperatures. Amplification rates of instability waves at inlet are evaluated by linear stability theory （LST）. It is found that increased velocity and core temperature would increase amplification rates substantially and such influence varies for different azimuthal wavenumbers. The most unstable modes in thin momentum thickness cases usually have higher frequencies and azimuthal wavenumbers. Mode switching is observed for low azimuthal wavenumbers, but it appears merely in high velocity cases. In addition, the results provided by linear parabolized stability equations show that the mean-flow divergence affects the spatial evolution of instability waves greatly. The most amplified instability waves globally are sometimes found to be different from that given by LST.

Preventing Pressure Oscillations Does Not Fix Local Linear Stability Issues of Entropy-Based Split-Form High-Order Schemes

Recently,it was discovered that the entropy-conserving/dissipative high-order split-form discontinuous Galerkin discretizations have robustness issues when trying to solve the sim-ple density wave propagation example for the compressible Euler equations.The issue is related to missing local linear stability,i.e.,the stability of the discretization towards per-turbations added to a stable base flow.This is strongly related to an anti-diffusion mech-anism,that is inherent in entropy-conserving two-point fluxes,which are a key ingredi-ent for the high-order discontinuous Galerkin extension.In this paper,we investigate if pressure equilibrium preservation is a remedy to these recently found local linear stability issues of entropy-conservative/dissipative high-order split-form discontinuous Galerkin methods for the compressible Euler equations.Pressure equilibrium preservation describes the property of a discretization to keep pressure and velocity constant for pure density wave propagation.We present the full theoretical derivation,analysis,and show corresponding numerical results to underline our findings.In addition,we characterize numerical fluxes for the Euler equations that are entropy-conservative,kinetic-energy-preserving,pressure-equilibrium-preserving,and have a density flux that does not depend on the pressure.The source code to reproduce all numerical experiments presented in this article is available online(https://doi.org/10.5281/zenodo.4054366).

Establishment of infinite dimensional Hamiltonian system of multilayer quasi-geostrophic flow & study on its linear stability

A multilayer flow is a stratified fluid composed of a finite number of layers with densities homogeneous within one layer but different from each other. It is an intermediate system between the single-layer barotropic model and the continuously stratified baroclinic model. Since this system can simulate the baroclinic effect simply, it is widely used to study the large-scale dynamic process in atmosphere and ocean. The present paper is concerned with the linear stability of the multilayer quasi-geostrophic flow, and the associated linear stability criteria are established. Firstly, the nonlinear model is turned into the form of a Hamiltonian system, and a basic flow is defined. But it cannot be an extreme point of the Hamiltonian function since the system is an infinite-dimensional one. Therefore, it is necessary to reconstruct a new Hamiltonian function so that the basic flow becomes an extreme point of it. Secondly, the linearized equations of disturbances in the multilayer quasi-geostrophic flow are derived by introducing infinitesimal disturbances superposed on the basic flows. Finally, the properties of the linearized system are discussed, and the linear stability criteria in the sense of Liapunov are derived under two different conditions with respect to certain norms.

Dynamic Stability Analysis of Linear Time-varying Systems via an Extended Modal Identification Approach

The problem of linear time-varying（LTV） system modal analysis is considered based on time-dependent state space representations, as classical modal analysis of linear time-invariant systems and current LTV system modal analysis under the ＂frozen-time＂ assumption are not able to determine the dynamic stability of LTV systems. Time-dependent state space representations of LTV systems are first introduced, and the corresponding modal analysis theories are subsequently presented via a stabilitypreserving state transformation. The time-varying modes of LTV systems are extended in terms of uniqueness, and are further interpreted to determine the system＇s stability. An extended modal identification is proposed to estimate the time-varying modes, consisting of the estimation of the state transition matrix via a subspace-based method and the extraction of the time-varying modes by the QR decomposition. The proposed approach is numerically validated by three numerical cases, and is experimentally validated by a coupled moving-mass simply supported beam exper- imental case. The proposed approach is capable of accurately estimating the time-varying modes, and provides anew way to determine the dynamic stability of LTV systems by using the estimated time-varying modes.

Linear global stability of a confined plume

A linear stability analysis is performed for a plume flow inside a cylinder of aspect ratio 1. The configu- ration is identical to that used by Lopez and Marques （2013） for their direct numerical simulation study, It is found that the first bifurcation, which leads to a periodic axisymmetric flow state, is accurately pre- dicted by linear analysis： both the critical Rayleigh number and the global frequency are consistent with the reported DNS results. It is further shown that pressure feedback drives the global mode, rather than absolute instability.

Stability of Three-Dimensional Interfacial Waves Under Subharmonic Disturbances

This study examines the stability regimes of three-dimensional interfacial gravity waves.The numerical results of the linear stability analysis extend the three-dimensional surface waves results of Ioualalen and Kharif(1994)to three-dimensional interfacial waves.An approach of the collocation type has been developed for this purpose.The equations of motion are reduced to an eigenvalue problem where the perturbations are spectrally decomposed into normal modes.The results obtained showed that the density ratio plays a stabilizing factor.In addition,the dominant instability is of three-dimensional structure,and it belongs to class I for all values of density ratio.