Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “differential evolution”

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

Evolutionary Computational Methods for Identifying Emergent Behavior in Autonomous Systems

A technique based on Evolutionary Computational Methods (ECMs) was developed that allows for the automated optimization of complex computationally modeled systems, such as autonomous systems. The primary technology, which enables the ECM to find optimal solutions in complex search spaces, derives from evolutionary algorithms such as the genetic algorithm and differential evolution. These methods are based on biological processes, particularly genetics, and define an iterative process that evolves parameter sets into an optimum. Evolutionary computation is a method that operates on a population of existing computational-based engineering models (or simulators) and competes them using biologically inspired genetic operators on large parallel cluster computers. The result is the ability to automatically find design optimizations and trades, and thereby greatly amplify the role of the system engineer.

Terrile, Richard J.↗

Experiments with a Parallel Multi-Objective Evolutionary Algorithm for Scheduling

Evolutionary multi-objective algorithms have great potential for scheduling in those situations where tradeoffs among competing objectives represent a key requirement. One challenge, however, is runtime performance, as a consequence of evolving not just a single schedule, but an entire population, while attempting to sample the Pareto frontier as accurately and uniformly as possible. The growing availability of multi-core processors in end user workstations, and even laptops, has raised the question of the extent to which such hardware can be used to speed up evolutionary algorithms. In this paper we report on early experiments in parallelizing a Generalized Differential Evolution (GDE) algorithm for scheduling long-range activities on NASA's Deep Space Network. Initial results show that significant speedups can be achieved, but that performance does not necessarily improve as more cores are utilized. We describe our preliminary results and some initial suggestions from parallelizing the GDE algorithm. Directions for future work are outlined.

scheduling↗

Evaluation of a Multizone Impedance Eduction Method

A computational study is used to evaluate the PyCHE impedance eduction method developed at the NASA Langley Research Center. This method combines an aeroacoustic duct propagation code based on numerical solution to the convected Helmholtz equation with a global optimizer that uses the Differential Evolution algorithm. The efficacy of this method is evaluated with acoustic pressure data simulated to represent that measured with one-zone, two-zone, and three-zone liners mounted in the NASA Langley Grazing Flow Impedance Tube. The PyCHE method has a normalized impedance error of approximately 0.2 for (uniform) one-zone liners with a length of at least 5”, and produces quite reasonable results for liners as short as 2”. Whereas the impedance of the liner has an effect on eduction accuracy, the amount of attenuation is shown to be the dominant parameter. Similar results are observed for two-zone liners, for which the impedance of each zone is unique. The two-zone results also indicate it is more difficult to accurately educe resistance than reactance, and a zone length of at least 6” (slightly longer than for uniform liners) is needed to limit the normalized error to 0.2. The PyCHE method is also demonstrated to successfully educe the impedances for each zone of a three-zone liner. These results are sufficiently encouraging to warrant the continued usage of the PyCHE impedance eduction method for single and multizone liners.

Jones, M. G.↗

SCExAO/CHARIS Near-infrared Integral Field Spectroscopy of the HD 15115 Debris Disk

We present new, near-infrared (1.1–2.4 μm) high-contrast imaging of the debris disk around HD 15115 with the Subaru Coronagraphic Extreme Adaptive Optics (SCExAO) system coupled with the Coronagraphic High Angular Resolution Imaging Spectrograph (CHARIS). The SCExAO/CHARIS resolves the disk down to ρ ∼ 0."2 (r(proj) ∼ 10 au), a factor of ∼3–5 smaller than previous recent studies. We derive a disk position angle of PA ∼ 279.°4–280.°5 and an inclination of i ∼ 85.°3–86.2.°. While recent SPHERE/IRDIS imagery of the system could suggest a significantly misaligned two-ring disk geometry, CHARIS imagery does not reveal conclusive evidence for this hypothesis. Moreover, optimizing models of both one- and two-ring geometries using differential evolution, we find that a single ring having a Hong-like scattering phase function matches the data equally well within the CHARIS field of view (ρ ≲ 1''). The disk's asymmetry, well evidenced at larger separations, is also recovered; the west side of the disk appears, on average, around 0.4 mag brighter across the CHARIS bandpass between 0farcs25 and 1''. Comparing Space Telescope Imaging Spectrograph (STIS) 50CCD optical photometry (2000–10500 Å) with CHARIS near-infrared photometry, we find a red (STIS/50CCD−CHARIS broadband) color for both sides of the disk throughout the 0farcs4–1'' region of overlap, in contrast to the blue color reported at similar wavelengths for regions exterior to ∼2''. Further, this color may suggest a smaller minimum grain size than previously estimated at larger separations. Finally, we provide constraints on planetary companions and discuss possible mechanisms for the observed inner disk flux asymmetry and color.

