Search NASA⌕ Search

SEARCH · Search NASA

Results for “ODE”

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 109 records · Page 6

Diagonally Implicit Runge-Kutta Methods for Ordinary Differential Equations. A Review

A review of diagonally implicit Runge-Kutta (DIRK) methods applied to rst-order ordinary di erential equations (ODEs) is undertaken. The goal of this review is to summarize the characteristics, assess the potential, and then design several nearly optimal, general purpose, DIRK-type methods. Over 20 important aspects of DIRKtype methods are reviewed. A design study is then conducted on DIRK-type methods having from two to seven implicit stages. From this, 15 schemes are selected for general purpose application. Testing of the 15 chosen methods is done on three singular perturbation problems. Based on the review of method characteristics, these methods focus on having a stage order of two, sti accuracy, L-stability, high quality embedded and dense-output methods, small magnitudes of the algebraic stability matrix eigenvalues, small values of aii, and small or vanishing values of the internal stability function for large eigenvalues of the Jacobian. Among the 15 new methods, ESDIRK4(3)6L[2]SA is recommended as a good default method for solving sti problems at moderate error tolerances.

Kennedy, Christopher A.↗

Machine Learning for the Prediction of Local Asteroid Damages

Risk assessment studies of local asteroid hazards traditionally simulate the physics of meteors with engineering models tailored to analyze tens-of-millions of scenarios. However, these simplified approaches still need to solve time-dependent ODEs to model the entry process and the resulting ground damage. With a computational cost of O(0.01 CPU.s) per scenario, simulating these large numbers of potential entry conditions in risk assessment studies can take several days on local computers. To improve computational efficiency, we propose in this paper an orthogonal approach based on machine learning models to predict the size of damaged areas given a list of entry parameters. We train 5 machine learning methods and compare the predictions to the outputs of the PAIR model, first only with primitive entry condition variables, and then with more advanced features. We find that complex models like neural networks are well-suited to estimate blast hazards, while simpler linear models can accurately assess thermal damage. For both types of hazards, the radii of damaged areas can be predicted with around 10% average errors and a coefficient of determination (R2) of 0.99. The CPU time is decreased by a factor O(10 3 ) compared to the PAIR model, which enables the simulation of millions of scenarios in minutes, on a local computer. We then use the same machine learning approaches for a classification task where the models are trained to predict if an asteroid will produce a given level of damage. Results show that complex models like the gradient boosting classifier and the neural network can perform this task with 98% accuracy. Beyond surrogate models, we finally incorporate the machine learning algorithms to the state-of-the-art Shapley sensitivity analysis and present a ranking of the entry parameters based on their contributions to ground damages.

SMD↗

Machine Learning Approach for Aircraft Performance Model Parameter Estimation for Trajectory Prediction Applications

Inaccurate prediction of aircraft trajectory by ground-based decision support tools (DST) is a major concern in air traffic management (ATM). Aircraft trajectory prediction tools rely on a simplified point-mass aircraft performance model (APM) to make their predictions. Even though the performance coefficients and weight of an aircraft are a vital part of the APM’s predictions and accuracy, these coefficients are proprietary in nature and therefore, unavailable to DSTs. Current ATM research focuses on improving the estimate of some APM parameters by freezing all other coefficients. This simplified approach introduces unwanted sources of bias and negatively impacts the accuracy of the performance model. In this paper, we apply machine learning (ML) techniques for the simultaneous prediction of three key APM parameters (two drag coefficients and the initial aircraft weight). To accomplish this, we employ an ordinary differential equation (ODE) fitting approach to generate optimized APM parameter labels customized to each individual flight record. Subsequently, we train ML models to capture the relationship between the historical data and the optimized APM parameters. Two different ML model solutions are applied and APM coefficients are predicted for unseen flights. The results indicate that the ML models are able to capture the relationship between APM parameters and flight-related features with good accuracy.

trajectory prediction, machine learning, aircraft ↗

