Search NASA⌕ Search

SEARCH · Search NASA

Results for “(non)linear 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 235 records · Page 13

Transient responses of phosphoric acid fuel cell power plant system

An analytical and computerized study of the steady state and transient response of a phosphoric acid fuel cell (PAFC) system was completed. Parametric studies and sensitivity analyses of the PAFC system's operation were accomplished. Four non-linear dynamic models of the fuel cell stack, reformer, shift converters, and heat exchangers were developed based on nonhomogeneous non-linear partial differential equations, which include the material, component, energy balance, and electrochemical kinetic features. Due to a lack of experimental data for the dynamic response of the components only the steady state results were compared with data from other sources, indicating reasonably good agreement. A steady state simulation of the entire system was developed using, nonlinear ordinary differential equations. The finite difference method and trial-and-error procedures were used to obtain a solution. Using the model, a PAFC system, that was developed under NASA Grant, NCC3-17, was improved through the optimization of the heat exchanger network. Three types of cooling configurations for cell plates were evaluated to obtain the best current density and temperature distributions. The steady state solutions were used as the initial conditions in the dynamic model. The transient response of a simplified PAFC system, which included all of the major components, subjected to a load change was obtained. Due to the length of the computation time for the transient response calculations, analysis on a real-time computer was not possible. A simulation of the real-time calculations was developed on a batch type computer. The transient response characteristics are needed for the optimization of the design and control of the whole PAFC system. All of the models, procedures and simulations were programmed in Fortran and run on IBM 370 computers at Cleveland State University and the NASA Lewis Research Center.

Lu, Cheng-Yi↗

Advances in Fluorescence Sensing Systems for the Remote Assessment of Nitrogen Supply in Field Corn

The studies described herein were conducted to better define changes in fluorescence properties of leaves from field grown corn (Zea mays L.) as they relate to varying levels of nitrogen (N) fertilization. This research was directed toward: 1) providing a remote non-destructive sensing technique to aid in the determination of optimal rates of N fertilization in corn crops and, 2) defining parameters for further development of fluorescence instrumentation to be operated remotely at field canopy levels. Fluorescence imaging bands centered in the blue (450 nm), green (525 nm), red (680 nm), and far-red (740 nm) and ratios of these bands were compared with the following plant parameters: rates of photosynthesis, N:C ratio, pigment concentrations, and grain yields. Both the fluorescence and physiological measures exhibited similar curvilinear responses to N fertilization level while significant linear correlations were obtained among fluorescence bands and band ratios to certain physiological measures of plant productivity. The red / blue, red / green, far-red / blue, far-red /green fluorescence ratios are well suited for remote observation and provided high correlations to grain yield, LAI, N:C, and chlorophyll contents. The results from this investigation indicate that fluorescence technology could aid in the determination of N fertilization requirements for corn. This discussion will also address design concepts and preliminary field trials of a mobile field-based Laser Induced Fluorescence Imaging System (LIFIS) capable of simultaneously acquiring images of four fluorescence emission bands from areas of plant canopies equaling 1 sq m and greater without interference of ambient solar radiation.

Corp, L. A.↗

Estimation of Sea Ice Thickness Distributions through the Combination of Snow Depth and Satellite Laser Altimetry Data

Combinations of sea ice freeboard and snow depth measurements from satellite data have the potential to provide a means to derive global sea ice thickness values. However, large differences in spatial coverage and resolution between the measurements lead to uncertainties when combining the data. High resolution airborne laser altimeter retrievals of snow-ice freeboard and passive microwave retrievals of snow depth taken in March 2006 provide insight into the spatial variability of these quantities as well as optimal methods for combining high resolution satellite altimeter measurements with low resolution snow depth data. The aircraft measurements show a relationship between freeboard and snow depth for thin ice allowing the development of a method for estimating sea ice thickness from satellite laser altimetry data at their full spatial resolution. This method is used to estimate snow and ice thicknesses for the Arctic basin through the combination of freeboard data from ICESat, snow depth data over first-year ice from AMSR-E, and snow depth over multiyear ice from climatological data. Due to the non-linear dependence of heat flux on ice thickness, the impact on heat flux calculations when maintaining the full resolution of the ICESat data for ice thickness estimates is explored for typical winter conditions. Calculations of the basin-wide mean heat flux and ice growth rate using snow and ice thickness values at the 70 m spatial resolution of ICESat are found to be approximately one-third higher than those calculated from 25 km mean ice thickness values.

