Search NASA⌕ Search

SEARCH · Search NASA

Results for “distributed optimization”

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

Influence of the Perception Constraint, Installation Effects, and Broadband SelfNoise on Low-Fidelity, Multi-point, UAM Vehicle Perception-Influenced-Design(PID) Optimization

UAM vehicles represent a unique challenge in acoustics because there are a limited number of full-scale datasets that can be used to calibrate and validate prediction models. In addition, acoustics is considered more often in the early design loop as a higher priority concern in UAM vehicle design, whereas with traditional vehicles, acoustics is typically addressed with post-design constraints. This paper will present an approach in design optimization of a proprotor in two flight conditions, hover, and forward flight, with several acoustic objective functions that simulate community response. An adjoint-capable vortex-lattice method (VLM), implemented in VSPAERO, will be used to compute the blade forces and inflow properties required to compute tonal and broadband noise using the second-generation Aircraft NOise Prediction Program (ANOPP2). In addition, a wing is placed near the proprotor to simulate installation effects on a vehicle’s wing. The influence of rotation rate, chord distribution, twist angle, and sweep on installation, broadband noise, and choice of perception constraint from the design optimization will be shown.

acoustics↗

Solid Modeling of Crew Exploration Vehicle Structure Concepts for Mass Optimization

Parametric solid and surface models of the crew exploration vehicle (CEV) command module (CM) structure concepts are developed for rapid finite element analyses, structural sizing and estimation of optimal structural mass. The effects of the structural configuration and critical design parameters on the stress distribution are visualized, examined to arrive at an efficient design. The CM structural components consisted of the outer heat shield, inner pressurized crew cabin, ring bulkhead and spars. For this study only the internal cabin pressure load case is considered. Component stress, deflection, margins of safety and mass are used as design goodness criteria. The design scenario is explored by changing the component thickness parameters and materials until an acceptable design is achieved. Aluminum alloy, titanium alloy and an advanced composite material properties are considered for the stress analysis and the results are compared as a part of lessons learned and to build up a structural component sizing knowledge base for the future CEV technology support. This independent structural analysis and the design scenario based optimization process may also facilitate better CM structural definition and rapid prototyping.

Mukhopadhyay, Vivek↗

Effects of Planetary Speed Reduction Ratio on Mean Service Life

Planetary gear transmissions are compact speed reductions which use parallel-load paths to transmit high power. The range of possible reduction ratios is bounded from below and above by limits on the relative size of the planet gears. For a single plane planetary transmission, the speed-reduction ratio must be greater than two for planet gears to exist. As the ratio increase, so does the size of the planets relative to the sizes of the sun and ring. Which ratio is best for a planetary reduction can be resolved by studying a series of optimal designs. In this series, each design is obtained by maximizing the service life for a spur-gear planetary with a fixed size, speed reduction ratio, input speed, power and materials. The planetary gear reduction service life is modeled as a function of the two-parameter Weibull distributed service lives of the ball bearings and spur gears in the reduction. Planet-bearing life strongly influences the optimal reduction lives which point to an optimal planetary reduction ratio in the neighborhood of four to five.

Savage, M.↗

Investigating Global Ion and Neutral Atom Populations with IBEX and Voyager

