Fast solar radiation pressure modelling with ray tracing and multiple reflections

Abstract 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 % , 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  nm / s 2 and 0.4  nm / 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.

[1]  O. Montenbruck,et al.  Enhanced solar radiation pressure modeling for Galileo satellites , 2015, Journal of Geodesy.

[2]  A. S. Ganeshan,et al.  GNSS Satellite Geometry and Attitude Models , 2015 .

[3]  M. Ziebart.  Generalized Analytical Solar Radiation Pressure Modeling Algorithm for Spacecraft of Complex Shape , 2004 .

[4]  Peter Shirley,et al.  An Efficient and Robust Ray-Box Intersection Algorithm , 2005, J. Graph. Tools.

[5]  Vlastimil Havran,et al.  Incremental BVH construction for ray tracing , 2015, Comput. Graph..

[6]  Marek Ziebart,et al.  Sunlight Illumination Models for Spacecraft Surface Charging , 2017, IEEE Transactions on Plasma Science.

[7]  Z. Kang,et al.  Precise orbit determination for the GRACE mission using only GPS data , 2006 .

[8]  H.-P. Seidel,et al.  Realtime Ray Tracing on GPU with BVH-based Packet Traversal , 2007, 2007 IEEE Symposium on Interactive Ray Tracing.

[9]  R. Dach,et al.  Geocenter coordinates estimated from GNSS data as viewed by perturbation theory , 2013 .

[10]  Kai Hormann,et al.  The point in polygon problem for arbitrary polygons , 2001, Comput. Geom..

[11]  E. F. Nichols,et al.  The Pressure Due To Radiation , 2015 .

[12]  Michiel Otten,et al.  Sentinel-1A - First precise orbit determination results , 2017 .

[13]  Dinesh Manocha,et al.  Fast BVH Construction on GPUs , 2009, Comput. Graph. Forum.

[14]  Rodríguez Solano,et al.  Impact of non-conservative force modeling on GNSS satellite orbits and global solutions , 2014 .

[15]  E. F. Nichols,et al.  Pressure Due to Light and Heat Radiation , 1902 .

[16]  Peter Lebedew,et al.  Untersuchungen über die Druckkräfte des Lichtes , 1901 .

[17]  R. Dach,et al.  CODE’s new solar radiation pressure model for GNSS orbit determination , 2015, Journal of Geodesy.

[18]  M. Ziebart,et al.  Combined radiation pressure and thermal modelling of complex satellites: Algorithms and on-orbit tests , 2004 .

[19]  M. Zhong,et al.  Analysis of HY2A precise orbit determination using DORIS , 2015 .

[20]  F. Gini,et al.  Fully Dynamic Approach for GOCE Precise Orbit Determination , 2013 .

[21]  J. Ries,et al.  Precision Orbit Determination Standards for the Jason Series of Altimeter Missions , 2010 .

[22]  H. Fliegel,et al.  Global Positioning System Radiation Force Model for geodetic applications , 1992 .

[23]  Baocheng Zhang,et al.  A new analytical solar radiation pressure model for current BeiDou satellites: IGGBSPM , 2016, Scientific reports.

[24]  G. Beutler,et al.  A New Solar Radiation Pressure Model for GPS Satellites , 1999, GPS Solutions.

[25]  Oliver Montenbruck,et al.  Satellite Orbits: Models, Methods and Applications , 2000 .

[26]  Paul Cross,et al.  GPS Block IIR Non-Conservative Force Modeling: Computation and Implications , 2003 .

[27]  Jan P. Weiss,et al.  An evaluation of solar radiation pressure strategies for the GPS constellation , 2011 .

[28]  M. Ziebart.  High Precision Analytical Solar Radiation Pressure Modelling for GNSS Spacecraft , 2001 .