Search NASA⌕ Search

SEARCH · Search NASA

Results for “derivative free 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Optimal Low Energy Earth-Moon Transfers

The optimality of a low-energy Earth-Moon transfer is examined for the first time using primer vector theory. An optimal control problem is formed with the following free variables: the location, time, and magnitude of the transfer insertion burn, and the transfer time. A constraint is placed on the initial state of the spacecraft to bind it to a given initial orbit around a first body, and on the final state of the spacecraft to limit its Keplerian energy with respect to a second body. Optimal transfers in the system are shown to meet certain conditions placed on the primer vector and its time derivative. A two point boundary value problem containing these necessary conditions is created for use in targeting optimal transfers. The two point boundary value problem is then applied to the ballistic lunar capture problem, and an optimal trajectory is shown. Additionally, the ballistic lunar capture trajectory is examined to determine whether one or more additional impulses may improve on the cost of the transfer.

Griesemer, Paul Ricord↗
NASA NTRS1972

A survey of methods of feasible directions for the solution of optimal control problems

Three methods of feasible directions for optimal control are reviewed. These methods are an extension of the Frank-Wolfe method, a dual method devised by Pironneau and Polack, and a Zontendijk method. The categories of continuous optimal control problems are shown as: (1) fixed time problems with fixed initial state, free terminal state, and simple constraints on the control; (2) fixed time problems with inequality constraints on both the initial and the terminal state and no control constraints; (3) free time problems with inequality constraints on the initial and terminal states and simple constraints on the control; and (4) fixed time problems with inequality state space contraints and constraints on the control. The nonlinear programming algorithms are derived for each of the methods in its associated category.

Polak, E.↗
NASA NTRS2016

Effect of Stagger on the Vibroacoustic Loads from Clustered Rockets

The effect of stagger startup on the vibro-acoustic loads that form during the end- effects-regime of clustered rockets is studied using both full-scale (hot-gas) and laboratory scale (cold gas) data. Both configurations comprise three nozzles with thrust optimized parabolic contours that undergo free shock separated flow and restricted shock separated flow as well as an end-effects regime prior to flowing full. Acoustic pressure waveforms recorded at the base of the nozzle clusters are analyzed using various statistical metrics as well as time-frequency analysis. The findings reveal a significant reduction in end- effects-regime loads when engine ignition is staggered. However, regardless of stagger, both the skewness and kurtosis of the acoustic pressure time derivative elevate to the same levels during the end-effects-regime event thereby demonstrating the intermittence and impulsiveness of the acoustic waveforms that form during engine startup.

Rojo, Raymundo↗
NASA NTRS1991

Sufficient conditions for a local minimum of the Bolza problem with multiple terminal point constraints

Sufficient conditions for a weak relative minimum are derived for a form of the Bolza problem of variational calculus. The derivation ties together first-order and second-order conditions in a unified consistent manner, addresses controllability considerations in some detail, and clarifies some small inconsistencies in earlier work. The resulting second-order conditions of optimality involve the integration of fewer backward-sweep matrix elements than the standard conditions in the literature for problems with an unspecified final time. As a result, the backward-sweep matrices have the same general structure and dimensionality whether the final time is specified or free, with the terminal values of two of three sweep matrices being more complicated in the latter case.

Wood, Lincoln J.↗
NASA NTRS1984

Explicit guidance of drag modulated aeroassisted transfer between elliptical orbits

This paper presents the complete analysis of the problem of minimum-fuel aeroassisted transfer between coplanar elliptical orbits in the case where the orientation of the final orbit is free for selection in the optimization process. The comparison between the optimal pure propulsive transfer and the idealized aeroassisted transfer, by several passages through the atmosphere, is made. In the case where aeroassisted transfer provides fuel saving, a practical scheme for its realization by one passage is proposed. The maneuver consists of three phases: a deorbit phase for nonzero entry angle, followed by an atmospheric fly-through with variable drag control and completed by a postatmospheric phase. An explicit guidance formula for drag control is derived and it is shown that the required exit speed for ascent to the final orbit can be obtained with a very high degree of accuracy.

Vinh, N. X.↗
NASA NTRSDate not supplied

