Stability analysis for a second-order continuous finite-time control system subject to a disturbance

In this paper, using finite-time control method, we consider the disturbance analysis of a second-order system with unknown but bounded disturbance. We show that the states of the second-order system will be stabilized to a region containing the origin. The radius of this region is determined by the control parameters and can be rendered as small as desired. The rigorous stability analysis is also given. Compared with the conventional PD control law, the finite-time control law yields a better disturbance rejection performance. Numerical simulation results show the effectiveness of the method.

Second-Order Adjoint Sensitivity Analysis Methodology for Computing Exactly Response Sensitivities to Uncertain Parameters and Boundaries of Linear Systems: Mathematical Framework

This work presents the “Second-Order Comprehensive Adjoint Sensitivity Analysis Methodology (2nd-CASAM)” for the efficient and exact computation of 1st- and 2nd-order response sensitivities to uncertain parameters and domain boundaries of linear systems. The model’s response (i.e., model result of interest) is a generic nonlinear function of the model’s forward and adjoint state functions, and also depends on the imprecisely known boundaries and model parameters. In the practically important particular case when the response is a scalar-valued functional of the forward and adjoint state functions characterizing a model comprising N parameters, the 2nd-CASAM requires a single large-scale computation using the First-Level Adjoint Sensitivity System (1st-LASS) for obtaining all of the first-order response sensitivities, and at most N large-scale computations using the Second-Level Adjoint Sensitivity System (2nd-LASS) for obtaining exactly all of the second-order response sensitivities. In contradistinction, forward other methods would require (N2/2 + 3 N/2) large-scale computations for obtaining all of the first- and second-order sensitivities. This work also shows that constructing and solving the 2nd-LASS requires very little additional effort beyond the construction of the 1st-LASS needed for computing the first-order sensitivities. Solving the equations underlying the 1st-LASS and 2nd-LASS requires the same computational solvers as needed for solving (i.e., “inverting”) either the forward or the adjoint linear operators underlying the initial model. Therefore, the same computer software and “solvers” used for solving the original system of equations can also be used for solving the 1st-LASS and the 2nd-LASS. Since neither the 1st-LASS nor the 2nd-LASS involves any differentials of the operators underlying the original system, the 1st-LASS is designated as a “first-level” (as opposed to a “first-order”) adjoint sensitivity system, while the 2nd-LASS is designated as a “second-level” (rather than a “second-order”) adjoint sensitivity system. Mixed second-order response sensitivities involving boundary parameters may arise from all source terms of the 2nd-LASS that involve the imprecisely known boundary parameters. Notably, the 2nd-LASS encompasses an automatic, inherent, and independent “solution verification” mechanism of the correctness and accuracy of the 2nd-level adjoint functions needed for the efficient and exact computation of the second-order sensitivities.

Static inelastic analysis of RC shear walls

A macro-model of a reinforced concrete (RC) shear wall is developed for static inelastic analysis.The model is composed of RC column elements and RC membrane elements.The column elements are used to model the boundary zone and the membrane elements are used to model the wall panel.Various types of constitutive relationships of concrete could be adopted for the two kinds of elements.To perform analysis,the wall is divided into layers along its height.Two adjacent layers are connected with a rigid beam.There are only three unknown displacement components for each layer.A method called single degree of freedom compensation is adopted to solve the peak value of the capacity curve.The post-peak stage analysis is performed using a forced iteration approach.The macro-model developed in the study and the complete process analysis methodology are verified by the experimental and static inelastic analytical results of four RC shear wall specimens.

Geometric Scaling Analysis of Deep Inelastic Scattering Data Including Heavy Quarks

An analytic massive total cross section of photon proton scattering is derived, which has geometric scaling. A geometric scaling is used to perform a global analysis of the deep inelastic scattering data on inclusive structure function F2 measured in lepton-hadron scattering experiments at small values of Bjorken x. It is shown that the descriptions of the inclusive structure function F2 and longitudinal structure function FL are improved with the massive analytic structure function, which may imply the gluon saturation effect dominating the parton evolution process at HERA. The inclusion of the heavy quarks prevent the divergence of the lepton-hadron cross section, which plays a significant role in the description of the photoproduction region.