Kellen Lawson↗

Progress on Inverse Estimation Technique of Non-Linear Pitch Damping Coefficient Curves Using Free-Flight CFD Generated Trajectories

Characterization of entry vehicle pitch damping coefficient curves is crucial to ensure appropriate re-entry and overall mission success. The pitch damping coefficient (C_(m_q )+C_(m_α ̇ )) is used to encapsulate the oscillatory growth or decay of a body during a trajectory. The inverse estimation technique utilizes an existing Free-Flight CFD (FF-CFD) dataset and wraps a reconstruction algorithm in an optimizer. The reconstruction integrates the planar equations of motion derived by Schoenenberger, Queen [1] using Python’s scipy.integrate.solve_ivp. The optimizer’s objective function is the normalized 𝐿2 residual of the angle of attack peaks between the reconstructed trajectory and the original data produced with FF-CFD. Inclusion of the peak times in this residual calculation allows for simultaneous optimization of the pitch moment coefficient, C_(m_α ). This residual equation is shown below in Eq. 1. The optimizer scipy.optimize.minimize was used with the gradient-based Powell method for the analysis presented, however the differential evolution method was investigated as means of comparison, and was found to produce marginally lower residual values with prohibitively longer run times. Further, the pitch damping curve is found by fitting a cubic interpolation function to a set of (α, (C_(m_q )+C_(m_α ̇ ))) control points, where the α points are held constant and the (C_(m_q )+C_(m_α ̇ )) values are the optimized parameters. The pitch moment curve uses a linear interpolation between the minimum and maximum α in the dataset. FF-CFD generated trajectories using the Dragonfly capsule geometry with the Genesis ballistic range model parameters were simulated and used for this analysis. These FF-CFD trajectories simulate planar motion, as restricted by the reconstructing the equations of motion, of three different cases: 1-DoF (free-to-pitch), 2-DoF (free-to-pitch and heave), and 3-DoF (free-to-pitch, heave, and decelerate). Pitch damping coefficient curves generated using this inverse estimation curve technique with FF-CFD 1-DoF Dragonfly data are found in Fig. 1. Preliminary results reconstructing ballistic range shots using these FF-CFD derived predictions of the pitch damping curve (Fig. 1) are shown in Fig. 2. It should be noted that the ballistic range shot used a Genesis model whereas the FF-CFD data used a Dragonfly geometry, however these geometries are similar.

entry↗

Implementation of Combinatorial Optimization Techniques for Automated Fiber Placement Through Thickness Defect Stack-Up Minimization

The Computer Aided Process Planning (CAPP) module was developed to facilitate and accelerate the process planning workflow for Automated Fiber Placement (AFP). CAPP assists process planners in identifying optimal starting point locations and layup strategies for each ply of a laminate. Ply optimization operates on measurement and scoring of geometry-based defects such as gaps, overlaps, angle deviation, and steering. This paper expands on the established framework for analyzing defect stack-up through thickness of a laminate. Four different combinatorial optimization algorithms are implemented and evaluated: (1) genetic algorithm, (2) differential evolution, (3) particle swarm, and (4) greedy search. The algorithms identify the optimal combination of ply-level layup strategies, by scoring potential laminates on defect stacking, using two different objective functions. A final optimization approach is also presented which trades some performance for a large gain in efficiency. These approaches are compared to a randomized combination using a complex tool surface in a virtual case study. The result is a streamlined methodology for comparing different laminate-level manufacturing strategies and minimizing the through thickness defect stack up.

CAPP↗

LUNA: LUT-Based Neural Architecture for Fast and Low-Cost Qubit Readout

