Search NASA⌕ Search

SEARCH · Search NASA

Results for “order reduction”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Jet Noise Prediction for Chevron Nozzle Flows with Wall-Modeled Large-Eddy Simulation

This paper presents results from ongoing research on jet noise prediction with wall-modeled large eddy simulations (WMLES) performed with the LAVA computational framework. In particular the focus of this study is on mixing enhancements from a single stream chevron nozzle at Reynolds number of 1×10(exp 6). Although the concept of chevron nozzles to reduce jet noise is not new, our understanding of its impact on the overall noise is still not well understood. As a first step towards predicting noise reduction due to mixing enhancement concepts from first principles with WMLES, we simulate the noise generated by a single stream chevron nozzle SMC001 as well as its equivalent axisymmetric round jet nozzle SMC000. Detailed comparisons are made with a dedicated experiment conducted at NASA’s Glenn Research Center and good agreement was achieved. Two different approaches to introduce a turbulent boundary layer were compared but show no major impact on the results. This is especially important given future work were multi-stream nozzles are considered and extended costs of resolving the internal BL would have a bigger cost impact. A permeable Ffowcs Williams Hawkings (FWH) surface enclosing the jet is used to predict far-field noise from the simulated flow-field and excellent comparison to microphone array measurements is achieved within the resolved frequency band. Sensitivity of far-field noise predictions to grid resolution is systematically documented. Near-field comparisons to PIV data shows great agreement for both velocity and normal stresses, however a systematic TKE overshoot at the nozzle exit is seen in the shear-layer. The paper shows a cost reduction of an order of magnitude compared to an earlier study of this configuration due to algorithmic and software improvements and demonstrates that WMLES can be used as a cost-competitive approach for jet noise predictions.

CST↗
NASA NTRSDate not supplied

How Well Satellite Remote Sensing Can Inform Surface Level Ozone Production Sensitivity to Precursor Trace Gas Emissions?

Surface-level ozone (O3) pollution mitigation strategies rely on a clear understanding of the atmospheric chemical processes involving O3, NOx (nitric oxide (NO) + nitrogen dioxide (NO2)) and volatile organic compounds (VOCs). A robust spatiotemporal classification of O3 photochemical regimes with relative chemical sensitivity of local O3 formation to emission reductions of NOx (NOx-limited regime) versus VOCs (radical-limited regime) is required to design anthropogenic emission reduction policies in order to improve surface air quality. Numerous studies found that the ratio of satellite-retrieved Vertical Column Densities (VCDs) of formaldehyde (HCHO) to NO2 can diagnose O3 production sensitivity. However, uncertainty still remains about the accuracy of such estimates using satellite-based retrievals. This study conducts an investigation to identify how well satellite data can inform surface-level O3 sensitivity by validating satellite-based estimates against aircraft-based retrievals taken during the Long Island Sound Tropospheric Ozone Study (LISTOS-2018) and Ozone Water-Land Environmental Transition Study (OWLETS-2) field campaigns. For this study, we use HCHO and NO2 retrievals from the Ozone Monitoring Instrument (OMI) onboard the Aura satellite, and high spatiotemporal resolution TROPOspheric Monitoring Instrument (TROPOMI) onboard the Sentinel-5 Precursor (S5P) satellite. The Community Multiscale Air Quality (CMAQ) modelling system is used to provide trace gas vertical profiles to recalculate air mass factor for deriving consistent satellite-based VCDs. To demonstrate the ability of satellites to diagnose O3 production regimes, we validate satellite-derived O3 sensitivity regimes, along with the OMI and TROPOMI products applied, using HCHO and NO2 retrievals from the aircraft-based instruments: Geostationary Trace gas and Aerosol Sensor Optimization (GeoTASO) spectrometer and Geostationary Coastal and Air Pollution Events Airborne Simulator (GCAS) instruments deployed during the LISTOS-2018 and OWLETS-2 campaigns. The results of this work will advance our fundamental knowledge of the spatiotemporal evolution of the complex O3-NOx-VOC chemistry, and will advance our understanding about the quality of remote-sensing, suborbital, and satellite data for accurately observing surface-level O3 production sensitivity.

Satellite↗
NASA NTRSDate not supplied

The Stratospheric Aerosol and Gas Experiment (SAGE) IV Pathfinder

