Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “monotonicity”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Hybrid Differential Dynamic Programming with Stochastic Search

Differential dynamic programming (DDP) has been demonstrated as a viable approach to low-thrust trajectory optimization, namely with the recent success of NASAs Dawn mission. The Dawn trajectory was designed with the DDP-based Static Dynamic Optimal Control algorithm used in the Mystic software. Another recently developed method, Hybrid Differential Dynamic Programming (HDDP) is a variant of the standard DDP formulation that leverages both first-order and second-order state transition matrices in addition to nonlinear programming (NLP) techniques. Areas of improvement over standard DDP include constraint handling, convergence properties, continuous dynamics, and multi-phase capability. DDP is a gradient based method and will converge to a solution nearby an initial guess. In this study, monotonic basin hopping (MBH) is employed as a stochastic search method to overcome this limitation, by augmenting the HDDP algorithm for a wider search of the solution space.

Aziz, Jonathan↗
NASA NTRS2016

Pilot Critical Incident Reports as a Means to Identify Human Factors of Remotely Piloted Aircraft

It has been estimated that aviation accidents are typically preceded by numerous minor incidents arising from the same causal factors that ultimately produced the accident. Accident databases provide in-depth information on a relatively small number of occurrences, however incident databases have the potential to provide insights into the human factors of Remotely Piloted Aircraft System (RPAS) operations based on a larger volume of less-detailed reports. Currently, there is a lack of incident data dealing with the human factors of unmanned aircraft systems. An exploratory study is being conducted to examine the feasibility of collecting voluntary critical incident reports from RPAS pilots. Twenty-three experienced RPAS pilots volunteered to participate in focus groups in which they described critical incidents from their own experience. Participants were asked to recall (1) incidents that revealed a system flaw, or (2) highlighted a case where the human operator contributed to system resilience or mission success. Participants were asked to only report incidents that could be included in a public document. During each focus group session, a note taker produced a de-identified written record of the incident narratives. At the end of the session, participants reviewed each written incident report, and made edits and corrections as necessary. The incidents were later analyzed to identify contributing factors, with a focus on design issues that either hindered or assisted the pilot during the events. A total of 90 incidents were reported. Human factor issues included the impact of reduced sensory cues, traffic separation in the absence of an out-the-window view, control latencies, vigilance during monotonous and ultra-long endurance flights, control station design considerations, transfer of control between control stations, the management of lost link procedures, and decision-making during emergencies. Pilots participated willingly and enthusiastically in the study, and generally had little difficulty recalling critical incidents. The results suggest that pilot interviews can be a productive method of gathering information on incidents that might not otherwise be reported. Some of the issues described in the reports have received significant attention in the literature, or are analogous to human factors of manned aircraft. In other cases, incident reports involved human factors that are poorly understood, and have not yet been the subject of extensive study. Although many of the reported incidents were related to pilot error, the participants also provided examples of the positive contribution that humans make to the operation of highly-automated systems.

unmanned aircraft systems↗
NASA NTRS2016

Dynamic Characterization of Galfenol (Fe81.6Ga18.4)

Galfenol has the potential to transform the smart materials industry by allowing for the development of multifunctional, load-bearing devices. One of the primary technical challenges faced by this development is the very limited experimental data on Galfenol's frequency-dependent response to dynamic stress, which is critically important for the design of such devices. This report details a novel and precise characterization of the constitutive behavior of polycrystalline Galfenol (Fe81.6Ga18.4) under quasi-static (1 Hz) and dynamic (4 to 1000 Hz) stress loadings. Mechanical loads are applied using a high-frequency load frame. Quasi-static minor and major hysteresis loop measurements of magnetic flux density and strain are presented for constant electromagnet currents (0 to 1.1 A) and constant magnetic fields 0 to 14 kA/m (0 to 180 Oe). The dynamic stress amplitude for minor and major loops is 2.88 and 31.4 MPa (418 and 4550 psi), respectively. Quasi-static material properties closely match published values for similar Galfenol materials. Quasi-static actuation responses are also measured and compared to quasi-static sensing responses; the high degree of reversibility seen in the comparison is consistent with published measurements and modeling results. Dynamic major and minor loops are measured for dynamic stresses up to 31 MPa (4496 psi) and 1 kHz, and the bias condition resulting in maximum, quasi-static sensitivity. Eddy current effects are quantified by considering solid and laminated Galfenol rods. Three key sources of error in the dynamic measurements are accounted for: (1) electromagnetic noise in strain signals due to Galfenol's magnetic response, (2) error in load signals due to the inertial force of fixturing, and (3) phase misalignment between signals due to conditioning electronics. For dynamic characterization, strain error is kept below 1.2 percent of full scale by wiring two collocated gauges in series (noise cancellation) and through leadwire weaving. Inertial force error is kept below 0.41 percent by measuring the dynamic force in the specimen using a nearly collocated piezoelectric load washer. The phase response of all conditioning electronics is explicitly measured and corrected for. In general, as frequency is increased, the sensing response becomes more linear because of an increase in eddy currents. As frequency increases above approximately 100 Hz, the elbow in the strain-versus-stress response disappears as the active (soft) regime stiffens toward the passive (hard) regime. Under constant-field conditions, the loss factors of the solid rod peak between 200 and 600 Hz, rather than exhibiting a monotonic increase. Compared to the solid rod, the laminated rod exhibits much slower increases in hysteresis with frequency, and its quasi-static behavior extends to higher frequencies. The elastic modulus of the laminated rod decreases between 100 and 300 Hz; this trend is currently unexplained.

