Design and Cryogenic Performance of a Hexapod Platform for a Large Ground-based Wide Field Survey Telescope

The thermal gradient is an important factor that causes degradation to the image quality of telescopes. In order to ensure the accurate alignment of the primary focus unit and the primary mirror, the hexapod platform(as a corrector) is investigated in this paper. First, a ground-based telescope with 2.5 m aperture and 3.5 deg field of view is described. The telescope is under construction, and it is expected to be finished in 2023. Second, the hexapod platform with flexure hinges utilized to adjust the primary focus unit is proposed, which is applied as a corrector.Then, the inverse kinematics of the platform is established and an open-loop control system is built based on it.Finally, the cryogenic performance test for the hexapod platform is performed. The experimental results show that the resolution and repeatability of the translation for the hexapod platform can be achieved at the micrometer level.The resolution and repeatability of the rotation can be achieved at the arc-second level. Therefore, the cryogenic performance of the hexapod platform can meet the optical imaging requirements of the wide-field ground-based telescope. The kinematic analysis and cryogenic performance tests in the paper provide a technical reference for the precise alignment of the primary focus unit and the primary mirror, which can improve the imaging quality of the telescope.

Disturbance Observer-based Pointing Control of Leighton Chajnantor Telescope

Leighton Chajnantor Telescope(LCT), i.e., the former Caltech Submillimeter Observatory telescope, will be refurbished at the new site in Chajnantor Plateau, Chile in 2023. The environment of LCT will change significantly after its relocation, and the telescope will be exposed to large wind disturbances directly because its enclosure will be completely open during observation. The wind disturbance is expected to be a challenge for LCT's pointing control since the existing control method cannot reject this disturbance very well. Therefore, it is very necessary to develop a new pointing control method with good capability of disturbance rejection. In this research, a disturbance observer—based composite position controller(DOB-CPC) is designed, in which an H∞feedback controller is employed to compress the disturbance, and a feedforward linear quadratic regulator is employed to compensate the disturbance precisely based on the estimated disturbance signal. Moreover, a controller switching policy is adopted, which applies the proportional controller to the transient process to achieve a quick response and applies the DOB-CPC to the steady state to achieve a small position error. Numerical experiments are conducted to verify the good performance of the proposed pointing controller(i.e., DOB-CPC) for rejecting the disturbance acting on LCT.

A class of fully third-order accurate projection methods for solving the incompressible Navier-Stokes equations

In this paper, a fully third-order accurate projection method for solving the incompressible Navier-Stokes equations is proposed. To construct the scheme, a continuous projection procedure is firstly presented. We then derive a sufficient condition for the continuous projection equations to be temporally third-order accurate approximations of the original Navier-Stokes equations by means of the localtruncation-error-analysis technique. The continuous projection equations are discretized temporally and spatially to third-order accuracy on the staggered grids, resulting in a fully third-order discrete projection scheme. The possibility to design higher-order projection methods is thus demonstrated in the present paper. A heuristic stability analysis is performed on this projection method showing the probability of its being stable. The stability of the present scheme is further verified through numerical tests. The third-order accuracy of the present projection method is validated by several numerical test cases.

PSRDP:A Parallel Processing Method for Pulsar Baseband Data

To address the problem of real-time processing of ultra-wide bandwidth pulsar baseband data,we designed and implemented a pulsar baseband data processing algorithm(PSRDP)based on GPU parallel computing technology.PSRDP can perform operations such as baseband data unpacking,channel separation,coherent dedispersion,Stokes detection,phase and folding period prediction,and folding integration in GPU clusters.We tested the algorithm using the J0437-4715 pulsar baseband data generated by the CASPSR and Medusa backends of the Parkes,and the J0332+5434 pulsar baseband data generated by the self-developed backend of the Nan Shan Radio Telescope.We obtained the pulse profiles of each baseband data.Through experimental analysis,we have found that the pulse profiles generated by the PSRDP algorithm in this paper are essentially consistent with the processing results of Digital Signal Processing Software for Pulsar Astronomy(DSPSR),which verified the effectiveness of the PSRDP algorithm.Furthermore,using the same baseband data,we compared the processing speed of PSRDP with DSPSR,and the results showed that PSRDP was not slower than DSPSR in terms of speed.The theoretical and technical experience gained from the PSRDP algorithm research in this article lays a technical foundation for the real-time processing of QTT(Qi Tai radio Telescope)ultra-wide bandwidth pulsar baseband data.

