Search NASA⌕ Search

SEARCH · Search NASA

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Accelerating Bilevel Optimization With Hierarchical Many-Threaded Parallel Differential Evolution

Bilevel optimization is encountered in many relevant real-world applications. The main feature of this type of problem is that an upper-level optimization problem is constrained by a nested lower-level optimization problem. Because of this nested structure, bilevel problems (BLPs) are usually computationally expensive to solve. Differential evolution (DE) has demonstrated promising results in solving BLPs of relatively small scales. As the problem scale increases, the decision space becomes intrinsically larger, requiring a growing number of function evaluations for the method to work properly. In this context, heavy parallelization and high-performance computing techniques are indispensable to enable the resolution of more complex and challenging optimization problems. Hence, we propose a hierarchical many-threaded parallel DE approach for BLPs, where both levels are parallelized. The computational experiments demonstrate that the parallel implementation achieved runtime speeds ranging from 44 to 2559 times faster than the sequential version on a well-known scalable SMD benchmark test problem when executed on an NVIDIA A100 GPU. The findings indicate that the algorithm’s convergence is strongly influenced by the number of both upper- and lower-level generations. Moreover, the success of experiments with large-scale problems is closely linked to the choice of small population sizes.

Dufek, Amanda S↗
DOE OSTI2017

Solving Inverse Transport Problems with Neutron Multiplication Measurements and Improved Differential Evolution

In the problem of inverse radiation transport, measurements of particle leakages from radioactive source/shield systems are used to infer unknown parameters within the systems. This reconstruction can be accomplished by finding the physical parameters of the unknown system that minimize the difference between calculated detector responses and measured detector responses. Recently, the differential evolution (DE) method was shown to be highly adept at solving inverse transport problems in which the measured quantities were leakages of unscattered gamma-ray lines. In this work, we present an extension of the DE method to problems in which both gamma-ray lines and neutron multiplication are measured. The neutron multiplication measurement is demonstrated to be an important addition when few gamma-ray lines are available. We also discuss the implementation of a variant of DE that significantly reduces the number of forward transport computations required to solve the inverse problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Identification of Stellar Flares Using Differential Evolution Template Optimization

In this paper, we explore methods for the identification of stellar flare events in irregularly sampled data of ground-based time domain surveys. In particular, we describe a new technique for identifying flaring stars, which we have implemented in a publicly available Python module called "PyVAN." The approach uses the Differential Evolution algorithm to optimize parameters of empirically derived light curve templates for different types of stars to fit a candidate light curve. The difference of the likelihoods that these best-fit templates produced the observed data is then used to delineate targets that are well-explained by a flare template but simultaneously poorly explained by templates of common contaminants. By testing on light curves of known identity and morphology, we show that our technique is capable of recovering flaring status in 69% of all light curves containing a flare event above thresholds drawn to include <1% of any contaminant population. By applying to Palomar Transient Factory data, we show consistency with prior samples of flaring stars, and identify a small selection of candidate flaring G-type stars for possible follow-up.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2006

Single- and Multiple-Objective Optimization with Differential Evolution and Neural Networks

Genetic and evolutionary algorithms have been applied to solve numerous problems in engineering design where they have been used primarily as optimization procedures. These methods have an advantage over conventional gradient-based search procedures became they are capable of finding global optima of multi-modal functions and searching design spaces with disjoint feasible regions. They are also robust in the presence of noisy data. Another desirable feature of these methods is that they can efficiently use distributed and parallel computing resources since multiple function evaluations (flow simulations in aerodynamics design) can be performed simultaneously and independently on ultiple processors. For these reasons genetic and evolutionary algorithms are being used more frequently in design optimization. Examples include airfoil and wing design and compressor and turbine airfoil design. They are also finding increasing use in multiple-objective and multidisciplinary optimization. This lecture will focus on an evolutionary method that is a relatively new member to the general class of evolutionary methods called differential evolution (DE). This method is easy to use and program and it requires relatively few user-specified constants. These constants are easily determined for a wide class of problems. Fine-tuning the constants will off course yield the solution to the optimization problem at hand more rapidly. DE can be efficiently implemented on parallel computers and can be used for continuous, discrete and mixed discrete/continuous optimization problems. It does not require the objective function to be continuous and is noise tolerant. DE and applications to single and multiple-objective optimization will be included in the presentation and lecture notes. A method for aerodynamic design optimization that is based on neural networks will also be included as a part of this lecture. The method offers advantages over traditional optimization methods. It is more flexible than other methods in dealing with design in the context of both steady and unsteady flows, partial and complete data sets, combined experimental and numerical data, inclusion of various constraints and rules of thumb, and other issues that characterize the aerodynamic design process. Neural networks provide a natural framework within which a succession of numerical solutions of increasing fidelity, incorporating more realistic flow physics, can be represented and utilized for optimization. Neural networks also offer an excellent framework for multiple-objective and multi-disciplinary design optimization. Simulation tools from various disciplines can be integrated within this framework and rapid trade-off studies involving one or many disciplines can be performed. The prospect of combining neural network based optimization methods and evolutionary algorithms to obtain a hybrid method with the best properties of both methods will be included in this presentation. Achieving solution diversity and accurate convergence to the exact Pareto front in multiple objective optimization usually requires a significant computational effort with evolutionary algorithms. In this lecture we will also explore the possibility of using neural networks to obtain estimates of the Pareto optimal front using non-dominated solutions generated by DE as training data. Neural network estimators have the potential advantage of reducing the number of function evaluations required to obtain solution accuracy and diversity, thus reducing cost to design.