Qubit readout is a critical operation in quantum computing systems, which maps the analog response of qubits into discrete classical states. Deep neural networks (DNNs) have recently emerged as a promising solution to improve readout accuracy . Prior hardware implementations of DNN-based readout are resource-intensive and suffer from high inference latency, limiting their practical use in low-latency decoding and quantum error correction (QEC) loops. This paper proposes LUNA, a fast and efficient superconducting qubit readout accelerator that combines low-cost integrator-based preprocessing with Look-Up Table (LUT) based neural networks for classification. The architecture uses simple integrators for dimensionality reduction with minimal hardware overhead, and employs LogicNets (DNNs synthesized into LUT logic) to drastically reduce resource usage while enabling ultra-low-latency inference. We integrate this with a differential evolution based exploration and optimization framework to identify high-quality design points. Our results show up to a 10.95x reduction in area and 30% lower latency with little to no loss in fidelity compared to the state-of-the-art. LUNA enables scalable, low-footprint, and high-speed qubit readout, supporting the development of larger and more reliable quantum computing systems.

Farooq, M. A. [Arizona State U., Tempe]↗

Measurement-informed Dynamic Aggregation of Distribution Systems

This paper proposes a measurement-informed dynamic aggregation methodology in order to create equivalent representations of distribution systems that are compatible with large-scale transmission analysis. By optimizing an equivalent feeder parameters using time-series measurements of active power, reactive power, and voltage at the Point of Interconnection (POI), the approach yields simplified yet dynamically accurate equivalents. Implemented in PSCAD with models of photovoltaic–battery systems, three-phase motors, and static loads, the method employs hybrid differential evolution and bounded least-squares optimization laying the foundation for for real-time state estimation and optimized sensor placement in distribution networks.

