Search NASA⌕ Search

SEARCH · Search NASA

Results for “projected gradient method”

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 73 records · Page 4

Investigation of Navier-Stokes Code Verification and Design Optimization

With rapid progress made in employing computational techniques for various complex Navier-Stokes fluid flow problems, design optimization problems traditionally based on empirical formulations and experiments are now being addressed with the aid of computational fluid dynamics (CFD). To be able to carry out an effective CFD-based optimization study, it is essential that the uncertainty and appropriate confidence limits of the CFD solutions be quantified over the chosen design space. The present dissertation investigates the issues related to code verification, surrogate model-based optimization and sensitivity evaluation. For Navier-Stokes (NS) CFD code verification a least square extrapolation (LSE) method is assessed. This method projects numerically computed NS solutions from multiple, coarser base grids onto a freer grid and improves solution accuracy by minimizing the residual of the discretized NS equations over the projected grid. In this dissertation, the finite volume (FV) formulation is focused on. The interplay between the xi concepts and the outcome of LSE, and the effects of solution gradients and singularities, nonlinear physics, and coupling of flow variables on the effectiveness of LSE are investigated. A CFD-based design optimization of a single element liquid rocket injector is conducted with surrogate models developed using response surface methodology (RSM) based on CFD solutions. The computational model consists of the NS equations, finite rate chemistry, and the k-6 turbulence closure. With the aid of these surrogate models, sensitivity and trade-off analyses are carried out for the injector design whose geometry (hydrogen flow angle, hydrogen and oxygen flow areas and oxygen post tip thickness) is optimized to attain desirable goals in performance (combustion length) and life/survivability (the maximum temperatures on the oxidizer post tip and injector face and a combustion chamber wall temperature). A preliminary multi-objective optimization study is carried out using a geometric mean approach. Following this, sensitivity analyses with the aid of variance-based non-parametric approach and partial correlation coefficients are conducted using data available from surrogate models of the objectives and the multi-objective optima to identify the contribution of the design variables to the objective variability and to analyze the variability of the design variables and the objectives. In summary the present dissertation offers insight into an improved coarse to fine grid extrapolation technique for Navier-Stokes computations and also suggests tools for a designer to conduct design optimization study and related sensitivity analyses for a given design problem.

Vaidyanathan, Rajkumar↗

Using Digital Image Correlation for Material Testing at Cryogenic Temperatures

Recent efforts were made to parametrize an aluminum alloy at cryogenic temperatures to calibrate finite element models simulating material deformation. Digital image correlation (DIC) was used to measure material properties using front/back strains of specimens during tensile tests. Data acquired from DIC was compared to physical extensometer measurements. Two camera pairs (front and back) were mounted via aluminum extrusion to the exterior of an environmental chamber of a hydraulic test stand. Images were taken through single-pane optical glass windows of the environment chamber. Positive pressure bags encased the cameras, attaching to the perimeter of the glass to keep out moisture and prevent windows from frosting. Figure 1 shows the environmental setup and the camera setup during cryogenic testing. Some of the challenges during testing at cryogenic temperatures include the gaseous flow of liquid nitrogen in the optical path of the mounted cameras, frost buildup on the windows as seen in Figure 1b, and frost on the test coupon. The temperature gradient present during the introduction of liquid nitrogen into the chamber can increase the projection error of images captured of the specimen. Figure 2 shows the estimated projection error from DIC as a function of temperature during the cool down process. This presentation addresses some of the concerns and observations when using DIC for testing at cryogenic temperatures using methods of isolated camera system and techniques for maintaining the glass panel clear from frost.

Joseph Cochran↗

Design of a lunar transportation system

