Fast solar radiation pressure modelling with ray tracing and multiple reflections

Physics based SRP (Solar Radiation Pressure) models using ray tracing methods are powerful tools when modelling the forces on complex real world space vehicles. Currently high resolution (1?mm) ray tracing with secondary intersections is done on high performance computers at UCL (University College London). This study introduces the BVH (Bounding Volume Hierarchy) into the ray tracing approach for physics based SRP modelling and makes it possible to run high resolution analysis on personal computers. The ray tracer is both general and efficient enough to cope with the complex shape of satellites and multiple reflections (three or more, with no upper limit). In this study, the traditional ray tracing technique is introduced in the first place and then the BVH is integrated into the ray tracing. Four aspects of the ray tracer were tested for investigating the performance including runtime, accuracy, the effects of multiple reflections and the effects of pixel array resolution.Test results in runtime on GPS IIR and Galileo IOV (In Orbit Validation) satellites show that the BVH can make the force model computation 30–50 times faster. The ray tracer has an absolute accuracy of several nanonewtons by comparing the test results for spheres and planes with the analytical computations. The multiple reflection effects are investigated both in the intersection number and acceleration on GPS IIR, Galileo IOV and Sentinel-1 spacecraft. Considering the number of intersections, the 3rd reflection can capture $99.12\%\text{,}99.14\%$, and $91.34\%$ of the total reflections for GPS IIR, Galileo IOV satellite bus and the Sentinel-1 spacecraft respectively. In terms of the multiple reflection effects on the acceleration, the secondary reflection effect for Galileo IOV satellite and Sentinel-1 can reach 0.2?$\text{nm}/{\text{s}}^2$ and 0.4?$\text{nm}/{\text{s}}^2$ respectively. The error percentage in the accelerations magnitude results show that the 3rd reflection should be considered in order to make it less than $0.035\%$. The pixel array resolution tests show that the dimensions of the components have to be considered when choosing the spacing of the pixel in order not to miss some components of the satellite in ray tracing. This paper presents the first systematic and quantitative study of the secondary and higher order intersection effects. It shows conclusively the effect is non-negligible for certain classes of misson.     

基于统计学特征的密集型气液雾化随机模型优化方法研究

在高气液动量比的空气雾化流场中，液滴和液丝从液核剥离的过程具有高自由度、分布密集的特点，传统理论模型难以对其准确预测。本文对结合大涡模拟方法的随机雾化模型进行优化，在初始雾化过程，提出使用液滴统计平均温度来表征液滴碰撞统计学特性的改进方法。液滴的统计平均温度分别采用气液相对动能模型的粒子追踪法和亚网格动能模型的粒子追踪法。研究表明，使用改进的气液雾化随机模型预测密集型空气雾化流场，大幅改善了传统雾化随机模型在初始雾化区域过预测的缺陷，平均动能的相对误差为15.5%，平均索特尔直径的相对误差为7.2%，与未改进前的模拟结果相比，误差降低了41.1%和15.0%。此外，本文还探究了喷雾张角模型对雾化流场预测结果的影响，分别将实验所得经验公式法、相界面气液动量平衡所得模拟法与亚网格动能模型的粒子追踪法结合。结果表明喷雾张角经验公式预测结果更为准确，在平均索特尔直径预测方面准确性提高了17.3%。     

First measurements of periodicities and anisotropies of cosmic ray flux observed with a water-Cherenkov detector at the Marambio Antarctic base

A new water-Cherenkov radiation detector, located at the Argentine Marambio Antarctic Base (64.24S-56.62 W), has been monitoring the variability of galactic cosmic ray (GCR) flux since 2019. One of the main aims is to provide experimental data necessary to study interplanetary transport of GCRs during transient events at different space/time scales. In this paper we present the detector and analyze observations made during one full year. After the analysis and correction of the GCR flux variability due to the atmospheric conditions (pressure and temperature), a study of the periodicities is performed in order to analyze modulations due to heliospheric phenomena. We can observe two periods: (a) 1 day, associated with the Earth’s rotation combined with the spatial anisotropy of the GCR flux; and (b) $～$ 30 days due to solar impact of stable solar structures combined with the rotation of the Sun. From a superposed epoch analysis, and considering the geomagnetic effects, the mean diurnal amplitude is $～$ 0.08% and the maximum flux is observed in $～$ 15 h local time (LT) direction in the interplanetary space. In such a way, we determine the capability of Neurus to observe anisotropies and other interplanetary modulations on the GCR flux arriving at the Earth.     

Real-time ray-based channel generation and emulation for UAV communications

A real-time ray-based hardware emulator for Unmanned Aerial Vehicle (UAV) communication channels which suits for the Three-Dimensional (3D) dynamic scenarios and considers the movements of both terminals is developed in this paper. The time-variant channel parameters, i.e., ray delay, ray gain, and ray Doppler frequency are precalculated in the host by using the Ray Tracing (RT) method. Meanwhile, RT simulation dramatically increases the number of valid rays. To address the problem of resource limitation and huge computational burden in the implementation, an efficient ray coefficients generation method based on iteration is proposed and implemented. With the advantages of low cost and high flexibility, a Software Defined Radio (SDR) hardware platform is used to emulate the ray-based UAV channels by utilizing the compact architecture including the Time-Division (TD) scheme and Tapped-Delay Line (TDL) for channel convolution. Finally, hardware measurement results demonstrate that the properties of emulated channel, i.e., Power Delay Profile (PDP) and Doppler Power Spectrum Density (DPSD) consist with the simulated ones, which verifies the correctness of hardware implementation. The proposed channel emulator provides an efficient way for optimization, verification, and evaluation of UAV communication systems.