SEARCH · Engineering Papers
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.
With the recent completion of the Global Positioning System (GPS) constellation and lower cost spaceborne receivers GPS usage is expanding rapidly in the area of space flight projects. Several of these are discussed, such as precise orbit determination, Earth gravity model improvement, and other global geodesy enhancements; high resolution 2D/3D ionospheric imaging; and atmospheric limb sounding.
We have developed a new method for solving low-thrust fuel-optimal orbit transfer problems in the vicinity of a large body (planet or asteroid), considering a high-fidelity spherical harmonic gravity model. The algorithm is formulated via the indirect optimization method, leading to a two-point boundary value problem (TPBVP). We make use of a hyperbolic tangent smoothing law for performing continuation on the thrust magnitude to reduce the sharpness of the control switches in early iterations and thus promote convergence. The TPBVP is solved using the method of particular solutions (MPS) shooting method and Picard-Chebyshev numerical integration. Application of Picard-Chebyshev integration affords an avenue for increased efficiency that is not available with step-by-step integrators. We demonstrate that computing the particular solutions with only a low-fidelity force model greatly increases the efficiency of the algorithm while ultimately achieving near machine precision accuracy. A salient feature of the MPS is that it is parallelizable, and thus further speedups are available. It is also shown that, for near Earth orbits and over a small number of en-route revolutions around the Earth, only the zonal perturbation terms are required in the costate equations to obtain a solution that is accurate to machine precision and optimal to engineering precision. The proposed framework can be used for trajectory design around small asteroids and also for orbit debris rendezvous and removal tasks.
A system of Lunar relay satellites has been proposed to address communication and navigation needs and ensure robustness for the variety of upcoming robotic and human exploration missions to the Moon. The relays are envisioned to estimate self position and time knowledge onboard, allowing the system to provide in-situ navigation services to missions in the Lunar and cis-Lunar environment. The quality and accuracy of those services are highly dependent on the navigation performance of the relay itself. To assess the Lunar relay navigation performance, a series of orbit determination (OD) Monte Carlo (MC) simulations are run using Lunar gravity modeling up to a degree and order of 250 and a variety of onboard clocks and measurement types including weak-signal GNSS, Ground Network (GN) pseudorange (PR) and Doppler, and optical navigation (OpNav) center-finding (CF). The estimated trajectories produced by these Lunar relay MC simulations, along with the associated errors, and transmitted navigation reference signal parameters, are used to evaluate the expected navigation performance of a user on a descent trajectory to the Lunar surface. The scenario features a lander system performing onboard navigation relying on one-way range and Doppler measurements from reference signals emitted by the Lunar relay. This paper can be used as a reference in determining the onboard clock and measurement types necessary to obtain acceptable navigation performance for the Lunar relay, and as a baseline for Lunar lander navigation performance using accurate measurements from a relay reference signal.
A system of Lunar relay satellites has been proposed to address communication and navigation needs and ensure robustness for the variety of upcoming robotic and human exploration missions to the Moon. The relays are envisioned to estimate self-position and time knowledge onboard, allowing the system to provide in-situ navigation services to missions in the Lunar and cis-Lunar environment. The quality and accuracy of those services are highly dependent on the navigation performance of the relay itself. To assess the Lunar relay navigation performance, a series of orbit determination (OD) Monte Carlo (MC) simulations are run using Lunar gravity modeling up to a degree and order of 250 and a variety of onboard clocks and measurement types including weak-signal GNSS, Ground Network (GN) pseudorange (PR) and Doppler, and optical navigation (OpNav) center-finding (CF). The estimated trajectories produced by these Lunar relay MC simulations, along with the associated errors, and transmitted navigation reference signal parameters, are used to evaluate the expected navigation performance of a user on a descent trajectory to the Lunar surface. The scenario features a lander system performing onboard navigation relying on one-way range and Doppler measurements from reference signals emitted by the Lunar relay. This paper can be used as a reference in determining the onboard clock and measurement types necessary to obtain acceptable navigation performance for the Lunar relay, and as a baseline for Lunar lander navigation performance using accurate measurements from a relay reference signal.
For the realization of the 2020 International Terrestrial Reference Frame (ITRF2020), the International DORIS Service delivered to the International Earth Rotation and Reference Systems Service (IERS) a set of 1456 weekly solution files from 1993.0 to 2021.0 including station coordinates and Earth orientation parameters (EOPs). The data come from fourteen DORIS satellites: TOPEX/Poseidon, SPOT-2, SPOT-3, SPOT-4, SPOT-5, Envisat, Jason-1, Jason-2, Cryosat-2, Saral, HY-2A, Jason-3, Sentinel-3A and Sentinel-3B. In their processing, the four analysis centers which contributed to the DORIS combined solution used the latest time variable gravity models, the new mean pole and diurnal-subdiurnal tidal EOP models recommended by IERS. In addition, all the analysis centers included in their processing precise SPOT-5 solar panel angle values and quaternions for, at least, the Jason satellites. Furthermore, a new Alcatel phase center variation model was implemented for the ITRF2020 processing. The main objective of this study is to present the combination process and to analyze the impact of the new modeling on the performance of the new combined solution. Comparisons with the IDS contribution to ITRF2014 show that i) the application of the new phase center variations for the Alcatel DORIS ground antennas in the data processing combined with the gradual replacement over time of the Alcatel by Starec antennas implies a scale drift from 1993.0 to 2002.5 and ii) thanks to a better modeling of the surface forces on the satellites, the new combined solution shows smaller annual and 118-day signals in the geocenter. A new DORIS terrestrial reference frame was computed to evaluate the intrinsic quality of the new combined solution. That evaluation shows that over almost the full time span the intrinsic IDS scale values lie in a range of ±5 mm. After mid-2008, the new DORIS reference frame has an internal position consistency in North-East-Up better than 7.5 mm.
Numerous vehicles will visit Gateway throughout the expected mission lifetime, and much like Gateway, these vehicles will operate in a Near Rectilinear Halo Orbit (NRHO). Short-term loitering strategies are explored for these vehicles. The objective of this analysis is to illuminate concepts that enable safe loitering, while considering the operational limitations to the design process. To highlight relative motion, a linear sensitivity analysis is applied. The linear approach offers an initial guess strategy for generating loiter architectures in the high-fidelity gravity model. Methods of constructing transfers between loitering regimes are explored.
Starlette Satellite Laser Ranging (SLR) data were used, along with several other satellite data sets, for the solution of a preliminary gravity field model for TOPEX, PTGF1. A further improvement in the earth gravity model was accomplished using data collected by 12 satellites to solve another preliminary gravity model for TOPEX, designated PTGF2. The solution for the Earth Rotation Parameter (ERP) was derived from the analysis of SLR data to Starlette during the MERIT Campaign. Starlette orbits in 1976 and 1983 were analyzed for the mapping of the tidal response of the earth. Publications and conference presentations pertinent to research are listed.
The venting of cylindrical containers partially filled with initially saturated liquids was conducted under zero gravity conditions and compared with an analytical model which determined the effect of interfacial mass transfer on the ullage pressure response during venting. A model is proposed to improve the estimation of the interfacial mass transfer. Duhammel's superposition integral is incorporated in this analysis to approximate the transient temperature response of the interface, treating the liquid as a semiinfinite solid with conduction heat transfer. This approach to estimating interfacial mass transfer gives improved response when compared to previous models. The model still predicts a pressure decrease greater than those in the experiments reported.
The objective is to predict the pressure response of a saturated liquid-vapor system when undergoing a venting or depressurization process in zero gravity at low vent rates. An experimental investigation of the venting of cylindrical containers partially filled with initially saturated liquids was previously conducted under zero-gravity conditions and compared with an analytical model which incorporated the effect of interfacial mass transfer on the ullage pressure response during venting. A new model is presented to improve the estimation of the interfacial mass transfer. Duhammel's superposition integral is incorporated to approximate the transient temperature response of the interface, treating the liquid as a semi-infinite solid with conduction heat transfer. Account is also taken of the condensation taking place within the bulk of a saturated vapor as isentropic expansion takes place. Computational results are presented for the venting of R-11 from a given vessel and initial state for five different venting rates over a period of three seconds, and compared to prior NASA experiments. An improvement in the prediction of the final pressure takes place, but is still considerably below the measurements.
The general theory of MST radar observations of gravity wave spectra is developed. This effort extends the previous results to include anisotropy and Doppler effects for the spectra, as well as the consequences for the multibeam configuration. The relationships between the observed one- or two-dimensional spectra for the line-of-sight velocity in the gravity wave spectra are derived. Expressions for cross spectra, as well as covariances between velocities observed on different beams, are computed. Using these results, studies are carried out to show how model parameters for gravity wave spectra can be estimated from the observed quantities. Model parameters include the variance, power law indices, anisotropy parameters, Doppler parameters, mean scale sizes, etc. Cases with different numbers of beams are investigated.
This paper presents a method to model the external gravity field and to estimate the internal density variation of a small-body. We first discuss the modeling problem, where we assume the polyhedral shape and internal density distribution are given, and model the body interior using finite elements definitions, such as cubes and spheres. The gravitational attractions computed from these approaches are compared with the true uniform-density polyhedral attraction and the level of accuracies are presented. We then discuss the inverse problem where we assume the body shape, radiometric measurements, and a priori density constraints are given, and estimate the internal density variation by estimating the density of each finite element. The result shows that the accuracy of the estimated density variation can be significantly improved depending on the orbit altitude, finite-element resolution, and measurement accuracy.
Flat-Earth modeling is a desirable alternative to the complex spherical-Earth modeling process. These methods were compared using 2 1/2 dimensional flat-earth and spherical modeling to compute gravity and scalar magnetic anomalies along profiles perpendicular to the strike of variably dimensioned rectangular prisms at altitudes of 150, 300, and 450 km. Comparison was achieved with percent error computations (spherical-flat/spherical) at critical anomaly points. At the peak gravity anomaly value, errors are less than + or - 5% for all prisms. At 1/2 and 1/10 of the peak, errors are generally less than 10% and 40% respectively, increasing to these values with longer and wider prisms at higher altitudes. For magnetics, the errors at critical anomaly points are less than -10% for all prisms, attaining these magnitudes with longer and wider prisms at higher altitudes. In general, in both gravity and magnetic modeling, errors increase greatly for prisms wider than 500 km, although gravity modeling is more sensitive than magnetic modeling to spherical-Earth effects. Preliminary modeling of both satellite gravity and magnetic anomalies using flat-Earth assumptions is justified considering the errors caused by uncertainties in isolating anomalies.
Sea surface height variability measured by TOPEX is analyzed in the tropical Pacific Ocean by way of assimilation into a wind-driven, reduced-gravity, shallow water model using an approximate Kalman filter and smoother. The analysis results in an optimal fit of the dynamic model to the observations, providing it dynamically consistent interpolation of sea level and estimation of the circulation. Nearly 80% of the expected signal variance is accounted for by the model within 20 deg of the equator, and estimation uncertainty is substantially reduced by the voluminous observation. Notable features resolved by the analysis include seasonal changes associated with the North Equatorial Countercurrent and equatorial Kelvin and Rossby waves. Significant discrepancies are also found between the estimate and TOPEX measurements, especially near the eastern boundary. Improvements in the estimate made by the assimilation are validated by comparisons with independent tide gauge and current meter observations. The employed filter and smoother are based on approximately computed estimation error covariance matrices, utilizing a spatial transformation and an symptotic approximation. The analysis demonstrates the practical utility of a quasi-optimal filter and smoother.
A Gravity Field Workshop was convened to review the actions which could be taken prior to a GRAVSAT mission to improve the Earth's gravity field model. This review focused on the potential improvements in the Earth's gravity field which could be obtained using the current satellite and surface gravity data base. In particular, actions to improve the quality of the gravity field determination through refined measurement corrections, selected data augmentation and a more accurate reprocessing of the data were considered. In addition, recommendations were formulated which define actions which NASA should take to develop the necessary theoretical and computation techniques for gravity model determination and to use these approaches to improve the accuracy of the Earth's gravity model.
The dominant error source in estimating the orbital position of a satellite from ground based tracking data is the modeling of the Earth's gravity field. The resulting orbit error due to gravity field model errors are predominantly long wavelength in nature. This results in an orbit error signature that is strongly correlated over distances on the size of ocean basins. Anderle and Hoskin (1977) have shown that the orbit error along a given ground track also is correlated to some degree with the orbit error along adjacent ground tracks. This cross track correlation is verified here and is found to be significant out to nearly 1000 kilometers in the case of TOPEX/POSEIDON when using the GEM-T1 gravity model. Finally, it was determined that even the orbit error at points where ascending and descending ground traces cross is somewhat correlated. The implication of these various correlations is that the orbit error due to gravity error is geographically correlated. Such correlations have direct implications when using altimetry to recover oceanographic signals.
Viewgraphs on the Middeck Active Control Experiment (MACE): Gravity and Suspension Effects are presented. Topics covered include: gravity and suspension influences; gravity and suspension effects on structure; gravity effects on sensors and actuators; modeling of gravity and suspension effects on structure; geometric stiffness theory and initial static deformation calculation; modeling gravity's effect on accelerometers and PMA's; application to MACE; MACE EM configuration study; and MACE DM configuration study.
What is the threshold gravity (minimum gravity level) required for the nominal functioning of the human system? What dosage is required? Do human cell lines behave differently in microgravity in response to an external stimulus? The critical need for such a gravity simulator is emphasized by recent experiments on human epithelial cells and lymphocytes on the Space Shuttle clearly showing that cell growth and function are markedly different from those observed terrestrially. Those differences are also dramatic between cells grown in space and those in Rotating Wall Vessels (RWV), or NASA bioreactor often used to simulate microgravity, indicating that although morphological growth patterns (three dimensional growth) can be successfully simulated using RWVs, cell function performance is not reproduced - a critical difference. If cell function is dramatically affected by gravity off-loading, then cell response to stimuli such as radiation, stress, etc. can be very different from terrestrial cell lines. Yet, we have no good gravity simulator for use in study of these phenomena. This represents a profound shortcoming for countermeasures research. We postulate that we can use magnetic levitation of cells and tissue, through the use of strong magnetic fields and field gradients, as a terrestrial microgravity model to study human cells. Specific objectives of the research are: 1. To develop a tried, tested and benchmarked terrestrial microgravity model for cell culture studies; 2. Gravity threshold determination; 3. Dosage (magnitude and duration) of g-level required for nominal functioning of cells; 4. Comparisons of magnetic levitation model to other models such as RWV, hind limb suspension, etc. and 5. Cellular response to reduced gravity levels of Moon and Mars. The paper will discuss experiments md modeling work to date in support of this project.
Numerical modeling has been performed of the fluid dynamics in a prototypical physical vapor transport crystal growing situation. Cases with and without gravity have been computed. Dependence of the flows upon the dimensionless parameters aspect ratio and Peclet, Rayleigh, and Schmidt numbers is demonstrated to a greater extent than in previous works. Most notably, it is shown that the effects of thermally-induced buoyant convection upon the mass flux on the growth interface crucially depend upon the temperature boundary conditions on the sidewall (e.g., whether adiabatic or of a fixed profile, and in the latter case the results depend upon the shape of the profile assumed).