The establishment of lunar bases is the next logical step in the exploration of space. Permanent lunar bases will support scientific investigation, the industrialization of space, and the development of self-sufficiency on the Moon. Scientific investigation and research and development would lead to applications utilizing lunar material resources. By utilizing these resources, the industrialization of space can become a reality. The above two factors coupled with the development of key and enabling technologies would lead to achievement of self-sufficiency of the lunar base. Attention was focused on specific design(s) to be pursued during subsequent stages in advanced courses. Some of the objectives in the project included: (1) minimizing the transportation of construction material and fuel from earth, or maximizing the use of the lunar material; (2) use of novel materials and light weight structures; (3) use of new manufacturing methods and technology such as magnetically levitated, or superconducting materials; and (4) innovative concepts of effectively utilizing the exotic lunar conditions, i.e. high thermal gradients, lack of atmosphere, zero wind forces, and lower gravity, etc.

Source record↗

Function-space quasi-Newton algorithms for optimal control problems with bounded controls and singular arcs

Two existing function-space quasi-Newton algorithms, the Davidon algorithm and the projected gradient algorithm, are modified so that they may handle directly control-variable inequality constraints. A third quasi-Newton-type algorithm, developed by Broyden, is extended to optimal control problems. The Broyden algorithm is further modified so that it may handle directly control-variable inequality constraints. From a computational viewpoint, dyadic operator implementation of quasi-Newton methods is shown to be superior to the integral kernel representation. The quasi-Newton methods, along with the steepest descent method and two conjugate gradient algorithms, are simulated on three relatively simple (yet representative) bounded control problems, two of which possess singular subarcs. Overall, the Broyden algorithm was found to be superior. The most notable result of the simulations was the clear superiority of the Broyden and Davidon algorithms in producing a sharp singular control subarc.

Edge, E. R.↗

Electromagnetic characterization of conformal antennas

The ultimate objective of this project is to develop a new technique which permits an accurate simulation of microstrip patch antennas or arrays with various feed, superstrate and/or substrate configurations residing in a recessed cavity whose aperture is planar, cylindrical or otherwise conformed to the substructure. The technique combines the finite element and boundary integral methods to formulate a system suitable for solution via the conjugate gradient method in conjunction with the fast Fourier transform. The final code is intended to compute both scattering and radiation patterns of the structure with an affordable memory demand. With upgraded capabilities, the four included papers examined the radar cross section (RCS), input impedance, gain, and resonant frequency of several rectangular configurations using different loading and substrate/superstrate configurations.

Volakis, John L.↗

Parallel Implicit Algorithms for CFD

The main goal of this project was efficient distributed parallel and workstation cluster implementations of Newton-Krylov-Schwarz (NKS) solvers for implicit Computational Fluid Dynamics (CFD.) "Newton" refers to a quadratically convergent nonlinear iteration using gradient information based on the true residual, "Krylov" to an inner linear iteration that accesses the Jacobian matrix only through highly parallelizable sparse matrix-vector products, and "Schwarz" to a domain decomposition form of preconditioning the inner Krylov iterations with primarily neighbor-only exchange of data between the processors. Prior experience has established that Newton-Krylov methods are competitive solvers in the CFD context and that Krylov-Schwarz methods port well to distributed memory computers. The combination of the techniques into Newton-Krylov-Schwarz was implemented on 2D and 3D unstructured Euler codes on the parallel testbeds that used to be at LaRC and on several other parallel computers operated by other agencies or made available by the vendors. Early implementations were made directly in Massively Parallel Integration (MPI) with parallel solvers we adapted from legacy NASA codes and enhanced for full NKS functionality. Later implementations were made in the framework of the PETSC library from Argonne National Laboratory, which now includes pseudo-transient continuation Newton-Krylov-Schwarz solver capability (as a result of demands we made upon PETSC during our early porting experiences). A secondary project pursued with funding from this contract was parallel implicit solvers in acoustics, specifically in the Helmholtz formulation. A 2D acoustic inverse problem has been solved in parallel within the PETSC framework.

Keyes, David E.↗

An Active Subspace Method for Accelerating Convergence in Delaunay-Based Optimization via Dimension Reduction