Considering Cell Proliferation to Optimize Detection of Radiation-induced 53BP1+ Foci in 15 Mouse Strains ex vivo

Due to high metabolic activity, proliferating cells continuously generate free radicals, which induce DNA double strand breaks (DSB). Fluorescently tagged nuclear foci of DNA repair protein 53 binding protein-1(53BP1) are used as a standard metric for measuring DSB formation at baseline and in response to environmental insults such as radiation. Here we demonstrate that the background level of spontaneous 53BP1+ foci formation can be modeled mathematically as a function of cell confluence, which is a metric of proliferation rate. This model was validated using spontaneous 53BP1+ foci data from 68 cultures of primary skin fibroblasts derived from 15 different strains of mice, showing a ~10 fold decrease from low to full confluence that is independent of mouse strain. We developed an online open access tool to correct for the impact of cell confluence on the detection of radiation-induced 53BP1+ foci (RIF). This tool provides guidelines for the number of cells required to reach statistical significance for the detection of excess foci induced by low doses of ionizing radiation as a function of confluence and time post-irradiation. We hope that this quantification tool will help the radiation biology community in the design of future experiments that utilize 53BP1+ foci-based quantification of radiation responses in vitro. Our “tool for enhanced results of RIF in cells” (terRIFic) can be found at: https://radbiolab.shinyapps.io/terrific/

Radiation, DNA damage, confluence↗
NASA NTRS2005

Vitamin E provides protection for bone in mature hindlimb unloaded male rats

The deleterious effects of skeletal unloading on bone mass and strength may, in part, result from increased production of oxygen-derived free radicals and proinflammatory cytokines. This study was designed to evaluate the ability of vitamin E (alpha-tocopherol), a free-radical scavenger with antiinflammatory properties, to protect against bone loss caused by skeletal unloading in mature male Sprague-Dawley rats. A 2 x 3 factorial design was used with either hindlimb unloading (HU) or normal loading (ambulatory; AMB), and low-dose (LD; 15 IU/kg diet), adequate-dose (AD; 75 IU/kg diet), or high-dose (HD; 500 IU/kg diet) vitamin E (DL-alpha-tocopherol acetate). To optimize the effects of vitamin E on bone, dietary treatments were initiated 9 weeks prior to unloading and continued during the 4-week unloading period, at which time animals were euthanized and blood and tissue samples were collected. Serum vitamin E was dose-dependently increased, confirming the vitamin E status of animals. The HD treatment improved oxidation parameters, as indicated by elevated serum ferric-reducing ability and a trend toward reducing tissue lipid peroxidation. Histomorphometric analysis of the distal femur revealed significant reductions in trabecular thickness (TbTh), double-labeled surface (dLS/BS), and rate of bone formation to bone volume (BFR/BV) due by HU. AMB animals on the HD diet and HU animals on the LD diet had reduced bone surface normalized to tissue volume (BS/TV) and trabecular number (TbN); however, the HD vitamin E protected against these changes in the HU animals. Our findings suggest that vitamin E supplementation provides modest bone protective effects during skeletal unloading.

NASA Discipline Musculoskeletal↗
NASA NTRS1999

Practical Aerodynamic Design Optimization Based on the Navier-Stokes Equations and a Discrete Adjoint Method

