Orbit determination of BDS-3 satellite based on regional ground tracking station and inter-satellite link observations

The BeiDou global navigation satellite system (BDS-3) has established the Ka-band inter-satellite link (ISL) to realize a two-way ranging function between satellites, which provides a new observation technology for the orbit determination of BDS-3 satellites. Therefore, this study presents a BDS satellite orbit determination model based on ground tracking station (GTS) observations and ISL ranging observations firstly to analyze the impact of the ISL ranging observations on the orbit determination of BDS-3 satellites. Subsequently, considering the data fusion processing, the variance component estimation (VCE) algorithm is applied to the parameter estimation process of the satellite orbit determination. Finally, using the measured data from China’s regional GTS observations and BDS-3 ISL ranging observations, the effects of ISL ranging observations on the orbit determination accuracy of BDS-3 satellites are analyzed. Moreover, the impact of the VCE algorithm on the fusion data processing is evaluated from the aspects of orbit determination accuracy, Ka-band hardware delay parameter stability, and ISL ranging observation residuals. The results show that for China’s regional GTSs, the addition of BDS-3 ISL ranging observations can significantly improve the orbit determination accuracy of BDS-3 satellites. The observed orbit determination accuracy of satellite radial component is improved from 48 cm to 4.1 cm. In addition, when the initial weight ratio between GTS observations and ISL ranging observations is not appropriate, the various indicators which include orbit determination accuracy, ISL hardware delay, and ISL observation residuals were observed to have improved after the adjustment of the VCE algorithm. These results validate the effectiveness of the VCE algorithm for the fusion data processing of the GTS observations and ISL ranging observations.     

Validation of non-linear split window algorithm for land surface temperature estimation using Sentinel-3 satellite imagery: Case study; Tehran Province,Iran

In recent years, land surface temperature (LST) has become critical in environmental studies and earth science. Remote sensing technology enables spatiotemporal monitoring of this parameter on large scales. This parameter can be estimated by satellite images with at least one thermal band. Sentinel-3 SLSTR data provide LST products with a spatial resolution of 1 km. In this research, direct and indirect validation procedures were employed to evaluate the Sentinel-3 SLSTR LST products over the study area in different seasons from 2018 to 2019. The validation method was based on the absolute (direct) evaluation of this product with field data and comparison (indirect) evaluation with the MODIS LST product and the estimated LST using the non-linear split-window (NSW) algorithm. Also, two emissivity estimation methods, (1) NDVI thresholding method (NDVI-THM) and (2) classification-based emissivity method (CBEM), were used to estimate the LST using the NSW method according to the two thermal bands of Sentinel-3 images. Then, the accuracy of these methods in estimating LST was evaluated using field data and temporal changes of vegetation, which the NDVI-THM method generated better results. For indirect evaluation between the Sentinel-3 LST product, MODIS LST product, and LST estimated using NSW, four filters based on spatial and temporal separates between pairs of pixels and pixel quality were used to ensure the accuracy and consistency of the compared pairs of a pixel. In general, the accuracy results of the LST products of MODIS and Sentinel-3, and LST estimated using NSW showed a similar trend for LST changes during the seasons. With respect to the two absolute and comparative validations for the Sentinel-3 LST products, summer with the highest values of bias (?1.24 K), standard deviation (StDv = 2.66 K), and RMSE (2.43 K), and winter with the lowest ones (bias of 0.14 K, StDv of 1.13 K, and RMSE of 1.12 K) provided the worst and best results for the seasons in the period of 2018–2019, respectively. According to both absolute and comparative evaluation results, the Sentinel-3 SLSTR LST products provided reliable results for all seasons on a large temporal and spatial scale over our studied area.     

Topology optimization in lightweight design of a 3D-printed flapping-wing micro aerial vehicle

Topology optimization is an effective method to obtain a lightweight structure that meets the requirements of structural strength. Whether the optimization results meet the actual needs mainly depends on the accuracy of the material properties and the boundary conditions, especially for a tiny Flapping-wing Micro Aerial Vehicle (FMAV) transmission system manufactured by 3D printing. In this paper, experimental and numerical computation efforts were undertaken to gain a reliable topology optimization method for the bottom of the transmission system. First, the constitutive behavior of the ultraviolet (UV) curable resin used in fabrication was evaluated. Second, a numerical computation model describing further verified via experiments. Topology optimization modeling considering nonlinear factors, e.g. contact, friction and collision, was presented, and the optimization results were verified by both dynamic simulation and experiments. Finally, detailed discussions on different load cases and constraints were presented to clarify their effect on the optimization. Our methods and results presented in this paper may shed light on the lightweight design of a FMAV.     

The adjusted optical properties for Galileo/BeiDou-2/QZS-1 satellites and initial results on BeiDou-3e and QZS-2 satellites