Jacobi-Sobolev Orthogonal Polynomialsand Spectral Methods for Elliptic Boundary Value Problems

Generalized Jacobi polynomials with indexes α,β∈ R are introduced and some basic properties are established. As examples of applications,the second- and fourth-order elliptic boundary value problems with Dirichlet or Robin boundary conditions are considered,and the generalized Jacobi spectral schemes are proposed. For the diagonalization of discrete systems,the Jacobi-Sobolev orthogonal basis functions are constructed,which allow the exact solutions and the approximate solutions to be represented in the forms of infinite and truncated Jacobi series. Error estimates are obtained and numerical results are provided to illustrate the effectiveness and the spectral accuracy.

Edge sensor based real-time calculation method for active panel position of QTT

QiTai Telescope(QTT),the world’s biggest full-steerable telescope,will be constructed in Xinjiang,China.Its extra high operating frequency(115 GHz)imposes strict requirements on the accuracy of the reflector,while its large aperture(110 m)increases the impact of antenna weight and environment on surface accuracy.However,the panels of reflector will deform under the influence of gravity and the environment.Therefore,to compensate for the performance degradation caused by this deformation,a technique called active surface adjustment has been proposed.The existing adjustment methods cannot detect the deformation caused by environment of the reflector surface in real time,which can result in a delay in the compensation.Consequently,it is difficult to achieve an optimal compensating result.To solve this problem,a real-time method to estimate the large-aperture reflector antenna surface by calculating the antenna panel position based on edge sensors is proposed in this paper.In the proposed method,a panel coordinate transfer matrix has been formulated and the data measured from these edge sensors can then be treated as input for the proposed transfer matrix to calculate the actual position of the antenna panel in real time.A numerical simulation has been carried out in the QTT model and the results obtained show that the proposed real-time method is a promising tool to estimate the large reflector surface position in real time.It is believed that this active surface adjustment method has laid the foundation for new methods that will be developed to compensate for the reflector electrical performance.

A hyperbolic Lindstedt-Poincare method for homoclinic motion of a kind of strongly nonlinear autonomous oscillators

A hyperbolic Lindstedt-Poincare method is presented to determine the homoclinic solutions of a kind of nonlinear oscillators, in which critical value of the homoclinic bifurcation parameter can be determined. The generalized Lienard oscillator is studied in detail, and the present method＇s predictions are compared with those of Runge-Kutta method to illustrate its accuracy.

Simulation method of Wolter type Ⅰ grazing incidence telescope imaging of an X-ray region

This paper presents a method for simulating Wolter type Ⅰ grazing incidence telescope imaging of an X-ray region. The imaging quality of a soft X-ray telescope is mainly affected by geometric aberrations,surface scatter and alignment errors. Using a Legendre-Fourier polynomial to fit the figure errors in the cylindrical coordinate system, the geometry of the mirror barrel with the figure errors can be modeled by the User Defined Object in ZEMAX. Based on the Harvey-Shack surface scatter theory, surface scatter is achieved by the Bidirectional Scatter Distribution Function(BSDF) scattering model. The alignment errors are calculated by a combination of experiment and simulation according to the shape of the spot resulting from the finite distance imaging experiment. By combining these factors that affect image quality,the telescope imaging of an infinity point is simulated by ZEMAX. According to the simulation results, the angular resolution of the telescope is calculated to be less than or equal to 1.86′ at a wavelength of 1.33 nm.The simulation method improves the ability of ZEMAX to simulate grazing incidence telescope imaging.This is of great significance for the research and manufacture of X-ray grazing incidence telescopes.