Dynamic sensing↗
NASA NTRS2016

Warming of the Global Ocean: Spatial Structure and Water-Mass Trends

This study investigates the multidecadal warming and interannual-to-decadal heat content changes in the upper ocean (0-700 m), focusing on vertical and horizontal patterns of variability. These results support a nearly monotonic warming over much of the World Ocean, with a shift toward Southern Hemisphere warming during the well-observed past decade. This is based on objectively analyzed gridded observational datasets and on a modeled state estimate. Besides the surface warming, a warming climate also has a subsurface effect manifesting as a strong deepening of the midthermocline isopycnals, which can be diagnosed directly from hydrographic data. This deepening appears to be a result of heat entering via subduction and spreading laterally from the high-latitude ventilation regions of subtropical mode waters. The basin-average multidecadal warming mainly expands the subtropical mode water volume, with weak changes in the temperature-salinity (u-S) relationship (known as ''spice'' variability). However, the spice contribution to the heat content can be locally large, for example in Southern Hemisphere. Multidecadal isopycnal sinking has been strongest over the southern basins and weaker elsewhere with the exception of the Gulf Stream/North Atlantic Current/subtropical recirculation gyre. At interannual to decadal time scales, wind-driven sinking and shoaling of density surfaces still dominate ocean heat content changes, while the contribution from temperature changes along density surfaces tends to decrease as time scales shorten.

Global Ocean↗
NASA NTRS2016

The Cognitive Challenges of Flying a Remotely Piloted Aircraft

A large variety of Remotely Piloted Aircraft (RPA) designs are currently in production or in development. These aircraft range from small electric quadcopters that are flown close to the ground within visual range of the operator, to larger systems capable of extended flight in airspace shared with conventional aircraft. Before RPA can operate routinely and safely in civilian airspace, we need to understand the unique human factors associated with these aircraft. The task of flying an RPA in civilian airspace involves challenges common to the operation of other highly-automated systems, but also introduces new considerations for pilot perception, decision-making, and action execution. RPA pilots participated in focus groups where they were asked to recall critical incidents that either presented a threat to safety, or highlighted a case where the pilot contributed to system resilience or mission success. Ninety incidents were gathered from focus-groups. Human factor issues included the impact of reduced sensory cues, traffic separation in the absence of an out-the-window view, control latencies, vigilance during monotonous and ultra-long endurance flights, control station design considerations, transfer of control between control stations, the management of lost link procedures, and decision-making during emergencies. Some of these concerns have received significant attention in the literature, or are analogous to human factors of manned aircraft. The presentation will focus on issues that are poorly understood, and have not yet been the subject of extensive human factors study. Although many of the reported incidents were related to pilot error, the participants also provided examples of the positive contribution that humans make to the operation of highly-automated systems.

remote pilot station↗
NASA NTRS2016

Human Factors of Remotely Piloted Aircraft Systems: Lessons from Incident Reports

An exploratory study is being conducted to examine the feasibility of collecting voluntary critical incident reports from RPAS pilots. Twenty-three experienced RPAS pilots volunteered to participate in focus groups in which they described critical incidents from their own experience. Participants were asked to recall (1) incidents that revealed a system flaw, or (2) highlighted a case where the human operator contributed to system resilience or mission success. Participants were asked to only report incidents that could be included in a public document. A total of 90 incidents were reported. Human factor issues included the impact of reduced sensory cues, traffic separation in the absence of an out-the-window view, control latencies, vigilance during monotonous and ultra-long endurance flights, control station design considerations, transfer of control between control stations, the management of lost link procedures, and decision-making during emergencies.

