Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed parallelization”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

A Tower-based Prototype VHF/UHF Radar for Subsurface Sensing: System Description and Data Inversion Results

Knowledge of subsurface characteristics such as permittivity variations and layering structure could provide a breakthrough in many terrestrial and planetary science disciplines. For Earth science, knowledge of subsurface and subcanopy soil moisture layers can enable the estimation of vertical flow in the soil column linking surface hydrologic processes with that in the subsurface. For planetary science, determining the existence of subsurface water and ice is regarded as one of the most critical information needs for the study of the origins of the solar system. The subsurface in general can be described as several near-parallel layers with rough interfaces. Each homogenous rough layer can be defined by its average thickness, permittivity, and rms interface roughness assuming a known surface spectral distribution. As the number and depth of layers increase, the number of measurements needed to invert for the layer unknowns also increases, and deeper penetration capability would be required. To nondestructively calculate the characteristics of the rough layers, a multifrequency polarimetric radar backscattering approach can be used. One such system is that we have developed for data prototyping of the Microwave Observatory of Subcanopy and Subsurface (MOSS) mission concept. A tower-mounted radar makes backscattering measurements at VHF, UHF, and L-band frequencies. The radar is a pulsed CW system, which uses the same wideband antenna to transmit and receive the signals at all three frequencies. To focus the beam at various incidence angles within the beamwidth of the antenna, the tower is moved vertically and measurements made at each position. The signals are coherently summed to achieve focusing and image formation in the subsurface. This requires an estimate of wave velocity profiles. To solve the inverse scattering problem for subsurface velocity profile simultaneously with radar focusing, we use an iterative technique based on a forward numerical solution of the layered rough surface problem. The layers are each defined in terms of a small number of unknown distributions as given above. An a priori estimate of the solution is first assumed, based on which the forward problem is solved for the backscattered measurements. This is compared with the measured data and using iterative techniques an update to the solution for the unknowns is calculated. The process continues until convergence is achieved. Numerical results will be shown using actual radar data acquired with the MOSS tower radar system in Arizona in Fall 2003, and compared with in-situ measurements.

Moghaddam, Mahta↗
NASA NTRS2008

Spaceborne Processor Array

A Spaceborne Processor Array in Multifunctional Structure (SPAMS) can lower the total mass of the electronic and structural overhead of spacecraft, resulting in reduced launch costs, while increasing the science return through dynamic onboard computing. SPAMS integrates the multifunctional structure (MFS) and the Gilgamesh Memory, Intelligence, and Network Device (MIND) multi-core in-memory computer architecture into a single-system super-architecture. This transforms every inch of a spacecraft into a sharable, interconnected, smart computing element to increase computing performance while simultaneously reducing mass. The MIND in-memory architecture provides a foundation for high-performance, low-power, and fault-tolerant computing. The MIND chip has an internal structure that includes memory, processing, and communication functionality. The Gilgamesh is a scalable system comprising multiple MIND chips interconnected to operate as a single, tightly coupled, parallel computer. The array of MIND components shares a global, virtual name space for program variables and tasks that are allocated at run time to the distributed physical memory and processing resources. Individual processor- memory nodes can be activated or powered down at run time to provide active power management and to configure around faults. A SPAMS system is comprised of a distributed Gilgamesh array built into MFS, interfaces into instrument and communication subsystems, a mass storage interface, and a radiation-hardened flight computer.

Chow, Edward T.↗
NASA NTRS1995

The Galactic Distribution of Aliphatic Hydrocarbons in the Diffuse Interstellar Medium

The infrared absorption feature near 2950/ cm (3.4 microns), characteristic of dust in the diffuse interstellar medium (ISM), is attributed to C-H stretching vibrations of aliphatic hydrocarbons. We show here that the strength of the band does not scale linearly with visual extinction everywhere, but instead increases more rapidly for objects near the center of the Galaxy, a behavior that parallels that of the Si-O stretching band due to silicate materials in the diffuse ISM. This implies that the grains responsible for the diffuse medium aliphatic C-H and silicate Si-O stretching bands are different from those responsible for much of the observed visual extinction. It also suggests that the distribution of the carbonaceous component of the diffuse ISM is not uniform throughout the Galaxy, but instead may increase in density toward the center of the Galaxy. The similar behavior of the C-H and Si-O stretching bands suggests that these two components may be coupled, perhaps in the form of silicate-core, organic-mantle grains. Several possible models of the distribution of this material are presented and it is demonstrated that the inner parts of the Galaxy has a carrier density that is 5 to 35 times higher than in the local ISM. Depending on the model used, the density of aliphatic material in the local ISM is found to be about 1 to 2 -CH3 groups /cc and about 2 to 5 -CH2- groups/cc. These densities are consistent with the strengths of the 2955 and 2925/ cm (3.38 and 3.42 microns) subfeatures (due to -CH3 and -CH2- groups, respectively) within the overall 2950/ cm (3.4 microns) band being described by the relations A(sub upsilon)/tau(sub 2925/cm) = 270 +/- 40 and A(sub upsilon)/tau(sub 2925/cm) = 250 +/- 40 in the local diffuse ISM.

