Comparing Physical Drag Coefficients Computed Using Different Gas-Surface Interaction Models

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,mostsatelliteshavecomplexgeometries,andamoresophisticatedmethodofcalculatingthedragcoefficient 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 threedimensional 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–Lordmodels.Resultsshowthatthetwogas–surfaceinteractionmodelscomparewellataltitudes below ∼500 km during solar maximum conditions and below ∼400 km during solar minimum conditions. The differenceindragcoefficientofasphereat ∼800 kmcalculated usingthetwogas–surfaceinteractionmodels 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.