The main objective of this project was to investigate pickup ion (PUI) production in the solar wind and heliosheath (the region between the termination shock and the heliopause) and compute the distributed energetic neutral atom fluxes throughout the helioshpere. The simulations were constrained by comparing the model output against observations from Ulysses, New Horizons, Voyager 1 and 2, and IBEX space probes. As evidenced by the number of peer reviewed journal publications resulting from the project (13 plus three submitted) and their citation rate (156 citations over three years), the project has made a lasting contribution to the field. The outcome is a significant improvement of our understanding of the pickup ion production and distribution in the distant heliosphere. The team has accomplished the entire set of tasks A-H set forth in the proposal. Namely, the transport modeling framework has been augmented with two populations of pickup ions (PUIs), the boundary conditions for the plasma and interstellar neutral hydrogen were verified against Ulysses and New Horizons PUI and an optimal set of velocity diffusion parameters established. The multi-component fluxes of PUIs were computed and isotropic velocity distributions generated for each cell in the computer simulation that covered the heliosphere from 1.5 AU to the heliopause. The distributions were carefully compared with in situ measurements at 3 AU (Ulysses), 12 AU (New Horizons), and 80-90 AU (Voyager 1 and 2) as well as those inferred from ENA fluxes measured by Cassini and IBEX (Wu et al., 2016). Some examples of modeldata comparison are shown in Figure 1. We have used coupled MHD-plasma and kinetic-neutral code to investigate the likely range of plasma and magnetic field parameters in the local interstellar medium (LISM), based on the assumption that the shape of the IBEX ribbon could be used to determine the orientation of the interstellar magnetic field. While the magnetic field is believed to be oriented toward the center of the ribbon, constraining its strength requires comparing the model-predicted angular diameter and circularity of the ribbon with the observations. The study, published in Heerikhuisen et al. (2014), found that the most likely range for the LISM magnetic field strength is between 0.2 and 0.3 nT, which is less than previously thought. Figure 2 shows the IBEX data (left) and compares it to the simulation with a 0.2 nT interstellar magnetic field (center) and a 0.4 nT (right).

Florinski, Vladimir↗

Communication Needs Assessment for Distributed Turbine Engine Control

Control system architecture is a major contributor to future propulsion engine performance enhancement and life cycle cost reduction. The control system architecture can be a means to effect net weight reduction in future engine systems, provide a streamlined approach to system design and implementation, and enable new opportunities for performance optimization and increased awareness about system health. The transition from a centralized, point-to-point analog control topology to a modular, networked, distributed system is paramount to extracting these system improvements. However, distributed engine control systems are only possible through the successful design and implementation of a suitable communication system. In a networked system, understanding the data flow between control elements is a fundamental requirement for specifying the communication architecture which, itself, is dependent on the functional capability of electronics in the engine environment. This paper presents an assessment of the communication needs for distributed control using strawman designs and relates how system design decisions relate to overall goals as we progress from the baseline centralized architecture, through partially distributed and fully distributed control systems.

Culley, Dennis E.↗

A Comparative Study on Computation of Cumulative Distribution Function in Predicting Time of Failure of Engineering Systems

Estimating accurate Time-of-Failure (ToF) of a system is key in making the decisions that impact operational safety and optimize cost. In this context, it is interesting to note that different approaches have been explored to tackle the problem of estimating ToF. The difference is in part characterized by different definitions of the hazard zones. The conventional definition for the cumulative distribution function (CDF) calculation is assumed to have well-defined hazard zones, that is, hazard zones defined as a function of the system state trajectory. An alternate method suggests the use of hazard zones defined as a function of the system state at time , instead of hazard zones defined as a function of system state up to and including time k (Acuna and Orchard 2018, 2017). This paper explores these differences and their impact on ToF estimation. Results for the conventional CDF definition indicated that, (i) the cumulative distribution function is always an increasing function of time, even when realizations of the degradation process are not monotonic, (ii) the sum of all probabilities is always 1 and does not need to be normalized, and (iii) all probabilities are positive and less than or equal to 1. Similar results are not observed for CDF calculation with hazard zones defined as a function only of the system state at time . Results for ToF estimation using Acuna's definition differ, suggesting that there is an underlying assumption of independence in the hazard zone definition. Therefore, we present an alternate definition of hazard zone which guarantees the properties of a well-defined CDF with a more straightforward ToF definition.

Hazard Zone↗

ISS-Crystal Growth of Photorefractive Materials (BSO): Critical Design Issues for Optimized Data Extraction from Space Experiments

