Trigonometric Regularization and Continuation Method Based Time-Optimal Control of Hypersonic Vehicles

Aiming at the time-optimal control problem of hypersonic vehicles(HSV)in ascending stage,a trigonometric regularization method(TRM)is introduced based on the indirect method of optimal control.This method avoids analyzing the switching function and distinguishing between singular control and bang-bang control,where the singular control problem is more complicated.While in bang-bang control,the costate variables are unsmooth due to the control jumping,resulting in difficulty in solving the two-point boundary value problem(TPBVP)induced by the indirect method.Aiming at the easy divergence when solving the TPBVP,the continuation method is introduced.This method uses the solution of the simplified problem as the initial value of the iteration.Then through solving a series of TPBVP,it approximates to the solution of the original complex problem.The calculation results show that through the above two methods,the time-optimal control problem of HSV in ascending stage under the complex model can be solved conveniently.

Non-cooperative Space Target Estimation Algorithm Without Prior Information Dependence Based on Temporal Line of Sight Constraint

Under single-satellite observation,the parameter estimation of the boost phase of high-precision space noncooperative targets requires prior information.To improve the accuracy without prior information,we propose a parameter estimation model of the boost phase based on trajectory plane parametric cutting.The use of the plane passing through the geo-center and the cutting sequence line of sight(LOS)generates the trajectory-cutting plane.With the coefficient of the trajectory cutting plane directly used as the parameter to be estimated,a motion parameter estimation model in space non-cooperative targets is established,and the Gauss-Newton iteration method is used to solve the flight parameters.The experimental results show that the estimation algorithm proposed in this paper weakly relies on prior information and has higher estimation accuracy,providing a practical new idea and method for the parameter estimation of space non-cooperative targets under single-satellite warning.

Relative orbit determination algorithm of space targets with passive observation

Angles-only relative orbit determination for space non-cooperative targets based on passive sensor is subject to weakly observable problem of the relative state between two spacecraft. Previously, the evidence for angles-only observability was found by using cylindrical dynamics, however, the solution of orbit determination is still not provided. This study develops a relative orbit determination algorithm with the cylindrical dynamics based on differential evolution. Firstly, the relative motion dynamics and line-of-sight measurement model for nearcircular orbit are established in cylindrical coordinate system.Secondly, the observability is qualitatively analyzed by using the dynamics and measurement model where the unobservable geometry is found. Then, the angles-only relative orbit determination problem is modeled into an optimal searching frame and an improved differential evolution algorithm is introduced to solve the problem. Finally, the proposed algorithm is verified and tested by a set of numerical simulations in the context of highEarth and low-Earth cases. The results show that initial relative orbit determination(IROD) solution with an appropriate accuracy in a relative short span is achieved, which can be used to initialize the navigation filter.

Experimental and numerical analyses of electrohydrodynamic force according to air pressure

Electrohydrodynamic(EHD) force produced by corona discharge is considered as a new thrust for solar-powered aircraft and stratosphere balloons in near space. However, its performance at low air pressures remains to be clarified.An experiment of measuring the EHD force at 0.02 atm–1.0 atm(1 atm = 1.01325×10^(5)Pa) is carried out with the wireto-cylinder geometric structure. The ion distribution is analyzed by using the drift-diffusion model with two-dimensional numerical simulation. The experimental result shows that the EHD force is not linearly related to the corona discharge current at low air pressures. Numerical simulation finds that the proportion of ions in the counter-direction electric field increases from approximately 0.36% to 30% when the pressure drops from 1.0 atm to 0.2 atm. As a result, the EHD force with a constant power supply drops faster than the previous theoretical prediction in the ground experiment, suggesting that the consideration of counter-direction EHD force is necessary to improve the ionic wind propulsion efficiency in near-space applications.

Vegetation scattering attenuation characteristics of terahertz wave

A terahertz(THz)wave transmitted through vegetation experiences both absorption and scattering caused by the air molecules and leaves.This paper presents the scattering attenuation characteristics of vegetation in a THz range.The theoretical path loss model near the vegetation yields the average attenuation of THz waves in a mixed channel composed of air and vegetation leaves.Furthermore,a simplified model of the vegetation structure is obtained for generic vegetation types based on a variety of parameters,such as leaf size,distribution,and moisture content.Finally,based on specific vegetation species and different levels of air humidity,the attenuation characteristics under different conditions are calculated,and the influence of different model parameters on the attenuation characteristics is obtained.

Distributed Dynamic Load in Structural Dynamics by the Impulse-Based Force Estimation Algorithm