Kurtz, Nathan T.↗

Issues in deep space radiation protection

The exposures in deep space are largely from the Galactic Cosmic Rays (GCR) for which there is as yet little biological experience. Mounting evidence indicates that conventional linear energy transfer (LET) defined protection quantities (quality factors) may not be appropriate for GCR ions. The available biological data indicates that aluminum alloy structures may generate inherently unhealthy internal spacecraft environments in the thickness range for space applications. Methods for optimization of spacecraft shielding and the associated role of materials selection are discussed. One material which may prove to be an important radiation protection material is hydrogenated carbon nanofibers. c 2001. Elsevier Science Ltd. All rights reserved.

short duration↗

Optimization of Orifice Geometry for Cross-Flow Mixing in a Cylindrical Duct

Mixing of gaseous jets in a cross-flow has significant applications in engineering, one example of which is the dilution zone of a gas turbine combustor. Despite years of study, the design of jet injection in combustors is largely based on practical experience. A series of experiments was undertaken to delineate the optimal mixer orifice geometry. A cross-flow to core-flow momentum-flux ratio of 40 and a mass flow ratio of 2.5 were selected as representative of an advanced design. An experimental test matrix was designed around three variables: the number of orifices, the orifice aspect ratio (long-to-short dimension), and the orifice angle. A regression analysis was performed on the data to arrive at an interpolating equation that predicted the mixing performance of orifice geometry combinations within the range of the test matrix parameters. Results indicate that mixture uniformity is a non-linear function of the number of orifices, the orifice aspect ratio, and the orifice angle. Optimum mixing occurs when the asymptotic mean jet trajectories are in the range of 0.35 less than r/R less than 0.5 (where r = 0 is at the mixer wall) at z/R = 1.0. At the optimum number of orifices, the difference between shallow-angled slots with large aspect ratios and round holes is minimal and either approach will lead to good mixing performance. At the optimum number of orifices, it appears possible to have two local optimums where one corresponds to an aspect ratio of 1.0 and the other to a high aspect ratio.

Sowa, W. A.↗

Multiple Mode Actuation of a Turbulent Jet

The effects of multiple mode periodic excitation on the evolution of a circular turbulent jet were studied experimentally. A short, wide-angle diffuser was attached to the jet exit. Streamwise and cross-stream excitations were introduced at the junction between the jet exit and the diffuser inlet on opposing sides of the jet. The introduction of high amplitude, periodic excitation in the streamwise direction enhances the mixing and promotes attachment of the jet shear-layer to the diffuser wall. Cross-stream excitation applied over a fraction of the jet circumference can deflect the jet away from the excitation slot. The two modes of excitation were combined using identical frequencies and varying the relative phase between the two actuators in search of an optimal response. It is shown that, for low and moderate periodic momentum input levels, the jet deflection angles depend strongly on the relative phase between the two actuators. Optimum performance is achieved when the phase difference is pi +/- pi/6. The lower effectiveness of the equal phase excitation is attributed to the generation of an azimuthally symmetric mode that does not produce the required non-axisymmetric vectoring. For high excitation levels, identical phase becomes more effective, while phase sensitivity decreases. An important finding was that with proper phase tuning, two unsteady actuators can be combined to obtain a non-linear response greater than the superposition of the individual effects.

Pack, LaTunia G.↗

Multiple Mode Actuation of a Turbulent Jet

The effects of multiple mode periodic excitation on the evolution of a circular turbulent jet were studied experimentally. A short, wide-angle diffuser was attached to the jet exit. Streamwise and cross-stream excitations were introduced at the junction between the jet exit and the diffuser inlet on opposing sides of the jet. The introduction of high amplitude, periodic excitation in the streamwise direction enhances the mixing and promotes attachment of the jet shear-layer to the diffuser wall. Cross-stream excitation applied over a fraction of the jet circumference can deflect the jet away from the excitation slot. The two modes of excitation were combined using identical frequencies and varying the relative phase between the two actuators in search of an optimal response. It is shown that, for low and moderate periodic momentum input levels, the jet deflection angles depend strongly on the relative phase between the two actuators. Optimum performance is achieved when the phase difference is pi +/- pi/6. The lower effectiveness of the equal phase excitation is attributed to the generation of an azimuthally symmetric mode that does not produce the required non-axisymmetric vectoring. For high excitation levels, identical phase becomes more effective, while phase sensitivity decreases. An important finding was that with proper phase tuning, two unsteady actuators can be combined to obtain a non-linear response greater than the superposition of the individual effects.