Optical fiber positioning based on four-quadrant detector with Gaussian fitting method

With the development of large-scale spectral surveys, fiber positioning technology has been developing rapidly. Because of the performance advantages of a four-quadrant（4Q） detector, a fiber positioning and real-time monitoring system based on the 4Q detector is proposed. The detection accuracy of this system is directly determined by the precision of the center of the spot. A Gaussian fitting algorithm based on the 4Q detector is studied and applied in the fiber positioning process to improve the calculated accuracy of the spot center. The relationship between the center position of the incident spot and the detector output signal is deduced. An experimental platform is built to complete the simulated experiment. Then we use the Gaussian fitting method to process experimental data, compare the fitting value with the theoretical one and calculate the corresponding error.

Surface and Texture Properties of Tb-Doped Ceria-Zirconia Solid Solution Prepared by Sol-Gel Method

The three-way catalysts （TWCs） promoters Ce0.6Zr0.4- x TbxO2-y were prepared by sol-gel method. BET surface areas analysis indicated that an increase of the dopant Tb content from x = 0.05 to x = 0.15 favors an increase of surface area from 66.8 to 80.4 m^2· g^-1 compared with the undoped sample Ce0 .6oZr0.40O2 65.1 m^2·g^- 1 after calcination at 650℃. Transmission electron microscopy （TEM） observation indicated that the doped samples have a higher thermal stability. The XRD and Raman spectra confirmed that the Ce0.6Zr0.4-xTbxO2-y cubic solid solution is formed. XPS analysis revealed that Ce and Tb mainly existed in the form of Ce^4＋ and Tb^3 ＋ , and Zr existed in the form of Zr^4＋ on the surface of the samples. The doped samples were homogenous in composition ; the introduction of Tb into the CeO2-ZrO2 promoters resuited in the formation of a solid solution, and the concentration of surface lattice oxygen was increased.

Construction of n-sided polygonal spline element using area coordinates and B-net method

In general, triangular and quadrilateral elements are commonly applied in two-dimensional finite element methods. If they are used to compute polycrystalline materials, the cost of computation can be quite significant. Polygonal elements can do well in simulation of the materials behavior and provide greater flexibility for the meshing of complex geometries. Hence, the study on the polygonal element is a very useful and necessary part in the finite element method. In this paper, an n-sided polygonal element based on quadratic spline interpolant, denoted by PS2 element, is presented using the triangular area coordinates and the B-net method. The PS2 element is conforming and can exactly model the quadratic field. It is valid for both convex and non-convex polygonal element, and insensitive to mesh distortions. In addition, no mapping or coordinate transformation is required and thus no Jacobian matrix and its inverse are evaluated. Some appropriate examples are employed to evaluate the performance of the proposed element.

Multi-train modeling and simulation integrated with traction power supply solver using simplified Newton–Raphson method

Multi-train modeling and simulation plays a vital role in railway electrification during operation and planning phase. Study of peak power demand and energy consumed by each traction substation needs to be deter- mined to verify that electrical energy flowing in its railway power feeding system is appropriate or not. Gauss-Seidel, conventional Newton-Raphson, and current injection methods are well-known and widely accepted as a tool for electrical power network solver in DC railway power supply study. In this paper, a simplified Newton-Raphson method has been proposed. The proposed method employs a set of current-balance equations at each electrical node instead of the conventional power-balance equation used in the conventional Newton-Raphson method. This concept can remarkably reduce execution time and computing complexity for multi-train simulation. To evaluate its use, Sukhumvit line of Bangkok transit system （BTS） of Thai- land with 21.6-km line length and 22 passenger stopping stations is set as a test system. The multi-train simulation integrated with the proposed power network solver is developed to simulate 1-h operation service of selected 5-min headway. From the obtained results, the proposed method is more efficient with approximately 18 % faster than the conventional Newton-Raphson method and just over 6 % faster than the current injection method.

