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.

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.

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.

Instability characteristics of a co-rotating wingtip vortex pair based on bi-global linear stability analysis

The Stereo Particle Image Velocimetry(SPIV)technology is applied to measure the wingtip vortices generated by the up-down symmetrical split winglet.Then,the temporal biglobal Linear Stability Analysis(bi-global LSA)is performed on this nearly equal-strength corotating vortex pair,which is composed of an upper vortex(vortex-u)and a down vortex(vortex-d).The results show that the instability eigenvalue spectrum illustrated by(ωr,ω_(i))contains two types of branches:discrete branch and continuous branch.The discrete branch contains the primary branches of vortex-u and vortex-d,the secondary branch of vortex-d and coupled branch,of which all of the eigenvalues are located in the unstable half-plane ofω_(i)>0,indicating that the wingtip vortex pair is temporally unstable.By contrast,the eigenvalues of the continuous branch are concentrated on the half-plane ofω_(i)<0 and the perturbation modes correspond to the freestream perturbation.In the primary branches of vortex-u and vortex-d,Mode P_(u) and Mode Pd are the primary perturbation modes,which exhibit the structures enclosed with azimuthal wavenumber m and radial wavenumber n,respectively.Besides,the results of stability curves for vortex-u and vortex-d demonstrate that the instability growth rates of vortex-u are larger than those of vortex-d,and the perturbation energy of Mode P_(u) is also larger than that of Mode Pd.Moreover,the perturbation energy of Mode P_(u) is up to 0.02650 and accounts for 33.56%percent in the corresponding branch,thereby indicating that the instability development of wingtip vortex is dominated by Mode P_(u).By further investigating the topological structures of Mode P_(u) and Mode Pd with streamwise wavenumbers,the most unstable perturbation mode with a large azimuthal wavenumber of m=5-6 is identified,which imposes on the entire core region of vortex-u.This large azimuthal wavenumber perturbation mode can suggest the potential physical-based flow control strategy by manipulating it.

Linear Stability Analysis on the Influences of the Spatial Variations in Thermal Conductivity and Expansivity on the Flow Patterns of Thermal Convection with Strongly Temperature-Dependent Viscosity

A series of linear stability analysis is carried out on the onset of thermal convection in the presence of spatial variations of viscosity, thermal conductivity and expansivity. We consider the temporal evolution of an infinitesimal perturbation superimposed to a static （motionless） and con- ductive state in a basally-heated planar layer. From the changes in flow patterns with increasing the amplitudes of temperature dependence of viscosity, we identified the transition into the ＂stagnant-lid＂ （ST） regime, where the convection occurs only beneath a thick and stagnant-lid of cold fluid at the top surface. Detailed analysis showed a significant increase of the aspect ratio of convection cells in ST regime induced by the spatial variations in thermal conductivity and/or expansivity： the horizon- tal length scale of ST convection can be enlarged by up to 50% with 10 times increase of thermal conductivity with depth. We further developed an analytical model of ST convection which success- fully reproduced the mechanism of increasing horizontal length scale of ST regime convection cells for given spatial variations in physical properties. Our findings may highlight the essential roles of the spatial variation of thermal conductivity on the convection patterns in the mantle.

Effect of multi-velocity-difference in traffic flow

Based on the optimal velocity models, an extended model is proposed, in which multi-veloclty-dllterence aheacl is taken into consideration. The damping effect of the multi-velocity-difference ahead has been investigated by means of analytical and numerical methods. Results indicate that the multi-velocity-difference leads to the enhancement of stability of traffic flow, suppression of the emergence of traffic jamming, and reduction of the energy consumption.

Linear Stability Analysis of Thermocapillary Flow in Rotating Shallow PoolsHeated from Inner Wall

Thermocapillary flow of silicon melt(Pr=0.011)in shallow annular pool heated from inner wall was simulated at the dimensionless rotation ratewranging from 0 to 7000.The effect of pool rotation on the stability of the thermocapillary flow was investigated.The steady axisymmetric basic state was solved by using the spectral element method;the critical stability parameters were determined by linear stability analysis;the mechanism of the flow instability was explored by the analysis of energy balance.A stability diagram,exhibiting the variation of the critical Marangoni number versus the dimensionless rotation ratewwas presented.The results reveal that only one Hopf bifurcation point appeared in the intervals ofω<3020 andω>3965,and the corresponding instability was caused by the shear energy,which was provided by the thermocapillary force and pool rotation,respectively.In addition,the competition between thermocapillary force and pool rotation leads to three Hopf bifurcation points in the range of 3020<ω<3965 with the increase of Marangoni number.

