Influence of time variable geopotential models on precise orbits of altimetry satellites,global and regional mean sea level trends

During the last decade a significant progress has been reached in the investigation of the gravity field of the Earth. Besides static, also time variable geopotential models have been recently created. In this paper we investigate the impact of the recent time variable geopotential models on altimetry satellite orbits and such altimetry products based on these orbits, as global and regional mean sea level trends. We show that the modeling of time variable gravity improves the orbit solutions, at least for the GRACE period where time variable gravity is sufficiently accurately observed by this mission. Our analysis includes six geopotential models jointly developed by GFZ German Research Centre for Geosciences and Space Geodesy Research Group (CNES/GRGS) Toulouse: the stationary model EIGEN-GL04S, a stationary version of EIGEN-6S (EIGEN-6S_stat), a corrected version of EIGEN-6S and three enhanced versions of EIGEN-6S called EIGEN-6S2, EIGEN-6S2A and EIGEN-6S2B. By “stationary” we mean “containing periodic parameters such as annual and semi-annual variations, but no secular (drift) terms”. We computed precise orbits for the radar altimetry satellites ERS-1, ERS-2, TOPEX/Poseidon, and Envisat over 20 years between 1991 and 2011. The orbit, single-mission and multi-mission altimetry crossover analyses show that the time variable models EIGEN-6S_corrected, EIGEN-6S2 and its two precursors EIGEN-6S2A/B perform notably better than the stationary models for the GRACE period from 2003 onwards. Thus, using EIGEN-6S2 and EIGEN-6S2A/B we have got 3.6% smaller root mean square fits of satellite laser ranging observations for Envisat, as when using EIGEN-GL04S. However, for the pre-GRACE period 1991–2003, the stationary geopotential models EIGEN-GL04S and EIGEN-6S_stat as well as EIGEN-6S2 having no drift terms for degree 3–50 at this time interval perform superior compared to EIGEN-6S_correct and EIGEN-6S2A/B which contain drifts for this period. We found, that the time variable geopotential models have a low (0.1–0.2 mm/yr) impact on our results for the global mean sea level trend. However, we found strong East/West differences up to 3 mm/yr in the regional mean sea level trends when using orbits of all four satellites based on time variable and stationary geopotential models. We show that these differences are related to the relative drifts of the centers-of-origin between the orbit solutions based on the time variable and stationary geopotential models. From the results of our detailed study, we conclude that the final version of the time variable gravity field model EIGEN-6S2 performs best for the four satellites tested. This model provides the most reliable and mission-consistent sea level estimates for the whole time period from 1992 to 2010. This model is of maximum spherical harmonic degree and order 260 and contains time series for drifts as well as annual and semiannual variations of the spherical harmonic coefficients for degree 2–50.     

Lunar geophysical parameters inversion based on gravity/topography admittance and particle swarm optimization

This paper presents a FORTRAN computer program. The program as code will be used for lunar parameter inversions based on gravity/topography admittance. This will be done by assuming that the lunar lithosphere is modeled as a thin elastic spherical shell. The parameters discussed here include; load ratio, crustal thickness, subsurface load depth, crustal density and elastic lithosphere thickness. The admittance of the best-fitting model can be found through automatically adjusting misfits between one theoretical admittance and an observed one. The results in this paper indicate that this research’s theoretical model is reasonable for exploring the best-fitting parameters. In addition, this code is not only able to automatically and simultaneously calculate the global optimum solution of the parameters studied, but also performs well in computational speed. The code can be easily modified to include more parameter inversions; such as the inversion for subsurface density anomaly and the case of considering infilling material in some lunar mare basins.     

Characteristics of the seasonal variation of the global tropopause revealed by COSMIC/GPS data