Analysis of the Diffractive Deep Inelastic Scattering Data with Running Coupling and Gluon Number Fluctuations

We study the effects of running coupling and gluon number fluctuations in the latest diffractive deep inelastic scattering data. It is found that the description of the data is improved once the running coupling and gluon number fluctuations are included with x2/d.o.f. = 0.867, x2/d.o.f. = 0.923 and x2/d.o.f. = 0.878 for three different groups of experimental data. The values of diffusive coefficient subtracted from the fit are smaller than the ones obtained by considering only the gluon number fluctuations in our previous studies. The smaller values of the diffusive coefficient are in agreement with the theoretical predictions, where the gluon number fluctuations are suppressed by the running coupling which leads to smaller values of the diffusive coefficient.

Pharmacokinetic Analysis of Four Bioactive Iridoid and Secoiridoid Glycoside Components of Radix Gentianae Macrophyllae and Their Synergistic Excretion by HPLC-DAD Combined with Second-Order Calibration

An HPLC-DAD method combined with second-order calibration based on the alternating trilinear decomposition(ATLD)algorithm with the aid of region selection was developed to simultaneously and quantitatively characterize the synergistic relationships and cumulative excretion of the four bioactive ingredients of Radix Gentianae Macrophyllae in vivo.Although the analytes spectra substantially overlapped with that of the biological matrix,the overlapping profiles between analytes and co-eluting interferences can be successfully separated and accurately quantified by the ATLD method on the basis of the strength of region selection.The proposed approach not only determined the content change but also revealed the synergistic relationships and the cumulative excretion in vivo of the four ingredients in urine and feces samples collected at different excretion time intervals.In addition,several statistical parameters were employed to evaluate the accuracy and precision of the method.Quantitative results were confirmed by HPLC-mass spectrometry.Satisfactory results indicated that the proposed approach can be utilized to investigate the pharmacokinetics of Radix Gentianae Macrophyllae excretion in vivo.

Evaluating the inelastic displacement ratios of moment-resisting steel frames designed according to the Egyptian code

Seismic codes estimate the maximum displacements of building structures under the design-basis earthquakes by amplifying the elastic displacements under the reduced seismic design forces with a deflection amplification factor(DAF). The value of DAF is often estimated as ρ× R, where R is the force reduction factor and ρ is the inelastic displacement ratio that accounts for the inelastic action of the structure according to the definition presented by FEMA P695. The purpose of this study is to estimate the ρ-ratio of moment resisting steel frames(MRSFs) designed according to the Egyptian code. This is achieved by conducting a series of elastic and inelastic time-history analyses by two sets of earthquakes on four MRSFs designed according to the Egyptian code and having 2, 4, 8 and 12 stories. The earthquakes are scaled to produce maximum story drift ratios(MSDRs) of 1.0%, 1.5%, 2.0% and 2.5%. The mean values of the ρ-ratio are calculated based on the displacement responses of the investigated frames. The results obtained in this study indicate that the consideration of ρ for both the roof drift ratios(RDRs) and the MSDRs equal to 1.0 is a reasonable estimation for MRSFs designed according to the Egyptian code.

Second-order consensus in networks of dynamic agents with communication time-delays

This paper proposes second-order consensus protocols with time-delays and gives the measure of the robustness of the protocols to the time-delay existing in the network of agents with second-order dynamics. By employing a frequency domain method, it is proven that the information states and their time derivatives of all the agents in the network achieve consensus asymptotically, respectively, for appropriate communication timedelay if the topology of weighted network is connected. Particularly, a tight upper bound on the communication time-delay that can be tolerated in the dynamic network is found. The consensus protocols are distributed in the sense that each agent only needs information from its neighboring agents, which reduces the complexity of connections between neighboring agents significantly. Numerical simulation results are provided to demonstrate the effectiveness and the sharpness of the theoretical results for second-order consensus in networks in the presence of communication time-delays.