Derivatives Across Adaptive-Step Integration Accurately and Efficiently

The numerical integration of ordinary differential equations (ODEs) is common in trajectory optimization. Previous studies have shown that the use of adaptive-step Runge-Kutta integration schemes can introduce errors in the derivatives due to non-differentiable behavior within the integration scheme itself. The authors previously demonstrated the use of an implicit adaptation to the Runge-Kutta methods to mitigate this issue, but at a significant computational cost. In this work, the authors demonstrate an explicit adaptation to the Runge-Kutta methods to address the issue that does not incur a significant computational cost.

Kaushik Ponnapalli↗

Minimum fuel trajectories for a low-thrust power-limited mission to the moon and to Lagrange points L4 and L5

Minimum fuel trajectories from a low earth parking orbit to Lagrange points L4 or L5 and to the moon are obtained for a low-thrust limited-power spacecraft, with thrust acceleration levels of the order of 0.001 G. The procedure to find a trajectory to the libration point starts from an analytical description of a slightly elliptical spiral, given by Breakwell and Rauch. The earth moon trajectory is found by matching an earth spiral to a moon spiral on the sphere of influence. Earth oblateness effect is considered.

Breakwell, John V.↗

Minimum fuel trajectory for the aerospace-plane

An overall trajectory for a single-stage-to-orbit vehicle with an initial weight of 234 tons is calculated, and four different propulsion models including turbojet, ramjet, scramjet, and rocket are considered. First, the atmospheric flight in the thicker atmosphere is discussed with emphasis on trajectory optimization, optimization problem, aerodynamic problem, propulsion model, and initial conditions. The performance of turbojet and ramjet-scramjet engines is analyzed; and then the flight to orbit is assessed from the optimization point of view. It is shown that roll modulation saves little during the trajectory, and the combined application of airbreathing propulsion and aerodynamic lift is suggested.

Breakwell, John V.↗

Minimum fuel lunar trajectories for a low-thrust power-limited spacecraft

Minimum-fuel trajectories from a low earth parking orbit to a low moon orbit are obtained for a low-thrust power-limited spacecraft with thrust acceleration levels of the order of 0.001 G. The trajectories are found by matching an earth spiral to a moon spiral at some intermediate distance. Results are given for the planar case and for the three-dimensional case where the moon orbit is polar.

Breakwell, John V.↗

Minimum fuel lunar trajectories for a low-thrust power-limited spacecraft

Minimum fuel trajectories from a low Earth parking orbit to a low Moon orbit are obtained for a low-thrust power-limited spacecraft with thrust acceleration levels of the order of 10(exp -3) G. The trajectories are found by matching an Earth spiral to a Moon spirit at some intermediate distance. Results are given for the planar case and for the three dimensional case where the Moon orbit is polar.

Golan, Oded M.↗

Communication by Fire (and Smoke) Signals in the Kingdom of Judah

This paper examines the use of ancient fire and smoke signals for communication in the Kingdom of Judah. Historical and biblical references are cited that discuss this communication system. The current physical and political landscape of Israel precludes testing of hypotheses using traditional techniques. The use of a GIS is enlisted to overcome these obstacles and predict line-of-sight patterns that are conducive for a fire signal communication system. Final demonstration of this predictive model will incorporate state-of-the-art technology and in-field data acquisition to provide the sufficient accuracy that is required for proof-of-concept. This research will provide insight into the technical capability of the ancient Israelites for communication across a mountainous, desert environment.

Borowski, Oded↗

Valley Network Morphology and Topographic Gradients on Mars

Data returned from the Mars Orbiter Laser Altimeter allows construction of a high precision digital elevation model. Quantitative investigations into the geomorphic properties of drainage features, similar to ones carried out on Earth, are now possible Additional information is contained in the original extended abstract.

Aharonson, Oded↗

Insights into the Earliest History of Mars: A New Synthesis