Compressible and incompressible versions of a three-dimensional unstructured mesh Reynolds-averaged Navier-Stokes flow solver have been differentiated and resulting derivatives have been verified by comparisons with finite differences and a complex-variable approach. In this implementation, the turbulence model is fully coupled with the flow equations in order to achieve this consistency. The accuracy demonstrated in the current work represents the first time that such an approach has been successfully implemented. The accuracy of a number of simplifying approximations to the linearizations of the residual have been examined. A first-order approximation to the dependent variables in both the adjoint and design equations has been investigated. The effects of a "frozen" eddy viscosity and the ramifications of neglecting some mesh sensitivity terms were also examined. It has been found that none of the approximations yielded derivatives of acceptable accuracy and were often of incorrect sign. However, numerical experiments indicate that an incomplete convergence of the adjoint system often yield sufficiently accurate derivatives, thereby significantly lowering the time required for computing sensitivity information. The convergence rate of the adjoint solver relative to the flow solver has been examined. Inviscid adjoint solutions typically require one to four times the cost of a flow solution, while for turbulent adjoint computations, this ratio can reach as high as eight to ten. Numerical experiments have shown that the adjoint solver can stall before converging the solution to machine accuracy, particularly for viscous cases. A possible remedy for this phenomenon would be to include the complete higher-order linearization in the preconditioning step, or to employ a simple form of mesh sequencing to obtain better approximations to the solution through the use of coarser meshes. An efficient surface parameterization based on a free-form deformation technique has been utilized and the resulting codes have been integrated with an optimization package. Lastly, sample optimizations have been shown for inviscid and turbulent flow over an ONERA M6 wing. Drag reductions have been demonstrated by reducing shock strengths across the span of the wing. In order for large scale optimization to become routine, the benefits of parallel architectures should be exploited. Although the flow solver has been parallelized using compiler directives. The parallel efficiency is under 50 percent. Clearly, parallel versions of the codes will have an immediate impact on the ability to design realistic configurations on fine meshes, and this effort is currently underway.

Grossman, Bernard↗
NASA NTRS1988

An improved photometric analysis of SX Aurigae

A new analysis of existing photometric data on SX Aurigae is presented which changes current conclusions in the literature concerning the star. It is concluded, based on fitting of the continuous spectrum to Kurucz model atmospheres and simplex optimization procedure utilized, that the system consists of two normal main-sequence stars of spectral types B2 and slightly earlier than B5. The components are very close to the inner Roche limiting surface. The secondary component is unevolved, and the primary component is extremely close to contact with its inner Roche lobe. The evidence indicates that mass flow from the primary to the secondary has not yet begun. Allowing the bolometric albedo of the secondary to be a free parameter produces a solution with an excellent fit to the B, V data. The derived value, A2 = 0.40, disagrees with the theoretical value for radiative atmospheres and agrees with the value for convective equilibrium.

Linnel, Albert P.↗
NASA NTRS1999

Practical Aerodynamic Design Optimization Based on the Navier-Stokes Equations and a Discrete Adjoint Method

The technical details are summarized below: Compressible and incompressible versions of a three-dimensional unstructured mesh Reynolds-averaged Navier-Stokes flow solver have been differentiated and resulting derivatives have been verified by comparisons with finite differences and a complex-variable approach. In this implementation, the turbulence model is fully coupled with the flow equations in order to achieve this consistency. The accuracy demonstrated in the current work represents the first time that such an approach has been successfully implemented. The accuracy of a number of simplifying approximations to the linearizations of the residual have been examined. A first-order approximation to the dependent variables in both the adjoint and design equations has been investigated. The effects of a "frozen" eddy viscosity and the ramifications of neglecting some mesh sensitivity terms were also examined. It has been found that none of the approximations yielded derivatives of acceptable accuracy and were often of incorrect sign. However, numerical experiments indicate that an incomplete convergence of the adjoint system often yield sufficiently accurate derivatives, thereby significantly lowering the time required for computing sensitivity information. The convergence rate of the adjoint solver relative to the flow solver has been examined. Inviscid adjoint solutions typically require one to four times the cost of a flow solution, while for turbulent adjoint computations, this ratio can reach as high as eight to ten. Numerical experiments have shown that the adjoint solver can stall before converging the solution to machine accuracy, particularly for viscous cases. A possible remedy for this phenomenon would be to include the complete higher-order linearization in the preconditioning step, or to employ a simple form of mesh sequencing to obtain better approximations to the solution through the use of coarser meshes. . An efficient surface parameterization based on a free-form deformation technique has been utilized and the resulting codes have been integrated with an optimization package. Lastly, sample optimizations have been shown for inviscid and turbulent flow over an ONERA M6 wing. Drag reductions have been demonstrated by reducing shock strengths across the span of the wing.

Grossman, Bernard↗
NASA NTRS2014

Use of MODIS Cloud Top Pressure to Improve Assimilation Yields of AIRS Radiances in GSI