Sandford, Scott A.↗
NASA NTRS1994

Analysis of coseismic surface displacement gradients using radar interferometry: New insights into the Landers earthquake

The map of the coseismic displacement field generated by interferometric processing of synthetic aperture radar (SAR) images taken before and after the June 28, 1992, Landers earthquake sequence brings new insights into the nature of deformation caused by these earthquakes. We use the interferometric map generated by Massonnet et al. (1993) to analyze the surface displacement field in the vicinity of the fault trace. Complexities in the fringe pattern near the fault reflect short-wavelength variations of the surface rupture and slip distribution, and attest to large displacement gradients. Along two sections of the fault, characteristic fringe patterns can be recognized, contrasting in density and direction with patterns observed away from the rupture. In order to understand the observed fringe patterns, we compute synthetic interferograms in three simple cases: (1) rigid-body rotations about a vertical axis, (2) about a horizontal axis (tilt), and (3) distributed, simple shear. The orientation and spatial separation of interferometric fringes predicted by these models help constrain near-field deformation and rupture parameters. Where the Kickapoo fault connects with the Homestead Valley fault, the interferogram shows a clear pattern of parallel N20 deg W fringes separated by about 160 m. This pattern and vertical offsets measured along the Kickapoo fault suggest that the block between this fault and the Johnson Valley fault may have been tilted, down to the west. A 5-km block lifted by 1 m on one side would be tilted by an angle of 0.01 deg (190 microrad), producing fringes separated by about 160 m, parallel to the tilt axis. Such a tilt, parallel to a N20 deg W direction, would account for the gradual, northward increase of the vertical slip component observed along the Kickapoo fault. This tilt may also explain the 1 m of reverse slip observed along the 'slip gap' section of the Homestead Valley break. Between the southern end of the Johnson Valley fault and the Eureka Peak fault, where no surface rupture has been mapped, the dense pattern of fringes implies distributed shear, probably resulting from fault slip at depth. The density and direction of the fringes in the gap are consistent with a right-lateral slip of 1.2-3.8 m on a blind fault locked above the depth of 1.5-2 km. Such observations of small wavelength features in the SAR interferogram bring new insights into the near-field displacement gradient and thus on response of the uppermost crust to seismic rupture.

Peltzer, Gilles↗
NASA NTRS1995

The galactic distribution of aliphatic hydrocarbons in the diffuse interstellar medium

The infrared absorption feature near 2950(exp -1) (3.4 micron), characteristic of dust in the diffuse interstellar medium (ISM), is attributed to C-H stretching vibrations of aliphatic hydrocarbons. We show here that the strength of the band does not scale linearly with visual extinction everywhere, but instead increases more rapidly for objects near the center of the Galaxy, a behavior that parallels that of the Si-O stretching band due to silicate materials in the diffuse ISM. This implies that the grains responsible for the diffuse medium aliphatic C-H and silicate Si-O stretching bands are different from those responsible for much of the observed visual extinction. It also suggests that the distribution of the carbonaceous component of the diffuse ISM is not uniform throughout the Galaxy, but instead may increase in density toward the center of the Galaxy. The similar behavior of the C-H and Si-O stretching bands suggests that these two components may be coupled, perhaps in the form of silicate-core, organic-mantle grains. Several possible models of the distribution of this material are presented and it is demonstrated that the inner parts of the Galaxy has a carrier density that is 5 to 35 times higher than in the local ISM. Depending on the model used, the density of aliphatic material in the local ISM is found to be about 1 to 2 -CH3 groups m(exp -3) and about 2 to 5 -CH2- groups m(exp -3). These densities are consistent with the strengths of the 2955 and 2925 cm(exp -1) (3.4 micron) band being described by the relations A(sub nu)/tau(sub 2955 cm(exp -1)) = 270 +/- 40 and A(sub nu)/tau(sub 2925 cm(exp -1)) = 250 +/- 40 in the local diffuse ISM.