A hybrid vertex-centered finite volume/element method for viscous incompressible flows on non-staggered unstructured meshes

This paper proposes a hybrid vertex-centered fi- nite volume/finite element method for solution of the two di- mensional （2D） incompressible Navier-Stokes equations on unstructured grids. An incremental pressure fractional step method is adopted to handle the velocity-pressure coupling. The velocity and the pressure are collocated at the node of the vertex-centered control volume which is formed by join- ing the centroid of cells sharing the common vertex. For the temporal integration of the momentum equations, an im- plicit second-order scheme is utilized to enhance the com- putational stability and eliminate the time step limit due to the diffusion term. The momentum equations are discretized by the vertex-centered finite volume method （FVM） and the pressure Poisson equation is solved by the Galerkin finite el- ement method （FEM）. The momentum interpolation is used to damp out the spurious pressure wiggles. The test case with analytical solutions demonstrates second-order accuracy of the current hybrid scheme in time and space for both veloc- ity and pressure. The classic test cases, the lid-driven cavity flow, the skew cavity flow and the backward-facing step flow, show that numerical results are in good agreement with the published benchmark solutions.

A meshfree-based local Galerkin method with condensation of degree of freedom for elastic dynamic analysis

Condensation technique of degree of freedom is first proposed to improve the computational efficiency of meshfree method with Galerkin weak form for elastic dy- namic analysis. In the present method, scattered nodes with- out connectivity are divided into several subsets by cells with arbitrary shape. Local discrete equation is established over each cell by using moving Kriging interpolation, in which the nodes that located in the cell are used for approxima- tion. Then local discrete equations can be simplified by con- densation of degree of freedom, which transfers equations of inner nodes to equations of boundary nodes based on cells. The global dynamic system equations are obtained by as- sembling all local discrete equations and are solved by using the standard implicit Newmark＇s time integration scheme. In the scheme of present method, the calculation of each cell is carried out by meshfree method, and local search is imple- mented in interpolation. Numerical examples show that the present method has high computational efficiency and good accuracy in solving elastic dynamic problems.

Short-term solar flare prediction using multi-model integration method

A multi-model integration method is proposed to develop a multi-source and heterogeneous model for short-term solar flare prediction. Different prediction models are constructed on the basis of extracted predictors from a pool of observation databases. The outputs of the base models are normal- ized first because these established models extract predictors from many data resources using different prediction methods. Then weighted integration of the base models is used to develop a multi-model integrated model （MIM）. The weight set that single models assign is optimized by a genetic algorithm. Seven base models and data from Solar and Heliospheric Observatory/Michelson Doppler Imager lon- gitudinal magnetograms are used to construct the MIM, and then its performance is evaluated by cross validation. Experimental results showed that the MIM outperforms any individual model in nearly every data group, and the richer the diversity of the base models, the better the performance of the MIM. Thus, integrating more diversified models, such as an expert system, a statistical model and a physical model, will greatly improve the performance of the MIM.

Pseudo-beam method for compressive buckling characteristics analysis of space inflatable load-carrying structures

This paper extends Le van＇s work to the case of nonlinear problem and the complicated configuration. The wrinkling stress distribution and the pressure effects are also included in our analysis. Pseudo-beam method is presented based on the inflatable beam theory to model the inflatable structures as a set of inflatable beam elements with a prestressed state. In this method, the discretized nonlinear equations are given based upon the virtual work principle with a 3-node Timoshenko＇s beam model. Finite element simulation is performed by using a 3-node BEAM189 element incorporating ANSYS nonlinear program. The pressure effect is equivalent included in our method by modifying beam element cross-section parameters related to pressure. A benchmark example, the bending case of an inflatable cantilever beam, is performed to verify the accuracy of our proposed method. The comparisons reveal that the numerical results obtained with our method are close to open published analytical and membrane finite element results. The method is then used to evaluate the whole buckling and the loadcarrying characteristics of an inflatable support frame subjected to a compression force. The wrinkling stress and region characteristics are also shown in the end. This method gives better convergence characteristics, and requires much less computation time. It is very effective to deal with the whole load-carrying ability analytical problems for large scale inflatable structures with complex configuration.

