Search NASA⌕ Search

SEARCH · Search NASA

Results for “atmospheric drag modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

High altitude atmospheric modeling

Five empirical models were compared with 13 data sets, including both atmospheric drag-based data and mass spectrometer data. The most recently published model, MSIS-86, was found to be the best model overall with an accuracy around 15 percent. The excellent overall agreement of the mass spectrometer-based MSIS models with the drag data, including both the older data from orbital decay and the newer accelerometer data, suggests that the absolute calibration of the (ensemble of) mass spectrometers and the assumed drag coefficient in the atomic oxygen regime are consistent to 5 percent. This study illustrates a number of reasons for the current accuracy limit such as calibration accuracy and unmodeled trends. Nevertheless, the largest variations in total density in the thermosphere are accounted for, to a very high degree, by existing models. The greatest potential for improvements is in areas where we still have insufficient data (like the lower thermosphere or exosphere), where there are disagreements in technique (such as the exosphere) which can be resolved, or wherever generally more accurate measurements become available.

Hedin, Alan E.↗

Atmospheric Models for Aerocapture Systems Studies

Aerocapture uses atmospheric drag to decelerate into captured orbit from interplanetary transfer orbit. This includes capture into Earth orbit from, for example, Lunar-return or Mars-return orbit. Eight Solar System destinations have sufficient atmosphere for aerocapture to be applicable - three of the rocky planets (Venus, Earth, and Mars), four gas giants (Jupiter, Saturn, Uranus, and Neptune), and Saturn's moon Titan. These destinations fall into two groups: (1) The rocky planets, which have warm surface temperatures (approx. 200 to 750 K) and rapid decrease of density with altitude, and (2) the gas giants and Titan, which have cold temperatures (approx. 70 to 170 K) at the surface or 1-bar pressure level, and slow rate of decrease of density with altitude. The height variation of average density with altitude above 1-bar pressure level for the gas giant planets is shown. The periapsis density required for aerocapture of spacecraft having typical values of ballistic coefficient (a measure of mass per unit cross-sectional area) is also shown. The aerocapture altitudes at the gas giants would typically range from approx. 150 to 300 km. Density profiles are compared for the rocky planets with those for Titan and Neptune. Aerocapture at the rocky planets would occur at heights of approx. 50 to 100 km. For comparison, typical density and altitudes for aerobraking operations (circularizing a highly elliptical capture orbit, using multiple atmospheric passes) are also indicated.

Justus, C. G.↗

Neutral Orbital Altitude Density Effects on the International Space Station

One of the design requirements of the International Space Station (ISS) is that each year accelerations of one micro-g cannot be exceeded at the ISS internal payload location for 6 periods of not less than 30 consecutive days. Although there are other causes, this study deals only with the accelerations caused by atmospheric drag. The critical ambient neutral density, computed using the Marshall Engineering Thermosphere Model, required to produce accelerations of one micro-g on the ISS, is estimated using an atmospheric drag acceleration equation. Results show that the design requirements may be difficult to meet during periods of extremely high solar activity; the planned reboost and altitude strategies for the ISS may have to be revised to allow for the uncertainty in the prediction of neutral atmospheric density within the 100-day period established for orbital decay before reboost.

Smith, O.E.↗

Contribution of the Lower Atmosphere to the Day-to-day Variation of Thermospheric Density

In this paper we carried out a numerical experiment using the Specified Dynamics mode of the Whole Atmosphere Community Climate Model with thermosphere and ionosphere eXtension (SD-WACCM-X). One SD-WACCM-X run was with realistic Kp and F10.7 and the other with constant Kp and F10.7. By comparing the day-to-day variability of thermosphere mass density at 300 km (low earth orbit, LEO) and 120 km (reentry level) in these two runs, we find that the density variation at 300 km is mainly driven by geomagnetic and solar forcing while at 120 km it is exclusively controlled by the lower atmosphere. At LEO altitudes, during solar minimum and geomagnetic quiet days, the impact from the lower atmosphere is much smaller than the effect of solar and geomagnetic variations but is not negligible (5–10% vs 20%).

Thermosphere Density↗

Recent improvements in our knowledge of neutral atmosphere structure from satellite drag measurements

Observational results on the density in the thermosphere and lower exosphere (i.e., within the altitude range from about 150 to about 1000 km) are discussed in this paper. Most observational results on total gas density were obtained from orbital drag and more recently also from in-situ drag analysis. The primary parameter measured is atmospheric density, with temperature as a secondary structural parameter deduced with the help of theory and/or atmospheric models. Both the merits and shortcomings of the drag analysis method are outlined in view of a comparison of temperature deduced from total density and kinetic gas temperature measured by incoherent scatter. Recent improvements of our knowledge of the known density variations are presented.

Roemer, M.↗