Motivated by the latest data from Mars Global Surveyor and from recent analyses of Martian meteorites, we offer a new synthesis of the relative timing of major events in the early geological history of Mars together with associated uncertainties. Additional information is contained in the original extended abstract.

Solomon, Sean C.↗

Localized Gravity/Topography Admittance and Correlation Spectra on Mars: Implications for Regional and Global Evolution

[i] From gravity and topography data collected by the Mars Global Surveyor spacecraft we calculate gravity/topography admittances and correlations in the spectral domain and compare them to those predicted from models of lithospheric flexure. On the basis of these comparisons we estimate the thickness of the Martian elastic lithosphere (T(sub e)) required to support the observed topographic load since the time of loading. We convert T(sub e) to estimates of heat flux and thermal gradient in the lithosphere through a consideration of the response of an elastic/plastic shell. In regions of high topography on Mars (e.g., the Tharsis rise and associated shield volcanoes), the mass-sheet (small-amplitude) approximation for the calculation of gravity from topography is inadequate. A correction that accounts for finite-amplitude topography tends to increase the amplitude of the predicted gravity signal at spacecraft altitudes. Proper implementation of this correction requires the use of radii from the center of mass (collectively known as the planetary shape ) in lieu of topography referenced to a gravitational equipotential. Anomalously dense surface layers or buried excess masses are not required to explain the observed admittances for the Tharsis Montes or Olympus Mons volcanoes when this correction is applied. Derived T, values generally decrease with increasing age of the lithospheric load, in a manner consistent with a rapid decline of mantle heat flux during the Noachian and more modest rates of decline during subsequent epochs.

McGovern, Patrick J.↗

Why are There So Few Magnetic Anomalies in Martian Lowlands and Basins?

The discovery of large areas of strongly magnetized crust on Mars [1,2] provides important information on the timing of both crustal and deeper interior processes on that planet. Critical to an understanding of that timing, as well as to the processes that have contributed to the magnetization and demagnetization of crustal materials, is the geographical distribution of magnetic anomalies discernible from spacecraft orbit. The paucity of resolved magnetic anomalies in the northern lowlands and within and surrounding the best-preserved major impact basins has been noted since the crustal field was first globally mapped [1], but no straightforward explanation of that full pattern has yet been offered. Here we suggest that ancient hydrothermal alteration of magnetic carriers in Martian lowlands and basins may have contributed to the magnetization distribution observed today.

Solomon, Sean C.↗

Depth and Distribution of CO2 Snow on Mars

The dynamic role of volatiles on the surface of Mars has been a subject of longstanding interest. In the pre-Viking era, much of the debate was necessarily addressed by theoretical considerations. A particularly influential treatment by Leighton and Murray put forth a simple model relying on solar energy balance, and led to the conclusion that the most prominent volatile exchanging with the atmosphere over seasonal cycles is carbon dioxide. Their model suggested that due to this exchange, atmospheric CO2 partial pressure is regulated by polar ice. While current thinking attributes a larger role to H2O ice than did the occasional thin polar coating this model predicted, the CO2 cycle appears to be essentially correct. There are a number of observational constraints on the seasonal exchange of surface volatiles with the atmosphere. The growth and retreat of polar CO2 frost is visible from Earth-based telescopes and from spacecraft in Mars orbit, both at visible wavelengths and in thermal IR properties of the surface. Recently, variations in Gamma ray and neutron fluxes have also been used to infer integrated changes in CO2 mass on the surface. Measurements made by Viking's Mars Atmospheric Water Detector experiment were sensitive to atmospheric H2O vapor abundance. Surface condensates and their transient nature were detected by the Viking landers. The study here is motivated by recent data collected by the Mars Global Surveyor, affording the opportunity to not only detect the lateral distribution of volatiles, but also to constrain the variable volumes of the reservoirs. We elaborate on a technique first employed by Smith et al. By examining averages of a large number of topographic measurements collected by the Mars Orbiter Laser Altimeter (MOLA), that study showed that the zonal pattern of deposition and sublimation of CO2 can be determined. In their first approach, reference surfaces were fit to all measurements in narrow latitude annuli, and the time dependent variations about those mean surfaces were examined. In their second approach, height measurements from pairs of tracks that cross on the surface were interpolated and differenced, forming a set of crossover residuals. These residuals were then examined as a function of time and latitude. The initial studies averaged over longitude to maximize signal and minimize noise in order to isolate the expected small signal. In this follow-up study we now attempt to extract the elevation change pattern also as a function of longitude, and we focus on the crossover approach.

