Engineering PapersSearch

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 37 records · Page 2

The nonlinear evolution of inviscid Goertler vortices in three-dimensional boundary layers

The nonlinear development of inviscid Gortler vortices in a three-dimensional boundary layer is considered. We do not follow the classical approach of weakly nonlinear stability problems and consider a mode which has just become unstable. Instead we extend the method of Blackaby, Dando, and Hall (1992), which considered the closely related nonlinear development of disturbances in stratified shear flows. The Gortler modes we consider are initially fast growing and we assume, following others, that boundary-layer spreading results in them evolving in a linear fashion until they reach a stage where their amplitudes are large enough and their growth rates have diminished sufficiently so that amplitude equations can be derived using weakly nonlinear and non-equilibrium critical-layer theories. From the work of Blackaby, Dando and Hall (1993) is apparent, given the range of parameters for the Gortler problem, that there are three possible nonlinear integro-differential evolution equations for the disturbance amplitude. These are a cubic due to viscous effects, a cubic which corresponds to the novel mechanism investigated in this previous paper, and a quintic. In this paper we shall concentrate on the two cubic integro-differential equations and in particular, on the one due to the novel mechanism as this will be the first to affect a disturbance. It is found that the consideration of a spatial evolution problem as opposed to temporal (as was considered in Blackaby, Dando, and Hall, 1992) causes a number of significant changes to the evolution equations.

Blackaby, Nicholas

The fully nonlinear development of Goertler vortices in growing boundary layers

The fully nonlinear development of small wavelength Goertler vortices in a growing boundary layer is investigated using a combination of asymptotic and numerical methods. The starting point for the analysis is the weakly nonlinear theory of Hall (1982b) who discussed the initial development of small amplitude vortices in a neighborhood of the location where they first become linearly unstable. That development is unusual in the context of nonlinear stability theory in that it is not described by the Stuart-Watson approach. In fact the development is governed by a pair of coupled nonlinear partial differential evolution equations for the vortex flow and the mean flow correction. Here the further development of this interaction is considered for vortices so large that the mean flow correction driven by them is as large as the basic state. Surprisingly it is found that such a nonlinear interaction can still be described by asymptotic means. It is shown that the vortices spread out across the boundary layer and effectively drive the boundary layer. In fact the system obtained by writing down the equations for the fundamental component of the vortex generate a differential equation for the basic state. Thus the mean flow adjusts so as to make these large amplitude vortices locally neutral. Moreover in the region where the vortices exist the mean flow has a square-root profile and the vortex velocity field can be written down in closed form. The upper and lower boundaries of the region of vortex activity are determined by a free-boundary problem involving the boundary layer equations. In general it is found that this region ultimately includes almost all of the original boundary layer and much of the free-stream. In this situation the mean flow has essentially no relationship to the flow which exists in the absence of the vortices.

Hall, Philip

Linking boundary-layer circulations and surface processes during FIFE89. Part 1: Observational analysis

Surface, aircraft, and satellite observations are analyzed for the 21-day 1989 intensive field campaign of the First ISLSCP Field Experiment (FIFE) to determine the effect of precipitation, vegetation, and soil moisture distributions on the thermal properties of the surface including the heat and moisture fluxes, and the corresponding response in the boundary-layer circulation. Mean and variance properties of the surface variables are first documented at various time and space scales. These calculations are designed to set the stage for Part 2, a modeling study that will focus on how time-space dependent rainfall distribution influences the intensity of the feedback between a vegetated surface and the atmospheric boundary layer. Further analysis shows strongly demarked vegetation and soil moisture gradients extending across the FIFE experimental site that were developed and maintained by the antecedent and ongoing spatial distribution of rainfall over the region. These gradients are shown to have a pronounced influence on the thermodynamic properties of the surface. Furthermore, perturbation surface wind analysis suggests for both short-term steady-state conditions and long-term averaged conditions that the gradient pattern maintained a diurnally oscillating local direct circulation with perturbation vertical velocities of the same order as developing cumulus clouds. Dynamical and scaling considerations suggest that the embedded perturbation circulation is driven by surface heating/cooling gradients and terrain ef fects rather than the manifestation of an inertial oscillation. The implication is that at even relatively small scales (less than 30 km), the differential evolution in vegetation density and soil moisture distribution over a relatively homogenous ecotone can give rise to preferential boundary-layer circulations capable of modifying local-scale horizontal and vertical motions.

Smith, Eric A.

Evolutionary optimization of a Genetically Refined Truss

Structural optimization is a field of research that has experienced noteworthy growth for many years. Researchers in this area have developed optimization tools to successfully design and model structures, typically minimizing mass while maintaining certain deflection and stress constraints. Numerous optimization studies have been performed to minimize mass, deflection and stress on a benchmark cantilever truss problem. Predominantly traditional optimization theory is applied to this problem. The cross-sectional area of each member is optimized to minimize the aforementioned objectives. This paper will present a structural optimization technique that has been previously applied to compliant mechanism design. This technique demonstrates a method that combines topology optimization, geometric refinement, finite element analysis, and two forms of evolutionary computation: Genetic Algorithms and Differential Evolution to successfully optimize a benchmark structural optimization problem. An non-traditional solution to the benchmark problem is presented in this paper, specifically a geometrically refined topological solution. The design process begins with an alternate control mesh formulation, multilevel geometric smoothing operation, and an elastostatic structural analysis. The design process is wrapped in an evolutionary computing optimization toolset.