Two 8-node quadrilateral spline elements by B-net method

Isoparametric quadrilateral elements are widely used in the finite element method, but the accuracy of the isoparametric quadrilateral elements will drop obviously deteriorate due to mesh distortions. Spline functions have some properties of simplicity and conformality. Two 8-node quadrilateral elements have been developed using the trian- gular area coordinates and the B-net method, which can ex- actly model the quadratic field for both convex and concave quadrangles. Some appropriate examples are employed to evaluate the performance of the proposed elements. The nu- merical results show that the two spline elements can obtain solutions which are highly accurate and insensitive to mesh distortions.

Comparison of two earthquake early warning location methods

According to earthquake catalog records of Fujian Seismic Network, the Tnow method and the fourstation continuous location method put forward by Jin Xing are inspected by using P-wave arrival information of the first four stations in each earthquake. It shows that the fourstation continuous location method can locate more seismic events than the Tnow method. By analyzing the results, it is concluded that the reason for this is that the Tnow method makes use of information from stations without being triggered, while some stations failed to be reflected in earthquake catalog because of discontinuous records or unclear records of seismic phases. For seismic events whose location results can be given, there is no obvious difference in location results of the two methods and positioning deviation of most seismic events is also not significant. For earthquakes outside the network, the positioning deviation may amplify as the epicentral distance enlarges, which may relate to the situation that the seismic stations are centered on one side of epicenter and the opening angle between seismic stations used for location and epicenter is small.

Preparation of BaCe_(0.5)Zr_(0.4)Y_(0.1)O_(3-α) by Sol-Gel Method and Its Electrical Properties

A precursor of BaCe0.5Zr0.4Y0. 1O3-α electrolytes was synthesized by the sol-gel method and sintered at temperature which were 150 - 250 ℃ lower than by solid state reaction. The AC impedance spectrums of electrolytes were measured by AUTOLAB PGSTA30 electrochemical measuring device at different temperatures. The conductivities of the electrolytes are 1.62×10^-4 - 6.43×10^-3, 2.52×10^-5 - 3.73×10^-3S·cm^-1 in the temperature range of 350-800℃. The activity energies are 0.54 and 0.84 eV. At the same time BaCe0.9Y0.1O3-α was prepared by direct solid state reaction. The conductivity of BaCe0.9Y0.1O3-α is 1 × 10^-4- 4×10^-3 S·cm^-1 and the activation energy is 0.50 eV at the same condition. The results show that conduction of electrolyte prepared by sol-gel method is higher than the one by solid state reaction. As far as BaCe0.9Y0.1O3-α concerned, its conductivity of the Zr-substituted specimens is decreased.

Simulation Study of Network Reconfiguration and Load-balancing Method for the Xinjiang Astronomical Observatory Data Center

The Xinjiang Astronomical Observatory Data Center faces issues related to delay-affected services. As a result, these services cannot be implemented in a timely manner due to the overloading of transmission links. In this paper, the software-defined network technology is applied to the Xinjiang Astronomical Observatory Data Center Network(XAODCN). Specifically, a novel reconfiguration method is proposed to realise the software-defined Xinjiang Astronomical Observatory Data Center Network(SDXAO-DCN), and a network model is constructed. To overcome the congestion problem, a traffic load-balancing algorithm is designed for fast transmission of the service traffic by combining three factors: network structure, congestion level and transmission service. The proposed algorithm is compared with current commonly load-balancing algorithms which are used in data center to verify its efficiency. Simulation experiments show that the algorithm improved transmission performance and transmission quality for the SDXAO-DCN.