Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical integration”

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 991 records · Page 55

Hyper-X Mach 10 Trajectory Reconstruction

This paper discusses the formulation and development of a trajectory reconstruction tool for the NASA X-43A/Hyper-X high speed research vehicle, and its implementation for the reconstruction and analysis of flight test data. Extended Kalman filtering techniques are employed to reconstruct the trajectory of the vehicle, based upon numerical integration of inertial measurement data along with redundant measurements of the vehicle state. The equations of motion are formulated in order to include the effects of several systematic error sources, whose values may also be estimated by the filtering routines. Additionally, smoothing algorithms have been implemented in which the final value of the state (or an augmented state that includes other systematic error parameters to be estimated) and covariance are propagated back to the initial time to generate the best-estimated trajectory, based upon all available data. The methods are applied to the problem of reconstructing the trajectory of the Hyper-X vehicle from data obtained during the Mach 10 test flight, which occurred on November 16th 2004.

Karlgaard, Christopher D.↗

The orbits of the outer Jovian satellites

This article describes a fit of numerically integrated orbits for the eight outer Jovian satellites to Earth-based astrometric observations.

Atmospheric observation Jupiter ephermerides↗

Upper bounds to error probabilities of coded systems over AWGN and fading channels

A family of upper bounds to error probabilities of coded systems on the additive white Gaussian noise channel was recently proposed by Divsalar. Their calculation depends only on the weight spectrum of the code words. We first elaborate upon these bounds to show how they can be further tightened by using numerical integration instead of a Chernoff bound, and by reducing the number of code words to be included in the bound. Next, we extend them to finding channels.

coded systems↗

The Orbit of Pheobe from Earthbased and Voyager Observations

This article presents the results of a fit of a numerically integrated orbit for the Saturnian satellite Pheobe to Earthbased astrometric observations (from 1904 to 1996) and imaging data acquired by the Voyager 2 spacecraft during its encounter with Saturn. The primary results are the epoch state vector used in the integration and a set of mean elements which provide a geometrical representation of the orbit. We also assess the quality of the fit and the accuracy of the orbit.

Ephemirides planet ephemirides satellite ephemirid↗

A Search for Low V Earth-to-Moon Trajectories

A search for low V Earth-to-Moon trajectories has been initiated. Numerical integration of the equations of motion from the circular restricted three-body problem has resulted in the computation of a trajectory that saves more that 100 m/s over a Hohmann transfer, although the flight time is almost ten months.

trajectory↗

A Continuous Square Root in Formation Filter-Swoother with Discrete Data Update

A differential equation for the square root information matrix is derived and adapted to the problems of filtering and smoothing. The resulting continuous square root information filter (SRIF) performs the mapping of state and process noise by numerical integration of the SRIF matrix and admits data via a discrete least square update.

SRIF mapping↗

Predicting Close Earth Approaches of Asteroids and Comets

The motions of all known Earth approaching asteroids and comets with reasonably secure orbits have been numerically integrated forward in time to A.D. 2200. Special care was taken to use the best available initial conditions including orbits based upon radar data.

Comets vaporizing ices Earth asteroid 2340 Hathor↗

Realistic Covariance Prediction For the Earth Science Constellations

Routine satellite operations for the Earth Science Constellations (ESC) include collision risk assessment between members of the constellations and other orbiting space objects. One component of the risk assessment process is computing the collision probability between two space objects. The collision probability is computed via Monte Carlo techniques as well as numerically integrating relative probability density functions. Each algorithm takes as inputs state vector and state vector uncertainty information for both objects. The state vector uncertainty information is expressed in terms of a covariance matrix. The collision probability computation is only as good as the inputs. Therefore, to obtain a collision calculation that is a useful decision-making metric, realistic covariance matrices must be used as inputs to the calculation. This paper describes the process used by NASA Goddard's Earth Science Mission Operations Project to generate realistic covariance predictions for three of the ESC satellites: Aqua, Aura, and Terra

Duncan, Matthew↗

Realistic Covariance Prediction for the Earth Science Constellation

Routine satellite operations for the Earth Science Constellation (ESC) include collision risk assessment between members of the constellation and other orbiting space objects. One component of the risk assessment process is computing the collision probability between two space objects. The collision probability is computed using Monte Carlo techniques as well as by numerically integrating relative state probability density functions. Each algorithm takes as inputs state vector and state vector uncertainty information for both objects. The state vector uncertainty information is expressed in terms of a covariance matrix. The collision probability computation is only as good as the inputs. Therefore, to obtain a collision calculation that is a useful decision-making metric, realistic covariance matrices must be used as inputs to the calculation. This paper describes the process used by the NASA/Goddard Space Flight Center's Earth Science Mission Operations Project to generate realistic covariance predictions for three of the Earth Science Constellation satellites: Aqua, Aura and Terra.