Ahmed, Kazi Ishrak [University of Tennessee, Knoxv↗

Evolution of a terrestrial magma ocean: Thermodynamics, kinetics, rheology, convection, differentiation

The evolution of an initially totally molten magma ocean is constrained on the basis of analysis of various physical problems in the magma ocean. First of all an equilibrium thermodynamics of the magma ocean is developed in the melting temperature range. The equilibrium thermodynamical parameters are found as functions only of temperature and pressure and are used in the subsequent models of kinetics and convection. Kinematic processes determine the crystal size and also determine a non-equilibrium thermodynamics of the system. Rheology controls all dynamical regimes of the magma ocean. The thermal convection models for different rheological laws are developed for both the laminar convection and for turbulent convection in the case of equilibrium thermodynamics of the multiphase system. The evolution is estimated on the basis of all the above analysis.

Solomatov, V. S.↗

The evolution of partially differentiated planetesimals Evidence from iron meteorite groups IAB and IIICD

Some of the properties of IAB and IIICD iron meteorites thought to be derived from partially differentiated planetesimals are summarized, and the physical aspects that may have controlled parent body differentiation and affected the composition of the sulfide melt are outlined. The chemical evolution of the parent body is then discussed, and observations supporting the partial differentiation model are examined. Finally, an attempt is made to reinterpret barometric and chronometric data in light of the partial differentiation model, and tentative conclusions are presented.

Kracher, A.↗

Thermal evolution of a differentiated Ganymede and implications for surface features

Thermodynamic models are developed for the processes which controlled the evolution of the surface Ganymede, an icy Jovian satellite assumed to have a rock-rich core surrounded by a water-ice mantle. Account is taken of a heat pulse which would have arisen from a Rayleigh-Taylor instability at a deep-seated liquid-solid water interface, rapid fracturing from global stresses imposed by warm ice diapiric upwelling, impacts by large meteorites, and resurfacing by ice flows (rather than core formation). Comparisons are made with existing models for the evolution of Callisto, and the difficulties in defining a mechanism which produced the groove terrain of Ganymede are discussed.

Kirk, R. L.↗

Sm-Nd constraints on early lunar differentiation and the evolution of KREEP

The internal Sm-Nd systematics of lunar KREEP basalt 15386 and its implications for the evolutionary history of KREEP on the moon are reported. The age of the 15386 rocklet indicated by the measured Sm-Nd mineral isochron is consistent with the reported Rb-Sr and (Ar-39)-(Ar-40) ages of other KREEP-rich samples, which is interpreted as the time 15386 crystallized from a liquid on or near the lunar surface. This suggests that the eruption of these incompatible element-rich liquids was related to deep impact events during the final bombardment phase of the surface of the moon. It is concluded that the extremely narrow spread of Sm-Nd model ages for these samples around 4.36 AE, and the compatibility of this age with that indicated by the U-Pb and Rb-Sr systems, indicate that the source of later KREEP volcanism was produced in the closing stages of an early global scale lunar differentiation episode.

Carlson, R. W.↗

Thermal evolution and chemical differentiation of the terrestrial magma ocean

The release of gravitational energy resulted in global melting and formation of a magma ocean during accretion of the Earth. Although it is believed that the formation of the magma ocean resulted in gravitational differentiation of melt and solid, the differentiation might be disturbed by the following processes: (1) convective mixing; (2) cooling and solidification; and (3) growth of the earth, which results in secular increase of pressure, and stirring by planetesimal impacts. The purpose of this study is to investigate the differentiation processes of the terrestrial magma ocean by taking into account various disturbing processes.

Abe, Y.↗

Analysis of planetary evolution with emphasis on differentiation and dynamics

In order to address the early stages of nebula evolution, a three-dimensional collapse code which includes not only hydrodynamics and radiative transfer, but also the effects of ionization and, possibly, magnetic fields is being addressed. As part of the examination of solar system evolution, an N-body code was developed which describes the latter stages of planet formation from the accretion of planetesimals. To test the code for accuracy and run-time efficiency, and to develop a stronger theoretical foundation, problems were studied in orbital dynamics. A regional analysis of the correlation in the gravity and topography fields of Venus was performed in order to determine the small and intermediate scale subsurface structure.

Kaula, William M.↗

Application of Aerosol Hygroscopicity Measured at the Atmospheric Radiation Measurement Program's Southern Great Plains Site to Examine Composition and Evolution

A Differential Mobility Analyzer/Tandem Differential Mobility Analyzer (DMA/TDMA) was used to measure submicron aerosol size distributions, hygroscopicity, and occasionally volatility during the May 2003 Aerosol Intensive Operational Period (IOP) at the Central Facility of the Atmospheric Radiation Measurement Program's Southern Great Plains (ARM SGP) site. Hygroscopic growth factor distributions for particles at eight dry diameters ranging from 0.012 micrometers to 0.600 micrometers were measured throughout the study. For a subset of particle sizes, more detailed measurements were occasionally made in which the relative humidity or temperature to which the aerosol was exposed was varied over a wide range. These measurements, in conjunction with backtrajectory clustering, were used to infer aerosol composition and to gain insight into the processes responsible for evolution. The hygroscopic growth of both the smallest and largest particles analyzed was typically less than that of particles with dry diameters of about 0.100 micrometers. It is speculated that condensation of secondary organic aerosol on nucleation mode particles is largely responsible for the minimal hygroscopic growth observed at the smallest sizes considered. Growth factor distributions of the largest particles characterized typically contained a nonhygroscopic mode believed to be composed primarily of dust. A model was developed to characterize the hygroscopic properties of particles within a size distribution mode through analysis of the fixed size hygroscopic growth measurements. The performance of this model was quantified through comparison of the measured fixed size hygroscopic growth factor distributions with those simulated through convolution of the size-resolved concentration contributed by each of the size modes and the mode-resolved hygroscopicity. This transformation from sizeresolved hygroscopicity to mode-resolved hygroscopicity facilitated examination of changes in the hygroscopic properties of particles within a size distribution mode that accompanied changes in the sizes of those particles. This model was used to examine three specific cases in which the sampled aerosol evolved slowly over a period of hours or days.

Gasparini, Roberto↗

Ulysses near-ecliptic observations of differential flow between protons and alphas in the solar wind

The evolution of differential streaming between protons and alpha particles in the solar wind was observed with the solar wind plasma experiment on the Ulysses spacecraft over the solar range of 1.15 to 5.40 AU between November 18, 1990, and May 5, 1992. The correlation of the difference in ion speeds, Delta V = the absolute value of V(sub alpha) - the absolute value of V(sub p), with the proton speed V(sub p) observed by other spacecraft at solar distances less than or equal to 1 AU disappeared at approximately 2 AU. At solar distances greater than or equal to 2.85 AU, the largest values of both V(sub alpha p) = the absolute value of V(sub alpha p) = the absolute value of V(sub alpha) - V(sub p) and the absolute value of Delta V were found in the interaction regions on the leading edges of high-speed streams. The differential streaming was typically enhanced just downstream of strong forward and reverse shocks, and large negative values of Delta V were frequently encountered in the interaction regions. A correlation between V(sub alpha p) and the ratio tau(sub zero)/tau(sub e) of Coulomb collision time to expansion time was observed at all distances, but it is suggested that at the larger values of tau(sub zero)/tau(sub e) observed correlation may arise from enhanced production of differential streaming by processes that also increase the entropy of the solar wind protons.

Neugebauer, M.↗