This paper proposes a novel approach for identifying distributed dynamic loads in the time domain.Using polynomial andmodal analysis,the load is transformed intomodal space for coefficient identification.This allows the distributed dynamic load with a two-dimensional form in terms of time and space to be simultaneously identified in the form of modal force,thereby achieving dimensionality reduction.The Impulse-based Force Estimation Algorithm is proposed to identify dynamic loads in the time domain.Firstly,the algorithm establishes a recursion scheme based on convolution integral,enabling it to identify loads with a long history and rapidly changing forms over time.Secondly,the algorithm introduces moving mean and polynomial fitting to detrend,enhancing its applicability in load estimation.The aforementioned methodology successfully accomplishes the reconstruction of distributed,instead of centralized,dynamic loads on the continuum in the time domain by utilizing acceleration response.To validate the effectiveness of the method,computational and experimental verification were conducted.

Modified filter for mean elements estimation with state jumping

To investigate the real-time mean orbital elements(MOEs)estimation problem under the influence of state jumping caused by non-fatal spacecraft collision or protective orbit trans-fer,a modified augmented square-root unscented Kalman filter(MASUKF)is proposed.The MASUKF is composed of sigma points calculation,time update,modified state jumping detec-tion,and measurement update.Compared with the filters used in the existing literature on MOEs estimation,it has three main characteristics.Firstly,the state vector is augmented from six to nine by the added thrust acceleration terms,which makes the fil-ter additionally give the state-jumping-thrust-acceleration esti-mation.Secondly,the normalized innovation is used for state jumping detection to set detection threshold concisely and make the filter detect various state jumping with low latency.Thirdly,when sate jumping is detected,the covariance matrix inflation will be done,and then an extra time update process will be con-ducted at this time instance before measurement update.In this way,the relatively large estimation error at the detection moment can significantly decrease.Finally,typical simulations are per-formed to illustrated the effectiveness of the method.

Temperature and Salinity Dual-parameter Sensing Based on Forward Brillouin Scattering in 1060-XP SMF

A novel temperature and salinity discriminative sensing method based on forward Brillouin scattering(FBS)in 1060-XP single-mode fiber(SMF)is proposed.The measured frequency shifts corresponding to different radial acoustic modes in 1060-XP SMF show different sensitivities to temperature and salinity.Based on the new phenomenon that different radial acoustic modes have different frequency shift-temperature and frequency shift-salinity coefficients,we propose a novel method for simultaneously measuring temperature and salinity by measuring the frequency shift changes of two FBS scattering peaks.In a proof-of-concept experiment,the temperature and salinity measurement errors are 0.12℃and 0.29%,respectively.The proposed method for simultaneously measuring temperature and salinity has the potential applications such as ocean surveying,food manufacturing and pharmaceutical engineering.

A Computational Framework for Parachute Inflation Based on Immersed Boundary/Finite Element Approach

A computational framework for parachute inflation is developed based on the immersed boundary/finite element approach within the open-source IBAMR library.The fluid motion is solved by Peskin's diffuse-interface immersed boundary(IB)method,which is attractive for simulating moving-boundary flows with large deformations.The adaptive mesh refinement technique is employed to reduce the computational cost while retain the desired resolution.The dynamic response of the parachute is solved with the finite element approach.The canopy and cables of the parachute system are modeled with the hyperelastic material.A tether force is introduced to impose rigidity constraints for the parachute system.The accuracy and reliability of the present framework is validated by simulating inflation of a constrained square plate.Application of the present framework on several canonical cases further demonstrates its versatility for simulation of parachute inflation.

Tradeoff Analysis of Factors Affecting Longitudinal Carrier Landing Performance for Small UAV Based on Backstepping Controller

Tradeoff analysis of the factors,including external environment and unmanned aerial vehicle(UAV)aerodynamic attributes,which affect longitudinal carrier landing performance,is important for small UAV.First,small UAV longitudinal carrier landing system is established,as well as the nonlinear dynamics and kinematics model,and then the longitudinal flight control system using backstepping technology with minimum information about the aerodynamic is designed.To assess the landing performance,a variety of influencing factors are considered,resulting in the constraints of aerodynamic attributes of carrier UAV.The simulation results show that the severe sea condition has the greatest influence on landing dispersion,while air wake is the primary factor on impact velocity.Among the longitudinal aerodynamic parameters,the lift curve slope is the most important factor affecting the landing performance,and increasing lift curve slope can improve the landing performance significantly.A better system performance will be achieved when the lift curve slope is larger than 2per radian.

Preparation and thermal conductivity of tungsten coated diamond/copper composites

Tungsten was plated on the surface of diamond by using thermal diffusion method.Different process parameters were employed to prepare the composites with tungsten,diamond and copper.The micro morphology of different samples was observed,and the thermal conductivity of samples was measured by laser flash method.The optimal process parameters for preparing diamond/copper composites with high thermal conductivity were investigated.The results indicated that plating tungsten on diamond could modify the interface bonding.When the diamond was plated for 60 min,the coating appeared intact,uniform and flat,and the thermal conductivity of the sample could reach as high as 486 W/(m·K).The integrity and uniformity were more important than thickness for the coating.When the tungsten-plated diamond was further annealed,the metallurgical bonding between the coating and the diamond was enhanced,and the thermal conductivity rose to 559 W/(m·K).

