Combination of terrestrial reference frames based on space geodetic techniques in SHAO:methodology and main issues

Based on years of input from the four geodetic techniques （SLR, GPS, VLBI and DORIS）, the strategies of the combination were studied in SHAO to generate a new global terrestrial reference frame as the material realization of the ITRS defined in IERS Conventions. The main input includes the time series of weekly solutions （or fortnightly for SLR 1983-1993） of observational data for satellite techniques and session-wise normal equations for VLBI. The set of estimated unknowns includes 3- dimensional Cartesian coordinates at the reference epoch 2005.0 of the stations distributed globally and their rates as well as the time series of consistent Earth Orientation Parameters （EOPs） at the same epochs as the input. Besides the final solution, namely SOL-2, generated by using all the inputs before 2015.0 obtained from short-term observation processing, another reference solution, namely SOL- 1, was also computed by using the input before 2009.0 based on the same combination of procedures for the purpose of comparison with ITRF2008 and DTRF2008 and for evaluating the effect of the latest six more years of data on the combined results. The estimated accuracy of the x-component and y-component of the SOL- 1 TRF-origin was better than 0.1 mm at epoch 2005.0 and better than 0.3 mm yr- 1 in time evolution, either compared with ITRF2008 or DTRF2008. However, the z-component of the translation parameters from SOL-1 to ITRF2008 and DTRF2008 were 3.4 mm and -1.0 ram, respectively. It seems that the z-component of the SOL-1 TRF-origin was much closer to the one in DTRF2008 than the one in ITRF2008. The translation parameters from SOL-2 to ITRF2014 were 2.2, -1.8 and 0.9 mm in the x-, y- and z-components respectively with rates smaller than 0.4 mmyr-1. Similarly, the scale factor transformed from SOL-1 to DTRF2008 was much smaller than that to ITRF2008. The scale parameter from SOL-2 to ITRF2014 was -0.31 ppb with a rate lower than 0.01 ppb yr-1. The external precision （WRMS） compared with IERS EOP 08 C04 of the combined EOP series was smaller than 0.06 mas for the polar motions, smaller than 0.01 ms for the UT1-UTC and smaller than 0.02 ms for the LODs. The precision of the EOPs in SOL-2 was slightly higher than that of SOL-1.

Hopping Parameter Expansion Technique in Noise Method for Disconnected Quark Loops

We investigate the effectiveness of the hopping parameter expansion（HPE） combined with the Z（2） noise method in the calculation of the trace of the inverse of Wilson＇s Dirac operator and some other disconnected contributions.A numerical comparison of the standard deviation for the Z（2） noise method and HPE with the Z（2） noise method is carried out. It is found that there are noise reductions in all the quantities we calculated using the HPE with the Z（2） noise method. For the trace of the inverse of Wilson＇s Dirac operator, the HPE can reduce the statistical error by about 60%.

Inexact Newton method via Lanczos decomposed technique for solving box-constrained nonlinear systems

This paper proposes an inexact Newton method via the Lanczos decomposed technique for solving the box-constrained nonlinear systems. An iterative direction is obtained by solving an affine scaling quadratic model with the Lanczos decomposed technique. By using the interior backtracking line search technique, an acceptable trial step length is found along this direction. The global convergence and the fast local convergence rate of the proposed algorithm are established under some reasonable conditions. Furthermore, the results of the numerical experiments show the effectiveness of the pro- posed algorithm.

Methods and Techniques of Electricity Thieving in Pakistan

The situation of electricity in Pakistan has been alarming from the last ten years. The deficiency in electricity has not only obstructed the business activities but also affected the domestic consumers, educational institutes and hospitals. Usually electrical power companies are liable for electricity shortfall and power interruption. However, electricity end consumers are also equally responsible behind strange shortfall and unusual power interruption. Frequently, the consumers use the heavy electrical equipment in their homes including heaters, geysers, irons and water motors which causes the more electricity consumption, load shedding and huge amount of bills. For escaping the huge amount of electricity bill, the consumers commit the illegal and unethical connections. The illegal usage of electrical power failed the power companies to plan schedule load shedding accordingly and the other side the damaged electricity wires or Pole Mount Transformer increased faults due to overburdening, which directly affected on extend power interruption. In addition that, responsible teams of electrical power companies cannot reach instantly to repair faults and prevent the theft. Electricity thieving is social crime committed by the consumers or meter readers which causes the electricity strange shortfall within country. This paper presents the practical demonstration about the common energy theft methods and techniques done by electricity consumers within their home and residential building. In Pakistan EPC (electrical power companies) deploy the traditional electromechanical meters for electricity consumption measurements, however, these meters do not have any real time communication. Therefore there are many easy ways to manipulate the meter reading as well as internal structural of metering system.