Sandford, Scott A.↗
NASA NTRS1991

Biofilm monitoring coupon system and method of use

An apparatus and method is disclosed for biofilm monitoring of a water distribution system which includes the mounting of at least one fitting in a wall port of a manifold in the water distribution system with a passage through the fitting in communication. The insertion of a biofilm sampling member is through the fitting with planar sampling surfaces of different surface treatment provided on linearly arrayed sample coupons of the sampling member disposed in the flow stream in edge-on parallel relation to the direction of the flow stream of the manifold under fluid-tight sealed conditions. The sampling member is adapted to be aseptically removed from or inserted in the fitting and manifold under a positive pressure condition and the fitting passage sealed immediately thereafter by appropriate closure means so as to preclude contamination of the water distribution system through the fitting. The apparatus includes means for clamping the sampling member and for establishing electrical continuity between the sampling surfaces and the system for minimizing electropotential effects. The apparatus may also include a plurality of fittings and sampling members mounted on the manifold to permit extraction of the sampling members in a timed sequence throughout the monitoring period.

Sauer, Richard L.↗
NASA NTRSDate not supplied

Adding GPU Support to the Markov Chain Monte Carlo Code Catmip

In geophysics, we are confronted with many under-determined inverse problems. For example, all of our observations of earthquakes are made at the Earth’s surface. So, when we try to infer how slip during an earthquake evolves in space and time, we find that there are many potential slip histories that are consistent with our limited observations and our understanding of earthquake physics. One way to approach these problems is with Bayesian analysis which allows us to infer the ensemble of all potential slip models that satisfy the observations and our prior knowledge of earthquake physics. In Bayesian analysis, our prior knowledge is known as the prior probability density function or prior PDF, the fit to the data is known as the data likelihood, and the target PDF that satisfies both the prior PDF and data likelihood is known as the posterior PDF. However, simulating the posterior PDF typically requires using Markov Chain Monte Carlo (MCMC) to draw tens of billions of random realizations of earthquake slip models, which may not be computationally feasible. To make this and similar geophysical inversions computationally tractable, we developed the Cascading Adaptive Transitional Metropolis In Parallel (CATMIP) algorithm. CATMIP is an efficient parallel Markov Chain Monte Carlo (MCMC) sampler that is used for model fitting and uncertainty quantification in geophysics. Example use cases are earthquake rupture modeling, determining mineral composition on Mars, reconstructing the history of ocean salinity, and historical earthquake relocation. CATMIP employs many parallel instances of the Metropolis algorithm for sampling in a transitioning framework. Transitioning is a process in which a set of random samples at equilibrium with a known probability density function (PDF) are used as seeds for the Markov chains to sample successive target PDFs that incrementally move the distribution from the starting seeds to the final desired PDF that describes the relative plausibility of potential values for the model parameters. The algorithm is implemented as a Master-Worker model employing MPI for communication. The worker processes are loosely coupled with global parameters periodically optimized by the master process. This provides a very high amount of parallelism with little communication between updates. During the presentation we will discuss the history of the algorithm and elaborate the earthquake rupture modeling use case for the CATMIP package. Our first step toward GPU optimization was to optimize the code for the CPU. CPU profiling revealed that most of the compute time is spent in calls to level 2 BLAS routines and calls to GSL random number generators. We revised the algorithm to employ level 3 BLAS routines instead. In our presentation we will describe how this was accomplished. Adding GPU support to CATMIP consisted mostly of replacing the calls to GSL with calls to GPU vendor-provided library routines. A small number of loops were directly implemented in CUDA. In the presentation will provide implementation details. Finally, we will discuss methods for profiling and opportunities for further optimizing GPU execution. By creating a code with the flexibility to run on either a CPU or GPU architecture, CATMIP can be used on systems ranging from large CPU-based HPC environments to single servers with GPU acceleration and everything in between.

HECC↗
NASA NTRS1982

The PKI collector