As evidenced by quadrennial holding of this symposium and publishing of the Scientific Assessment of Ozone Depletion, the continued monitoring of and reporting on the stratospheric ozone layer is paramount to determining the efficacy of the Montreal Protocol. Measurements from satellites provided data on the initial decline of ozone in the late 1970s and early 1980s that supported the adoption of the Montreal Protocol, with current observations hinting at recovery. Adequate determination of that recovery requires continuous and, in the case of multiple instruments, overlapping data records. However, most current satellite systems are well beyond their expected lifetimes and are large and expensive. A new measurement paradigm is needed to enable cost-effective, sustainable measurements of atmospheric ozone into the coming decades when ozone is expected to recover. The Stratospheric Aerosol and Gas Experiment (SAGE) IV is a solar occultation imager capable of measuring ozone, aerosol, and other trace gas species with the same quality as previous SAGE instruments (including SAGE III/ISS). Taking advantage of recent technological advancements, a single sensorcraft occupies a 6U CubeSat form-factor providing an order of magnitude reduction in cost over traditional missions. SAGE IV is an example of an innovative and affordable mission that can sustain a crucial science measurement well into the future.

Robert Damadeo↗
NASA NTRS2021

Expected Thrust Fraction: Resilient Trajectory Design Applied to the Earth Return Orbiter

Electric propulsion is an enabling technology for NASA’s Mars Sample Return mission that affords an order of magnitude reduction in fuel. Unfortunately, elec- tric propulsion’s low-thrust capabilities require long thrusting arcs that make space- craft susceptible to disturbances resulting in missed thrust events. One method for designing trajectories resilient to such missed thrust events is expected thrust frac- tion, which embeds their stochastic nature into a time-varying duty cycle. This paper investigates the application of expected thrust fraction to the Earth Return Orbiter’s outbound Earth-to-Mars trajectory and compares it with trajectories de- signed using traditional methods to mitigate missed thrust events. Through the use of expected thrust fraction, a trajectory with a baseline 65.4% success rate is im- proved to an 86.6% success rate at the cost of only a 12 kg reduction in delivered mass. Additionally, when expected thrust fraction is used in conjunction with a terminal coast, a 96.0% success rate is achieved, which is higher than any single mitigation technique.

Laipert, Frank↗
NASA NTRS2023

Sidelobe Modeling and Mitigation for A Three Mirror Anastigmat Cosmic Microwave Background Telescope

Telescopes measuring cosmic microwave background (CMB) polarization on large angular scales require exquisite control of systematic errors to ensure the fidelity of the cosmological results. In particular, far-sidelobe contamination from wide angle scattering is a potentially prominent source of systematic error for large aperture microwave telescopes. Here we describe and demonstrate a ray-tracing-based modeling technique to predict far sidelobes for a three mirror anastigmat telescope designed to observe the CMB from the South Pole. Those sidelobes are produced by light scattered in the receiver optics subsequently interacting with the walls of the surrounding telescope enclosure. After comparing simulated sidelobe maps and angular power spectra for different enclosure wall treatments, we propose a highly scattering surface that would provide more than an order of magnitude reduction in the degree-scale far-sidelobe contrast compared to a typical reflective surface. We conclude by discussing the fabrication of a prototype scattering wall panel and presenting measurements of its angular scattering profile.

Ian Gullett↗
NASA NTRSDate not supplied

AstroDOME: Standardized, Mass-Produceable, Small NASA Astrophysics Observatories for Big Missions

NASA’s Goddard Space Flight Center is investigating the use of small, commercial off the shelf (COTS) spacecraft to enable breakthrough astrophysics research at radically reduced cost over previous generations of missions. This paper will focus on the Mission Systems perspective of this investigation. It will cover the spacecraft to instrument interface, and lessons learned in early concept payload development. Low-cost COTS spacecraft have proliferated in the last ~5 years. They represent an order of magnitude reduction in the cost of spacecraft bus development and acquisition. They are enabled by the advent of low-cost smallsat launch options in the immediately preceding years, and are seeded by recent distributed and disaggregated missions in other sectors such as telecommunications megaconstellations (e.g. Starlink, OneWeb, Project Kuiper) and ongoing Space Development Agency (SDA) developments. Our study was tasked with developing a standardized, mass-produceable, astrophysics payload which is reconfigurable to different wavelength bands while maintaining a single optics bench design, and a standardized mechanical, electrical, data, and thermal interface to the spacecraft bus, and which may be reconfigured to multiple COTS spacecraft offerings. In this paper we will discuss lessons learned from our case study. We will also present identified keys to mission success for performing cutting edge astrophysics observations while taking advantage of economies of scale with a disaggregated, standardized architecture.