The Effect of Drag Model on Heatshield Recession due to Particle Impacts for Martian Spacecraft

A spacecraft entering the Martian atmosphere during a dust storm may experience recession to the heatshield due to dust particle impacts. Aerodynamic drag is the primary force that determines the trajectory of the dust particles through the shock layer. This paper examines the effect of particle drag model on the heatshield recession. Three particle drag models are assessed including two that are intended to be applicable over a wide range of particle flow conditions. Particle trajectories are computed in conditions measured during the 2007 major global dust storm. It was found that accounting for Knudsen number and compressibility effects made a large difference in the estimated particle impact velocity. The two drag models that were valid for transitional, compressible particle flow environments predicted only slightly-different amounts of heatshield recession due to dust particle impacts. A brief description of the effects of non-spherical particles on drag coefficients is provided.

Heatshield↗

The Effect of Drag Model on Heatshield Recession due to Particle Impacts for Martian Spacecraft

A spacecraft entering the Martian atmosphere during a dust storm may experience recession to the heatshield due to dust particle impacts. Aerodynamic drag is the primary force that determines the trajectory of the dust particles through the shock layer. This paper examines the effect of particle drag model on the heatshield recession. Three particle drag models are assessed including two that are intended to be applicable over a wide range of particle flow conditions. Particle trajectories are computed in conditions measured during the 2007 major global dust storm. It was found that accounting for Knudsen number and compressibility effects made a large difference in the estimated particle impact velocity. The two drag models that were valid for transitional, compressible particle flow environments predicted only slightly-different amounts of heatshield recession due to dust particle impacts. A brief description of the effects of non-spherical particles on drag coefficients is provided.

Heatshield↗

Techniques for the analysis of geodynamic effects using laser data

New orbit computation techniques have been developed to realize the full precision of laser ranging measurements from a single tracking station used to accurately determine the orbital inclination of a satellite. In order to evaluate earth and ocean tidal effects on the satellite and polar motion effects on the station latitude, improved computational techniques are described for perturbations significantly influencing the satellite's inclination, such as solar radiation pressure and geopotential resonance. By using the time independent value of maximum latitude reached by the satellite as the experimental variable, orbit tracking errors caused by imprecise modelling of the gravity field and atmospheric drag have been largely overcome and made possible long term analysis of osculating elements. With these techniques, quarter day spans of laser data have been employed to monitor the inclination of the satellite to the order 0.01 arcseconds precision over a period of seventeen months.

Dunn, P. J.↗

The JPL orbit determination software system

The Jet Propulsion Laboratory (JPL) orbit determination software system is described. Consisting of a set of computer programs, this system has been developed at JPL for the primary purpose of determining the flight path of deep-space mission spacecraft in NASA's planetary program. Secondary uses include science studies in celestial mechanics and extraterrestrial atmospherics. Attention is given to such features as the capability for processing 30 data types from up to 15 tracking stations, and for precise navigation to natural satellites of outer planets. In addition, the system contains models for such conditions as mascons, atmospheric drag, attitude control, and finite motor burns, to provide for extensive solve-on capabilities.

Ekelund, J. E.↗

Modelling of dispersal and deposition of impact glass spherules from the Cretaceous-Tertiary boundary deposit

The dispersal of glass spherules or tektites from a bolide impact with the Earth is modelled as ballistic trajectories in standard atmosphere. Ballistic dispersal of Cretaceous-Tertiary boundary impact glass spherules found in Haiti and Mimbral, Mexico requires a fireball radius in excess of 50 km but less than 100 km to account for the observed distribution. Glass spherules from 1 and up to 8 mm in diameter have been found at the KT boundary at Beloc in Haiti, at Mimbral, Mexico, and at DSDP Sites 536 and 540 in the Gulf of Mexico corresponding to paleodistances of 600 to 1000 km from the Chicxulub crater. In Haiti the basal and major glass-bearing unit at the KT boundary is attributed to fallout on basis of sedimentologic features. When compared with theoretical and observed dispersal of volcanic ejecta, the grain size versus distance relationship of the KT boundary tektite fallout is extreme, and rules out a volcanic fallout origin. At a comparable distance from source, the KT impact glass spherules are more than an order of mangitude coarser than ejecta of the largest known volcanic events. We model the dispersal of KT boundary impact glass spherules as ballistic ejecta from a fireball generated by the impact of a 10 km diameter bolide. Mass of ejecta in the fireball is taken as twice the bolide mass. Melt droplets are accelerated by gas flow in the fireball cloud, and leave the fireball on ballistic trajectories within the atmosphere, subject to drag, depending on angle of ejection and altitude. The model for ballistic dispersal is based on equations of motion, drag and ablation for silicate spheres in standard atmosphere.

Espindola, J. M.↗