Solar Radiation Pressure (SRP) is the dominant non-gravitational perturbation for GNSS (Global Navigation Satellite System) satellites. In the absence of precise surface models, the Empirical CODE Orbit Models (ECOM, ECOM2) are widely used in GNSS satellite orbit determination. Based on previous studies, the use of an a priori box-wing model enhances the ECOM model, especially if the spacecraft is a stretched body satellite. However, so far not all the GNSS system providers have published their metadata. To ensure a precise use of the a priori box-wing model, we estimate the optical parameters of all the Galileo, BeiDou-2, and QZS-1 (Quasi Zenith Satellite System) satellites based on the physical processes from SRP to acceleration. Validation using orbit prediction proves that the adjusted parameters of Galileo and QZS-1 satellites exhibit almost the same performance as the corresponding published and “best guess” values. Whereas, the estimated parameters of BeiDou-2 satellites demonstrate an improvement of more than 60% over the initial “guess” values. The resulting optical parameters of all the satellites are introduced into an a priori box-wing model, which is jointly used with ECOM and ECOM2 model in the orbit determination. Results show that the pure ECOM2 model exhibits better performance than the pure ECOM model for Galileo, BeiDou-2 GEO and QZS-1 orbits. Combined with the a priori box-wing model the ECOM model (ECOM+BW) results in the best Galileo, BeiDou-2 GEO and QZS-1 orbits. The standard deviation (STD) of satellite laser ranging residuals reduce by about 20% and 5% with respect to the pure ECOM2 model for Galileo and BeiDou-2 GEO orbits, while the reductions are about 40% and 60% for QZS-1 orbits in yaw-steering and orbit-normal mode respectively. BeiDou-2 IGSO and MEO satellite orbits do not benefit much from the a priori box-wing model. In summary, we suggest setting up a unified SRP model of ECOM+BW for Galileo, QZS-1, and BeiDou-2 orbits based on the adjusted metadata. In addition, we estimate the optical parameters of BeiDou-3e and QZS-2 satellites using a limited number of tracking stations. Results regarding the unified SRP model indicate the same advantages, the STD of satellite laser ranging residuals reduces by about 30% and 20% for QZS-2 and BeiDou-3e orbits respectively over orbit products without a priori model. The estimation procedure is effective and easy to apply to the new emerging satellites in the future.     

对转压气机外伸激波对叶顶泄漏流影响

摘要:为了揭示对转压气机下游转子外伸激波对上游转子泄漏流的影响规律,针对上游转子叶顶间隙分别为0.2、0.5、0.8 mm的对转压气机开展了非定常数值模拟研究。研究发现:受下游转子外伸激波掠扫影响,上游转子尾缘附近压力面会形成弱压缩波,且随上游转子泄漏流增强而逐渐减弱;而该外伸激波在上游转子尾缘附近吸力面,会形成与型线切向相垂直的较强压缩波,且其位置基本不受叶顶间隙大小影响;外伸激波使上游转子尾缘附近吸、压力面压差增大,叶顶泄漏流增强,进而导致其损失增大;随着叶顶间隙增大,上游转子叶尖区弦长前半段压力波动的频率,由通道激波转为叶顶泄漏流主导,且呈现减小的趋势,而弦长后半段压力波动的频率主要由外伸激波主导,且基本不变。      

某航天器输氢管道系统结构完整性评估

为保障某航天器在服役年限内不发生事故，针对航天器上采用插入焊接技术的精密管道系统进行结构完整性评估。在考虑到系统存在整圈的环形未焊透缺陷基础上，结合工程实践，假设在焊接区域存在更危险的平面缺陷。首先，开展高压气相热充氢试验下的标准拉伸和三点弯试验得到管道系统材料"J-75"在临氢环境下应力应变曲线和断裂韧性。其次，采用有限元程序计算系统在内压、惯性载荷与离面位移作用下的应力分布，考虑到间隙配合中的不确定性，改变模型中的最大间隙多次模拟计算，发现最大应力位置基本不变，以此区域作为假想缺陷的具体萌生区域。最后，采用失效评估图（FAD）方法对存在4种缺陷模式的系统分别进行结构完整性评估。通过引入安全裕度的概念，发现系统在含有初始尺寸（0.15 mm）的4种平面缺陷时，均处于安全状态，并且安全裕度都在2倍以上；通过工程临界分析（ECA）发现随着缺陷尺寸的增大，系统中含有缺陷2（位于管道焊接区域最大应力处的轴对称平面且沿焊深方向扩展的缺陷）时最为危险，且当缺陷深度超过0.61 mm时将断裂失效；考虑到测试数据的弥散性，针对系统最危险的模式，以安全裕度为可靠性指标，将失效评估图方法与可靠性评价相结合，得到系统在含有尺寸为0.15 mm的缺陷2情况下保有足够安全裕度的可靠性下限值超过0.999 5。     

