Search

Satellite drag modeling remains the largest source of uncertainty affecting space operations in low Earth orbit. The uncertainty stems from inaccurate models for mass density and drag coefficient. Drag coefficient modeling also impacts scientific knowledge on the physics and dynamics of the upper atmosphere through the estimation of high-fidelity mass density from measurements of acceleration on-board satellites. Efforts over the last decade have pushed drag coefficient modeling in the right direction, however, have resulted in multiple methods and tools. We provide a comprehensive review of the drag coefficient modeling methods and tools. Current scale differences between thermospheric data sets mostly originate from errors in the aerodynamic modeling, specifically in the modeling of the satellite outer surface geometry and the gas-surface interactions. Enhancing these models' accuracy is intrinsically connected to the satellite drag fidelity for science and operations. A team of invested scientists recently met under the community-driven International Space Weather Action Teams (ISWAT) initiative to discuss and consolidate the efforts towards a drag coefficient modeling baseline for consistency in science and operations. In this paper, we compare the available methods for drag coefficient modeling, their impact on the derived density estimates, and make recommendations for adoption of baseline methods for science and operations. Results show that the differences in derived densities estimates can reach tens of percent at altitudes above 4̃50 km during solar minimum conditions resulting mainly from differences in the modeling of gas-surface interactions. As a result, we conclude and recommend that robust uncertainty quantification be an integral part of any modeling efforts that employ the high-fidelity accelerometer derived density estimates. We also conclude and recommend that gas-surface interaction models that account for impact of altitude and solar variations be employed moving forward. Finally, we recommend future science missions to improve our understanding of gas-surface interactions and eventually the upper thermosphere variability. [ABSTRACT FROM AUTHOR]

Drag coefficient is a major source of uncertainty in calculating the aerodynamic forces on satellites in low Earth orbit. Closed-form solutions are available for simple geometries under the assumption of free molecular flow; however, most satellites have complex geometries, and a more sophisticated method of calculating the drag coefficient is needed. This work builds toward modeling physical drag coefficients using the direct simulation Monte Carlo method capable of accurately modeling flow shadowing and concave geometries. The direct simulation three-dimensional visual program and the direct simulation Monte Carlo analysis code are used to compare the effects of two separate gas-surface interaction models: diffuse reflection with incomplete accommodation and quasi-specular Cercignani-Lampis-Lord models. Results show that the two gas-surface interaction models compare well at altitudes below ~500 km during solar maximum conditions and below ~400 km during solar minimum conditions. The difference in drag coefficient of a sphere at ~800 km calculated using the two gas-surface interaction models is ~6% during solar maximum and increases to ~10% during solar minimum. The difference in drag coefficient of the GRACE satellite computed using the two gas-surface interaction models at ~500 km differs by ~15% during solar minimum conditions and by ~2-3 % during solar maximum conditions. [ABSTRACT FROM AUTHOR]

Abstract: Drag coefficient is a major source of uncertainty in predicting the orbit of a satellite in low Earth orbit (LEO). Computational methods like the Test Particle Monte Carlo (TPMC) and Direct Simulation Monte Carlo (DSMC) are important tools in accurately computing physical drag coefficients. However, the methods are computationally expensive and cannot be employed real time. Therefore, modeling of the physical drag coefficient is required. This work presents a technique of developing parameterized drag coefficients models using the DSMC method. The technique is validated by developing a model for the Gravity Recovery and Climate Experiment (GRACE) satellite. Results show that drag coefficients computed using the developed model for GRACE agree to within 1% with those computed using DSMC. [Copyright &y& Elsevier]

Ambitious satellite constellation projects by commercial entities and the ease of access to space in recent times have led to a dramatic proliferation of low-Earth space traffic. It jeopardizes space safety and long-term sustainability, necessitating better space domain awareness (SDA). Correct modeling of uncertainties in force models and orbital states, among other things, is an essential part of SDA. For objects in the low-Earth orbit (LEO) region, the uncertainty in the orbital dynamics mainly emanate from limited knowledge of the atmospheric drag-related parameters and variables. In this paper, which extends the work by Paul et al. (2021), we develop a feed-forward deep neural network model for the prediction of the satellite drag coefficient for the full range of satellite attitude (i.e., satellite pitch ∈ (- 90 ° , + 90 °) and satellite yaw ∈ (0 ° , + 360 °)). The model simultaneously predicts the mean and the standard deviation and is well-calibrated. We use numerically simulated physical drag coefficient data for training our neural network. The numerical simulations are carried out using the test particle Monte Carlo method using the diffuse reflection with incomplete accommodation gas-surface interaction model. Modeling is carried out for the well-known CHAllenging Minisatellite Payload (CHAMP) satellite. Finally, we use the Monte Carlo approach to propagate CHAMP over a three-day period under various modeling scenarios to investigate the distribution of radial, along-track, and cross-track orbital errors caused by drag coefficient uncertainty. The key takeaways of this paper are - (a) a constant drag coefficient cannot be used for reliable SDA purposes, and (b) stochastic machine learning models allow for the computation of drag coefficients in a timely manner while providing reliable uncertainty estimates. [ABSTRACT FROM AUTHOR]

For satellites in the Low Earth Orbit (LEO) region, the drag coefficient is a primary source of uncertainty for orbit determination and prediction. Researchers at the Los Alamos National Laboratory (LANL) have created the so-called Response Surface Modeling (RSM) toolkit to provide the community with a resource for simulating and modeling satellite drag coefficients for satellites with complex geometries (modeled using triangulated facets) in the free molecular flow (FMF) regime. The toolkit fits an interpolation surface using non-parametric Gaussian Process Regression (GPR) over drag coefficient data computed using the numerical Test Particle Monte Carlo (TPMC) method. The fitted response surface provides a substantial computational benefit over numerical approaches for calculating drag coefficients. In this work, the RSM toolkit is further developed into a versatile software with extended capabilities. The capabilities are specifically expanded to include uncertainty quantification and adaptation for automatic development of regression models for satellites with non-stationary components (e.g. rotating solar panels). Furthermore, the toolkit uses Python 3.x and C programming languages to provide an open source software package with a OSI approved GPL license. To assist the end user, the new RSM toolkit has been developed to have a user-friendly installation process and is provided with extensive documentation. The analysis of two different conceptual satellites is performed during this work: a simple cube and a CubeSat consisting of a simple cube body with 2 rotating solar panels. During the creation of the regression model for each satellite for different atmospheric species, it is found that the cube's minimum Root Mean Squared Error (RMSE) is 0.00211 and the maximum RMSE is 0.00350. The CubeSat has a minimum RMSE of 0.00304 and the maximum is 0.00498. These results are overall conducive of a well performing regression model. [ABSTRACT FROM AUTHOR]