Using the Global Navigation Satellite System (GNSS) radio occultation observations from Formosa Satellite mission-3/Constellation Observing System for Meteorology, Ionosphere, and Climate (FORMOSAT-3/COSMIC) from 2007 to 2012, the climatological characteristics of the global tropopause was studied, with the following features identified. The overall results generally agree with previous studies. The tropopause has an obvious zonal structure, with more zonal characteristics in the Southern Hemisphere than the Northern Hemisphere. The vertical shape of the tropopause is sharp in the tropics and broad in the sub-tropical latitudes, with the sharpest latitudinal gradient in the mid-latitudes of both hemispheres. The global tropopause exists in a large range between 8 km and 17 km (or between 100 hPa and 340 hPa). The highest tropopause is over the South Asian monsoon regions for the entire year. The spatial structure of the tropopause in the polar region is of concentric structure, with an altitude between 7.5 km and 10 km. It is more symmetric in the Antarctic than the Arctic. Differing from other places, the height of the tropopause in the Antarctic is higher in winter as opposed to summer. The tropopause has distinct seasonal variability, especially in polar regions.     

Friction and Wear Properties of Amorphous Carbon Nitride Coatings in Water Lubrication

The friction and wear properties of amorphous carbon nitride(a-CN x)coatings in water lubrication were reviewed.The influences of mating materials and tribological variables such as normal load(W)and sliding speed(V)on the friction and wear properties of the a-CN x coatings were analyzed.It was indicated that the specific wear rate of the a-CN x coatings was related to the hydration reaction of mating materials with water.If the mating materials were easily hydrated,the specific wear rate of a-CN x coatings was low.The water-lubricated properties of the a-CN x coatings were better in comparison to the a-C coatings.The a-CN x/Si-based non-oxide ceramics tribo-pairs exhibited the lowest friction coefficient and wear rate.To describe their friction and wear properties at the normal loads of 3—15Nand the sliding speeds of 0.05—0.5m/s,the wear-mechanism maps for the a-CN x/SiC(Si3N4)tribo-pairs in water were developed.     

基于华氏宏观经济数学模型的高等职业教育动态投入产出模型

伴随着经济的发展和产业结构的变化,高等职业教育越来越受重视,高职教育的投入产出关系也引起了人们的关注。通过分析高职教育的投入产出关系,应用华氏宏观经济数学模型建立了高等职业教育动态投入产出模型。它对分析高职教育各项投入的效率具有实际意义。     

Design for aircraft engine multi-objective controllers with switching characteristics

The aircraft engine multi-loop control system is described and the switching control theory is introduced to solve the regulating and protecting control problems in this paper. The aircraft engine multi-loop control system is firstly described and the control problems are formulated. Secondly, the theory of the smooth switching control is devoted and a new extended scheme for the smooth switching of a switched control system is introduced. Then, for the key technologies of aero-engines switching control, a design algorithm is presented which can determine which candidate controller should be put in feedback with the plant to achieve a desired performance and the procedure to design the aircraft engine multi-loop control system is detailed. The switching performance objectives and the switching scheme are given and a family of PID controllers and compensators is designed. The simulation shows that using the switching control design method can not only improve the dynamic performance of the aircraft engine control system and reduce the switching times, but also guarantee the stability in some peculiar occasions.     

Numerical investigation on flow field characteristics of dual synthetic cold/hot jets using POD and DMD methods

In order to improve the control ability of synthetic jets in compressible boundary layer, a novel control method based on dual synthetic cold/hot jets coupled control of velocity profile and temperature profile was proposed. As fundamental investigations on the effects of synthetic jet temperature on the jet behavior and flow field characteristics were essentially necessary, preliminary numerical simulations were conducted to study the influence of temperature (200 K and 400 K) on the flow field characteristics of synthetic jets using Large Eddy Simulations (LES) model. Time-averaged flow fields showed that different temperatures led to variable behavior of two strands of jets. For dual synthetic cold jets, a potential-core arose apparently with its height ranging from 0.01 to 0.03 m, while for dual synthetic hot jets, two strands of jets emerged downstream. The modal decomposition of instantaneous flow fields had been done using both Proper Orthogonal Decomposition (POD) and Dynamic Mode Decomposition (DMD). Various modes showed different characteristics of the flow fields. As the POD method focuses on the energy of flow while the DMD method focuses on the frequency, the first two modes had many similarities, but the third and fourth modes demonstrated completely different vortex structures. The current researches play a role of preliminary investigations for further and comprehensive exploration of novel flow control measures in global velocity field.     