Direct Ion Beam Figuring Process and Rotational Measurement Method for Ultra-smooth Aspherical Surfaces of a 46.5 nm Telescope

This paper describes a fabrication process for the hyperboloidal concave mirror of a 46.5 nm telescope. The180 mm aperture hyperboloidal concave mirror and 70 mm aperture compensator are machined directly from chemical mechanical polishing of a spherical surface to a high-accuracy aspherical surface by ion beam figuring.The aspherical measurement method is the Dall null test. To minimize system errors in the measurement process,the rotational measurement method with six rotations is used in the null test. The results of the analysis for the ME(first solve the machined surface profile, then solve the system errors) and EM(first solve the system errors, then solve the machined surface profile) methods of calculation in the measurement are given. The ME method is a more accurate rotational test method, and the six rotations are appropriate for rotational measurements. After the figuring process, the hyperboloidal concave mirror surface profile reached 8.27 nm rms and the compensator surface profile is approximately 4 nm rms. The roughness of the hyperboloidal concave mirror is smooth to0.160 nm rms.

The Velocity Aberration Effect of the CSST Main Survey Camera

In this study,we conducted simulations to find the geometric aberrations expected for images taken by the Main Survey Camera of the Chinese Space Station Telescope(CSST)due to its motion.As anticipated by previous work,our findings indicate that the geometric distortion of light impacts the focal plane's apparent scale,with a more pronounced influence as the size of the focal plane increases.Our models suggest that the effect consistently influences the pixel scale in both the vertical and parallel directions.The apparent scale variation follows a sinusoidal distribution throughout one orbital period.Simulations reveal that the effect is particularly pronounced in the center of the Galaxy and gradually diminishes along the direction of ecliptic latitude.At low ecliptic latitudes,the total aberration leads to about a 0.94 pixel offset(a 20 minute exposure)and a 0.26 pixel offset(a 300 s exposure)at the edge of the field of view.Appropriate processings for the geometric effect during the CSST pre-and post-observation phases are presented.

Digital Beamforming Demonstration of a Microstrip Antenna Phased Array Feeds(PAF)at 1.25 GHz

We introduce the structure of a radio astronomy phased array feeds(PAF)beamforming demonstrator.In a laboratory environment,we have demonstrated beamforming on a received 1.25 GHz sinusoidal signal and used digital weighting techniques to plot the 2D pattern of the PAF.The radio frequency part of the demonstrator includes a 4×4 linearly polarized microstrip antenna array,all of which is connected in series with a low-noise amplifier.The signals from the central 4×2 array elements are injected into a radio frequency system-on-chip digital board,which can receive eight inputs with a bandwidth of 512 MHz.Combining the principle of undersampling,the beamforming is completed at a frequency of 1.25 GHz for the offline data,and a 2D image of the beam is plotted using beam scanning technology.

A Wideband Microwave Holography Methodology for Reflector Surface Measurement of Large Radio Telescopes

Most reflector surface holographic measurements of a large radio telescope utilize a geostationary satellite as the signal source. The shortcoming is that those measurements could only be done at a limited elevation angle due to the satellite’s relatively stationary state. This paper proposed a new wideband microwave holographic measurement method based on radio sources to achieve full-elevation-angle measurement with small size reference antenna. In theoretical derivation, the time delay and phase change due to path length and device difference between the antenna under test and reference antenna are compensated first. Then the correct method of wideband holography effect, which is because of antenna pattern differing under different wavelengths when receiving a wideband signal, is presented. To verify the proposed methodology, a wideband microwave holographic measurement system is established, the data processing procedure is illustrated, and the reflector surface measurement experiments on a 40 m radio telescope at different elevation angles are conducted. The result shows that the primary reflector surface root-mean-square at around elevation angles of 28°, 44°, 49°, and 75° are respectively 0.213 mm, 0.170 mm, 0.188 mm, and 0.199 mm. It is basically consistent with the real data, indicating that the proposed wideband microwave holography methodology is feasible.

A method and results of color calibration for the Chang'e-3 terrain camera and panoramic camera

The terrain camera （TCAM） and panoramic camera （PCAM） are two of the major scientific payloads installed on the lander and rover of the Chang＇e 3 mission re- spectively. They both use a Bayer color filter array covering CMOS sensor to capture color images of the Moon＇s surface. RGB values of the original images are related to these two kinds of cameras. There is an obvious color difference compared with human visual perception. This paper follows standards published by the International Commission on Illumination to establish a color correction model, designs the ground calibration experiment and obtains the color correction coefficient. The image qual- ity has been significantly improved and there is no obvious color difference in the corrected images. Ground experimental results show that： （1） Compared with uncor- rected images, the average color difference of TCAM is 4.30, which has been reduced by 62.1%. （2） The average color differences of the left and right cameras in PCAM are 4.14 and 4.16, which have been reduced by 68.3% and 67.6% respectively.