Realization of the full potential of photorefractive materials in device technology is seriously impeded by our inability to achieve controlled formation of critical defects during single crystal growth and by difficulties in meeting the required degree of compositional uniformity on a micro-scale over macroscopic dimensions. The exact nature and origin of the critical defects which control photorefractivity could not as yet be identified because of gravitational interference. There exists, however, strong evidence that the density of defect formation and their spatial distribution are adversely affected by gravitational interference which precludes the establishment of quantifiable and controllable heat and mass transfer conditions during crystal growth. The current, NASA sponsored research at MIT is directed at establishing a basis for the development of a comprehensive approach to the optimization of property control during melt growth of photorefractive materials, making use of the m-g environment, provided in the International Space Station. The objectives to be pursued in m-g research on photorefractive BSO (Bi12SiO20) are: (a) identification of the x-level(s) responsible for photorefractivity in undoped BSO; (b) development of approaches leading to the control of x-level formation at uniform spatial distribution; (c) development of doping and processing procedures for optimization of the critical, application specific parameters, spectral response, sensitivity, response time and matrix stability. The presentation will focus on: the rationale for the justification of the space experiment, ground-based development efforts, design considerations for the space experiments, strategic plan of the space experiments, and approaches to the quantitative analysis of the space experiments.

Hyers, Robert W.↗

Effects of Planetary Speed-Reduction Ratio on Mean Service Life

Planetary-gear transmissions are compact speed reductions which use parallel-load paths to transmit high power. The range of possible reduction ratios is bounded from below and above by limits on the relative size of the planet gears. For a single-plane planetary transmission, the speed-reduction ratio must be greater than two for planet gears to exist. As the ratio increases, so does the size of the planets relative to the sizes of the sun and ring. Which ratio is best for a planetary reduction can be resolved by studying a series of optimal designs. In this series, each design is obtained by maximizing the service life for a spur-gear planetary with a fixed size, speed-reduction ratio, input speed, power and materials. The planetary-gear reduction service life is modeled as a function of the two-parameter Weibull distributed service lives of the ball bearings and spur gears in the reduction. Planet-bearing life strongly influences the optimal reduction lives which point to an optimal planetary reduction ratio in the neighborhood of four to five.

Savage, M.↗

Multidisciplinary High-Fidelity Analysis and Optimization of Aerospace Vehicles: Preliminary Results - Part 2

An objective of the High Performance Computing and Communication Program at the NASA Langley Research Center is to demonstrate multidisciplinary shape and sizing optimization of a complete aerospace vehicle configuration by using high-fidelity finite-element structural analysis and computational fluid dynamics aerodynamic analysis in a distributed, heterogeneous computing environment that includes high performance parallel computing. A software system has been designed and implemented to integrate a set of existing discipline analysis codes, some of them computationally intensive, into a distributed computational environment for the design of a high-speed civil transport configuration. The paper describes both the preliminary results from implementing and validating the multidisciplinary analysis and the results from an aerodynamic optimization. The discipline codes are integrated by using the Java programming language and a Common Object Request Broker Architecture compliant software product. A companion paper describes the formulation of the multidisciplinary analysis and optimization system.

Walsh, J. L.↗

Optimization of Training Sets for Neural-Net Processing of Characteristic Patterns from Vibrating Solids

Artificial neural networks have been used for a number of years to process holography-generated characteristic patterns of vibrating structures. This technology depends critically on the selection and the conditioning of the training sets. A scaling operation called folding is discussed for conditioning training sets optimally for training feed-forward neural networks to process characteristic fringe patterns. Folding allows feed-forward nets to be trained easily to detect damage-induced vibration-displacement-distribution changes as small as 10 nm. A specific application to aerospace of neural-net processing of characteristic patterns is presented to motivate the conditioning and optimization effort.

Decker, Arthur J.↗

Three-Dimensional Venturi Sensor for Measuring Extreme Winds