Pack, LaTunia G.↗

Massively Parallel Dantzig-Wolfe Decomposition Applied to Traffic Flow Scheduling

Optimal scheduling of air traffic over the entire National Airspace System is a computationally difficult task. To speed computation, Dantzig-Wolfe decomposition is applied to a known linear integer programming approach for assigning delays to flights. The optimization model is proven to have the block-angular structure necessary for Dantzig-Wolfe decomposition. The subproblems for this decomposition are solved in parallel via independent computation threads. Experimental evidence suggests that as the number of subproblems/threads increases (and their respective sizes decrease), the solution quality, convergence, and runtime improve. A demonstration of this is provided by using one flight per subproblem, which is the finest possible decomposition. This results in thousands of subproblems and associated computation threads. This massively parallel approach is compared to one with few threads and to standard (non-decomposed) approaches in terms of solution quality and runtime. Since this method generally provides a non-integral (relaxed) solution to the original optimization problem, two heuristics are developed to generate an integral solution. Dantzig-Wolfe followed by these heuristics can provide a near-optimal (sometimes optimal) solution to the original problem hundreds of times faster than standard (non-decomposed) approaches. In addition, when massive decomposition is employed, the solution is shown to be more likely integral, which obviates the need for an integerization step. These results indicate that nationwide, real-time, high fidelity, optimal traffic flow scheduling is achievable for (at least) 3 hour planning horizons.

Rios, Joseph Lucio↗

Thermal traction contact performance evaluation under fully flooded and starved conditions

Ultra high speed traction tests were performed on two traction fluids commonly employed. Traction data on these fluids is required for purposes of traction drive design optimization techniques. To obtain the traction data, an existing twin disc traction test machine was employed. This machine was modified to accommodate the range of test variables. All the data reported was obtained under conditions of side slip, a technique whereby only low power levels are required to simulate real traction drive contacts. Theoretical traction predictions were performed for a representative number of curves that showed the influence of rolling velocity, of contact pressure and of aspect ratio. To establish the accuracy of the thermal model the predictions were performed ith increasing levels of independence of experimentally determined parameters. In the final resulting prediction only two non linear thermal parameters were used for the prediction of 15 different traction curves covering the entire range of variables as used in the investigation, with the exception of the influence of asperity traction. Comparison of these theoretical curves and corresponding experimental traces show very good agreement.

Tevaarwerk, J. L.↗

Effects of walking velocity on vertical head and body movements during locomotion

Trunk and head movements were characterized over a wide range of walking speeds to determine the relationship between stride length, stepping frequency, vertical head translation, pitch rotation of the head, and pitch trunk rotation as a function of gait velocity. Subjects (26-44 years old) walked on a linear treadmill at velocities of 0.6-2.2 m/s. The head and trunk were modeled as rigid bodies, and rotation and translation were determined using a video-based motion analysis system. At walking speeds up to 1.2 m/s there was little head pitch movement in space, and the head pitch relative to the trunk was compensatory for trunk pitch. As walking velocity increased, trunk pitch remained approximately invariant, but a significant head translation developed. This head translation induced compensatory head pitch in space, which tended to point the head at a fixed point in front of the subject that remained approximately invariant with regard to walking speed. The predominant frequency of head translation and rotation was restricted to a narrow range from 1.4 Hz at 0.6 m/s to 2.5 Hz at 2.2 m/s. Within the range of 0.8-1.8 m/s, subjects tended to increase their stride length rather than step frequency to walk faster, maintaining the predominant frequency of head movement at close to 2.0 Hz. At walking speeds above 1.2 m/s, head pitch in space was highly coherent with, and compensatory for, vertical head translation. In the range 1.2-1.8 m/s, the power spectrum of vertical head translation was the most highly tuned, and the relationship between walking speed and head and trunk movements was the most linear. We define this as an optimal range of walking velocity with regard to head-trunk coordination. The coordination of head and trunk movement was less coherent at walking velocities below 1.2 m/s and above 1.8 m/s. These results suggest that two mechanisms are utilized to maintain a stable head fixation distance over the optimal range of walking velocities. The relative contribution of each mechanism to head orientation depends on the frequency of head movement and consequently on walking velocity. From consideration of the frequency characteristics of the compensatory head pitch, we infer that compensatory head pitch movements may be produced predominantly by the angular vestibulocollic reflex (aVCR) at low walking speeds and by the linear vestibulocollic reflex (1VCR) at the higher speeds.