Method for array gain and phase uncertainties calibration based on ISM and ESPRIT

A new method for array calibration of array gain and phase uncertainties, which severely degrade the performance of spatial spectrum estimation, is presented. The method is based on the idea of the instrumental sensors method （ISM）, two well-calibrated sensors are added into the original array. By applying the principle of estimation of signal parameters via rotational invariance techniques （ESPRIT）, the direction-of-arrivals （DOAs） and uncertainties can be estimated simultaneously through eigen-decomposition. Compared with the conventional ones, this new method has less computational complexity while has higher estimation precision, what＇s more, it can overcome the problem of ambiguity. Both theoretical analysis and computer simulations show the effectiveness of the proposed method.

Simulation of oxygen transfer in liquid lead under influence of nanoparticles by using lattice Boltzmann method

Oxygen transfer presents a serious challenge in the application of liquid lead as a nuclear coolant in advanced reactors. To mitigate corrosion by liquid lead in contact with steel, carefully controlling the oxygen concentration has been used as an effective way. Oxygen needs to mix in liquid lead uniformly and quickly. To enhance oxygen transport in liquid lead, nanoparticles are added to the liquid metal. In the current study, a lattice Boltzmann method is applied to investigate natural convection of copper/lead and aluminum oxide/lead in two-dimensional simplified container. Two thermal boundary cases are evaluated in order to check the effect of different natural convection flow patterns on oxygen transport. Some useful information are obtained such as improvement in natural convection and reduction in oxygen equilibrium time.

GPU Accelerated Marine Data Visualization Method

The study of marine data visualization is of great value. Marine data, due to its large scale, random variation and multiresolution in nature, are hard to be visualized and analyzed. Nowadays, constructing an ocean model and visualizing model results have become some of the most important research topics of ‘Digital Ocean'. In this paper, a spherical ray casting method is developed to improve the traditional ray-casting algorithm and to make efficient use of GPUs. Aiming at the ocean current data, a 3D view-dependent line integral convolution method is used, in which the spatial frequency is adapted according to the distance from a camera. The study is based on a 3D virtual reality and visualization engine, namely the VV-Ocean. Some interactive operations are also provided to highlight the interesting structures and the characteristics of volumetric data. Finally, the marine data gathered in the East China Sea are displayed and analyzed. The results show that the method meets the requirements of real-time and interactive rendering.

Mode Ⅲ fracture analysis by Trefftz boundary element method

This paper presents a hybrid Trefftz （HT） boundary element method （BEM） by using two indirect techniques for mode III fracture problems. Two Trefftz complete functions of Laplace equation for normal elements and a special purpose Trefftz function for crack elements are proposed in deriving the Galerkin and the collocation techniques of HT BEM. Then two auxiliary functions are introduced to improve the accuracy of the displacement field near the crack tips, and stress intensity factor （SIF） is evaluated by local crack elements as well. Furthermore, numerical examples are given, including comparisons of the present results with the analytical solution and the other numerical methods, to demonstrate the efficiency for different boundary conditions and to illustrate the convergence influenced by several parameters. It shows that HT BEM by usingthe Galerkin and the collocation techniques is effective for mode III fracture problems.

Noise reduction methods in the analysis of near infrared lunar occultation light curves for high angular resolution measurements

A lunar occultation （LO） technique in the near-infrared （NIR） provides angular resolution down to milliarcseconds for an occulted source, even with ground- based 1 m class telescopes. LO observations are limited to brighter objects because they require a high signal-to-noise ratio （S/N ~40） for proper extraction of angular diameter values. Hence, methods to improve the S/N ratio by reducing noise using Fourier and wavelet transforms have been explored in this study. A sample of 54 NIR LO light curves observed with the IR camera at Mt Abu Observatory has been used. It is seen that both Fourier and wavelet methods have shown an improvement in S/N compared to the original data. However, the application of wavelet transforms causes a slight smoothing of the fringes and results in a higher value for angular diameter. Fourier transforms which reduce discrete noise frequencies do not distort the fringe. The Fourier transform method seems to be effective in improving the S/N, as well as improving the model fit, particularly in the fainter regime of our sample. These methods also provide a better model fit for brighter sources in some cases, though there may not be a significant imorovement in S/N.

How to co-add images? I. A new iterative method for image reconstruction of dithered observations

By employing the previous Voronoi approach and replacing its nearest neighbor approx- imation with Drizzle in iterative signal extraction, we develop a fast iterative Drizzle algorithm, namedfiDrizzle, to reconstruct the underlying band-limited image from undersampled dithered frames. Compared with the existing iDrizzle, the new algorithm improves rate of convergence and accelerates the computational speed. Moreover, under the same conditions （e.g. the same number of dithers and iterations）, fiDrizzle can make a better quality reconstruction than iDrizzle, due to the newly discov- ered High Sampling caused Decelerating Convergence （HSDC） effect in the iterative signal extraction process.fiDrizzle demonstrates its powerful ability to perform image deconvolution from undersampled dithers.