Delaunay-based derivative-free optimization, ∆DOGS, is an efficient and provably-convergent global optimization method for the problems which has computationally expensive objection function and the analytical expression for the objective function is not available. ∆-DOGS is a novel optimization scheme in the family of response surface methods (RSMs); however, it suffers from the curse of dimensionality since the computational cost increases dramatically as the number of design parameters increases. As a result, the number of design parameters in ∆-DOGS algorithm is relatively low (n.10). To avoid such problems, this paper proposes a combination of derivative-free optimization, seeking the global minimizer of an expensive and nonconvex objective function f(x) and active subspace method, detecting the directions of the most variability using evaluations of the gradient. The contribution of other directions to the objective function is bounded by a sufficiently small constant. This new algorithm iteratively applied Delaunay-based derivative-free optimization to seek the minimizer on the d-dimensional active subspace that has most function variation. Inverse mapping is needed to project data from active subspace to full-model for evaluating function values. This task is overcome by solving an inequality constrained problem that curves the response surface of the objective function. The test results show that this strategy is effective on a handful of optimization problems.

Bewley, Thomas R.↗

Advancement of the General Aviation Synthesis Program Using Python to Enable Optimization-Based Hybrid-Propulsion Aircraft Design

In support of the Electrified Powertrain Flight Demonstrator and Advanced Air Transport Technologies projects at NASA, a new tool has been developed at NASA’s Ames and Glenn Research Centers to enable coupled engine and airframe optimization and analysis. The new tool combines the engineering-level analysis methods and empirical models of the FORTRAN General Aviation Synthesis Program (GASP) with the Python-based OpenMDAO framework to provide a modular framework for efficient gradient-based optimization with the aim of incorporating new subsystem models for unconventional configurations. The tool has been verified against GASP analyses of several aircraft models and mission formulations. Preliminary efforts have been made to integrate pyCycle, a thermodynamic cycle analysis tool, to enable simultaneous optimization of hybrid propulsion system and vehicle parameters while taking full mission performance and constraints into account. This will improve current capabilities to assess impacts of electrified powertrain technologies on future aircraft designs.

Kenneth R. Lyons↗

Advancement of the General Aviation Synthesis Program Using Python to Enable Optimization-Based Hybrid-Propulsion Aircraft Design

In support of the Electrified Powertrain Flight Demonstrator and Advanced Air Transport Technologies projects at NASA, a new tool has been developed at NASA's Ames and Glenn Research Centers to enable coupled engine and airframe optimization and analysis. The new tool combines the engineering-level analysis methods and empirical models of the FORTRAN General Aviation Synthesis Program (GASP) with the Python-based OpenMDAO framework to provide a modular framework for efficient gradient-based optimization with the aim of incorporating new subsystem models for unconventional configurations. The tool has been verified against GASP analyses of several aircraft models and mission formulations. Preliminary efforts have been made to integrate pyCycle, a thermodynamic cycle analysis tool, to enable simultaneous optimization of hybrid propulsion system and vehicle parameters while taking full mission performance and constraints into account. This will improve current capabilities to assess impacts of electrified powertrain technologies on future aircraft designs.

Kenneth R Lyons↗

Early results of the ionospheric experiment of the Apollo-Soyuz Test Project

A description is presented of a spacecraft-to-spacecraft Doppler-tracking experiment which was performed by the Smithsonian Astrophysical Observatory on the occasion of the Apollo-Soyuz Test Project (ASTP). The experiment involved the measurement of the relative velocity between the ASTP docking module and the Apollo command service module by a Doppler-tracking method. The objectives of the ionospheric experiment include the measurement of the time changes of the columnar electron content between the two spacecraft. The obtained data can provide a basis for the determination of the horizontal gradients of electron density at the height of 220 km.

Grossi, M. D.↗

Relaxation and Preconditioning for High Order Discontinuous Galerkin Methods with Applications to Aeroacoustics and High Speed Flows