NASA Discipline Neuroscience↗

An optimized index of human cardiovascular adaptation to simulated weightlessness

Prolonged exposure to weightlessness is known to produce a variety of cardiovascular changes, some of which may influence the astronaut's performance during a mission. In order to find a reliable indicator of cardiovascular adaptation to weightlessness, we analyzed data from nine male subjects after a 24-hour period of normal activity and after a period of simulated weightlessness produced by two hours in a launch position followed by 20 hours of 6 degrees head-down tilt plus pharmacologically induced diuresis (furosemide). Heart rate, arterial pressure, thoracic fluid index, and radial flow were analyzed. Autoregressive spectral estimation and decomposition were used to obtain the spectral components of each variable from the subjects in the supine position during pre- and post-simulated weightlessness. We found a significant decrease in heart rate power and an increase in thoracic fluid index power in the high frequency region (0.2-0.45 Hz) and significant increases in radial flow and arterial pressure powers in the low frequency region (<0.2 Hz) in response to simulated weightlessness. However, due to the variability among subjects, any single variable appeared limited as a dependable index of cardiovascular adaptation to weightlessness. The backward elimination algorithm was then used to select the best discriminatory features from these spectral components. Fisher's linear discriminant and Bayes' quadratic discriminant were used to combine the selected features to obtain an optimal index of adaptation to simulated weightlessness. Results showed that both techniques provided improved discriminant performance over any single variable and thus have the potential for use as an index to track adaptation and prescribe countermeasures to the effects of weightlessness.

NASA Discipline Number 14-10↗

Optimizations on Designing High-Resolution Finite-Difference Schemes

We describe a general optimization procedure for both maximizing the resolution characteristics of existing finite differencing schemes as well as designing finite difference schemes that will meet the error tolerance requirements of numerical solutions. The procedure is based on an optimization process. This is a generalization of the compact scheme introduced by Lele in which the resolution is improved for single, one-dimensional spatial derivative, whereas in the present approach the complete scheme, after spatial and temporal discretizations, is optimized on a range of parameters of the scheme and the governing equations. The approach is to linearize and Fourier analyze the discretized equations to check the resolving power of the scheme for various wave number ranges in the solution and optimize the resolution to satisfy the requirements of the problem. This represents a constrained nonlinear optimization problem which can be solved to obtain the nodal weights of discretization. An objective function is defined in the parametric space of wave numbers, Courant number, Mach number and other quantities of interest. Typical criterion for defining the objective function include the maximization of the resolution of high wave numbers for acoustic and electromagnetic wave propagations and turbulence calculations. The procedure is being tested on off-design conditions of non-uniform mesh, non-periodic boundary conditions, and non-constant wave speeds for scalar and system of equations. This includes the solution of wave equations and Euler equations using a conventional scheme with and without optimization and the design of an optimum scheme for the specified error tolerance.

Liu, Yen↗

Detection and Tracking of Aircraft in the Far-Field from Small Unmanned Aerial Systems

Onboard far-field aircraft detection is needed for safe non-cooperative traffic mitigation in autonomous small Unmanned Aerial System (sUAS) operations. Machine vision systems, based on standard optics and visible light detectors, possess the ideal size, weight, and power (SWaP) requirements for sUAS. This work presents the design and analysis of a novel aircraft detection and tracking pipeline based on optical sensing alone. Key contributions of the work include a refined range inequality model based on sensing and detection with FAA well-clear separation assurance distances between aircraft in mind, a detector fusion method to maximize the benefit of two image detectors, and a comparative analysis of Linear Kalman-filtering and Extended Kalman-filtering to seek optimal tracking performance. The pipeline is evaluated offline against multiple intruder platforms, using two types of flight encounters: multirotor sUAS vs. fixed-wing sUAS and multirotor sUAS vs. general aviation(GA)plane. Analysis is restricted to the rate-limiting head-on and departing collision volume cases vertically separated for safety. Results indicate that it is feasible to use the proposed optical spatial-temporal tracking algorithm to provide adequate alerting time to prevent penetration of well-clear separation volumes for both sUAS and GA aircraft.

Unmanned Aerial System↗

DSP Implementation of the Retinex Image Enhancement Algorithm

