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.     

Estimations of the Seismic Pressure Noise on Mars Determined from Large Eddy Simulations and Demonstration of Pressure Decorrelation Techniques for the Insight Mission

The atmospheric pressure fluctuations on Mars induce an elastic response in the ground that creates a ground tilt, detectable as a seismic signal on the InSight seismometer SEIS. The seismic pressure noise is modeled using Large Eddy Simulations (LES) of the wind and surface pressure at the InSight landing site and a Green’s function ground deformation approach that is subsequently validated via a detailed comparison with two other methods: a spectral approach, and an approach based on Sorrells’ theory (Sorrells, Geophys. J. Int. 26:71–82, 1971; Sorrells et al., Nat. Phys. Sci. 229:14–16, 1971). The horizontal accelerations as a result of the ground tilt due to the LES turbulence-induced pressure fluctuations are found to be typically \(\sim 2 \mbox{--} 40~\mbox{nm}/\mbox{s}^{2}\) in amplitude, whereas the direct horizontal acceleration is two orders of magnitude smaller and is thus negligible in comparison. The vertical accelerations are found to be \(\sim 0.1\mbox{--}6~\mbox{nm}/\mbox{s}^{2}\) in amplitude. These are expected to be worst-case estimates for the seismic noise as we use a half-space approximation; the presence at some (shallow) depth of a harder layer would significantly reduce quasi-static displacement and tilt effects.We show that under calm conditions, a single-pressure measurement is representative of the large-scale pressure field (to a distance of several kilometers), particularly in the prevailing wind direction. However, during windy conditions, small-scale turbulence results in a reduced correlation between the pressure signals, and the single-pressure measurement becomes less representative of the pressure field. The correlation between the seismic signal and the pressure signal is found to be higher for the windiest period because the seismic pressure noise reflects the atmospheric structure close to the seismometer.In the same way that we reduce the atmospheric seismic signal by making use of a pressure sensor that is part of the InSight Auxiliary Payload Sensor Suite, we also the use the synthetic noise data obtained from the LES pressure field to demonstrate a decorrelation strategy. We show that our decorrelation approach is efficient, resulting in a reduction by a factor of \(\sim 5\) in the observed horizontal tilt noise (in the wind direction) and the vertical noise. This technique can, therefore, be used to remove the pressure signal from the seismic data obtained on Mars during the InSight mission.