A three-dimensional (3D) Venturi sensor is being developed as a compact, rugged means of measuring wind vectors having magnitudes of as much as 300 mph (134 m/s). This sensor also incorporates auxiliary sensors for measuring temperature from -40 to +120 F (-40 to +49 C), relative humidity from 0 to 100 percent, and atmospheric pressure from 846 to 1,084 millibar (85 to 108 kPa). Conventional cup-and-vane anemometers are highly susceptible to damage by both high wind forces and debris, due to their moving parts and large profiles. In addition, they exhibit slow recovery times contributing to an inaccurately high average-speed reading. Ultrasonic and hot-wire anemometers overcome some of the disadvantages of the cup and-vane anemometers, but they have other disadvantageous features, including limited dynamic range and susceptibility to errors caused by external acoustic noise and rain. In contrast, the novel 3D Venturi sensor is less vulnerable to wind damage because of its smaller profile and ruggedness. Since the sensor has no moving parts, it provides increased reliability and lower maintenance costs. It has faster response and recovery times to changing wind conditions than traditional systems. In addition, it offers wide dynamic range and is expected to be relatively insensitive to rain and acoustic energy. The Venturi effect in this sensor is achieved by the mirrored double-inflection curve, which is then rotated 360 to create the desired detection surfaces. The curve is optimized to provide a good balance of pressure difference between sensor ports and overall maximum fluid velocity while in the shape. Four posts are used to separate the two shapes, and their size and location were chosen to minimize effects on the pressure measurements. The 3D Venturi sensor has smart software algorithms to map the wind pressure exerted on the surfaces of the design. Using Bernoulli's equation, the speed of the wind is calculated from the differences among the pressure readings at the various ports. The direction of the wind is calculated from the spatial distribution and magnitude of the pressure readings. All of the pressure port sizes and locations have been optimized to minimize measurement errors and to reside in areas demonstrating a stable pressure reading proportional to the velocity range.

Zysko, Jan A.↗

Coupled Aeropropulsive Design Optimization of a Podded Electric Propulsor

New aircraft concepts are increasingly relying on non-traditional propulsion systems to achieve lower energy consumption. These non-traditional technologies, such as boundary layer ingestion or distributed electric propulsion, require a tight integration of the propulsion sys- tem to the airframe, and therefore, a traditional design approach where the aerodynamics and propulsion disciplines are considered separately is likely to result in sub-optimal designs. As a result, we have to rely on coupled aeropropulsive design optimization, in which fully coupled aeropropulsive models are optimized to maximize the advantages of these tightly integrated propulsion systems. Despite its advantages, aeropropulsive design optimization is a relatively new field, and significant advancements are required for wide adoption of this approach, especially in the context of robustness. In this work, we introduce two fully coupled aeropropulsive design optimization approaches that are compatible with CFD models with body-force terms and boundary conditions to model the effects of the propulsion system in the flow field. To test these new approaches, we developed a podded electric fan model based on the aft-propulsor of NASA’s STARC-ABL concept. This simple design problem enables us to rapidly develop and test different aeropropulsive coupling approaches, while including the challenging physical interactions that arise from coupling CFD and propulsion models. The coupled aeropropulsive model is implemented using the MPhys library, which is built with NASA’s OpenMDAO frame- work. Using this benchmark design problem, we demonstrate the robustness of the approaches and our aeropropulsive design optimization framework by performing a sweep of optimizations for a total of 50 CFD-based aeropropulsive design optimizations. Furthermore, the new aero- propulsive coupling approaches enable multi-point design optimizations. We demonstrate this capability in a multi-point design optimization problem where we optimize cruise performance subject to a fan-face distortion constraint at rolling take-off conditions. The aeropropulsive modeling approaches we present in this work represent a significant milestone in the field of aeropropulsive design optimization. The framework we developed is extremely robust and flexible, and these developments will be crucial for future aeropropulsive design optimizations of complete turbofan engines.

Optimization propulsion aerodynamics↗

An Optimal Estimation Aerosol Retrieval Algorithm for Lidar-only and Combined Lidar Polarimeter Remote Sensors