D-region ionospheric disturbances associated with the Extremely Severe Cyclone Fani over North Indian Ocean as observed from two tropical VLF stations

The D-region ionospheric disturbances due to the tropical cyclone Fani over the Indian Ocean have been analysed using Very Low Frequency (VLF) radio communication signals from three transmitters (VTX, NWC and JJI) received at two low latitude stations (Kolkata-CUB and Cooch Behar-CHB). The cyclone Fani formed from a depression on 26th April, 2019 over the Bay of Bengal (Northeastern part of the Indian Ocean) and turned into an extremely severe cyclone with maximum 1-min sustained winds of 250 km/h on 2 May, 2019 which made landfall on 3 May, 2019. Out of six propagation paths, five propagation paths, except the JJI-CHB which was far away from the cyclone track, showed strong perturbations beyond $3\sigma$ level compared to unperturbed signals. Consistent good correlations of VLF signal perturbations with the wind speed and cyclone pressure have been seen for both the receiving stations. Computations of radio signal perturbations at CUB and CHB using the Long Wave Propagation Capability (LWPC) code revealed a Gaussian perturbation in the D-region ionosphere. Analysis of atmospheric temperature at different layers from the NASA’s TIMED satellite revealed a cooling effect near the tropopause and warming effects near the stratopause and upper mesosphere regions on 3 May, 2019. This study shows that the cyclone Fani perturbed the whole atmosphere, from troposphere to ionosphere and the VLF waves responded to the disturbances in the conductivity profiles of the lower ionosphere.     

Interaction mechanisms of shock waves with the boundary layer and wakes in a highly-loaded NGV using hybrid RANS/LES

Accurate predictions of Shock Waves and Boundary Layer Interaction (SWBLI) and strong Shock Waves and Wake Vortices Interaction (SWWVI) in a highly-loaded turbine propose challenges to the currently widely used Reynolds-Averaged Navier-Stokes (RANS) model. In this work, the SWBLI and the SWWVI in a highly-loaded Nozzle Guide Vane (NGV) are studied using a hybrid RANS/LES strategy. The Turbulence Kinetic Energy (TKE) budget and the Proper Orthogonal Decomposition (POD) method are used to analyze flow mechanisms. Results show that this hybrid RANS/LES method can obtain detailed flow structures for flow mechanisms analysis. Strong shock waves induce boundary layer separation, while the presence of a separation bubble can in turn lead to a Mach reflection phenomenon. The shock waves cause trailing-edge vortices to break clearly, and the wakes, in turn, can change the shocks intensity and direction. Furthermore, the Entropy Generation Rate (EGR) is used to analyze the irreversible loss. It turns out that the SWWVI can reduce the flow field loss. There are several weak shock waves in the NGV flow field, which can increase the irreversible loss. This work offers flow mechanisms analysis and presents the EGR distribution in SWBLI and SWWVI areas in a transonic turbine blade.     

不同喷铝参数对复合材料雷击防护性能模拟

为了研究复合材料不同的雷击防护(lightning strike protection,LSP)系统在雷电流作用下的损伤规律,对雷击损伤过程和烧蚀机理进行分析,建立复合材料层合板雷击防护的能量平衡数学模型。在此基础上,在ABAQUS中建立复合材料基准件、全喷铝和局部喷铝防护系统的碳纤维增强复合材料(carbon fiber reinforced polymer,CFRP)层合板电-热耦合有限元模型,和实验结果对比验证模型的有效性,对雷击烧蚀损伤特征进行分析,并引入雷击烧蚀损伤指数DI,得出三种不同模型在不同铝涂层厚度、不同峰值雷电流作用下的烧蚀损伤规律,并对两种不同的局部喷铝防护系统下复合材料的损伤特征进行对比分析。结果表明:雷击防护系统下复合材料的烧蚀损伤面积和铝涂层厚度关系进行函数拟合,两者均满足幂函数关系。