The design and manufacturing of a solar thermal collector is discussed. The collector has three primary subsystems: concentrator, receiver/fluid loop, and controls. Identical curved reflective columns are utilized in a faceted Fresnel design to support 864 one foot square flat inexpensive second-surface, silvered glass mirrors. The columns are ganged together and rotated through their centers of gravity to provide elevation tracking. The concentrator is supported by a lightweight spaceframe structure which distributes all wind and gravity loads to the base supports. The base of the structure is a track which rotates on wheels mounted on concrete piers. A parallel tube steel heat exchanger is mounted at the concentrator focal area in a well insulated, galvanized steel housing. Two rows of vertical close-packed, staggered tubes connect a mud header and a steam header. Automatic two axis tracking and operational control is provided with a microprocessor based package. Concentrator-mounted shadowbands are the basis for active tracking. A software program provides azimuthal tracking during cloudy periods.

Rice, M. P.↗
NASA NTRS1987

Parallel discrete event simulation: A shared memory approach

With traditional event list techniques, evaluating a detailed discrete event simulation model can often require hours or even days of computation time. Parallel simulation mimics the interacting servers and queues of a real system by assigning each simulated entity to a processor. By eliminating the event list and maintaining only sufficient synchronization to insure causality, parallel simulation can potentially provide speedups that are linear in the number of processors. A set of shared memory experiments is presented using the Chandy-Misra distributed simulation algorithm to simulate networks of queues. Parameters include queueing network topology and routing probabilities, number of processors, and assignment of network nodes to processors. These experiments show that Chandy-Misra distributed simulation is a questionable alternative to sequential simulation of most queueing network models.

Reed, Daniel A.↗
NASA NTRS1988

Parallel discrete event simulation using shared memory

With traditional event-list techniques, evaluating a detailed discrete-event simulation-model can often require hours or even days of computation time. By eliminating the event list and maintaining only sufficient synchronization to ensure causality, parallel simulation can potentially provide speedups that are linear in the numbers of processors. A set of shared-memory experiments, using the Chandy-Misra distributed-simulation algorithm, to simulate networks of queues is presented. Parameters of the study include queueing network topology and routing probabilities, number of processors, and assignment of network nodes to processors. These experiments show that Chandy-Misra distributed simulation is a questionable alternative to sequential-simulation of most queueing network models.

Reed, Daniel A.↗
NASA NTRS1993

Optoelectronic Control Of Phased-Array Antenna

Radiation patterns switched within microseconds. According to concept, digital beam-steering information destined for number of phased antenna elements transmitted in serial multiplex along optical fiber to optoelectronic integrated circuit (OEIC), demultiplexing information into parallel electronic channels to phase shifters associated with various elements. Because of small distances between elements of such arrays, compactness of optical fibers makes them attractive for distribution of beam-steering (phase-shifting) control signals. Also offers advantages of light weight, low attenuation, mechanical flexibility, large bandwidth, and immunity to both crosstalk and electromagnetic interference.

Bhasin, K. B.↗
NASA NTRS1995

Mississippi State University Center for Air Sea Technology FY95 Research Program

The Mississippi State University (MSU) Center for Air Sea Technology (CAST) evolved from the Institute for Naval Oceanography's (INO) Experimental Center for Mesoscale Ocean Prediction (ECMOP) which was started in 1989. MSU CAST subsequently began operation on 1 October 1992 under an Office of Naval Research (ONR) two-year grant which ended on 30 September 1994. In FY95 MSU CAST was successful in obtaining five additional research grants from ONR, as well as several other research contracts from the Naval Oceanographic Office via NASA, the Naval Research Laboratory, the Army Corps of Engineers, and private industry. In the past, MSU CAST technical research and development has produced tools, systems, techniques, and procedures that improve efficiency and overcome deficiency for both the operational and research communities residing with the Department of Defense, private industry, and university ocean modeling community. We continued this effort with the following thrust areas: to develop advanced methodologies and tools for model evaluation, validation and visualization, both oceanographic and atmospheric; to develop a system-level capability for conducting temporally and ; spatially scaled ocean simulations driven by or are responsive to ocean models, and take into consideration coupling to atmospheric models; to continue the existing oceanographic/atmospheric data management task with emphasis on distributed databases in a network environment, with database optimization and standardization, including use of Mosaic and World Wide Web (WWW) access; and to implement a high performance parallel computing technology for CAST ocean models

Yeske, Lanny↗
NASA NTRS1993

Viewpoint 9--molecular structure of aqueous interfaces