Hull, Patrick V.

Towards Robust Designs Via Multiple-Objective Optimization Methods

Fabricating and operating complex systems involves dealing with uncertainty in the relevant variables. In the case of aircraft, flow conditions are subject to change during operation. Efficiency and engine noise may be different from the expected values because of manufacturing tolerances and normal wear and tear. Engine components may have a shorter life than expected because of manufacturing tolerances. In spite of the important effect of operating- and manufacturing-uncertainty on the performance and expected life of the component or system, traditional aerodynamic shape optimization has focused on obtaining the best design given a set of deterministic flow conditions. Clearly it is important to both maintain near-optimal performance levels at off-design operating conditions, and, ensure that performance does not degrade appreciably when the component shape differs from the optimal shape due to manufacturing tolerances and normal wear and tear. These requirements naturally lead to the idea of robust optimal design wherein the concept of robustness to various perturbations is built into the design optimization procedure. The basic ideas involved in robust optimal design will be included in this lecture. The imposition of the additional requirement of robustness results in a multiple-objective optimization problem requiring appropriate solution procedures. Typically the costs associated with multiple-objective optimization are substantial. Therefore efficient multiple-objective optimization procedures are crucial to the rapid deployment of the principles of robust design in industry. Hence the companion set of lecture notes (Single- and Multiple-Objective Optimization with Differential Evolution and Neural Networks ) deals with methodology for solving multiple-objective Optimization problems efficiently, reliably and with little user intervention. Applications of the methodologies presented in the companion lecture to robust design will be included here. The evolutionary method (DE) is first used to solve a relatively difficult problem in extended surface heat transfer wherein optimal fin geometries are obtained for different safe operating base temperatures. The objective of maximizing the safe operating base temperature range is in direct conflict with the objective of maximizing fin heat transfer. This problem is a good example of achieving robustness in the context of changing operating conditions. The evolutionary method is then used to design a turbine airfoil; the two objectives being reduced sensitivity of the pressure distribution to small changes in the airfoil shape and the maximization of the trailing edge wedge angle with the consequent increase in airfoil thickness and strength. This is a relevant example of achieving robustness to manufacturing tolerances and wear and tear in the presence of other objectives.

Man Mohan, Rai

Evolutionary Optimization of a Geometrically Refined Truss

Structural optimization is a field of research that has experienced noteworthy growth for many years. Researchers in this area have developed optimization tools to successfully design and model structures, typically minimizing mass while maintaining certain deflection and stress constraints. Numerous optimization studies have been performed to minimize mass, deflection, and stress on a benchmark cantilever truss problem. Predominantly traditional optimization theory is applied to this problem. The cross-sectional area of each member is optimized to minimize the aforementioned objectives. This Technical Publication (TP) presents a structural optimization technique that has been previously applied to compliant mechanism design. This technique demonstrates a method that combines topology optimization, geometric refinement, finite element analysis, and two forms of evolutionary computation: genetic algorithms and differential evolution to successfully optimize a benchmark structural optimization problem. A nontraditional solution to the benchmark problem is presented in this TP, specifically a geometrically refined topological solution. The design process begins with an alternate control mesh formulation, multilevel geometric smoothing operation, and an elastostatic structural analysis. The design process is wrapped in an evolutionary computing optimization toolset.

Hull, P. V.

Steady-State ALPS for Real-Valued Problems

The two objectives of this paper are to describe a steady-state version of the Age-Layered Population Structure (ALPS) Evolutionary Algorithm (EA) and to compare it against other GAs on real-valued problems. Motivation for this work comes from our previous success in demonstrating that a generational version of ALPS greatly improves search performance on a Genetic Programming problem. In making steady-state ALPS some modifications were made to the method for calculating age and the method for moving individuals up layers. To demonstrate that ALPS works well on real-valued problems we compare it against CMA-ES and Differential Evolution (DE) on five challenging, real-valued functions and on one real-world problem. While CMA-ES and DE outperform ALPS on the two unimodal test functions, ALPS is much better on the three multimodal test problems and on the real-world problem. Further examination shows that, unlike the other GAs, ALPS maintains a genotypically diverse population throughout the entire search process. These findings strongly suggest that the ALPS paradigm is better able to avoid premature convergence then the other GAs.

Hornby, Gregory S.

Reducing the Volume of NASA Earth-Science Data

A computer program reduces data generated by NASA Earth-science missions into representative clusters characterized by centroids and membership information, thereby reducing the large volume of data to a level more amenable to analysis. The program effects an autonomous data-reduction/clustering process to produce a representative distribution and joint relationships of the data, without assuming a specific type of distribution and relationship and without resorting to domain-specific knowledge about the data. The program implements a combination of a data-reduction algorithm known as the entropy-constrained vector quantization (ECVQ) and an optimization algorithm known as the differential evolution (DE). The combination of algorithms generates the Pareto front of clustering solutions that presents the compromise between the quality of the reduced data and the degree of reduction. Similar prior data-reduction computer programs utilize only a clustering algorithm, the parameters of which are tuned manually by users. In the present program, autonomous optimization of the parameters by means of the DE supplants the manual tuning of the parameters. Thus, the program determines the best set of clustering solutions without human intervention.

Lee, Seungwon

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

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.