Heat transfer of nanofluids considering nanoparticle migration and second-order slip velocity

The heat transfer of a magnetohydrodynamics nanofluid inside an annulus considering the second-order slip condition and nanoparticle migration is theoret-ically investigated. A second-order slip condition, which appropriately represents the non-equilibrium region near the interface, is prescribed rather than the no-slip condition and the linear Navier slip condition. To impose different temperature gradients, the outer wall is subjected to q2, the inner wall is subjected to q1, and q1 〉 q2. A modified two-component four-equation non-homogeneous equilibrium model is employed for the nanofiuid, which have been reduced to two-point ordinary boundary value differential equations in the consideration of the thermally and hydrodynamically fully developed flow. The homotopy analysis method （HAM） is employed to solve the equations, and the h-curves are plotted to verify the accuracy and efficiency of the solutions. Moreover, the effects of the physical factors on the flow and heat transfer are discussed in detail, and the semi-analytical relation between NUB and NBT is obtained.

Two new predictor-corrector algorithms for second-order cone programming

Based on the ideas of infeasible interior-point methods and predictor-corrector algorithms, two interior-point predictor-corrector algorithms for the second-order cone programming （SOCP） are presented. The two algorithms use the Newton direction and the Euler direction as the predictor directions, respectively. The corrector directions belong to the category of the Alizadeh-Haeberly-Overton （AHO） directions. These algorithms are suitable to the cases of feasible and infeasible interior iterative points. A simpler neighborhood of the central path for the SOCP is proposed, which is the pivotal difference from other interior-point predictor-corrector algorithms. Under some assumptions, the algorithms possess the global, linear, and quadratic convergence. The complexity bound O（rln（εo/ε）） is obtained, where r denotes the number of the second-order cones in the SOCP problem. The numerical results show that the proposed algorithms are effective.

Robust Blind Separation for MIMO Systems against Channel Mismatch Using Second-Order Cone Programming

To improve the deteriorated capacity gain and source recovery performance due to channel mismatch problem,this paper reports a research about blind separation method against channel mismatch in multiple-input multiple-output(MIMO) systems.The channel mismatch problem can be described as a channel with bounded fluctuant errors due to channel distortion or channel estimation errors.The problem of blind signal separation/extraction with channel mismatch is formulated as a cost function of blind source separation(BSS) subject to the second-order cone constraint,which can be called as second-order cone programing optimization problem.Then the resulting cost function is solved by approximate negentropy maximization using quasi-Newton iterative methods for blind separation/extraction source signals.Theoretical analysis demonstrates that the proposed algorithm has low computational complexity and improved performance advantages.Simulation results verify that the capacity gain and bit error rate(BER) performance of the proposed blind separation method is superior to those of the existing methods in MIMO systems with channel mismatch problem.

A Multiscale Method for Damage Analysis of Quasi-Brittle Heterogeneous Materials

A novel multiscale algorithm based on the higher-order continuum at both micro-and macrostructural level is proposed for the consideration of the quasi-brittle damage response of heterogeneous materials.Herein,the microlevel damage is modelled by the degradation of the homogenized stress and tangent stiffness tensors,which are then upscaled to govern the localization at the macrolevel.The C^1 continuity finite element employing a modified case of Mindlin’s form II strain energy density is derived for the softening analysis.To the authors’knowledge,the finite element discretization based on the strain gradient theory is applied for the modeling of damage evolution at the microstructural level for heterogeneous materials for the first time.The advantage of the novel C1 finite element formulation in comparison with the standard finite element discretization in terms of the regularization efficiency as well as the objectivity has been shown.An isotropic damage law is used for the reduction of the constitutive and nonlocal material behaviour,which is necessary for the physically correct description of the localization formation in quasi-brittle materials.The capabilities of the derived finite element to capture the fully developed localization zones are tested on a random representative volume element(RVE)for several different loading cases.By employing the conventional second-order computational homogenization,the microstructural material constitutive response is averaged over the whole RVE area.In order to model the loss of structural integrity when sharp localization is formed across RVE,the specific conditions which detect a completely formed localization zone are developed.A new failure criterion at the microstructural level has been proposed.The derived finite element formulation,as well as the multiscale damage algorithm,are implemented into the finite element program ABAQUS.The capabilities of the presented multiscale scheme to capture the effects of the deformation localization are demonstrated by few benchmark numerical examples.