This project is about the investigation of the development of the discontinuous Galerkin finite element methods, for general geometry and triangulations, for solving convection dominated problems, with applications to aeroacoustics. Other related issues in high order WENO finite difference and finite volume methods have also been investigated. methods are two classes of high order, high resolution methods suitable for convection dominated simulations with possible discontinuous or sharp gradient solutions. In [18], we first review these two classes of methods, pointing out their similarities and differences in algorithm formulation, theoretical properties, implementation issues, applicability, and relative advantages. We then present some quantitative comparisons of the third order finite volume WENO methods and discontinuous Galerkin methods for a series of test problems to assess their relative merits in accuracy and CPU timing. In [3], we review the development of the Runge-Kutta discontinuous Galerkin (RKDG) methods for non-linear convection-dominated problems. These robust and accurate methods have made their way into the main stream of computational fluid dynamics and are quickly finding use in a wide variety of applications. They combine a special class of Runge-Kutta time discretizations, that allows the method to be non-linearly stable regardless of its accuracy, with a finite element space discretization by discontinuous approximations, that incorporates the ideas of numerical fluxes and slope limiters coined during the remarkable development of the high-resolution finite difference and finite volume schemes. The resulting RKDG methods are stable, high-order accurate, and highly parallelizable schemes that can easily handle complicated geometries and boundary conditions. We review the theoretical and algorithmic aspects of these methods and show several applications including nonlinear conservation laws, the compressible and incompressible Navier-Stokes equations, and Hamilton-Jacobi-like equations.

Shu, Chi-Wang↗

The control of float zone interfaces by the use of selected boundary conditions

The main goal of the float zone crystal growth project of NASA's Materials Processing in Space Program is to thoroughly understand the molten zone/freezing crystal system and all the mechanisms that govern this system. The surface boundary conditions required to give flat float zone solid melt interfaces were studied and computed. The results provide float zone furnace designers with better methods for controlling solid melt interface shapes and for computing thermal profiles and gradients. Documentation and a user's guide were provided for the computer software.

Foster, L. M.↗

Robot acting on moving bodies (RAMBO): Preliminary results

A robot system called RAMBO is being developed. It is equipped with a camera, which, given a sequence of simple tasks, can perform these tasks on a moving object. RAMBO is given a complete geometric model of the object. A low level vision module extracts and groups characteristic features in images of the object. The positions of the object are determined in a sequence of images, and a motion estimate of the object is obtained. This motion estimate is used to plan trajectories of the robot tool to relative locations nearby the object sufficient for achieving the tasks. More specifically, low level vision uses parallel algorithms for image enchancement by symmetric nearest neighbor filtering, edge detection by local gradient operators, and corner extraction by sector filtering. The object pose estimation is a Hough transform method accumulating position hypotheses obtained by matching triples of image features (corners) to triples of model features. To maximize computing speed, the estimate of the position in space of a triple of features is obtained by decomposing its perspective view into a product of rotations and a scaled orthographic projection. This allows the use of 2-D lookup tables at each stage of the decomposition. The position hypotheses for each possible match of model feature triples and image feature triples are calculated in parallel. Trajectory planning combines heuristic and dynamic programming techniques. Then trajectories are created using parametric cubic splines between initial and goal trajectories. All the parallel algorithms run on a Connection Machine CM-2 with 16K processors.

Davis, Larry S.↗

Robot Acting on Moving Bodies (RAMBO): Interaction with tumbling objects

Interaction with tumbling objects will become more common as human activities in space expand. Attempting to interact with a large complex object translating and rotating in space, a human operator using only his visual and mental capacities may not be able to estimate the object motion, plan actions or control those actions. A robot system (RAMBO) equipped with a camera, which, given a sequence of simple tasks, can perform these tasks on a tumbling object, is being developed. RAMBO is given a complete geometric model of the object. A low level vision module extracts and groups characteristic features in images of the object. The positions of the object are determined in a sequence of images, and a motion estimate of the object is obtained. This motion estimate is used to plan trajectories of the robot tool to relative locations rearby the object sufficient for achieving the tasks. More specifically, low level vision uses parallel algorithms for image enhancement by symmetric nearest neighbor filtering, edge detection by local gradient operators, and corner extraction by sector filtering. The object pose estimation is a Hough transform method accumulating position hypotheses obtained by matching triples of image features (corners) to triples of model features. To maximize computing speed, the estimate of the position in space of a triple of features is obtained by decomposing its perspective view into a product of rotations and a scaled orthographic projection. This allows use of 2-D lookup tables at each stage of the decomposition. The position hypotheses for each possible match of model feature triples and image feature triples are calculated in parallel. Trajectory planning combines heuristic and dynamic programming techniques. Then trajectories are created using dynamic interpolations between initial and goal trajectories. All the parallel algorithms run on a Connection Machine CM-2 with 16K processors.