Radiances from hyperspectral sounders such as the Atmospheric Infrared Sounder (AIRS) are routinely assimilated both globally and regionally in operational numerical weather prediction (NWP) systems using the Gridpoint Statistical Interpolation (GSI) data assimilation system. However, only thinned, cloud-free radiances from a 281-channel subset are used, so the overall percentage of these observations that are assimilated is somewhere on the order of 5%. Cloud checks are performed within GSI to determine which channels peak above cloud top; inaccuracies may lead to less assimilated radiances or introduction of biases from cloud-contaminated radiances.Relatively large footprint from AIRS may not optimally represent small-scale cloud features that might be better resolved by higher-resolution imagers like the Moderate Resolution Imaging Spectroradiometer (MODIS). Objective of this project is to "swap" the MODIS-derived cloud top pressure (CTP) for that designated by the AIRS-only quality control within GSI to test the hypothesis that better representation of cloud features will result in higher assimilated radiance yields and improved forecasts.

Zavodsky, Bradley↗
NASA NTRS2004

Aerodynamic Shape Optimization Based on Free-form Deformation

This paper presents a free-form deformation technique suitable for aerodynamic shape optimization. Because the proposed technique is independent of grid topology, we can treat structured and unstructured computational fluid dynamics grids in the same manner. The proposed technique is an alternative shape parameterization technique to a trivariate volume technique. It retains the flexibility and freedom of trivariate volumes for CFD shape optimization, but it uses a bivariate surface representation. This reduces the number of design variables by an order of magnitude, and it provides much better control for surface shape changes. The proposed technique is simple, compact, and efficient. The analytical sensitivity derivatives are independent of the design variables and are easily computed for use in a gradient-based optimization. The paper includes the complete formulation and aerodynamics shape optimization results.

Samareh, Jamshid A.↗
NASA NTRS1981

Optimum climb and descent trajectories for airline missions

The characteristics of optimum fixed-range trajectories whose structure is constrained to climb, steady cruise, and descent segments are derived by application of optimal control theory. The performance function consists of the sum of fuel and time costs, referred to as direct operating cost (DOC). The state variable is range to go and the independent variable is energy. In this formulation a cruise segment always occurs at the optimum cruise energy for sufficiently large range. At short ranges (400 n. mi. and less), a cruise segment may also occur below the optimum cruise energy. The existence of such a cruise segment depends primarily on the fuel flow vs thrust characteristics and on thrust constraints. If thrust is a free control variable along with airspeed, it is shown that such cruise segments will not generally occur. If thrust is constrained to some maximum value in climb and to some minimum in descent, such cruise segments generally will occur.

Erzberger, H.↗
NASA NTRS1980

Constrained optimum trajectories with specified range

The characteristics of optimum fixed-range trajectories whose structure is constrained to climb, steady cruise, and descent segments are derived by application of optimal control theory. The performance function consists of the sum of fuel and time costs, referred to as direct operating costs (DOC). The state variable is range-to-go and the independent variable is energy. In this formulation a cruise segment always occurs at the optimum cruise energy for sufficiently large range. At short ranges (500 n. mi. and less) a cruise segment may also occur below the optimum cruise energy. The existence of such a cruise segment depends primarily on the fuel flow vs thrust characteristics and on thrust constraints. If thrust is a free control variable along with airspeed, it is shown that such cruise segments will not generally occur. If thrust is constrained to some maximum value in climb and to some minimum in descent, such cruise segments generally will occur. The performance difference between free thrust and constrained thrust trajectories has been determined in computer calculations for an example transport aircraft.

Erzberger, H.↗
NASA NTRS1995

Damage-mitigating control of a reusable rocket engine for high performance and extended life