Global measurements of aerosol vertical profile, composition, concentration, and size distribution are very important due to aerosol impacts on air quality, climate, clouds, and ocean ecosystems. We have developed a lidar-only and combined lidar and polarimeter algorithm to retrieve vertically-resolved profiles of aerosol microphysical properties. The new retrieval system is modular in design. It consists of three modules, a vertically-resolved aerosol profile retrieval module for lidar data, an aerosol and cloud/ocean retrieval module for polarimeter data, and a combined retrieval module for both lidar and polarimeter measurements. In addition to performing optimal estimation retrievals on various data sources, we have designed the system so that it can be used to carry out aerosol retrieval performance trade studies for various lidar and polarimeter configurations. For example, the retrieval system can take inputs from attenuated backscatter lidar at two wavelengths (e.g. CALIPSO like instrument), aerosol backscattering measured at two wavelengths and and aerosol extinction measured at one wavelength (HSRL-1), or aerosol backscattering measured at three wavelengths and aerosol extinction measured at two wavelengths (HSRL-1) (HSRL-2). We have applied the lidar-only algorithm to both simulated data and various field campaign data (DISCOVER-AQ, CHARMS, TCAP, and ORACLES). For polarimeter-only modules, we have done the same for SABOR, TCAP, NAAMES, and ORACLES data. We have applied the combined retrieval algorithm to ORACLES data. The simulated retrieval studies show that the combined lidar+polarimeter retrieval provides much higher information content relative to their individual counterparts for retrieving effective radius, particle concentrations, and absorption properties.

Xu Liu↗

Development and Application of the Collaborative Optimization Architecture in a Multidisciplinary Design Environment

Collaborative optimization is a design architecture applicable in any multidisciplinary analysis environment but specifically intended for large-scale distributed analysis applications. In this approach, a complex problem is hierarchically de- composed along disciplinary boundaries into a number of subproblems which are brought into multidisciplinary agreement by a system-level coordination process. When applied to problems in a multidisciplinary design environment, this scheme has several advantages over traditional solution strategies. These advantageous features include reducing the amount of information transferred between disciplines, the removal of large iteration-loops, allowing the use of different subspace optimizers among the various analysis groups, an analysis framework which is easily parallelized and can operate on heterogenous equipment, and a structural framework that is well-suited for conventional disciplinary organizations. In this article, the collaborative architecture is developed and its mathematical foundation is presented. An example application is also presented which highlights the potential of this method for use in large-scale design applications.

Braun, R. D.↗

Tuning Monotonic Basin Hopping: Improving the Efficiency of Stochastic Search as Applied to Low-Thrust Trajectory Optimization