A fast Stokes inversion technique based on quadratic regression

Stokes inversion calculation is a key process in resolving polarization information on radiation from the Sun and obtaining the associated vector magnetic fields. Even in the cases of simple local thermo- dynamic equilibrium （LTE） and where the Milne-Eddington approximation is valid, the inversion problem may not be easy to solve. The initial values for the iterations are important in handling the case with mul- tiple minima. In this paper, we develop a fast inversion technique without iterations. The time taken for computation is only 1/100 the time that the iterative algorithm takes. In addition, it can provide available initial values even in cases with lower spectral resolutions. This strategy is useful for a filter-type Stokes spectrograph, such as SDO/HMI and the developed two-dimensional real-time spectrograph （2DS）.

Narrowband EUV Sc/Si Multilayer for the Solar Upper Transition Region Imager at 46.5 nm

The transition region is the key region between the lower solar atmosphere and the corona, which has been limitedly understood by human beings. Therefore, the Solar Upper Transition Region Imager(SUTRI) was proposed by Chinese scientists and launched in 2022 July. Right now, the first imaging observation of the upper transition region around 46.5 nm has been carried out by SUTRI. To ensure the spectral and temporal resolution of the SUTRI telescope, we have developed a narrowband Sc/Si multilayer. Based on the extreme ultraviolet(EUV)reflectivity measurements, the multilayer structure has been modified for ensuring the peak position of reflectivity was at 46.5 nm. Finally, the narrowband Sc/Si multilayer was successfully deposited on the hyperboloid primary mirror and secondary mirrors. The deviation of multilayer thickness uniformity was below than 1%, and the average EUV reflectivity at 46.1 nm was 27.8% with a near-normal incident angle of 5°. The calculated bandwidth of the reflectivity curve after primary and secondary mirrors was 2.82 nm, which could ensure the requirements of spectral resolution and reflectivity of SUTRI telescope to achieve its scientific goals.

The Edge Sensor of Segmented Mirror Based on Fringes of Equal Thickness

Co-phase and co-focus detection is one of the key technologies for large-aperture segmented mirror telescopes.In this paper,a new edge sensor based on fringes of equal thickness is developed,which can detect each segment's relative piston,tilt,and tip errors from the interferograms.Based on the co-focus demand for many ground-based seeing limited segmented mirror telescopes,an edge sensor prototype based on such a principle is built and applied in the indoor segmented mirror experiment system in the lab.According to the co-focus requirement of the Large Sky Area Multi-Object Fiber Spectroscopic Telescope,many simulations and experiments are carried out for cofocus error detection of the segmented mirror system.Experiment results show that the co-focus accuracy is better than 002 rms,which can meet the co-focus requirements of most large or extremely large segmented mirror astronomical telescopes.

The Chromatic Point-spread Function of Weak Lensing Measurement in the Chinese Space Station Survey Telescope

Weak gravitational lensing is a powerful tool in modern cosmology.To accurately measure the weak lensing signal,one has to control the systematic bias on a small level.One of the most difficult problems is how to correct the smearing effect of the Point-Spread Function(PSF)on the shape of the galaxies.The chromaticity of PSF for a broad-band observation can lead to new subtle effects.Since the PSF is wavelength-dependent and the spectrum energy distributions between stars and galaxies are different,the effective PSF measured from the star images will be different from those that smear the galaxies.Such a bias is called color bias.We estimate it in the optical bands of the Chinese Space Station Survey Telescope from simulated PSFs,and show the dependence on the color and redshift of the galaxies.Moreover,due to the spatial variation of spectra over the galaxy image,another higher-order bias exists:color gradient bias.Our results show that both color bias and color gradient bias are generally below 0.1%in CSST.Only for small-size galaxies,one needs to be careful about the color gradient bias in the weak lensing analysis using CSST data.

Distributed Control Software for the Active Surface System of Tian-ma Radio Telescope

The Tian-ma Radio Telescope(TMRT) applies an Active Surface System(ASFS), which corrects for large-scale deformations due to gravity and thermal on the primary reflector. The centralized and automated management of the ASFS using software has become a challenge, for which we have developed the TMRT Active Surface System Control Software(TASCS). This paper describes the design and implementation of TASCS for device control,status monitoring, human-computer interaction, and data management functionalities. TASCS adopts the opensource Tango Controls framework and distributes middleware technology to realize real-time automated adjustment of the primary reflector through remote centralized control of a large number of actuators. At present, it has been successfully deployed on the TMRT and has played an important role in Event Horizon Telescope observations.