Illustrative Application of the 2^nd-Order Adjoint Sensitivity Analysis Methodology to a Paradigm Linear Evolution/Transmission Model: Reaction-Rate Detector Response

This work continues the illustrative application of the “Second Order Comprehensive Adjoint Sensitivity Analysis Methodology” (2nd-CASAM) to a benchmark mathematical model that can simulate the evolution and/or transmission of particles in a heterogeneous medium. The model response considered in this work is a reaction-rate detector response, which provides the average interactions of particles with the respective detector or, alternatively, the time-average of the concentration of a mixture of substances in a medium. The definition of this model response includes both uncertain boundary points of the benchmark, thereby providing both direct and indirect contributions to the response sensitivities stemming from the boundaries. The exact expressions for the 1st- and 2nd-order response sensitivities to the boundary and model parameters obtained in this work can serve as stringent benchmarks for inter-comparing the performances of all (deterministic and statistical) sensitivity analysis methods.

Illustrative Application of the 2^nd-Order Adjoint Sensitivity Analysis Methodology to a Paradigm Linear Evolution/Transmission Model: Point-Detector Response

This work illustrates the application of the “Second Order Comprehensive Adjoint Sensitivity Analysis Methodology” (2nd-CASAM) to a mathematical model that can simulate the evolution and/or transmission of particles in a heterogeneous medium. The model response is the value of the model’s state function (particle concentration or particle flux) at a point in phase-space, which would simulate a pointwise measurement of the respective state function. This paradigm model admits exact closed-form expressions for all of the 1st- and 2nd-order response sensitivities to the model’s uncertain parameters and domain boundaries. These closed-form expressions can be used to verify the numerical results of production and/or commercial software, e.g., particle transport codes. Furthermore, this paradigm model comprises many uncertain parameters which have relative sensitivities of identical magnitudes. Therefore, this paradigm model could serve as a stringent benchmark for inter-comparing the performances of all deterministic and statistical sensitivity analysis methods, including the 2nd-CASAM.

Progress of neutron induced prompt gamma analysis technique in 1988～2003

This paper described new development of the neutron induced prompt gamma-ray analysis(NIPGA) technology in 1988-2003.The pulse fast-thermal neutron activation analysis method,which utilized the inelastic reaction and capture reaction jointly,was employed to measure the elemental contents more efficiently.Lifetime of the neutron generator was more than 10000h and the performance of detector and MCA reached a high level.At the same time,Monte Carlo library least-square method was used to solve the nonlinearity problem in the NIPGA.

Dynamic Analysis of a Turret-Moored Tanker

In this paper, the second-order perturbation method in frequency domain is used to calculate RAO and spectra of motion and mooring line tension of a turret-moored tanker in ballast condition. The calculated results are compared with corresponding experiment results. In the experiment the wave heading is 180 degrees, and the wave spectra is the P-M spectrum and white noise spectrum. In the theoretical calculations, the damping coefficient of slow oscillation of the tanker is determined on the basis of the damping obtained from a test of irregular waves where the mooring system is replaced by a nonlinear spring with nonlinear stiffness similar to that of the mooring system. From the comparison between theoretical calculations and experiment al results, it can be found that the theoretical results obtained by the second-order perturbation method in frequency domain are in good agreement with the experimental results, indicating that the damping coefficient of slow oscillation of the tanker required in frequency domain calculation can be determined by reference to the damping coefficient of the tanker moored by a spring system in irregular waves, and the second-order perturbation method can be used to analyze the dynamic response of a turret-moored tanker.