Matthew Marcus↗
DOE OSTI2024

Promoting the regulatory acceptance of combined ion and neutron irradiation for material degradation in nuclear reactors

The Advanced Materials and Manufacturing Technologies (AMMT) program within the Department of Energy (DOE) Office of Nuclear Energy has developed its current recommendation for promoting the use of combined ion irradiation and neutron irradiation for the accelerated qualification of materials to be deployed in nuclear reactors. This plan is intended to provide a collaborative path forward that can be adopted by academia, national laboratories, and industry, and has been developed with input from the regulatory research arm of the U.S. Nuclear Regulatory Commission (NRC). To deploy new materials or materials manufactured with new technologies, such as additive manufacturing, materials must be evaluated for reactor-induced degradation from the combination of harsh temperatures, corrosive environments, and radiation fields. However, rapid deployment of materials necessitates accelerated testing methods rather than relying on years of neutron irradiation in a material test reactor. Ion irradiation has demonstrated success in reproducing material microstructure and select property evolution resulting from neutron irradiation with three to four orders of magnitude reduction in time and cost, making it an ideal candidate for accelerated irradiation testing. This presentation provides context governing both the scientific and regulatory aspects of the proposed goal. The discussion is aimed at a broad audience including researchers from industry, national laboratories, and academia. The recommended path forward is presented as a conceptual framework of specific steps. In brief, the strategy entails developing an integrated ion and neutron irradiation test plan for the material property of interest based on the fundamental tenet of the linkage of microstructure and properties in materials. Physics-based modeling interprets ion irradiation data and predicts neutron irradiation microstructure and properties with uncertainty bounds. The first round of testing is sufficient for an initial licensing application using a risk-informed approach, while a minimum required neutron irradiation test plan reduces cost and time requirements. A surveillance program with witness specimens in-reactor provides additional data over time to improve model predictions to higher damage levels and further reduce uncertainty bounds, which can be used for license extensions or longer lifetimes in new license applications.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Build-To-Replace Strategy to Reduce O&M Costs of Advanced Reactors

This paper targets the goal of reducing fixed operations and maintenance (O&M) costs to $2/MWh for advanced reactor (AR) designs, an order of magnitude reduction from current nuclear fleets. Such goal will be severely constrained if the underlying assumptions of O&M approaches and practices are not questioned and reexamined, especially today when changes can be implemented effectively and efficiently. Achieving a 90% reduction in AR O&M costs requires a completely new way of thinking – a paradigm shift – not through incremental, technology-focused approaches alone. Here, we address this challenge by evaluating the impact of moving to shorter design life for major structures, systems and components (SSCs) and shorter, more predictable refurbishment cycles as modeled by the commercial airline industry. This paper provides a brief overview of this different mindset to O&M applied to ARs: the build-to-replace approach. We provide a direct example of analysis of a build-to-replace scenario by identifying and evaluating scenarios for reduced system and component lifetimes and associated replacement and refurbishment schedules to evaluate impacts on O&M costs and other lifecycle elements such as SSC reliability.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Modeling of a Four-Stage Linear Ionization Cooling Channel for a Muon Collider in G4Beamline

A previous study of an eight-stage rectilinear ionization cooling channel in the ICOOL software demonstrated a five-order-of-magnitude reduction in a muon beam’s 6D emittance. In this study, we look to compare the ways ICOOL and Muons, Inc.’s g4Beamline software model ionization cooling by comparing their modeling of the first four stages to this optimized cooling channel constructed in ICOOL. We begin by identifying the parameters used to construct the optimized ionization cooling channel in ICOOL. We then reconstruct this beam in g4Beamline with identical parameter specifications and simulate the cooling of an identical input beam. Finally, we compare the two simulations based on their beam transmission, longitudinal emittance, and transverse emittance along the channel length. Through this process, we demonstrate that G4Beamline accurately reproduces transverse cooling results but predicts systematically different longitudinal emittance evolution while maintaining similar overall cooling performance, reproducing a 97.9% reduction in 6D emittance over four stages.