Interpretation of gas-film cooling against aero-thermal heating for high-speed vehicles

The possible application of the film-cooling technique against aero-thermal heating for surfaces of high-speed flying vehicles is discussed. The technique has been widely used in the heat protection of gas turbine blades. It is shown in this paper that, by applying this technique to high-speed flying vehicles, the working principle is fundamentally different. Numerical simulations for two model problems axe performed to support the argument. Besides the heat protection, the appreciable drag reduction is found to be another favorable effect. For the second model problem, i.e., the gas cooling for an optical window on a sphere cone, the hydrodynamic instability of the film is studied by the linear stability analysis to observe possible occurrence of laminar-turbulent transition.

Cross-diffusive effects on the onset of double-diffusive convection in a horizontal saturated porous fluid layer heated and salted from above

Double-diffusive stationary and oscillatory instabilities at the marginal state in a saturated porous horizontal fluid layer heated and salted from above are investigated theoretically under the Darcy＇s framework for a porous medium. The contributions of Soret and Dufour coefficients are taken into account in the analysis. Linear stability analysis shows that the critical value of the Darcy-Rayleigh number depends on cross-diffusive parameters at marginally stationary convec- tion, while the marginal state characterized by oscillatory convection does not depend on the cross-diffusion terms even if the condition and frequency of oscillatory convection depends on the cross-diffusive parameters. The critical value of the Darcy-Rayleigh number increases with increasing value of the solutal Darcy-Rayleigh number in the absence of cross- diffusive parameters. The critical Darcy-Rayleigh number decreases with increasing Soret number, resulting in destabiliza- tion of the system, while its value increases with increasing Dufour number, resulting in stabilization of the system at the marginal state characterized by stationary convection. The analysis reveals that the Dufour and Soret parameters as well as the porosity parameter play an important role in deciding the type of instability at the onset. This analysis also indicates that the stationary convection is followed by the oscillatory convection for certain fluid mixtures. It is interesting to note that the roles of cross-diffusive parameters on the double-diffusive system heated and salted from above are reciprocal to the double-diffusive system heated and salted from below.

Newly modified method and its application to the coupled Boussinesq equation in ocean engineering with its linear stability analysis

Investigating the dynamic characteristics of nonlinear models that appear in ocean science plays an important role in our lifetime.In this research,we study some features of the paired Boussinesq equation that appears for two-layered fluid flow in the shallow water waves.We extend the modified expansion function method(MEFM)to obtain abundant solutions,as well as to find new solutions.By using this newly modified method one can obtain novel and more analytic solutions comparing to MEFM.Also,numerical solutions via the Adomian decomposition scheme are discussed and favorable comparisons with analytical solutions have been done with an outstanding agreement.Besides,the instability modulation of the governing equations are explored through the linear stability analysis function.All new solutions satisfy the main coupled equation after they have been put into the governing equations.

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.

Effect of vortex dynamics and instability characteristics on the induced drag of trailing vortices

In this paper, trailing vortices generated by three wingtip configurations, namely the M6wing and the M6 wing with a blended or split winglet, are experimentally investigated using the Stereo Particle Image Velocimetry(SPIV) technology. Then, linear stability analysis is performed to investigate instability characteristics. Three corresponding trailing vortex patterns, including the isolated trailing vortex without wake(pattern v) and with wake(pattern v-w), co-rotating vortex pair(pattern v-v), are observed in experiments. The strength of trailing vortices, characterized by circulation, is reduced after installing winglets as expected, and the strength of pattern v-v can be further suppressed compared with pattern v-w. Moreover, instability characteristics, such as the eigenvalue spectrum and perturbation mode, are distinctive among these three vortex patterns.The distribution of eigenvalue spectrums indicates that pattern v and pattern v-w are temporally“marginally stable”, while pattern v-v is temporally “unstable.” The primary perturbation mode of pattern v and pattern v-w is the m =-1 helical mode, while |m|>1 for the case of pattern v-v.The effect of vortex dynamics and instability characteristics can be concluded in two aspects.Firstly, the value of induced drag is polluted by about 3% from vortex wandering since vortex wandering affects the tangential velocity and streamwise vorticity of trailing vortices. Secondly, the growth rate and penetration depth perturbation mode affect trailing vortex evolution and further affect induced drag. Specifically, the larger the growth rate and penetration depth are, the more turbulence injects inside the vortex core, thus leading to a quicker and more intense attenuation of trailing vortex, as well as a smaller induced drag. This finding can guide us to manipulate the induced drag in flow control.