Duncan, Matthew↗

Low-thrust Orbit Transfers Using Candidate Lyapunov Functions with a Mechanism for Coasting

We consider low-thrust orbit transfers around a central body, where specified changes are sought in orbit elements except true anomaly. The desired changes in the remaining five elements can be arbitrarily large. Candidate Lyapunov functions are created based on analytic expressions for maximum rates of change of the orbit elements and the desired changes in the elements. These functions may be thought of as proximity quotients because they provide R measure of the proximity to the target orbit. The direction of thrust needed for steepest descent to the target orbit is also available analytically. The thrust is shutoff if the effectivity of the thrust at the current location on the osculating orbit is below some threshhold value. Thus, the equations of motion can be numerically integrated to obtain quickly and simply a transfer to the target orbit. A series of transfers can be easily computed to assess the trade-off between propellant mass and flight time. Preliminary comparisons to optimal solutions show that the method, while sub-optimal, performs well.

Candidate Lyapunoc functions↗

A Semi-Empirical Model for Forecasting Relativistic Electrons at Geostationary Orbit

We developed a new prediction model for forecasting relativistic (>2MeV) electrons, which provides a VERY HIGH correlation between predicted and actually measured electron fluxes at geostationary orbit. This model implies the multi-step particle acceleration and is based on numerical integrating two linked continuity equations for primarily accelerated particles and relativistic electrons. The model includes a source and losses, and used solar wind data as only input parameters. We used the coupling function which is a best-fit combination of solar wind/Interplanetary Magnetic Field parameters, responsible for the generation of geomagnetic activity, as a source. The loss function was derived from experimental data. We tested the model for four year period 2004-2007. The correlation coefficient between predicted and actual values of the electron fluxes for whole four year period as well as for each of these years is about 0.9. The high and stable correlation between the computed and actual electron fluxes shows that the reliable forecasting these electrons at geostationary orbit is possible. The correlation coefficient between predicted and actual electron fluxes is stable and incredibly high.

Lyatsky, Wladislaw↗

Relativistic Electrons at Geostationary Orbit: Modeling Results

We developed a new prediction model for forecasting relativistic (greater than 2MeV) electrons, which provides a VERY HIGH correlation between predicted and actually measured electron fluxes at geostationary orbit. This model implies the multi-step particle acceleration and is based on numerical integrating two linked continuity equations for primarily accelerated particles and relativistic electrons. The model includes a source and losses, and used solar wind data as only input parameters. We used the coupling function which is a best-fit combination of solar wind/interplanetary magnetic field parameters, responsible for the generation of geomagnetic activity, as a source. The loss function was derived from experimental data. We tested the model for four year period 2004-2007. The correlation coefficient between predicted and actual values of the electron fluxes for whole four year period as well as for each of these years is stable and incredibly high (about 0.9). The high and stable correlation between the computed and actual electron fluxes shows that the reliable forecasting these electrons at geostationary orbit is possible.

Khazanov, George V.↗

Prediction Model for Relativistic Electrons at Geostationary Orbit

We developed a new prediction model for forecasting relativistic (greater than 2MeV) electrons, which provides a VERY HIGH correlation between predicted and actually measured electron fluxes at geostationary orbit. This model implies the multi-step particle acceleration and is based on numerical integrating two linked continuity equations for primarily accelerated particles and relativistic electrons. The model includes a source and losses, and used solar wind data as only input parameters. We used the coupling function which is a best-fit combination of solar wind/interplanetary magnetic field parameters, responsible for the generation of geomagnetic activity, as a source. The loss function was derived from experimental data. We tested the model for four year period 2004-2007. The correlation coefficient between predicted and actual values of the electron fluxes for whole four year period as well as for each of these years is stable and incredibly high (about 0.9). The high and stable correlation between the computed and actual electron fluxes shows that the reliable forecasting these electrons at geostationary orbit is possible.

Khazanov, George V.↗

Newton-Euler Dynamic Equations of Motion for a Multi-body Spacecraft