In this review we summarize recent progress in our understanding of the structure of aqueous interfaces emerging from molecular level computer simulations. It is emphasized that the presence of the interface induces specific structural effects which, in turn, influence a wide variety of phenomena occurring near the phase boundaries. At the liquid-vapor interface, the most probable orientations of a water molecule is such that its dipole moment lies parallel to the interface, one O-H bond points toward the vapor and the other O-H bond is directed toward the liquid. The orientational distributions are broad and slightly asymmetric, resulting in an excess dipole moment pointing toward the liquid. These structural preferences persist at interfaces between water and nonpolar liquids, indicating that the interactions between the two liquids in contact are weak. It was found that liquid-liquid interfaces are locally sharp but broadened by capillary waves. One consequence of anisotropic orientations of interfacial water molecules is asymmetric interactions, with respect to the sign of the charge, of ions with the water surface. It was found that even very close to the surface ions retain their hydration shells. New features of aqueous interfaces have been revealed in studies of water-membrane and water-monolayer systems. In particular, water molecules are strongly oriented by the polar head groups of the amphiphilic phase, and they penetrate the hydrophilic head-group region, but not the hydrophobic core. At infinite dilution near interfaces, amphiphilic molecules exhibit behavior different from that in the gas phase or in bulk water. This result sheds new light on the nature of hydrophobic effect in the interfacial regions. The presence of interfaces was also shown to affect both equilibrium and dynamic components of rates of chemical reactions. Applications of continuum models to interfacial problems have been, so far, unsuccessful. This, again, underscores the importance of molecular-level information about interfaces.

Review, Tutorial↗
NASA NTRS2006

Phased-Array Study of Dual-Flow Jet Noise: Effect of Nozzles and Mixers

A 16-microphone linear phased-array installed parallel to the jet axis and a 32-microphone azimuthal phased-array installed in the nozzle exit plane have been applied to identify the noise source distributions of nozzle exhaust systems with various internal mixers (lobed and axisymmetric) and nozzles (three different lengths). Measurements of velocity were also obtained using cross-stream stereo particle image velocimetry (PIV). Among the three nozzle lengths tested, the medium length nozzle was the quietest for all mixers at high frequency on the highest speed flow condition. Large differences in source strength distributions between nozzles and mixers occurred at or near the nozzle exit for this flow condition. The beamforming analyses from the azimuthal array for the 12-lobed mixer on the highest flow condition showed that the core flow and the lobe area were strong noise sources for the long and short nozzles. The 12 noisy spots associated with the lobe locations of the 12-lobed mixer with the long nozzle were very well detected for the frequencies 5 KHz and higher. Meanwhile, maps of the source strength of the axisymmetric splitter show that the outer shear layer was the most important noise source at most flow conditions. In general, there was a good correlation between the high turbulence regions from the PIV tests and the high noise source regions from the phased-array measurements.

Soo Lee, Sang↗
NASA NTRS2007

The Effects of Rainfall Inhomogeneity on Climate Variability of Rainfall Estimated from Passive Microwave Sensors

Passive microwave rainfall estimates that exploit the emission signal of raindrops in the atmosphere are sensitive to the inhomogeneity of rainfall within the satellite field of view (FOV). In particular, the concave nature of the brightness temperature (T(sub b)) versus rainfall relations at frequencies capable of detecting the blackbody emission of raindrops cause retrieval algorithms to systematically underestimate precipitation unless the rainfall is homogeneous within a radiometer FOV, or the inhomogeneity is accounted for explicitly. This problem has a long history in the passive microwave community and has been termed the beam-filling error. While not a true error, correcting for it requires a priori knowledge about the actual distribution of the rainfall within the satellite FOV, or at least a statistical representation of this inhomogeneity. This study first examines the magnitude of this beam-filling correction when slant-path radiative transfer calculations are used to account for the oblique incidence of current radiometers. Because of the horizontal averaging that occurs away from the nadir direction, the beam-filling error is found to be only a fraction of what has been reported previously in the literature based upon plane-parallel calculations. For a FOV representative of the 19-GHz radiometer channel (18 km X 28 km) aboard the Tropical Rainfall Measuring Mission (TRMM), the mean beam-filling correction computed in this study for tropical atmospheres is 1.26 instead of 1.52 computed from plane-parallel techniques. The slant-path solution is also less sensitive to finescale rainfall inhomogeneity and is, thus, able to make use of 4-km radar data from the TRMM Precipitation Radar (PR) in order to map regional and seasonal distributions of observed rainfall inhomogeneity in the Tropics. The data are examined to assess the expected errors introduced into climate rainfall records by unresolved changes in rainfall inhomogeneity. Results show that global mean monthly errors introduced by not explicitly accounting for rainfall inhomogeneity do not exceed 0.5% if the beam-filling error is allowed to be a function of rainfall rate and freezing level and does not exceed 2% if a universal beam-filling correction is applied that depends only upon the freezing level. Monthly regional errors can be significantly larger. Over the Indian Ocean, errors as large as 8% were found if the beam-filling correction is allowed to vary with rainfall rate and freezing level while errors of 15% were found if a universal correction is used.