Rai, Man Mohan↗
DOE OSTI2020

Drive-pressure optimization in ramp-wave compression experiments through differential evolution

Ramp-wave dynamic-compression experiments are used to examine quasi-isentropic loading paths in materials. The gradual and continuous increase in pressure created by ramp waves make these types of experiments ideal for studying nonequilibrium material behavior, such as solidification kinetics. In ramp-wave compression experiments, the input drive pressure to the experimental setup may be exerted through one of a number of different mechanisms (e.g., magnetic fields, gas-gun-driven impactors, or high-energy lasers) and is generally required for simulating such experiments. Yet, regardless of the specific mechanism, this drive pressure cannot be measured directly (measurements are generally taken at a location near the back of the experimental setup through a transparent window), leading to an inverse problem where one must determine the drive pressure at the front of the experimental setup (i.e., the input) that corresponds to the particle velocity (the output) measured near the back of the experimental setup. Furthermore, we solve this inverse problem using a heuristic optimization algorithm, known as differential evolution, coupled with a multiphysics, hydrodynamics code that simulates the compression of the experimental setup. By running many rounds of forward simulations of the experimental setup, our optimization process iteratively searches for a drive pressure that is optimized to closely reproduce the experimentally measured particle velocity near the back of the experimental setup. While our optimization methodology requires a significant number of hydrodynamics simulations to be conducted, many of these can be performed in parallel, which greatly reduces the time cost of our methodology. One novel aspect of our method for determining the drive pressure is that it does not require physical modeling of the drive mechanism and can thus be broadly applied to many types of ramp-compression experiments, regardless of the drive mechanism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2012

The Ground Flash Fraction Retrieval Algorithm Employing Differential Evolution: Simulations and Applications

The ability to estimate the fraction of ground flashes in a set of flashes observed by a satellite lightning imager, such as the future GOES-R Geostationary Lightning Mapper (GLM), would likely improve operational and scientific applications (e.g., severe weather warnings, lightning nitrogen oxides studies, and global electric circuit analyses). A Bayesian inversion method, called the Ground Flash Fraction Retrieval Algorithm (GoFFRA), was recently developed for estimating the ground flash fraction. The method uses a constrained mixed exponential distribution model to describe a particular lightning optical measurement called the Maximum Group Area (MGA). To obtain the optimum model parameters (one of which is the desired ground flash fraction), a scalar function must be minimized. This minimization is difficult because of two problems: (1) Label Switching (LS), and (2) Parameter Identity Theft (PIT). The LS problem is well known in the literature on mixed exponential distributions, and the PIT problem was discovered in this study. Each problem occurs when one allows the numerical minimizer to freely roam through the parameter search space; this allows certain solution parameters to interchange roles which leads to fundamental ambiguities, and solution error. A major accomplishment of this study is that we have employed a state-of-the-art genetic-based global optimization algorithm called Differential Evolution (DE) that constrains the parameter search in such a way as to remove both the LS and PIT problems. To test the performance of the GoFFRA when DE is employed, we applied it to analyze simulated MGA datasets that we generated from known mixed exponential distributions. Moreover, we evaluated the GoFFRA/DE method by applying it to analyze actual MGAs derived from low-Earth orbiting lightning imaging sensor data; the actual MGA data were classified as either ground or cloud flash MGAs using National Lightning Detection Network[TM] (NLDN) data. Solution error plots are provided for both the simulations and actual data analyses.

Koshak, William↗
NASA NTRS1978

Lunar initial Nd-143/Nd-144 - Differential evolution of the lunar crust and mantle