The Magnetospheric MultiScale (MMS) mission employs a formation of spinning spacecraft with several flexible appendages and thruster-based control. To understand the complex dynamic interaction of thruster actuation, appendage motion, and spin dynamics, each spacecraft is modeled as a tree of rigid bodies connected by spherical or gimballed joints. The method presented facilitates assembling by inspection the exact, nonlinear dynamic equations of motion for a multibody spacecraft suitable for solution by numerical integration. The building block equations are derived by applying Newton's and Euler's equations of motion to an "element" consisting of two bodies and one joint (spherical and gimballed joints are considered separately). Patterns in the "mass" and L'force" matrices guide assembly by inspection of a general N-body tree-topology system. Straightforward linear algebra operations are employed to eliminate extraneous constraint equations, resulting in a minimum-dimension system of equations to solve. This method thus combines a straightforward, easily-extendable, easily-mechanized formulation with an efficient computer implementation.

Stoneking, Eric↗

Variations of Strahl Properties with Fast and Slow Solar Wind

The interplanetary solar wind electron velocity distribution function generally shows three different populations. Two of the components, the core and halo, have been the most intensively analyzed and modeled populations using different theoretical models. The third component, the strahl, is usually seen at higher energies, is confined in pitch-angle, is highly field-aligned and skew. This population has been more difficult to identify and to model in the solar wind. In this work we make use of the high angular, energy and time resolution and three-dimensional data of the Cluster/PEACE electron spectrometer to identify and analyze this component in the ambient solar wind during high and slow speed solar wind. The moment density and fluid velocity have been computed by a semi-numerical integration method. The variations of solar wind density and drift velocity with the general build solar wind speed could provide some insight into the source, origin, and evolution of the strahl.

Figueroa-Vinas, Adolfo↗

MHD Modeling of the Solar Wind with Turbulence Transport and Heating

We have developed a magnetohydrodynamic model that describes the global axisymmetric steady-state structure of the solar wind near solar minimum with account for transport of small-scale turbulence associated heating. The Reynolds-averaged mass, momentum, induction, and energy equations for the large-scale solar wind flow are solved simultaneously with the turbulence transport equations in the region from 0.3 to 100 AU. The large-scale equations include subgrid-scale terms due to turbulence and the turbulence (small-scale) equations describe the effects of transport and (phenomenologically) dissipation of the MHD turbulence based on a few statistical parameters (turbulence energy, normalized cross-helicity, and correlation scale). The coupled set of equations is integrated numerically for a source dipole field on the Sun by a time-relaxation method in the corotating frame of reference. We present results on the plasma, magnetic field, and turbulence distributions throughout the heliosphere and on the role of the turbulence in the large-scale structure and temperature distribution in the solar wind.

Goldstein, M. L.↗

Thermal Transport Model for Heat Sink Design

A document discusses the development of a finite element model for describing thermal transport through microcalorimeter arrays in order to assist in heat-sinking design. A fabricated multi-absorber transition edge sensor (PoST) was designed in order to reduce device wiring density by a factor of four. The finite element model consists of breaking the microcalorimeter array into separate elements, including the transition edge sensor (TES) and the silicon substrate on which the sensor is deposited. Each element is then broken up into subelements, whose surface area subtends 10 10 microns. The heat capacity per unit temperature, thermal conductance, and thermal diffusivity of each subelement are the model inputs, as are the temperatures of each subelement. Numerical integration using the Finite in Time Centered in Space algorithm of the thermal diffusion equation is then performed in order to obtain a temporal evolution of the subelement temperature. Thermal transport across interfaces is modeled using a thermal boundary resistance obtained using the acoustic mismatch model. The document concludes with a discussion of the PoST fabrication. PoSTs are novel because they enable incident x-ray position sensitivity with good energy resolution and low wiring density.

Chervenak, James A.↗

Multiple Exhaust Nozzle Effects on J-2X Gas Generator Outlet Impedance

The current test setup of the J-2X gas generator system uses a multiple nozzle configuration to exhaust hot gases to drive the propellant supply turbines. Combustion stability assessment of this gas generator design requires knowledge of the impedance effects the multiple nozzle configuration creates on the combustion chamber acoustic modes. Parallel work between NASA and Sierra Engineering is presented, showing two methods used to calculate the effective end impedance resulting from multiple nozzle configurations. The NASA method is a simple estimate of the effective impedance using the long wavelength approximation. Sierra Engineering has developed a more robust numerical integration method implemented in ROCCID to accommodate for multiple nozzles. Analysis using both methods are compared to J-2X gas generator test data collected over the past year.

Kenny, R. Jeremy↗