The goal of damage mitigating control in reusable rocket engines is to achieve high performance with increased durability of mechanical structures such that functional lives of the critical components are increased. The major benefit is an increase in structural durability with no significant loss of performance. This report investigates the feasibility of damage mitigating control of reusable rocket engines. Phenomenological models of creep and thermo-mechanical fatigue damage have been formulated in the state-variable setting such that these models can be combined with the plant model of a reusable rocket engine, such as the Space Shuttle Main Engine (SSME), for synthesizing an optimal control policy. Specifically, a creep damage model of the main thrust chamber wall is analytically derived based on the theories of sandwich beam and viscoplasticity. This model characterizes progressive bulging-out and incremental thinning of the coolant channel ligament leading to its eventual failure by tensile rupture. The objective is to generate a closed form solution of the wall thin-out phenomenon in real time where the ligament geometry is continuously updated to account for the resulting deformation. The results are in agreement with those obtained from the finite element analyses and experimental observation for both Oxygen Free High Conductivity (OFHC) copper and a copper-zerconium-silver alloy called NARloy-Z. Due to its computational efficiency, this damage model is suitable for on-line applications of life prediction and damage mitigating control, and also permits parametric studies for off-line synthesis of damage mitigating control systems. The results are presented to demonstrate the potential of life extension of reusable rocket engines via damage mitigating control. The control system has also been simulated on a testbed to observe how the damage at different critical points can be traded off without any significant loss of engine performance. The research work reported here is built upon concepts derived from the disciplines of Controls, Thermo-fluids, Structures, and Materials. The concept of damage mitigation, as presented in this report, is not restricted to control of rocket engines. It can be applied to any system where structural durability is an important issue.

Ray, Asok↗
NASA NTRS1998

Geometry Modeling and Grid Generation for Design and Optimization

Geometry modeling and grid generation (GMGG) have played and will continue to play an important role in computational aerosciences. During the past two decades, tremendous progress has occurred in GMGG; however, GMGG is still the biggest bottleneck to routine applications for complicated Computational Fluid Dynamics (CFD) and Computational Structures Mechanics (CSM) models for analysis, design, and optimization. We are still far from incorporating GMGG tools in a design and optimization environment for complicated configurations. It is still a challenging task to parameterize an existing model in today's Computer-Aided Design (CAD) systems, and the models created are not always good enough for automatic grid generation tools. Designers may believe their models are complete and accurate, but unseen imperfections (e.g., gaps, unwanted wiggles, free edges, slivers, and transition cracks) often cause problems in gridding for CSM and CFD. Despite many advances in grid generation, the process is still the most labor-intensive and time-consuming part of the computational aerosciences for analysis, design, and optimization. In an ideal design environment, a design engineer would use a parametric model to evaluate alternative designs effortlessly and optimize an existing design for a new set of design objectives and constraints. For this ideal environment to be realized, the GMGG tools must have the following characteristics: (1) be automated, (2) provide consistent geometry across all disciplines, (3) be parametric, and (4) provide sensitivity derivatives. This paper will review the status of GMGG for analysis, design, and optimization processes, and it will focus on some emerging ideas that will advance the GMGG toward the ideal design environment.

Samareh, Jamshid A.↗
NASA NTRS1978

Characteristics of constrained optimum trajectories with specified range

Necessary conditions of optimality are derived for trajectories whose structure is limited to climb, steady cruise, and descent segments. The performance function consists of the sum of fuel and time costs, referred to as direct operating cost (DOC). The state variable is range to go and the independent variable is energy. In this formulation a cruise segment always occurs at the optimum cruise energy for sufficiently large range. At short ranges (500 n. mi. and less) a cruise segment may also occur below the optimum cruise energy. The existence of such a cruise segment depends primarily on the fuel flow vs thrust characteristics and on thrust constraints. If thrust is a free control variable along with airspeed, it is shown that such cruise segments will not generally occur. If thrust is constrained to some maximum value in climb and to some minimum in descent, such cruise segments generally will occur. Computer calculations of typical short-range trajectories obtained about a 1% cost penalty for constraining the thrust.

Erzberger, H.↗
NASA NTRS1972

Optimal braking studies

To brake in minimum distance, the tire slip must be controlled to ride the peak of the mu-slip curve so that maximum ground force is developed between tire and pavement. The resulting control system differs from antiskid systems which react to impending wheel lockup. A simplified model is presented to permit development of a sound control strategy. Liapunov techniques are used to derive a peak riding adaptive controller applicable to each wheel of a breaking vehicle. The controller is applied to a more sophisticated model of a braking airplane with strut bending dynamics included. Simulation results verify the peak riding property of the controller and the rapid adaption of the controller to extreme runway conditions. The effect of actuator dynamics, perturbation frequency, type and location of sensors, absence of a free wheel, and a method in which the pilot's braking commands can be interfaced with the peak riding system are also considered.

Pazdera, J. S.↗