Keeler, Dominic [Purdue U., West Lafayette] (ORCID↗
NASA NTRS1990

Order (n) DISCOS for multibody systems with gear reduction

Recent developments in O(n) algorithms (where n is the number of bodies in the system) and parallel processing have drastically reduced the computer time needed to simulate systems involving many bodies. This paper presents a gear-reduction model for the O(n) version of DISCOS, a standard software package for simulation and analysis of flexible multibody systems. The gear-reduction model allows the accurate modeling of harmonic drives, which are commonly used in robot joints. The formulation has been implemented and validated with known results. The gear model can also be used for gear-train, rack-and-pinion, and screw joints.

Chun, Hon M.↗
NASA NTRS2008

Electromagnetically Clean Solar Arrays

The term 'electromagnetically clean solar array' ('EMCSA') refers to a panel that contains a planar array of solar photovoltaic cells and that, in comparison with a functionally equivalent solar-array panel of a type heretofore used on spacecraft, (1) exhibits less electromagnetic interferences to and from other nearby electrical and electronic equipment and (2) can be manufactured at lower cost. The reduction of electromagnetic interferences is effected through a combination of (1) electrically conductive, electrically grounded shielding and (2) reduction of areas of current loops (in order to reduce magnetic moments). The reduction of cost is effected by designing the array to be fabricated as a more nearly unitary structure, using fewer components and fewer process steps. Although EMCSAs were conceived primarily for use on spacecraft they are also potentially advantageous for terrestrial applications in which there are requirements to limit electromagnetic interference. In a conventional solar panel of the type meant to be supplanted by an EMCSA panel, the wiring is normally located on the back side, separated from the cells, thereby giving rise to current loops having significant areas and, consequently, significant magnetic moments. Current-loop geometries are chosen in an effort to balance opposing magnetic moments to limit far-0field magnetic interactions, but the relatively large distances separating current loops makes full cancellation of magnetic fields problematic. The panel is assembled from bare photovoltaic cells by means of multiple sensitive process steps that contribute significantly to cost, especially if electomagnetic cleanliness is desired. The steps include applying a cover glass and electrical-interconnect-cell (CIC) sub-assemble, connecting the CIC subassemblies into strings of series-connected cells, laying down and adhesively bonding the strings onto a panel structure that has been made in a separate multi-step process, and mounting the wiring on the back of the panel. Each step increases the potential for occurrence of latent defects, loss of process control, and attrition of components. An EMCSA panel includes an integral cover made from a transparent material. The silicone cover supplants the individual cover glasses on the cells and serves as an additional unitary structural support that offers the advantage, relative to glass, of the robust, forgiving nature of the silcone material. The cover contains pockets that hold the solar cells in place during the lamination process. The cover is coated with indium tin oxide to make its surface electrically conductive, so that it serves as a contiguous, electrically grounded shield over the entire panel surface. The cells are mounted in proximity to metallic printed wiring. The painted-wiring layer comprises metal-film traces on a sheet of Kapton (or equivalent) polyimide. The traces include contact pads on one side of the sheet for interconnecting the cells. Return leads are on the opposite side of the sheet, positioned to form the return currents substantially as mirror images of, and in proximity to, the cell sheet currents, thereby minimizing magnetic moments. The printed-wiring arrangement mimics the back-wiring arrangement of conventional solar arrays, but the current-loop areas and the resulting magnetic moments are much smaller because the return-current paths are much closer to the solar-cell sheet currents. The contact pads are prepared with solder fo electrical and mechanical bonding to the cells. The pocketed cover/shield, the solar cells, the printed-wiring layer, an electrical bonding agent, a mechanical-bonding agent, a composite structural front-side face sheet, an aluminum honeycomb core, and a composite back-side face sheet are all assembled, then contact pads are soldered to the cells and the agents are cured in a single lamination process.

Stem, Theodore G.↗
NASA NTRS1983

Applications of ethylene vinyl acetate as an encapsulation material for terrestrial photovoltaic modules

Terrestrial photovoltaic modules must undergo substantial reductions in cost in order to become economically attractive as practical devices for large scale production of electricity. Part of the cost reductions must be realized by the encapsulation materials that are used to package, protect, and support the solar cells, electrical interconnects, and other ancillary components. As many of the encapsulation materials are polymeric, cost reductions necessitate the use of low cost polymers. The performance and status of ethylene vinyl acetate, a low cost polymer that is being investigated as an encapsulation material for terrestrial photovoltaic modules, are described.

Cuddihy, E. F.↗
DOE OSTI2024

Medium-Range Structural Order as the Driver of Activated Dynamics and Complexity Reduction in Glass-Forming Liquids

Here, we analyze in depth the Elastically Collective Nonlinear Langevin Equation theory of activated dynamics in metastable liquids to establish that the predicted inter-relationships between the alpha relaxation time, local cage and collective elastic barriers, dynamic localization length, and shear modulus are causally related within the theory to the medium range order (MRO) static correlation length. The latter grows exponentially with density for metastable hard sphere fluids and as a nonuniversal inverse power law with temperature for supercooled liquids under isobaric conditions. The physical origin of predicted connections between the alpha time and other metrics of cage order and the thermodynamic inverse dimensionless compressibility is fully established. It is discovered that although kinetic constraints from the real space first coordination shell are important for the alpha time, they are of secondary importance compared to the consequences of the more universal MRO correlations in both the modestly and deeply metastable regimes. This understanding sheds new light on the theoretical basis for, and prior successes of, the predictive mapping of chemically complex thermal liquids to effective hard sphere fluids based on matching their dimensionless compressibilities, a scheme we call “complexity reduction”. In essence, the latter is equivalent to the physical requirement that the thermal liquid MRO correlation equals that of its effective hard sphere analog. The mapping alone is shown to provide a remarkable level of quantitative predictive power for the glass transition temperature T g of 21 molecular and polymer liquids. Predictions for the chemically specific absolute magnitude and growth with cooling of the MRO correlation length are obtained and lie in the window of 2–6 nm at T g . Dynamic heterogeneity, elastic facilitation, and beyond pair structure issues are briefly discussed. Future opportunities to theoretically analyze the equilibrated deep glass regime are outlined.

cancer↗
NASA NTRS1986

Turbulent boundary-layer wall pressure fluctuations downstream of a tandem LEBU

An initial measurement of turbulent wall pressure fluctuations downstream of tandem large eddy breakup devices (LEBUs) indicates a significant reduction by comparison with the reference case of a flat plate; the average magnitude of the reduction is 12.5 percent. Peak reduction is at 7-8 kHz, and is of the order of the C(f) reduction due to the tandem LEBUs. These data indicate a secondary benefit derivable from LEBUs, in addition to their skin friction reduction effect.

Beeler, G. B.↗
NASA NTRS1987

An electrochemical series of redox couples in silicate melts - A review and applications to geochemistry

An electrochemical series for redox couples in a glass-forming oxide melt is developed. This series is a quantitative numerical scale of reference reduction potentials of the redox couples in a silicate melt that is a model for basaltic magmas. The redox couples are ordered in terms of their reference reduction potentials; the order appears to be relatively independent of the exact melt composition and temperature. Thus, upon calibration to a desired composition, oxygen fugacity, and temperature, this electrochemical series can provide estimates of redox state proportions in basaltic magmas on different planetary bodies. The geochemical electrochemical series can also be used to understand the interrelationship of the redox state of the magma and the presence of volatile species such as oxygen, water, sulfur gases, and carbon gases.

Schreiber, Henry D.↗
NASA NTRS1998

Method and apparatus for reducing the drag of flows over surfaces

An apparatus, and its accompanying method, for reducing the drag of flows over a surface includes arrays of small disks and sensors. The arrays are embedded in the surface and may extend above, or be depressed below, the surface, provided they remain hydraulically smooth either when operating or when inactive. The disks are arranged in arrays of various shapes, and spaced according to the cruising speed of the vehicle on which the arrays are installed. For drag reduction at speeds of the order of 30 meters/second, preferred embodiments include disks that are 0.2 millimeter in diameter and spaced 0.4 millimeter apart. For drag reduction at speeds of the order of 300 meters/second, preferred embodiments include disks that are 0.045 millimeter in diameter and spaced 0.09 millimeter apart. Smaller and larger dimensions for diameter and spacing are also possible. The disks rotate in the plane of the surface, with their rotation axis substantially perpendicular to the surface. The rotating disks produce velocity perturbations parallel to the surface in the overlying boundary layer. The sensors sense the flow at the surface and connect to control circuitry that adjusts the rotation rates and duty cycles of the disks accordingly. Suction and blowing holes can be interspersed among, or made coaxial with, the disks for creating general three-component velocity perturbations in the near-surface region. The surface can be a flat, planar surface or a nonplanar surface, such as a triangular riblet surface. The present apparatus and method have potential applications in the field of aeronautics for improving performance and efficiency of commercial and military aircraft, and in other industries where drag is an obstacle, including gas and oil delivery through long-haul pipelines.

Keefe, Laurence R.↗
NASA NTRS1989

The role of topography in geodetic gravity field modelling

Masses associated with the topography, bathymetry, and its isostatic compensation are a dominant source of gravity field variations, especially at shorter wavelengths. On global scales the topographic/isostatic effects are also significant, except for the lowest harmonics. In practice, though, global effects need not be taken into account as such effects are included in the coefficients of the geopotential reference fields. On local scales, the short-wavelength gravity variations due to the topography may, in rugged terrain, be an order of magnitude larger than other effects. In such cases, explicit or implicit terrain reduction procedures are mandatory in order to obtain good prediction results. Such effects may be computed by space-domain integration or by fast Fourier transformation (FFT) methods. Numerical examples are given for areas of the Canadian Rockies. In principle, good knowledge of the topographic densities is required to produce the smoothest residual field. Densities may be determined from sample measurements or by gravimetric means, but both are somewhat troublesome methods in practice. The use of a standard density, e.g., 2.67 g/cu cm, may often yield satisfactory results and may be put within a consistent theoretical framework. The independence of density assumptions is the key point of the classical Molodensky approach to the geodetic boundary value problem. The Molodensky solutions take into account that land gravity field observations are done on a non-level surface. Molodensky's problem may be solved by integral expansions or more effective FFT methods, but the solution should not be intermixed with the use of terrain reductions. The methods are actually complimentary and may both be required in order to obtain the smoothest possible signal, least prone to aliasing and other effects coming from sparse data coverage, typical of rugged topography.

Forsberg, R.↗
NASA NTRS2011

Rapid Calculation of Spacecraft Trajectories Using Efficient Taylor Series Integration

A variable-order, variable-step Taylor series integration algorithm was implemented in NASA Glenn's SNAP (Spacecraft N-body Analysis Program) code. SNAP is a high-fidelity trajectory propagation program that can propagate the trajectory of a spacecraft about virtually any body in the solar system. The Taylor series algorithm's very high order accuracy and excellent stability properties lead to large reductions in computer time relative to the code's existing 8th order Runge-Kutta scheme. Head-to-head comparison on near-Earth, lunar, Mars, and Europa missions showed that Taylor series integration is 15.8 times faster than Runge- Kutta on average, and is more accurate. These speedups were obtained for calculations involving central body, other body, thrust, and drag forces. Similar speedups have been obtained for calculations that include J2 spherical harmonic for central body gravitation. The algorithm includes a step size selection method that directly calculates the step size and never requires a repeat step. High-order Taylor series integration algorithms have been shown to provide major reductions in computer time over conventional integration methods in numerous scientific applications. The objective here was to directly implement Taylor series integration in an existing trajectory analysis code and demonstrate that large reductions in computer time (order of magnitude) could be achieved while simultaneously maintaining high accuracy. This software greatly accelerates the calculation of spacecraft trajectories. At each time level, the spacecraft position, velocity, and mass are expanded in a high-order Taylor series whose coefficients are obtained through efficient differentiation arithmetic. This makes it possible to take very large time steps at minimal cost, resulting in large savings in computer time. The Taylor series algorithm is implemented primarily through three subroutines: (1) a driver routine that automatically introduces auxiliary variables and sets up initial conditions and integrates; (2) a routine that calculates system reduced derivatives using recurrence relations for quotients and products; and (3) a routine that determines the step size and sums the series. The order of accuracy used in a trajectory calculation is arbitrary and can be set by the user. The algorithm directly calculates the motion of other planetary bodies and does not require ephemeris files (except to start the calculation). The code also runs with Taylor series and Runge-Kutta used interchangeably for different phases of a mission.

Scott, James R.↗