human factors↗
NASA NTRS2016

Stair-Step Particle Flux Spectra on the Lunar Surface: Evidence for Nonmonotonic Potentials?

We present examples of unusual "stair‐step" differential flux spectra observed by the Apollo 14 Suprathermal Ion Detector Experiment on the lunar dayside surface in Earth's magnetotail. These spectra exhibit a relatively constant differential flux below some cutoff energy and then drop off precipitously, by about an order of magnitude or more, at higher energies. We propose that these spectra result from photoions accelerated on the lunar dayside by nonmonotonic potentials (i.e.,potentials that do not decay to zero monotonically) and present a model for the expected differential flux. The energy of the cutoff and the magnitude of the differential flux are related to the properties of the local space environment and are consistent with the observed flux spectra. If this interpretation is correct, these surface‐based ion observations provide a unique perspective that both complements and enhances the conclusions obtained by remote‐sensing orbiter observations on the Moon's exospheric and electrostatic properties.

Stair-step↗
NASA NTRS2016

Formaldehyde Column Density Measurements as a Suitable Pathway to Estimate Near-Surface Ozone Tendencies from Space

In support of future satellite missions that aim to address the current shortcomings in measuring air quality from space, NASA's Deriving Information on Surface Conditions from Column and Vertically Resolved Observations Relevant to Air Quality (DISCOVER-AQ) field campaign was designed to enable exploration of relationships between column measurements of trace species relevant to air quality at high spatial and temporal resolution. In the DISCOVER-AQ data set, a modest correlation (r2 = 0.45) between ozone (O3) and formaldehyde (CH2O) column densities was observed. Further analysis revealed regional variability in the O3-CH2O relationship, with Maryland having a strong relationship when data were viewed temporally and Houston having a strong relationship when data were viewed spatially. These differences in regional behavior are attributed to differences in volatile organic compound (VOC) emissions. In Maryland, biogenic VOCs were responsible for approx.28% of CH2O formation within the boundary layer column, causing CH2O to, in general, increase monotonically throughout the day. In Houston, persistent anthropogenic emissions dominated the local hydrocarbon environment, and no discernable diurnal trend in CH2O was observed. Box model simulations suggested that ambient CH2O mixing ratios have a weak diurnal trend (+/-20% throughout the day) due to photochemical effects, and that larger diurnal trends are associated with changes in hydrocarbon precursors. Finally, mathematical relationships were developed from first principles and were able to replicate the different behaviors seen in Maryland and Houston. While studies would be necessary to validate these results and determine the regional applicability of the O3-CH2O relationship, the results presented here provide compelling insight into the ability of future satellite missions to aid in monitoring near-surface air quality.

Surface Conditions from Column and Vertically Reso↗
NASA NTRS2016

K-Shell Photoabsorption and Photoionisation of Trace Elements I. Isoelectronic Sequences With Electron Number 3< or = N < or = 11

Context. With the recent launching of the Hitomi X-ray space observatory, K lines and edges of chemical elements with low cosmic abundances, namely F, Na, P, Cl, K, Sc, Ti, V, Cr, Mn, Co, Cu and Zn, can be resolved and used to determine important properties of supernova remnants, galaxy clusters and accreting black holes and neutron stars.Aims. The second stage of the present ongoing project involves the computation of the accurate photoabsorption and photoionisation cross sections required to interpret the X-ray spectra of such trace elements.Methods. Depending on target complexity and computer tractability, ground-state cross sections are computed either with the close-coupling Breit-Pauli R-matrix method or with the autostructure atomic structure code in the isolated-resonance approximation. The intermediate-coupling scheme is used whenever possible. In order to determine a realistic K-edge behaviour for each species, both radiative and Auger dampings are taken into account, the latter being included in the R-matrix formalism by means of an optical potential.Results. Photoabsorption and total and partial photoionisation cross sections are reported for isoelectronic sequences with electron numbers 3< or = N< or = 11. The Na sequence (N=11) is used to estimate the contributions from configurations with a 2s hole (i.e. [2s]) and those containing 3d orbitals, which will be crucial when considering sequences with N 11.Conclusions. It is found that the [2s/u] configurations must be included in the target representations of species with N> 11 as they contribute significantly to the monotonic background of the cross section between the L and K edges. Configurations with 3d orbitals are important in rendering an accurate L edge, but they can be practically neglected in the K-edge region.