The Sm-Nd evolution of Apollo 15 green glass is discussed. The ICE age (intercept with chondritic evolution) of 3.8 + or - 0.4 eons overlaps the range of reported (Ar-39)-(Ar-40) ages and implies a distinct source region for green glass, characterized by very low and unfractionated REE abundances. Evidence is presented that LINd (lunar initial Nd) is compatible with a 'chondritic'-type Nd isotopic evolution as observed in the Juvinas meteorite. This normalization is used to study the Sm-Nd system of various lunar rock types. The results obtained from a limited number of rocks clearly indicate differential Sm-Nd evolution for the lunar crust and mantle. High-Ti basalts returned by the Apollo 11 and 17 missions were derived from distinct source regions. The Nd-143 evolution in KREEP requires a source region which is clearly distinct from any mantle reservoir.

Lugmair, G. W.↗
NASA NTRS1995

Starspot evolution, differential rotation, and magnetic cycles in the chromospherically active binaries lambda andromedae, sigma Geminorum, II Pegasi, and V711 Tauri

We have analyzed 15-19 yr of photoelectric photometry, obtained manually and with automated telescopes, of the chromospherically active binaries lambda And, sigma Gem, II Peg, and V711 Tau. These observations let us identify individual dark starspots on the stellar surfaces from periodic dimming of the starlight, follow the evolution of these spots, and search for long-term cyclic changes in the properties of these starspots that might reveal magnetic cycles analogous to the Sun's 11 yr sunspot cycle. We developed a computer code to fit a simple two-spot model to our observed light curves that allows us to extract the most easily determinable and most reliable spot parameters from the light curves, i.e., spot longitudes and radii. We then used these measured properties to identify individual spots and to chart their life histories by constructing migration and amplitude curves. We identified and followed 11 spots in lambda And, 16 in sigma Gem, 12 in II Peg, and 15 in V711 Tau. Lifetimes of individual spots ranged from a few months to longer than 6 yr. Differential rotation coefficients, estimated from the observed range of spot rotation periods for each star and defined by equation (2), were 0.04 for lambda And, 0.038 for sigma Gem, 0.005 for II Peg, and 0.006 for V711 Tau, versus 0.19 for the Sun. We searched for cyclic changes in mean brightness, B-V color index, and spot rotation period as evidence for long-term cycles. Of these, long-term variability in mean brightness appears to offer the best evidence for such cycles in these four stars. Cycles of 11.1 yr for lambda And, 8.5 yr for sigma Gem, 11 yr for II Peg, and 16 yr V711 Tau are implied by these mean brightness changes. Cyclic changes in spot rotation period were found in lambda And and possibly II Peg. Errors in B-V were too large for any long-term changes to be detectable.

Henry, Gregory W.↗
DOE OSTI2013

WITNESSING THE DIFFERENTIAL EVOLUTION OF DISK GALAXIES IN LUMINOSITY AND SIZE VIA GRAVITATIONAL LENSING

We take advantage of the magnification in size and flux of a galaxy provided by gravitational lensing to analyze the properties of 62 strongly lensed galaxies from the Sloan Lens ACS (SLACS) Survey. The sample of lensed galaxies spans a redshift range of 0.20 ≤ z ≤ 1.20 with a median redshift of z = 0.61. We use the lens modeling code LENSFIT to derive the luminosities, sizes, and Sérsic indices of the lensed galaxies. The measured properties of the lensed galaxies show a primarily compact, {sup d}isk{sup -}like population with the peaks of the size and Sérsic index distributions corresponding to ∼1.50 kpc and n ∼ 1, respectively. Comparison of the SLACS galaxies to a non-lensing, broadband imaging survey shows that a lensing survey allows us to probe a galaxy population that reaches ∼2 mag fainter. Our analysis allows us to compare the (z) = 0.61 disk galaxy sample (n ≤ 2.5) to an unprecedented local galaxy sample of ∼670, 000 SDSS galaxies at z ∼ 0.1; this analysis indicates that the evolution of the luminosity-size relation since z ∼ 1 may not be fully explained by a pure-size or pure-luminosity evolution but may instead require a combination of both. Our observations are also in agreement with recent numerical simulations of disk galaxies that show evidence of a mass-dependent evolution since z ∼ 1, where high-mass disk galaxies (M{sub *} > 10{sup 9} M{sub ☉}) evolve more in size and low-mass disk galaxies (M{sub *} ≤ 10{sup 9} M{sub ☉}) evolve more in luminosity.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1995