Davis, Larry S.↗

NASA-UVA Light Aerospace Alloy and Structures Technology program (LA2ST)

The objective of the LA2ST Program is to conduct interdisciplinary graduate student research on the performance of next generation, light-weight aerospace alloys, composites, and thermal gradient structures in collaboration with NASA-Langley researchers. The general aim is to produce relevant data and basic understanding of material mechanical response, environment/corrosion behavior, and microstructure; new monolithic and composite alloys; advanced processing methods; new solid and fluid mechanics analyses; measurement and modeling advances; and a pool of educated students for aerospace technologies. Specific technical objectives are presented for each of the following research projects: time-temperature dependent fracture in advanced wrought ingot metallurgy, and spray deposited aluminum alloys; cryogenic temperature effects on the deformation and fracture of Al-Li-Cu-In alloys; effects of aging and temperature on the ductile fracture of AA2095 and AA2195; mechanisms of localized corrosion in alloys 2090 and 2095; hydrogen interactions in aluminum-lithium alloys 2090 and selected model alloys; mechanisms of deformation and fracture in high strength titanium alloys (effects of temperature and hydrogen and effects of temperature and microstructure); evaluations of wide-panel aluminum alloy extrusions; Al-Si-Ge alloy development; effects of texture and precipitates on mechanical property anisotropy of Al-Cu-Mg-X alloys; damage evolution in polymeric composites; and environmental effects in fatigue life prediction - modeling crack propagation in light aerospace alloys.

Starke, Edgar A., Jr.↗

Salton Sea Project, Phase 1

A feasibility study was made for a salt gradient solar pond power plant in or near the Salton Sea of California. The conclusions support continuance 5-MWe proof-of-concept experiment, and ultimate construction by an electric utility company of a 600-MWe plant. The Solar Pond concept would be an environmental benefit to the Salton Sea by reversing the increasing salinity trend. The greatest cost drivers are the lake dike construction and pond sealing. Problems to be resolved include method of brine production from Salton Sea water for the first unit (which requires evaporation pond area and time), the high turbidity and color content of the Salton Sea water (which requires pretreatment), and other questions related to pond permeability, bio-activity and soil/brine chemical reactions. All technical and environmental problems appear solvable and/or manageable if care is taken in mitigating impacts.

Peelgren, M. L.↗

An update on remote measurement of soil moisture over vegetation using infrared temperature measurements: A FIFE perspective

Using model development, image analysis and micrometeorological measurements, the object is to push beyond the present limitations of using the infrared temperature method for remotely determining surface energy fluxes and soil moisture over vegetation. Model development consists of three aspects: (1) a more complex vegetation formulation which is more flexible and realistic; (2) a method for modeling the fluxes over patchy vegetation cover; and (3) a method for inferring a two-layer soil vertical moisture gradient from analyses of horizontal variations in surface temperatures. HAPEX and FIFE satellite data will be used along with aircraft thermal infrared and solar images as input for the models. To test the models, moisture availability and bulk canopy resistances will be calculated from data collected locally at the Rock Springs experimental field site and, eventually, from the FIFE project.

Carlson, Toby N.↗

Shuttle ascent trajectory optimization with function space quasi-Newton techniques

A Space Shuttle ascent trajectory optimization problem from lift-off to orbital insertion is solved with a function space version of a quasi-Newton parameter optimization method developed by Broyden. The problem includes five parameter and one bounded-function controls, two state-variable constraints, and four terminal conditions. The bounded controls are treated directly, while the remaining constraints are adjoined to the performance index (maximum payload) with penalty functions. The problem is formulated as a four-phase variational problem (liftoff, pitch-over, gravity-turn, linear tangent steering), and the appropriate gradients are developed by first variation theory. A projection operator is introduced to aid in the interpretation of the algorithm with mixed parameter and function controls.

Edge, E. R.↗