Robust adaptive control of hypersonic vehicle considering inlet unstart

In this paper,a model reference adaptive control(MRAC)augmentation method of a linear controller is proposed for air-breathing hypersonic vehicle(AHV)during inlet unstart.With the development of hypersonic flight technology,hypersonic vehicles have been gradually moving to the stage of weaponization.During the maneuvers,changes of attitude,Mach number and the back pressure can cause the inlet unstart phenomenon of scramjet.Inlet unstart causes significant changes in the aerodynamics of AHV,which may lead to deterioration of the tracking performance or instability of the control system.Therefore,we firstly establish the model of hypersonic vehicle considering inlet unstart,in which the changes of aerodynamics caused by inlet unstart is described as nonlinear uncertainty.Then,an MRAC augmentation method of a linear controller is proposed and the radial basis function(RBF)neural network is used to schedule the adaptive parameters of MRAC.Furthermore,the Lyapunov function is constructed to prove the stability of the proposed method.Finally,numerical simulations show that compared with the linear control method,the proposed method can stabilize the attitude of the hypersonic vehicle more quickly after the inlet unstart,which provides favorable conditions for inlet restart,thus verifying the effectiveness of the augmentation method proposed in the paper.

Simulation Analysis and Experimental Verification of Spiral-tube-type Valveless Piezoelectric Pump with Gyroscopic Effect

The current research of the valveless piezoelectric pump focuses on increasing the flow rate and pressure differential. Compared with the valve piezoelectric pump, the valveless one has excellent performances in simple structure, low cost, and easy miniaturization. So, their important development trend is the mitigation of their weakness, and the multi-function integration. The flow in a spiral tube element is sensitive to the element attitude caused by the Coriolis force, and that a valveless piezoelectric pump is designed by applying this phenomenon. The pump has gyroscopic effect, and has both the actuator function of fluid transfer and the sensor function, which can obtain the angular velocity when its attitude changes. First, the present paper analyzes the flow characteristics in the tube, obtains the calculation formula for the pump flow, and identifies the relationship between pump attitude and flow, which clarifies the impact of flow and driving voltage, frequency, spiral line type and element attitude, and verifies the gyroscopic effect of the pump. Then, the finite element simulation is used to verify the theory. Finally, a pump is fabricated for experimental testing of the relationship between pump attitude and pressure differential. Experimental results show that when Archimedes spiral θ=4π is selected for the tube design, and the rotation speed of the plate is 70 r/min, the pressure differential is 88.2 Pa, which is 1.5 times that of 0 r/min rotation speed. The spiral-tube-type valveless piezoelectric pump proposed can turn the element attitude into a form of pressure output, which is important for the multi-function integration of the valveless piezoelectric pump and for the development of civil gyroscope in the future.

Pseudo-noise preamble based joint frame and frequency synchronization algorithm in OFDM communication systems

Frame and frequency synchronization are essential for orthogonal frequency division multiplexing （OFDM） systems. The frame offset owing to incorrect start point position of the fast Fourier transform （FFT） window, and the carrier frequency offset （CFO） due to Doppler frequency shift or the frequency mismatch between the transmitter and receiver oscil ators, can bring severe inter-symbol interference （ISI） and inter-carrier interference （ICI） for the OFDM system. Relying on the relatively good correlation charac-teristic of the pseudo-noise （PN） sequence, a joint frame offset and normalized CFO estimation algorithm based on PN preamble in time domain is developed to realize the frame and frequency synchronization in the OFDM system. By comparison, the perfor-mances of the traditional algorithm and the improved algorithm are simulated under different conditions. The results indicate that the PN preamble based algorithm both in frame offset estimation and CFO estimation is more accurate, resource-saving and robust even under poor channel condition, such as low signal-to-noise ratio （SNR） and large normalized CFO.

Determination of thickness of wax deposition in oil pipelines using gamma-ray transmission method

Wax deposition in pipelines is a crucial problem in the oil industry.An approach that combines the gammaray transmission method with scanning technology is proposed to detect the thickness of wax deposition.The performance of the method is validated through simulations with MCNP code.An experiment is also carried out with a 300 mCi ^(137)Cs source and a LaBr_3 detector.A good correspondence is observed between the simulation and experimental results.The results indicate that the approach is efficient for detecting the thickness of wax deposition in oil pipelines.

ON RESIDUAL COMPRESSIVE STRENGTH PREDICTION OF COMPOSITE SANDWICH PANELS AFTER LOW-VELOCITY IMPACT DAMAGE