Improved Capacity Spectrum Method Using Structural Inelastic Response

Existing problems of capacity-demand diagram were considered, and the conversion method between the results of the capacitydemand curves method and dynamic time-history analysis is presented. The factors causing the relative error between the different capacity curves, such as the distribution type for lateral load, assumptions of member bar model , surface irregularity of the structures and P-A effects, were analyzed to build the statistical correlation between the maximum base shear and the top displacements of the components under the different pushover analysis. ATC-40 （Procedure-B） method is verified to be more accurate in calculating the roof displacements of buildings. Meanwhile, the improved method is refined in three aspects and proved effective by a real concrete-freme building.

Three-dimensional mixed convection stagnation-point fow past a vertical surface with second-order slip velocity

This study is concerned with the three-dimensional(3D)stagnation-point for the mixed convection flow past a vertical surface considering the first-order and secondorder velocity slips.To the authors’knowledge,this is the first study presenting this very interesting analysis.Nonlinear partial differential equations for the flow problem are transformed into nonlinear ordinary differential equations(ODEs)by using appropriate similarity transformation.These ODEs with the corresponding boundary conditions are numerically solved by utilizing the bvp4c solver in MATLAB programming language.The effects of the governing parameters on the non-dimensional velocity profiles,temperature profiles,skin friction coefficients,and the local Nusselt number are presented in detail through a series of graphs and tables.Interestingly,it is reported that the reduced skin friction coefficient decreases for the assisting flow situation and increases for the opposing flow situation.The numerical computations of the present work are compared with those from other research available in specific situations,and an excellent consensus is observed.Another exciting feature for this work is the existence of dual solutions.An important remark is that the dual solutions exist for both assisting and opposing flows.A linear stability analysis is performed showing that one solution is stable and the other solution is not stable.We notice that the mixed convection and velocity slip parameters have strong effects on the flow characteristics.These effects are depicted in graphs and discussed in this paper.The obtained results show that the first-order and second-order slip parameters have a considerable effect on the flow,as well as on the heat transfer characteristics.

Elastic and Inelastic Response of Structural Systems in Seismic Pounding

The present paper addresses the comparative study of three adjacent single-degree-of freedom structures for elastic and inelastic system with and without pounding under seismic excitations. For the gap between three adjacent structures, the simulation is done by using linear spring element without damping. The entire numerical simulation is done in time domain by considering the inputs of four real ground motions. The results of the study show that the response of elastic system is much different to that of response of inelastic system in the absence and presence of pounding, especially in lighter or more flexible structures. Elastic structures show much severe pounding response than inelastic structures. Modeling of colliding structures behaving inelastically is really needed in order to obtain the accurate structural pounding involved response under seismic excitation.

Double stabilizations and convergence analysis of a second-order linear numerical scheme for the nonlocal Cahn-Hilliard equation

In this paper,we study a second-order accurate and linear numerical scheme for the nonlocal CahnHilliard equation.The scheme is established by combining a modified Crank-Nicolson approximation and the Adams-Bashforth extrapolation for the temporal discretization,and by applying the Fourier spectral collocation to the spatial discretization.In addition,two stabilization terms in different forms are added for the sake of the numerical stability.We conduct a complete convergence analysis by using the higher-order consistency estimate for the numerical scheme,combined with the rough error estimate and the refined estimate.By regarding the numerical solution as a small perturbation of the exact solution,we are able to justify the discrete?^(∞)bound of the numerical solution,as a result of the rough error estimate.Subsequently,the refined error estimate is derived to obtain the optimal rate of convergence,following the established?∞bound of the numerical solution.Moreover,the energy stability is also rigorously proved with respect to a modified energy.The proposed scheme can be viewed as the generalization of the second-order scheme presented in an earlier work,and the energy stability estimate has greatly improved the corresponding result therein.