Predicting Space Weather Effects on Close Approach Events

The NASA Robotic Conjunction Assessment Risk Analysis (CARA) team sends ephemeris data to the Joint Space Operations Center (JSpOC) for conjunction assessment screening against the JSpOC high accuracy catalog and then assesses risk posed to protected assets from predicted close approaches. Since most spacecraft supported by the CARA team are located in LEO orbits, atmospheric drag is the primary source of state estimate uncertainty. Drag magnitude and uncertainty is directly governed by atmospheric density and thus space weather. At present the actual effect of space weather on atmospheric density cannot be accurately predicted because most atmospheric density models are empirical in nature, which do not perform well in prediction. The Jacchia-Bowman-HASDM 2009 (JBH09) atmospheric density model used at the JSpOC employs a solar storm active compensation feature that predicts storm sizes and arrival times and thus the resulting neutral density alterations. With this feature, estimation errors can occur in either direction (i.e., over- or under-estimation of density and thus drag). Although the exact effect of a solar storm on atmospheric drag cannot be determined, one can explore the effects of JBH09 model error on conjuncting objects' trajectories to determine if a conjunction is likely to become riskier, less risky, or pass unaffected. The CARA team has constructed a Space Weather Trade-Space tool that systematically alters the drag situation for the conjuncting objects and recalculates the probability of collision for each case to determine the range of possible effects on the collision risk. In addition to a review of the theory and the particulars of the tool, the different types of observed output will be explained, along with statistics of their frequency.

Space Weather↗

Precision orbit determination for the TOPEX/Poseidon mission

Computation of precise orbits for the TOPEX/Poseidon (T/P) spacecraft is analyzed focusing on gravity field modeling, nonconservative force modeling, satellite tracking technologies, and orbit determination software. It was found that the radial orbit error budget for T/P allows 10 cm rms error due to gravity field mismodeling, 3 cm due to solid earth and ocean tides, and 6 cm due to radiative forces, and 3 cm due to atmospheric drag. It is concluded that the current models are capable of achieving the radial orbit error requirements.

Nerem, R. S.↗

Investigation of dynamic noise affecting geodynamics information in a tethered subsatellite

Measurement of the gradient of the gravitational acceleration from a satellite platform is likely to provide the next improvement in knowledge of the Earth's gravity field after the upcoming Geopotential Research Mission. Observations from the subsatellite of a tethered satellite system (TSS) would increase sensitivity and resolution due to the low altitude possible. However, the TSS is a dynamically noisy system and would be perturbed by atmospheric drag fluctuations. The dynamic noise is being modeled in order to evaluate the feasibility of TSS gradiometry and to design methods of abating the error caused by this noise. The demonstration flights of the TSS are to provide an opportunity to directly observe the dynamical environment and refine modeling techniques. Random vibration analysis as a technique for modeling the TSS under atmospheric perturbation was studied.

Gullahorn, G. E.↗

Techniques of orbital decay and long-term ephemeris prediction for satellites in earth orbit

In the special perturbation method, Cowell and variation-of-parameters formulations of the motion equations are implemented and numerically integrated. Variations in the orbital elements due to drag are computed using the 1970 Jacchia atmospheric density model, which includes the effects of semiannual variations, diurnal bulge, solar activity, and geomagnetic activity. In the general perturbation method, two-variable asymptotic series and automated manipulation capabilities are used to obtain analytical solutions to the variation-of-parameters equations. Solutions are obtained considering the effect of oblateness only and the combined effects of oblateness and drag. These solutions are then numerically evaluated by means of a FORTRAN program in which an updating scheme is used to maintain accurate epoch values of the elements. The atmospheric density function is approximated by a Fourier series in true anomaly, and the 1970 Jacchia model is used to periodically update the Fourier coefficients. The accuracy of both methods is demonstrated by comparing computed orbital elements to actual elements over time spans of up to 8 days for the special perturbation method and up to 356 days for the general perturbation method.

Barry, B. F.↗

An Investigation Into Transecting Satellites in Future Space Traffic Management Scenarios