Differential evolution of radio quiet and radio loud quasars

We have observed 40 optically selected high-redshift quasars with the Very Large Array (VLA) at 1515 MHz (20 cm). Three quasars are detected above our 3 sigma limit of 0.2 mJy. We review all available evidence on radio detections of optically selected quasars in terms of their ratio R of radio to optical luminosity. We find that the quasars ranked in the top 5-20 percent in R show strong evolution in their R distribution. We present two interpretations of this phenomenon. The first one poses that the R values decrease by a factor of around 100 from redshift 0.3 to 1.3. In the second interpretation, we consider that there are two populations of quasars, which we call radio silent and radio active, separated by R approximately 0.1. The radio-silent fraction of quasars increases from approximately 16% at redshift 0.3, to approximately 75% at redshift 1.3.

Schmidt, Maarten↗
NASA NTRS2016

Differential Evolution Optimization for Targeting Spacecraft Maneuver Plans

Previous analysis identified specific orbital parameters as being safer for conjunction avoidance for the TDRS fleet. With TDRS-9 being considered an at-risk spacecraft, a potential conjunction concern was raised should TDRS-9 fail while at a longitude of 12W. This document summarizes the analysis performed to identify if these specific orbital parameters could be targeted using the remaining drift-termination maneuvers for the relocation of TDRS-9 from 41W longitude to 12W longitude.

Manuever↗
DOE OSTI2020

LA17-V-InverseRadTransportProblem-PD2Kb. Final Project Report

In the inverse radiation transport problem, a radiation signal emanating from an unknown object is used to reconstruct geometric and material parameters of the object , e.g., the mass of special nuclear material (SNM), the isotopic composition of the SNM (e.g., the enrichment of 235U), the thickness of shielding, etc. The inverse radiation transport problem is one of the most fundamental and important in the emergency response mission, and has been addressed in previous NA-22 and NA-42/84 investments by using derivative-free optimization methods. Our project was aimed at (i) consolidating previous work by applying polynomial chaos expansion (PCE) for reducing significantly the number of forward model simulations required by derivative-free optimization methods; (ii) developing a general-purpose framework and stand-alone software for solving the inverse radiation transport problem by applying a new “predictive modeling” methodology (predictive modeling of coupled multi-physics systems, PM_CMPS) which enables the fusing of experimental and computational information; and (iii) developing several benchmarks (for unscattered gamma rays) for demonstrating the quantification of non-Gaussian features (e.g., skewness and kurtosis, which characterize asymmetries and long tails) of response distributions by deterministically and efficiently propagating parameter uncertainties (covariances and higher-order moments, where available) using high-order sensitivities computed by the recently developed comprehensive adjoint sensitivity analysis methodology (C-ASAM). With regard to (i), application of PCE, the most significant accomplishment of this project was an implementation of PCE into INVERSE that significantly reduced computation time for the DiffeRential Evolution Adaptive Metropolis (DREAM) and Differential Evolution (DE) methods. For inverse transport problems with measurements of neutron leakage, neutron multiplication, and neutron-induced gamma-rays, the solution times for the DREAM and DE methods were reduced by factors of 30-300, reducing run times from hours to minutes while maintaining the same accuracy as the standard (non-PCE) DREAM and DE. Speedups were further improved by exploring sparse grid quadrature sets using the Toolkit for Adaptive Stochastic Modeling and Non-Intrusive ApproximatioN (Tasmanian) toolset. An initial exploration of sparse grids indicated that further speedups factors of 5-10 are obtainable. The implementation of PCE in INVERSE is not in a form that users can easily access. With regard to (ii), solving the inverse problem using a predictive modeling paradigm, the most significant accomplishment of this project were the derivation and implementation of adjoint-based second derivatives in neutron and photon transport theory. There are two codes that compute these derivatives. SENSPG computes them for uncollided gamma rays and SENSMG computes them for neutron leakage. SENSMG is open-source code available on github. Other significant accomplishments included a stand-alone code, MULTI-PRED, that implements the predictive modeling paradigm, along with a detailed manual. The code was implemented in the existing INVERSE code for inverse problems, but a head-to-head comparison of inverse methods was not achieved. We also derived and implemented the capability to compute sensitivities for multiplicity measurements into the SENSMG code, and this capability was used in a data assimilation application, but it was not used for solving inverse problems. With regard to (iii), developing benchmarks to demonstrate the importance of high-order sensitivity and uncertainty analysis, the most significant accomplishment of this project were two test problems that showed how important second-order sensitivities are.

07 ISOTOPE AND RADIATION SOURCES↗