The Retinex is a general-purpose image enhancement algorithm that is used to produce good visual representations of scenes. It performs a non-linear spatial/spectral transform that synthesizes strong local contrast enhancement and color constancy. A real-time, video frame rate implementation of the Retinex is required to meet the needs of various potential users. Retinex processing contains a relatively large number of complex computations, thus to achieve real-time performance using current technologies requires specialized hardware and software. In this paper we discuss the design and development of a digital signal processor (DSP) implementation of the Retinex. The target processor is a Texas Instruments TMS320C6711 floating point DSP. NTSC video is captured using a dedicated frame-grabber card, Retinex processed, and displayed on a standard monitor. We discuss the optimizations used to achieve real-time performance of the Retinex and also describe our future plans on using alternative architectures.

Hines, Glenn↗

DSP Implementation of the Multiscale Retinex Image Enhancement Algorithm

The Retinex is a general-purpose image enhancement algorithm that is used to produce good visual representations of scenes. It performs a non-linear spatial/ spectral transform that synthesizes strong local contrast enhancement and color constancy. A real-time, video frame rate implementation of the Retinex is required to meet the needs of various potential users. Retinex processing contains a relatively large number of complex computations, thus to achieve real-time performance using current technologies requires specialized hardware and software. In this paper we discuss the design and development of a digital signal processor (DSP) implementation of the Retinex. The target processor is a Texas Instruments TMS320C6711 floating point DSP. NTSC video is captured using a dedicated frame grabber card, Retinex processed, and displayed on a standard monitor. We discuss the optimizations used to achieve real-time performance of the Retinex and also describe our future plans on using alternative architectures.

Hines, Glenn D.↗

Attitude Reconstruction of Free-Flight CFD Generated Trajectories Using Non-Linear Pitch Damping Coefficient Curves

Attitude history reconstruction of Free-flight CFD generated trajectories with non-linear pitch damping coefficient curves is investigated. Free-flight CFD simulations of the capsule shape used for the Genesis sample return mission and the upcoming Dragonfly mission to Titan are conducted for 1-, 2-, and 3-degree-of-freedom cases. Two different data reduction methodologies are employed to derive a pitch damping curve as a function of instantaneous angle of attack. These curves are then used to reconstruct the attitude history of the body which is compared to the raw simulation results. While both data reduction methods produce pitch damping curves that can generally reconstruct the motion seen in the Free-flight simulations, it is found that optimization of the pitch damping curve using an inverse estimation process yields superior and more generalizable results. Further refinement of this technique could allow pitch damping curves derived using CFD to serve as a capability complementary to existing techniques for dynamic stability characterization.

entry↗

Attitude Reconstruction of Free-Flight CFD Generated Trajectories Using Non-Linear Pitch Damping Coefficient Curves

Attitude history reconstruction of Free-flight CFD generated trajectories with non-linear pitch damping coefficient curves is investigated. Free-flight CFD simulations of the capsule shape used for the Genesis sample return mission and the upcoming Dragonfly mission to Titan are conducted for 1-, 2-, and 3-degree-of-freedom cases. Two different data reduction methodologies are employed to derive a pitch damping curve as a function of instantaneous angle of attack. These curves are then used to reconstruct the attitude history of the body which is compared to the raw simulation results. While both data reduction methods produce pitch damping curves that can generally reconstruct the motion seen in the Free-flight simulations, it is found that optimization of the pitch damping curve using an inverse estimation process yields superior and more generalizable results. Further refinement of this technique could allow pitch damping curves derived using CFD to serve as a capability complementary to existing techniques for dynamic stability characterization.

entry↗

Large scale nonlinear numerical optimal control for finite element models of flexible structures

This paper discusses the development of large scale numerical optimal control algorithms for nonlinear systems and their application to finite element models of structures. This work is based on our expansion of the optimal control algorithm (DDP) in the following steps: improvement of convergence for initial policies in non-convex regions, development of a numerically accurate penalty function method approach for constrained DDP problems, and parallel processing on supercomputers. The expanded constrained DDP algorithm was applied to the control of a four-bay, two dimensional truss with 12 soft members, which generates geometric nonlinearities. Using an explicit finite element model to describe the structural system requires 32 state variables and 10,000 time steps. Our numerical results indicate that for constrained or unconstrained structural problems with nonlinear dynamics, the results obtained by our expanded constrained DDP are significantly better than those obtained using linear-quadratic feedback control.

Shoemaker, Christine A.↗