The number of satellites expected to populate the near-Earth space environment is set to dramatically increase in the coming decade as new large constellations are approved and deployed. Current strategies for deploying new batches of these satellites often involve launching into an initial orbit, and then performing apogee raising maneuvers to reach a target altitude. Similarly, end-of-life planning for these constellation satellites can consist of de-orbit burns that lower perigee to permit disposal via re-entry. Both the raising and de-orbiting maneuvers can result in the individual satellites traveling in transecting orbits that have the potential to cross other spacecraft trajectories. While individual large constellations may be able to coexist in separate altitude and inclination bands, having thousands of satellites moving between these bands as new satellites are added and old satellites are removed could pose additional collision risks. Similar concerns have been raised regarding the impact that large numbers of university-class CubeSats might have in terms of their overall collision risk, especially as these satellites typically do not have propulsion systems for active maneuvering. To assess the impact that transecting satellites might have on future space traffic management strategies, this study explored a variety of realistic future scenarios using a high-fidelity simulation tool. The model can simulate the orbit of tens of thousands of resident space objects (RSOs) simultaneously, to include active satellites, debris, rocket bodies, or even future hypothetical satellite constellations, using a realistic force model that incorporates non-spherical gravity, atmospheric drag, and solar radiation pressure, as well as station-keeping. The simulation can be customized to accommodate different methods of calculating the probability of collision, as well as the process for determining probability ellipsoids and screening volumes. This makes it possible to replicate, and compare, different processes used by different spacecraft operators and space situational awareness (SSA) providers. As the model is run forward in time, each conjunction event is recorded, allowing for the analysis of statistics and meta-data related to these events, providing insight into the nature and frequency of potential collisions, such as whether are they active or passive objects, what size are the two objects, and who owns the objects (if known). This information makes it possible to characterize how changes to the status quo affect the number and type of conjunctions that occur, as well as the distributional effects on various types of satellite operators. To assess the general risks that transecting satellites might pose for hypothetical future space object environments, approximately 60,000 new large constellation satellites were considered, in addition to the existing catalog of approximately 7800 known resident space objects (RSOs), over a simulation period of one year. The results indicate that the future space environment will introduce a non-linear increase in conjunction events as the number of RSOs also increase. This will require adjustments to spacecraft fuel budgets in order to conduct the avoidance maneuvers necessary to minimize collision risk, both for existing and new satellites. This increase is due in large part to the higher density of RSOs and the overlap between some constellation orbits. Current catalog objects were shown to require three times more ∆V for collision avoidance (CA) maneuvers in the simulated future environment, and some constellation spacecraft were estimated to devote the majority of their annual ∆V to CA. The impact of small satellites was found to be proportional for the current space environment, and actually decreased in terms of percentage for the future scenario, suggesting that small satellites do not pose an outsized collision risk. Lastly, transecting satellites were found to contribute thousands of additional conjunctions outside of their operational orbit, and may require up to an additional 5% in CA maneuver fuel allocation.

conjunction assessment↗

A dynamical model of upper-atmospheric helium

This numerical study of the global distribution of helium in the terrestrial thermosphere is based on the extension of a three-dimensional single-fluid numerical model of the thermosphere previously developed by the authors to treat the dynamics of a minor gas imbedded in a background gas made up of N2, O2, and O. The Cira (1972) model of the upper atmosphere, based on satellite drag data, is used to specify the background gas density and temperature as functions of altitude, latitude, and local time for a given day of the year. Specific emphasis is placed on a study of the helium distribution at solstice, i.e., the 'winter helium bulge'. Effects of solar activity, eddy diffusion, and exospheric transport on the global distribution of helium are investigated. The model reproduces the general features of the winter helium bulge, including the variation of its amplitude with solar activity, as determined by comparison of data from several recent satellite-borne mass spectrometers.

Straus, J. M.↗

Asteroid Fragmentation Approaches for Modeling Atmospheric Energy Deposition

During asteroid entry, energy is deposited in the atmosphere through thermal ablation and momentum-loss due to aerodynamic drag. Analytic models of asteroid entry and breakup physics are used to compute the energy deposition, which can then be compared against measured light curves and used to estimate ground damage due to airburst events. This work assesses and compares energy deposition results from four existing approaches to asteroid breakup modeling, and presents a new model that combines key elements of those approaches. The existing approaches considered include a liquid drop or "pancake" model where the object is treated as a single deforming body, and a set of discrete fragment models where the object breaks progressively into individual fragments. The new model incorporates both independent fragments and aggregate debris clouds to represent a broader range of fragmentation behaviors and reproduce more detailed light curve features. All five models are used to estimate the energy deposition rate versus altitude for the Chelyabinsk meteor impact, and results are compared with an observationally derived energy deposition curve. Comparisons show that four of the five approaches are able to match the overall observed energy deposition profile, but the features of the combined model are needed to better replicate both the primary and secondary peaks of the Chelyabinsk curve.

Modeling↗

Modeling

A prediction of the future population of satellites, satellite fragments, and assorted spacecraft debris in Earth orbit can be reliably made only after three conditions are satisfied: (1) the size and spatial distributions of these Earth-orbiting objects are established at some present-day time; (2) the processes of orbital evolution, explosions, hypervelocity impact fragmentation, and atmospheric drag are understood; and (3) a reasonable traffic model for the future launch rate of Earth-orbiting objects is assumed. The theoretician will then take these three quantities as input data and will carry through the necessary mathematica and numerical analyses to project the present-day orbital population into the future.

Zook, H. A.↗