Analysis and Computation for a Fluid Mixture Model

A fluid mixture model of tissue deformations has been studied in this paper.The model is a mixed system of nonlinear hyperbolic and elliptic partial differential equations.Both theoretical linear stability and numerical analysis are presented.Comparisons between standard numerical methods that utilize Runge-Kutta methods coupled with the WENO scheme and the immersed interface methods are given.Numerical examples are also presented.

Thermal Convection in a Spherical Shell with Melting/Freezing at either or both of Its Boundaries

In a number of geophysical or planetological settings, including Earth＇s inner core, a sili- cate mantle crystallizing from a magma ocean, or an ice shell surrounding a deep water ocean--a situa- tion possibly encountered in a number of Jupiter and Saturn＇s icy satellites--a convecting crystalline layer is in contact with a layer of its melt. Allowing for melting/freezing at one or both of the boundaries of the solid layer is likely to affect the pattern of convection in the layer. We study here the onset of thermal convection in a viscous spherical shell with dynamically induced melting/freezing at either or both of its boundaries. It is shown that the behavior of each interface---permeable or impermeable-- depends on the value of a dimensionless number P （one for each boundary）, which is the ratio of a melting/freezing timescale over a viscous relaxation timescale. A small value of P corresponds to perme- able boundary conditions, while a large value of P corresponds to impermeable boundary conditions. Linear stability analysis predicts a significant effect of semi-permeable boundaries when the number P characterizing either of the boundary is small enough： allowing for melting/freezing at either of the boundary allows the emergence of larger scale convective modes. The effect is particularly drastic when the outer boundary is permeable, since the degree 1 mode remains the most unstable even in the case of thin spherical shells. In the case of a spherical shell with permeable inner and outer boundaries, the most unstable mode consists in a global translation of the solid shell, with no deformation. In the limit of a full sphere with permeable outer boundary, this corresponds to the ＂convective translation＂ mode recently proposed for Earth＇s inner core. As another example of possible application, we discuss the case of thermal convection in Enceladus＇ ice shell assuming the presence of a global subsurface ocean, and found that melting/freezing could have an important effect on the pattern of convection in the ice shell.

Numerical research on airfoil transition delay by alternative current dielectric barrier discharge actuation

A Dielectric Barrier Discharge(DBD) plasma actuator can create a body force which locally accelerates the base flow leading to an attenuation of broadband disturbance to delay the transition. In this study, numerical simulation on an NLF0416 airfoil is conducted to investigate transition delay and drag reduction by a DBD plasma actuator. To simulate plasma’s effect more accurately, boundary-layer data is acquired from Reynolds Averaged Navier Stocks(RANS) equations instead of laminar boundary layer equations, although RANS equations need a much finer boundary-layer grid, and the linear stability analysis method is used to analyze the boundary layer and get the transition point. In this study, the influences of different actuation intensities and positions are investigated, and results show that if the actuation intensity is stronger and the actuation position is closer to the base transition point, more drag reduction can be obtained. However, the efficiency of plasma transition delay is really low. For example, when the actuation voltage is 16 k V,the actuation frequency is 1 k Hz, and the main Mach number is 0.1, the saved power due to drag reduction is about 5.09 W, but the power consumed is about 32.61 W, and the efficiency is just15.6%.

Modeling Magma Dynamics with a Mixed Fourier Collocation-Discontinuous Galerkin Method

A high-order discretization consisting of a tensor product of the Fourier collocation and discontinuous Galerkin methods is presented for numerical modeling of magma dynamics.The physical model is an advection-reaction type system consisting of two hyperbolic equations and one elliptic equation.The high-order solution basis allows for accurate and efficient representation of compaction-dissolution waves that are predicted from linear theory.The discontinuous Galerkin method provides a robust and efficient solution to the eigenvalue problem formed by linear stability analysis of the physical system.New insights into the processes of melt generation and segregation,such as melt channel bifurcation,are revealed from two-dimensional time-dependent simulations.