Aharonson, Oded↗

Spatial and Temporal Patterns of Tidal Dissipation in Synchronous Satellites

Tidal heating is an important energy source for several solar system bodies, and there is a wide-spread perception that the pattern of surface heat flow is diagnostic of internal structure. We wish to clarify that situation. Our analysis depends upon two important assumptions: First, that heat transport is dominated by conduction. Second, that the body can be modeled by a sequence of spherically symmetric layers, each with a linear visco-elastic rheology. Under these assumptions, surface heat flow patterns in tidally dominated satellites will reflect radially integrated dissipation patterns. For synchronously rotating satellites with zero obliquity, this pattern depends quite strongly on orbital eccentricity but relatively little on purely radial variations in internal structure. The total amount of heat generated within the body does depend sensitively on internal structure, but the spatial pattern is rather insensitive to structure, especially at low orbital eccentricities.

Bills, Bruce G.↗

Constraints on the Volatile Distribution Within Shackleton Crater at the Lunar South Pole

Shackleton crater is nearly coincident with the Moon's south pole. Its interior receives almost no direct sunlight and is a perennial cold trap, making Shackleton a promising candidate location in which to seek sequestered volatiles. However, previous orbital and Earth-based radar mapping and orbital optical imaging have yielded conflicting interpretations about the existence of volatiles. Here we present observations from the Lunar Orbiter Laser Altimeter on board the Lunar Reconnaissance Orbiter, revealing Shackleton to be an ancient, unusually well-preserved simple crater whose interior walls are fresher than its floor and rim. Shackleton floor deposits are nearly the same age as the rim, suggesting that little floor deposition has occurred since the crater formed more than three billion years ago. At a wavelength of 1,064 nanometres, the floor of Shackleton is brighter than the surrounding terrain and the interiors of nearby craters, but not as bright as the interior walls. The combined observations are explicable primarily by downslope movement of regolith on the walls exposing fresher underlying material. The relatively brighter crater floor is most simply explained by decreased space weathering due to shadowing, but a one-micrometre-thick layer containing about 20 per cent surficial ice is an alternative possibility.

Zuber, Maria T.↗

Summary of the Results from the Lunar Orbiter Laser Altimeter after Seven Years in Lunar Orbit

In June 2009 the Lunar Reconnaissance Orbiter (LRO) spacecraft was launched to the Moon. The payload consists of 7 science instruments selected to characterize sites for future robotic and human missions. Among them, the Lunar Orbiter Laser Altimeter (LOLA) was designed to obtain altimetry, surface roughness, and reflectance measurements. The primary phase of lunar exploration lasted one year, following a 3-month commissioning phase. On completion of its exploration objectives, the LRO mission transitioned to a science mission. After 7 years in lunar orbit, the LOLA instrument continues to map the lunar surface. The LOLA dataset is one of the foundational datasets acquired by the various LRO instruments. LOLA provided a high-accuracy global geodetic reference frame to which past, present and future lunar observations can be referenced. It also obtained high-resolution and accurate global topography that were used to determine regions in permanent shadow at the lunar poles. LOLA further contributed to the study of polar volatiles through its unique measurement of surface brightness at zero phase, which revealed anomalies in several polar craters that may indicate the presence of water ice. In this paper, we describe the many LOLA accomplishments to date and its contribution to lunar and planetary science.

LRO↗