This paper introduces a nonlinear finite element analysis on damage propagation behavior of composite sandwich panels under in-plane uniaxial quasi-static compression after a low velocity impact. The major damage modes due to the impact, including the residual indentation on the impacted facesheet, the initially crushed core under the impacted area, and the delamination are incorporated into the model. A consequential core crushing mechanism is incorporated into the analysis by using an element deactivation technique. Damage propagation behavior, which corresponds to those observed in sandwich compression after impact （SCAI） tests, has been successfully captured in the numerical simulation. The critical far field stress corresponding to the onset of damage propagation at specified critical locations near the damage zone are captured successfully. They show a good correlation with experimental data. These values can be used to effectively predict the residual compressive strength of low-velocity impact damaged composite sandwich panels.

Diagnosing the Fine Structure of Electron Energy Within the ECRIT Ion Source

The ion source of the electron cyclotron resonance ion thruster （ECRIT） extracts ions from its ECR plasma to generate thrust, and has the property of low gas consumption （2 seem, standard-state cubic centimeter per minute） and high durability. Due to the indispensable effects of the primary electron in gas discharge, it is important to experimentally clarify the electron energy structure within the ion source of the ECRIT through analyzing the electron energy distribution function （EEDF） of the plasma inside the thruster. In this article the Langmuir probe diagnosing method was used to diagnose the EEDF, from which the effective electron temperature, plasma density and the electron energy probability function （EEPF） were deduced. The experimental results show that the magnetic field influences the curves of EEDF and EEPF and make the effective plasma parameter nonuniform. The diagnosed electron temperature and density from sample points increased from 4 eV/2 ×10^16 m-3 to 10 eV/4×10^16 m-3 with increasing distances from both the axis and the screen grid of the ion source. Electron temperature and density peaking near the wall coincided with the discharge process. However, a double Maxwellian electron distribution was unexpectedly observed at the position near the axis of the ion source and about 30 mm from the screen grid. Besides, the double Maxwellian electron distribution was more likely to emerge at high power and a low gas flow rate. These phenomena were believed to relate to the arrangements of the gas inlets and the magnetic field where the double Maxwellian electron distribution exits. The results of this research may enhance the understanding of the plasma generation process in the ion source of this type and help to improve its performance.

A novel weighed cooperative bandwidth spectrum sensing for spectrum occupancy of cognitive radio network

In cognitive radio network(CRN), a secondary user(SU) may utilize the spectrum resource of the primary user(PU) and avoid causing harmful interference to the primary network(PN) via spectrum sensing. In the traditional time spectrum sensing, the SU cannot detect the PU's presence during its transmission, thus increasing interference to the PN. In this work, a novel weighed cooperative bandwidth spectrum sensing method is proposed, which allows multiple SUs to use part of the bandwidth to perform cooperative spectrum sensing throughout the whole frame in order to detect the PU's reappearance in time. The SU's spectrum efficiency is maximized by jointly optimizing sensing bandwidth proportion, number of cooperative SUs and detection probability, subject to the constraints on the SU's interference and the false alarm probability. Simulation results show significant decrease on the interference and improvement on the spectrum efficiency using the proposed weighed cooperative bandwidth spectrum sensing method.

Time dependence of entropy flux and entropy production for a dynamical system driven by noises with coloured cross-correlation

This paper shows the Fokker-Planck equation of a dynamical system driven by coloured cross-correlated white noises in the absence and presence of a small external force. Based on the Fokker-Planck equation and the definition of Shannon＇s information entropy, the time dependence of entropy flux and entropy production can be calculated. The present results can be used to explain the extremal behaviour of time dependence of entropy flux and entropy production in view of the dissipative parameter γ of the system, coloured cross-correlation time τ and coloured cross-correlation strength λ.

Influence of Refractive Index Distribution on Brillouin Gain Spectrum in GeO_(2)-Doped Optical Fibers

GeO_(2) is commonly used as dopant to adjust the refractive index profile(RIP)and the acoustic velocity profile(AVP)in the fiber,thereby forming different Brillouin gain spectrum(BGS)characteristics such as Brillouin gain,acoustic mode number and peak intensity difference.When an optical fiber is used in optical fiber sensing or communication system,its BGS characteristics may play an important role in determining the performance of the system.In this paper,finite element analysis(FEA)method is used to study the influence of refractive index distribution and its corresponding AVP on the BGS in step-index,graded-index,and complex-index optical fibers.A new method has also been proposed to efficiently discriminate acoustic mode solution and obtain the new and full images of total Brillouin gain and acoustic modes number of the fiber as a function of the refractive index distribution,considering the influence of changing the refractive index difference and the geometric size simultaneously.For each type of optical fiber,the recommended parameter range is provided for optical fiber sensing and optical fiber communication.Moreover,the suitable optical fiber with close peak intensity in its multi-peak BGS is explored and achieved,which can be used in Brillouin beat spectrum detection systems to improve sensing accuracy.