Kummerow, Christian↗
NASA NTRS1993

Plains tectonism on Venus: The deformation belts of Lavinia Planitia

High-resolution radar images from the Magellan spacecraft have revealed the first details of the morphology of the Lavinia Planitia region of Venus. A number of geologic units can be distinguished, including volcanic plains units with a range of ages. Transecting these plains over much of the Lavinia region are two types of generally orthogonal features that we interpret to be compressional wrinkle ridges and extensional grooves. The dominant tectonic features of Lavinia are broad elevated belts of intense deformation that transect the plains with complex geometry. They are many tens to a few hundred kilometers wide, as much as 1000 km long, and elevated hundreds of meters above the surrounding plains. Two classes of deformation belts are seen in the Lavinia region. 'Ridge belts' are composed of parallel ridges, each a few hundred meters in elevation, that we interpret to be folds. Typical fold spacings are 5-10 km. 'Fracture belts' are dominated instead by intense faulting, with faults in some instances paired to form narrow grabens. There is also some evidence for modest amounts of horizontal shear distributed across both ridge and fracture belts. Crosscutting relationships among the belts show there to be a range in belt ages. In western Lavinia, in particular, many ridge and fracture belts appear to bear a relationship to the much smaller wrinkle ridges and grooves on the surrounding plains: ridge morphology tends to dominate belts that lie more nearly parallel to local plains wrinkle ridges, and fracture morphology tends to dominate belts that lie more nearly parallel to local plains grooves. We use simple models to explore the formation of ridge and fracture belts. We show that convective motions in the mantle can couple to the crust to cause horizontal stresses of a magnitude sufficient to induce the formation of deformation belts like those observed in Lavinia. We also use the small-scale wavelengths of deformation observed within individual ridge belts to place an approximate lower limit on the venusian thermal gradient in the Lavinia region at the time of deformation.

Squyres, Steven W.↗
NASA NTRS1992

Distributed project scheduling at NASA: Requirements for manual protocols and computer-based support

The increasing complexity of space operations and the inclusion of interorganizational and international groups in the planning and control of space missions lead to requirements for greater communication, coordination, and cooperation among mission schedulers. These schedulers must jointly allocate scarce shared resources among the various operational and mission oriented activities while adhering to all constraints. This scheduling environment is complicated by such factors as the presence of varying perspectives and conflicting objectives among the schedulers, the need for different schedulers to work in parallel, and limited communication among schedulers. Smooth interaction among schedulers requires the use of protocols that govern such issues as resource sharing, authority to update the schedule, and communication of updates. This paper addresses the development and characteristics of such protocols and their use in a distributed scheduling environment that incorporates computer-aided scheduling tools. An example problem is drawn from the domain of Space Shuttle mission planning.

Richards, Stephen F.↗
NASA NTRS2005

Evolutionary Computing for Low-thrust Navigation

The development of new mission concepts requires efficient methodologies to analyze, design and simulate the concepts before implementation. New mission concepts are increasingly considering the use of ion thrusters for fuel-efficient navigation in deep space. This paper presents parallel, evolutionary computing methods to design trajectories of spacecraft propelled by ion thrusters and to assess the trade-off between delivered payload mass and required flight time. The developed methods utilize a distributed computing environment in order to speed up computation, and use evolutionary algorithms to find globally Pareto-optimal solutions. The methods are coupled with two main traditional trajectory design approaches, which are called direct and indirect. In the direct approach, thrust control is discretized in either arc time or arc length, and the resulting discrete thrust vectors are optimized. In the indirect approach, a thrust control problem is transformed into a costate control problem, and the initial values of the costate vector are optimized. The developed methods are applied to two problems: 1) an orbit transfer around the Earth and 2) a transfer between two distance retrograde orbits around Europa, the closest to Jupiter of the icy Galilean moons. The optimal solutions found with the present methods are comparable to other state-of-the-art trajectory optimizers and to analytical approximations for optimal transfers, while the required computational time is several orders of magnitude shorter than other optimizers thanks to an intelligent design of control vector discretization, advanced algorithmic parameterization, and parallel computing.

optimization↗