过氧化氢/紫外线/臭氧氧化技术处理肼类污水的应用

采用放大试验规模（处理量1 m³/h）的H₂O₂（过氧化氢）/UV（紫外线）/O₃（臭氧）氧化技术处理肼类推进剂污水，在30.0±1.6 ℃，pH=9.0±0.2的条件下，研究肼类的过氧化氢增强光解臭氧化降解。对比不同氧化技术的协同效应以及对污水降解的影响，重点考察过氧化氢、紫外线、臭氧和初始浓度等工艺参数对降解效果的影响，对氧化技术的应用进行了优化。研究结果表明：该技术可使COD去除率提高27.66%，肼类的降解速率随着过氧化氢投加量、紫外线辐射强度、O₃投加速率和水质地提高而升高，随着初始质量浓度的提高而下降。在最佳工艺下，处理5 000 mg/L质量浓度的废水，处理60 min时，COD去除率分别为98.62%（偏二甲肼）、99.17%（甲基肼）、99.94%（肼）和93.25%（单推-3）。     

双模块下颌式内转进气道/圆锥前体一体化布局研究

为实现三维内转进气道的内收缩流场与圆锥前体的外压缩流场的良好匹配，提出了一种双模块下颌式内转进气道/圆锥前体(Double-modules Chin Inward-turning Inlet and Conical Forebody，DCII/CF)一体化设计方法，获得一种新颖的双发并置、侧向安装的DCII/CF一体化布局。针对该布局形式，开展了DCII/CF一体化构型与传统的单模块内转进气道/圆锥前体(Single-module Inward-turning Inlet and Conical Forebody，SII/CF)一体化构型的数值对比研究。结果表明：DCII/CF一体化布局不仅为内转进气道提供了优秀的前体附面层排移效果，还有效避免了传统SII/CF布局中前体附面层与进气道内部流场之间的相互干扰。在Ma∞=6.0设计状态，DCII/CF一体化布局的进气道总压恢复系数相较传统的SII/CF布局有了显著提高，从0.403提高至0.482；但由于前体附面层的排移，该布局的捕获流量略有降低， SII/CF的流量系数为0.956，该布局则为0.917。而在非设计状态，该布局形式同样具备较好的总压恢复性能，在Ma∞=5.0与Ma∞=4.0的总压恢复系数分别达到了0.586和0.682，明显高于SII/CF的总压恢复系数0.507和0.619。     

Daytime VHF amplitude scintillations recorded at an Indian low-latitude station,Waltair (17.7°N, 83.3°E) during 1997–2003

In this research, it is presented the daytime amplitude scintillations recorded at VHF frequency (244 MHz) at an Indian low-latitude station, Waltair (17.7°N, 83.3°E) during seven continuous years (1997–2003). Contrary to the nighttime scintillation seasonal trends, the occurrence of daytime scintillations maximizes during summer followed by winter and the equinox seasons. The fade depths, scintillation indices and the patch durations of daytime scintillations are meager when compared with their nighttime counterparts. A co-located digital high frequency (HF) ionosonde radar confirms the presence of sporadic (Es) layers when daytime scintillations are observed. The presence of daytime scintillations is evident when the critical frequency of the Es-layer (foEs) is ≥4 MHz and Es-layers are characterized by a highly diffuse range spread Es echoes as can be seen on ionograms. It is surmised that the gradient drift instability (GDI) seems to be the possible mechanism for the generation of these daytime scintillations. It is quite likely that the spread Es-F-layer coupling is done through polarization electric fields (Ep) that develop inside the destabilized patches of sporadic E layers, which are mapped up to the F region along the field lines as to initiate the daytime scintillations through the GDI mechanism. Further, the presence of additional stratification of ionosphere F-layer, popularly known as the F3-layer, is observed on ionograms once the Es-layers and daytime scintillations are ceased.     

扩压叶栅叶尖非定常流动的涡动力学机理

为了揭示压气机叶尖区旋涡结构与流动非定常性之间的关联，采用URANS对一亚声速平面扩压叶栅在不同攻角下的流场进行了求解，并借助Q判据提取了叶尖瞬态涡系结构。结果表明：泄漏涡的破碎现象能够通过诱导新的涡结构间接作用于相邻通道的叶尖流动，是导致叶尖流场失稳的关键因素。在-0.3°和+0.7°攻角下，叶尖泄漏涡发生了螺旋破碎，并伴随有非定常诱导涡的出现，诱导涡对相邻叶片载荷的影响使得叶尖泄漏涡发生周期性摆动；在+1.7°攻角下，泄漏涡破碎会导致反流涡的形成，反流涡的输运会给叶片载荷和来流攻角带来非定常扰动，反过来又会作用于泄漏涡的破碎和反流涡的生成，最终表现为一种自维持的非定常流动现象。