atomic data / X-rays: general↗
NASA NTRS2016

Observations of the Minor Species Al and Fe in Mercury's Exosphere

We report here on the first observational evidence of Al and Fe in the exosphere of Mercury, based on measurements of resolved emission lines of these metals with Keck-1/HIRES. Al emission was observed on two separate runs, in 2008 and 2013, with tangent column densities of 3.1 +/- 1.0 and 4.0 +/-1.5 x 10(exp 7) Al atoms cm(exp - 2) at altitudes of 1185 and 1870 km (1.5 and 1.75 R(sub M). The Al radiative intensity was seen to increase where the slit crossed the planetary penumbral shadow, and then decrease monotonically with altitude. Fe emission has been observed once, in 2009, indicating an extended source. We also present observed 3- Sigma Ca(+) upper limits near Mercury's equatorial anti-solar limb, from which an abundance limit of 4.0 x 10(exp 6) cm(exp -2) at 1650 km altitude is derived for the Ca ion. A simple model for zenith column abundances of the neutral species yields 1.9 -5.2 x 10(exp 7) Al cm(exp -2) , and 8.2 x 10(exp 8) Fe cm(exp -2) . The observations appear to be consistent with production of these species by impact vaporization, with a large fraction of the Al ejecta in molecular form, and that for Fe in mixed atomic and molecular forms. The scale height of the Al gas is consistent with a kinetic temperature of 6100-8000 K. The apparent high temperature and low density of the Al gas would suggest that it may be produced by dissociation of molecules.

Mercury↗
NASA NTRS2016

Compact, Lightweight, Ceramic Matrix Composite (CMC) Based Acoustic Liners for Reducing Subsonic Jet Aircraft Engine Noise

Recent developments have reduced fan and jet noise contributions to overall subsonic aircraft jet-engine noise. Now, aircraft designers are turning their attention toward reducing engine core noise. The NASA Glenn Research Center and NASA Langley Research Center have teamed to investigate the development of a compact, lightweight acoustic liner based on oxide/oxide ceramic matrix composite (CMC) materials. The NASA team has built upon an existing oxide/oxide CMC sandwich structure concept that provides monotonal noise reduction. Oxide/oxide composites have good high temperature strength and oxidation resistance, which could allow them to perform as core liners at temperatures up to 1000C (1832F), and even higher depending on the selection of the composite constituents. NASA has initiated the evaluation of CMC-based liners that use cells of different lengths (variable-depth channels) or effective lengths to achieve broadband noise reduction. Reducing the overall liner thickness is also a major goal, to minimize the volume occupied by the liner. As a first step toward demonstrating the feasibility of our concepts, an oxide/oxide CMC acoustic testing article with different channel lengths was tested. Our approach, summary of test results, current status, and goals for the future are reported.

testing↗
NASA NTRS2017

On High-Order Upwind Methods for Advection

In the fourth installment of the celebrated series of five papers entitled "Towards the ultimate conservative difference scheme", Van Leer (1977) introduced five schemes for advection, the first three are piecewise linear, and the last two, piecewise parabolic. Among the five, scheme I, which is the least accurate, extends with relative ease to systems of equations in multiple dimensions. As a result, it became the most popular and is widely known as the MUSCL scheme (monotone upstream-centered schemes for conservation laws). Schemes III and V have the same accuracy, are the most accurate, and are closely related to current high-order methods. Scheme III uses a piecewise linear approximation that is discontinuous across cells, and can be considered as a precursor of the discontinuous Galerkin methods. Scheme V employs a piecewise quadratic approximation that is, as opposed to the case of scheme III, continuous across cells. This method is the basis for the on-going "active flux scheme" developed by Roe and collaborators. Here, schemes III and V are shown to be equivalent in the sense that they yield identical (reconstructed) solutions, provided the initial condition for scheme III is defined from that of scheme V in a manner dependent on the CFL number. This equivalence is counter intuitive since it is generally believed that piecewise linear and piecewise parabolic methods cannot produce the same solutions due to their different degrees of approximation. The finding also shows a key connection between the approaches of discontinuous and continuous polynomial approximations. In addition to the discussed equivalence, a framework using both projection and interpolation that extends schemes III and V into a single family of high-order schemes is introduced. For these high-order extensions, it is demonstrated via Fourier analysis that schemes with the same number of degrees of freedom 𝐾 per cell, in spite of the different piecewise polynomial degrees, share the same sets of eigenvalues and thus, have the same stability and accuracy. Moreover, these schemes are accurate to order 2𝐾−1, which is higher than the expected order of 𝐾.

high-order methods↗
NASA NTRS2017

Mechanism of Supercooled Water Droplet Breakup near the Leading Edge of an Airfoil

This work presents the results of an experimental study on supercooled droplet deformation and breakup near the leading edge of an airfoil. The results are compared to prior room temperature droplet deformation results to explore the effects of droplet supercooling. The experiments were conducted in the Adverse Environment Rotor Test Stand (AERTS) at The Pennsylvania State University. An airfoil model placed at the end of the rotor blades mounted onto the hub in the AERTS chamber was moved at speeds ranging between 50 and 80 m/sec. The temperature of the chamber was set at -20˚C. A monotonic droplet generator was used to produce droplets that fell from above, perpendicular to the path of the airfoil. The supercooled state of the droplets was determined by measurement of the temperature of the drops at various locations below the droplet generator exit. A temperature prediction code was also used to estimate the temperature of the droplets based on vertical velocity and the distance traveled by droplets from the droplet generator to the airfoil stagnation line. High speed imaging was employed to observe the interaction between the droplets and the airfoil. The high speed imaging provided droplet deformation information as the droplet approached the airfoil near the stagnation line. A tracking software program was used to measure the horizontal and vertical displacement of the droplet against time. It was demonstrated that to compare the effects of water supercooling on droplet deformation, the ratio of the slip velocity and the initial droplet velocity must be equal. A case with equal slip velocity to initial velocity ratios was selected for room temperature and supercooled droplet conditions. The airfoil velocity was 60 m/s and the slip velocity for both sets of data was 40 m/s. In these cases, the deformation of the weakly supercooled and warm droplets did not present different trends. The similar behavior for both environmental conditions indicates that water supercooling has no effect on particle deformation for the limited range of the weak supercooling of the droplets tested and the selected impact velocity. The assumption of a constant surface tension value was further supported by the equal trend of the Bond number obtained for supercooled and room temperature droplets.

in-flight icing↗
NASA NTRS2017

Two-Pendulum Model of Propellant Slosh in Europa Clipper PMD Tank

The objective of this fluids analysis is to model propellant slosh for the Europa Clipper mission using a two-pendulum model, such that controls engineers can predict slosh behavior during the mission. Propellant slosh causes shifts in center of mass and exerts forces and torques on the spacecraft which, if not adequately controlled, can lead to mission failure. The two-pendulum model provides a computationally simple model that can be used to predict slosh for the Europa Clipper tank geometry. The Europa Clipper tank is cylindrical with a domed top and bottom and includes a propellant management device (PMD). Due to the lack of experimental data in low gravity environments, computational fluid dynamics (CFD) simulation results were used as 'real' slosh behavior for two propellants at three fill fractions. Key pendulum parameters were derived that allow the pendulum model's center of mass, forces, and moments to closely match the CFD data. The parameter trends were examined as a function of tank fill fraction and compared with solutions to analytic equations that describe the frequency of slosh in tanks with simple geometries. The trends were monotonic as expected, and parameters resembled analytical predictions; any differences could be explained by the specific differences in the geometry of the tank. This paper summarizes the new method developed at Goddard Space Flight Center (GSFC) for deriving pendulum parameters for two-pendulum equivalent sloshing models. It presents the results of this method and discusses the validity of the results. This analysis is at a completed stage and will be applied in the immediate future to the evolving tank geometry as Europa Clipper moves past its preliminary design review (PDR) phase.

Propellant Slosh↗
NASA NTRS2017

Flexural Fatigue Behavior of an EBC CMC Composite System In Air and Steam at High Temperature

Both coated and uncoated SiCSiC ceramic matrix composite (CMC) samples were tested in flexure under sustained peak low cycle fatigue (SPLCF) conditions in air or steam at elevated temperatures. The SiCSiC composites were reinforced with 2-D plies of boron nitride coated Hi-Nicalon Type-S SiC fibers which were woven as 5 harness satin (5HS) cloth. The composites were densified by chemical vapor infiltration (CVI) followed by slurry melt infiltration (SMI). A multilayer barium strontium aluminosilicate (BSAS) coating was applied to the samples by a plasma spray method. Fatigue loading limits were determined from monotonic flexure tests at room temperature and 1200oC. Stress levels under the proportional limit of the composite material were selected for the SPLCF tests. After cyclic testing, the composites were evaluated to determine crack propagation and failure modes in the coated and uncoated composites. Microstructural examination was used to identify coating degradation and failure modes of the EBCCMC system.

Steam testing↗
NASA NTRS2017

The Nanograv Nine-Year Data Set: Measurement and Analysis of Variations in Dispersion Measures

We analyze dispersion measure(DM) variations of 37 millisecond pulsars in the nine-year North American Nanohertz Observatory for Gravitational Waves (NANOGrav) data release and constrain the sources of these variations. DM variations can result from a changing distance between Earth and the pulsar, inhomogeneities in the interstellar medium, and solar effects. Variations are significant for nearly all pulsars, with characteristic timescales comparable to or even shorter than the average spacing between observations. Five pulsars have periodic annual variations, 14 pulsars have monotonically increasing or decreasing trends, and 14 pulsars show both effects. Of the four pulsars with linear trends that have line-of-sight velocity measurements, three are consistent with a changing distance and require an overdensity of free electrons local to the pulsar. Several pulsars show correlations between DM excesses and lines of sight that pass close to the Sun. Mapping of the DM variations as a function of the pulsar trajectory can identify localized interstellar medium features and, in one case, an upper limit to the size of the dispersing region of 4 au. Four pulsars show roughly Kolmogorov structure functions (SFs), and another four show SFs less steep than Kolmogorov. One pulsar has too large an uncertainty to allow comparisons. We discuss explanations for apparent departures from a Kolmogorov-like spectrum, and we show that the presence of other trends and localized features or gradients in the interstellar medium is the most likely cause.

pulsars: general↗
NASA NTRS2017

In-Situ Observations of Longitudinal Compression Damage in Carbon-Epoxy Cross Ply Laminates Using Fast Synchrotron Radiation Computed Tomography

The role of longitudinal compressive failure mechanisms in notched cross-ply laminates is studied experimentally with in-situ synchrotron radiation based computed tomography. Carbon/epoxy specimens loaded monotonically in uniaxial compression exhibited a quasi-stable failure process, which was captured with computed tomography scans recorded continuously with a temporal resolutions of 2.4 seconds and a spatial resolution of 1.1 microns per voxel. A detailed chronology of the initiation and propagation of longitudinal matrix splitting cracks, in-plane and out-of-plane kink bands, shear-driven fiber failure, delamination, and transverse matrix cracks is provided with a focus on kink bands as the dominant failure mechanism. An automatic segmentation procedure is developed to identify the boundary surfaces of a kink band. The segmentation procedure enables 3-dimensional visualization of the kink band and conveys the orientation, inclination, and spatial variation of the kink band. The kink band inclination and length are examined using the segmented data revealing tunneling and spatial variations not apparent from studying the 2-dimensional section data.

Bergan, Andrew C.↗
NASA NTRS2017

Gravity-Assist Trajectories to the Ice Giants: An Automated Method to Catalog Mass- Or Time-Optimal Solutions

This work presents an automated method of calculating mass (or time) optimal gravity-assist trajectories without a priori knowledge of the flyby-body combination. Since gravity assists are particularly crucial for reaching the outer Solar System, we use the Ice Giants, Uranus and Neptune, as example destinations for this work. Catalogs are also provided that list the most attractive trajectories found over launch dates ranging from 2024 to 2038. The tool developed to implement this method, called the Python EMTG Automated Trade Study Application (PEATSA), iteratively runs the Evolutionary Mission Trajectory Generator (EMTG), a NASA Goddard Space Flight Center in-house trajectory optimization tool. EMTG finds gravity-assist trajectories with impulsive maneuvers using a multiple-shooting structure along with stochastic methods (such as monotonic basin hopping) and may be run with or without an initial guess provided. PEATSA runs instances of EMTG in parallel over a grid of launch dates. After each set of runs completes, the best results within a neighborhood of launch dates are used to seed all other cases in that neighborhood-allowing the solutions across the range of launch dates to improve over each iteration. The results here are compared against trajectories found using a grid-search technique, and PEATSA is found to outperform the grid-search results for most launch years considered.

Hughes, Kyle M.↗