Trajectory optimization methods using MBH have become well developed during the past decade. An essential component of MBH is a controlled random search through the multi-dimensional space of possible solutions. Historically, the randomness has been generated by drawing RVs from a uniform probability distribution. Here, we investigate the generating the randomness by drawing the RVs from Cauchy and Pareto distributions, chosen because of their characteristic long tails. We demonstrate that using Cauchy distributions (as first suggested by Englander significantly improves MBH performance, and that Pareto distributions provide even greater improvements. Improved performance is defined in terms of efficiency and robustness, where efficiency is finding better solutions in less time, and robustness is efficiency that is undiminished by (a) the boundary conditions and internal constraints of the optimization problem being solved, and (b) by variations in the parameters of the probability distribution. Robustness is important for achieving performance improvements that are not problem specific. In this work we show that the performance improvements are the result of how these long-tailed distributions enable MBH to search the solution space faster and more thoroughly. In developing this explanation, we use the concepts of sub-diffusive, normally-diffusive, and super-diffusive RWs originally developed in the field of statistical physics.

mission design↗

Ames Optimized TCA Configuration

Configuration design at Ames was carried out with the SYN87-SB (single block) Euler code using a 193 x 49 x 65 C-H grid. The Euler solver is coupled to the constrained (NPSOL) and the unconstrained (QNMDIF) optimization packages. Since the single block grid is able to model only wing-body configurations, the nacelle/diverter effects were included in the optimization process by SYN87's option to superimpose the nacelle/diverter interference pressures on the wing. These interference pressures were calculated using the AIRPLANE code. AIRPLANE is an Euler solver that uses a unstructured tetrahedral mesh and is capable of computations about arbitrary complete configurations. In addition, the buoyancy effects of the nacelle/diverters were also included in the design process by imposing the pressure field obtained during the design process onto the triangulated surfaces of the nacelle/diverter mesh generated by AIRPLANE. The interference pressures and nacelle buoyancy effects are added to the final forces after each flow field calculation. Full details of the (recently enhanced) ghost nacelle capability are given in a related talk. The pseudo nacelle corrections were greatly improved during this design cycle. During the Ref H and Cycle 1 design activities, the nacelles were only translated and pitched. In the cycle 2 design effort the nacelles can translate vertically, and pitch to accommodate the changes in the lower surface geometry. The diverter heights (between their leading and trailing edges) were modified during design as the shape of the lower wing changed, with the drag of the diverter changing accordingly. Both adjoint and finite difference gradients were used during optimization. The adjoint-based gradients were found to give good direction in the design space for configurations near the starting point, but as the design approached a minimum, the finite difference gradients were found to be more accurate. Use of finite difference gradients was limited by the CPU time limit available on the Cray machines. A typical optimization run using finite difference gradients can use only 30 to 40 design variables and one optimization iteration within the 8 hour queue limit for the chosen grid size and convergence level. The efficiency afforded by the adjoint method allowed for 50-120 design variables and 5-10 optimization iterations in the 8 hour queue. Geometric perturbations to the wing and fuselage were made using the Hicks/Henne (HH) shape functions. The HH functions were distributed uniformly along the chords of the wing defining sections and lofted linearly. During single-surface design, constraints on thickness and volume at selected wing stations were imposed. Both fuselage camber and cross-sectional area distributions were permitted to change during design. The major disadvantage to the use of these functions is the inherent surface waviness produced by repeated use of such functions. Many smoothing operations were required following optimization runs to produce a configuration with reasonable smoothness. Wagner functions were also used on the wing sections but were never used on the fuselage. The Wagner functions are a family of increasingly oscillatory functions that have also been used extensively in airfoil design. The leading and trailing edge regions of the wing were designed by use of polynomial and monomial functions respectively. Twist was attempted but was abandoned because of little performance improvement available from changing the baseline twist.

Cliff, Susan E.↗

Approach for Uncertainty Propagation and Robust Design in CFD Using Sensitivity Derivatives

This paper presents an implementation of the approximate statistical moment method for uncertainty propagation and robust optimization for a quasi 1-D Euler CFD (computational fluid dynamics) code. Given uncertainties in statistically independent, random, normally distributed input variables, a first- and second-order statistical moment matching procedure is performed to approximate the uncertainty in the CFD output. Efficient calculation of both first- and second-order sensitivity derivatives is required. In order to assess the validity of the approximations, the moments are compared with statistical moments generated through Monte Carlo simulations. The uncertainties in the CFD input variables are also incorporated into a robust optimization procedure. For this optimization, statistical moments involving first-order sensitivity derivatives appear in the objective function and system constraints. Second-order sensitivity derivatives are used in a gradient-based search to successfully execute a robust optimization. The approximate methods used throughout the analyses are found to be valid when considering robustness about input parameter mean values.

Putko, Michele M.↗

Real time thermal imaging of high temperature semiconductor melts

A real time thermal imaging system with temperature resolution better than + or - 1 C and spatial resolution of better than 0.5 mm was developed and applied to the analysis of melt surface thermal field distributions in both Czochralski and liquid encapsulated Czochralski (LEC) growth configurations. The melt is viewed in near normal incidence by a high resolution charge coupled device camera to which is attached a very narrow bandpass filter. The resulting image is digitized and processed using a pipelined pixel processor operating at an effective 40 million operations per second thus permitting real time high frequency spatial and temporal filtering of the high temperature scene. A multi-pixel averaging algorithm was developed which permits localized, low noise sensing of temperature variations at any location in the hot zone as a function of time. This signial is used to implement initial elements of a feedforward growth control scheme which is aimed at reducing disturbances to the melt caused by the batch nature of the growth process. The effect of magnetic melt stabilization on radial melt temperature distributions was measured using this technique. Problems associated with residual internal reflections and non-optimized path geometry are discussed.

Wargo, Michael J.↗