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

A Variational Formulation for the Finite Element Analysis of Sound Wave Propagation in a Spherical Shell

Development of design tools to furnish optimal acoustic environments for lightweight aircraft demands the ability to simulate the acoustic system on a workstation. In order to form an effective mathematical model of the phenomena at hand, we have begun by studying the propagation of acoustic waves inside closed spherical shells. Using a fully-coupled fluid-structure interaction model based upon variational principles, we have written a finite element analysis program and are in the process of examining several test cases. Future investigations are planned to increase model accuracy by incorporating non-linear and viscous effects.

Lebiedzik, Catherine↗

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↗

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.↗

Space Shuttle propulsion parameter estimation using optimal estimation techniques, volume 1

The mathematical developments and their computer program implementation for the Space Shuttle propulsion parameter estimation project are summarized. The estimation approach chosen is the extended Kalman filtering with a modified Bryson-Frazier smoother. Its use here is motivated by the objective of obtaining better estimates than those available from filtering and to eliminate the lag associated with filtering. The estimation technique uses as the dynamical process the six degree equations-of-motion resulting in twelve state vector elements. In addition to these are mass and solid propellant burn depth as the ""system'' state elements. The ""parameter'' state elements can include aerodynamic coefficient, inertia, center-of-gravity, atmospheric wind, etc. deviations from referenced values. Propulsion parameter state elements have been included not as options just discussed but as the main parameter states to be estimated. The mathematical developments were completed for all these parameters. Since the systems dynamics and measurement processes are non-linear functions of the states, the mathematical developments are taken up almost entirely by the linearization of these equations as required by the estimation algorithms.

Source record↗

Modeling Orifice Pulse Tube Coolers

We have developed a calculational model that treats all the components of an orifice pulse tube cooler. We base our analysis on 1-dimensional thermodynamic equations for the regenerator and we assume that all mass flows, pressure oscillations and temperature oscillations are small and sinusoidal. Non-linear pressure drop effects are included in the regenerator to account for finite pressure amplitude effects. The resulting mass flows and pressures are matched at the boundaries with the other components of the cooler: compressor, aftercooler, cold heat exchanger, pulse tube, hot heat exchanger, orifice and reservoir. The results of the calculation are oscillating pressures, mass flows and enthalpy flows in the main components of the cooler. By comparing with the calculations of other available models, we show that our model is very similar to REGEN 3 from NIST and DeltaE from Los Alamos National Lab. Our model is much easier to use than other available models because of its simple graphical interface and the fact that no guesses are required for the operating pressures or mass flows. In addition, the model only requires a few minutes of running time allowing many parameters to be optimized in a reasonable time. A version of the model is available for use over the World Wide Web at http://irtek.arc.nasa.gov. Future enhancements include adding a bypass orifice and including second order terms in steady mass streaming and steady heat transfer. A two-dimensional anelastic approximation of the fluid equations will be used as the basis for the latter analysis. Preliminary results are given in dimensionless numbers appropriate for oscillating compressible flows. The model shows how transverse heat transfer reduces enthalpy flow, particularly for small pulse tubes. The model also clearly shows mass recirculation in the open tube on the order of the tube length. They result from the higher order Reynolds stresses. An interesting result of the linearized approach is that the steady mass streaming does not affect the enthalpy flow at second order. The major effect of recirculating mass streaming is to increase transverse temperature gradients, which leads to higher entropy production and reduced efficiency.

Roach, Kittel P.↗

Efficient Low-Speed Flight in a Wind Field

A new software tool was needed for flight planning of a high altitude, low speed unmanned aerial vehicle which would be flying in winds close to the actual airspeed of the vehicle. An energy modeled NLP (non-linear programming) formulation was used to obtain results for a variety of missions and wind profiles. The energy constraint derived included terms due to the wind field and the performance index was a weighted combination of the amount of fuel used and the final time. With no emphasis on time and with no winds the vehicle was found to fly at maximum lift to drag velocity, V(sub md). When flying in tail winds the velocity was less than V(sub md), while flying in head winds the velocity was higher than V(sub md). A family of solutions was found with varying times of flight and varying fuel amounts consumed which will aid the operator in choosing a flight plan depending on a desired landing time. At certain parts of the flight, the turning terms in the energy constraint equation were found to be significant. An analysis of a simpler vertical plane cruise optimal control problem was used to explain some of the characteristics of the vertical plane NLP results.

Feldman, Michael A.↗

Dynamics and control of coherent structures in the turbulent wall layer: An overview

We expand the velocity field in the vicinity of the wall in empirical eigenfunctions obtained from experiment. Truncating our system and using Galerkin projection, we obtain a closed set of non-linear ordinary differential equations with ten degrees of freedom. We find a rich dynamical behavior, including in particular a heteroclinic attracting orbit giving rise to intermittency. The intermittent jump from one attracting point to the other resembles in many respects the bursts observed in experiments. Specifically, the time between jumps and the duration of the jumps, is approximately that observed in a burst; the jump begins with the formation of a narrowed and intensified updraft, like the ejection phase of a burst, and is followed by a gentle, diffuse downdraft like the sweep phase of a burst. The magnitude of the Reynolds stress spike produced during a burst is limited by our truncation. The behavior is quite robust, much of it being due to the symmetries present (Aubry's group has examined dimensions up to 128 with persistence of the global behavior). We have examined eigenvalues and coefficients obtained from experiment, and from exact simulation, which differ in magnitude. Similar behavior is obtained in both cases; in the latter case, the heteroclinic orbits connect limit cycles instead of fixed points, corresponding to cross-stream waving of the streamwise rolls. The bifurcation diagram remains structurally similar, but somewhat distorted. The role of the pressure term is made clear - it triggers the intermittent jumps, which otherwise would occur at longer and longer intervals, as the system trajectory is attracted closer and closer to the heteroclinic cycle. The pressure term results in the jumps occurring at essentially random times, and the magnitude of the signal determines the average timing. Stretching of the wall region shows that the model is consistent with observations of polymer drag reduction. Change of the third order coefficients, corresponding to acceleration or deceleration of the mean flow, changes the heteroclinic cycles from attracting to repelling, increasing or decreasing the stability, in agreement with observations. The existence of fixed points is an artifact introduced by the projection; however, a decoupled model still displays the rich dynamics. Numerous assumptions made in Aubry et al. (1988) can now be proved exactly. Feeding back eigenfuncitons with the proper phase can delay the bursting, (the heteroclinic jump to the other fixed point), decreasing the drag. It is also possible to speed up the bursting, increasing mixing to control separation. Our approach is optimal for short time tracking in control.

Berkooz, Gal↗

Cloud Ozone Dust Imager (CODI): Investigation and Technical Plan - Volume 1

The Cloud Ozone Dust Imager (CODI) is proposed to investigate the current climatic balance of the Mars atmosphere, with particular emphasis on the important but poorly understood roles which dust and water ice aerosols play in this balance. The large atmospheric heating (20-50 K) resulting from global dust storms around Mars perihelion is well recognized. However, groundbased observations of Mars atmospheric temperatures, water vapor, and clouds since the Viking missions have identified a much colder, cloudier atmosphere around Mars aphelion that may prove as important as global dust storms in determining the interannual and long-term behavior of the Mars climate. The key climate issues CODI is designed to investigate are: 1) the degree to which non-linear interactions between atmospheric dust heating, water vapor saturation, and cloud nucleation influence the seasonal and interannual variability of the Mars atmosphere, and 2) whether the strong orbital forcing of atmospheric dust loading, temperatures and water vapor saturation determines the long-term balance of Mars water, as reflected in the north-south hemispheric asymmetries of atmospheric water vapor and polar water ice abundances. The CODI experiment will measure the daily, seasonal and (potentially) interannual variability of atmospheric dust and cloud opacities, and the key physical properties of these aerosols which determine their role in the climate cycles of Mars. CODI is a small (1.2 kg), fixed pointing camera, in which four wide-angle (+/- 70 deg) lenses illuminate fixed filters and CCD arrays. Simultaneous sky/surface imaging of Mars is obtained at an angular resolution of 0.28 deg/pixel for wavelengths of 255, 336, 502, and 673 nm (similar to Hubble Space Telescope filters). These wavelengths serve to measure atmospheric ozone (255 and 336 nm), discriminate ice and dust aerosols (336 and 673 nm), and construct color images (336, 502, and 673 nm). The CODI images are detected on four 512 x 512 pixel arrays, as partitioned on two 1024 x 1024 CCD's operated in frame transfer mode. The center of the CODI field-of-view is canted 40 deg from the zenith direction to obtain sky brightness measurements and a 20 deg surface field-of-view. Daily image observations will be conducted when the Sun is greater than or equal to 5 deg outside the edge of the CODI field-of-view, and twilight and nighttime imaging will obtained on a weekly basis. The 673 nm channel includes a polarizer wheel to obtain sky/surface polarimetry. A dust cover protects the entire lens assemblies of all four CODI channels. This opaque dust cover, which is normally opened for CODI imaging, includes a small fixed mirror and transparent window positioned above the 673 nm lens, to redirect the 673 nm field-of-view to the surface for descent imaging. Fixed pointing, internal data buffering, low operating power (2-4 W for less than or equal to 30 seconds), selective data transmission, and simple operational characteristics of the CODI experiment place minimum resource and operational demands on the Mars Surveyor 1998 lander. The CODI science goals are optimized for, but not restricted to, a low-latitude landing site (20 deg S-30 deg N). The primary CODI measurement objectives are the opacities, wave forms, particle properties (size, shape, and alignment), and heights of clouds; the opacities, particle properties, and vertical distribution of dust; and the opacity and vertical distribution of ozone. The variability of cloud, ozone, and dust opacities will be determined on diurnal, daily, and seasonal timescales. Wind velocities will be determined from cloud motions and wave characteristics; and the temporal variability of atmospheric water vapor, with limited altitude information, will be inferred from the CODI ozone observations. Secondary measurement objectives include limited descent imaging capability, surface uv-visible photometry and polarimetry, photochemistry, and meteorite infall rates.

Clancy, R. Todd↗

Synthesis of a non-linear feedback system with significant plant-ignorance for prescribed time-domain tolerances

The problem considered is the design of a feedback system containing a linear, time invariant, minimum phase plant, whose parameters are known only within given bounds, such that the time response of the system remains within specified limits. A quasi-optimal design, for given design constraints, is one which minimizes the effect of white sensor noise on the input to the plant. An investigation was conducted on the use of the non linear device known as the Clegg integrator in the design of such a system. The describing function of the Clegg integrator has the same magnitude characteristic, apart from a scale factor, as the linear integrator, but has 52 deg less phase-lag, at all frequencies, than the linear integrator; thus, when used in a feedback system, it provides a larger stability margin than the linear integrator. This property allows the nonlinear feedback system to be designed so that the sensor noise is attenuated more than in the linear design.

Krishnan, K. R.↗

Dual characterizations of optimal systems.

The complementary variational principle developed in a Hilbert space setting provides a duality principle in the calculus of variations with dynamic constraints. This concept is adopted in this paper to investigate dual characterizations of optimal control systems. Systems under consideration include those with dynamics governed by linear ordinary differential equations, linear partial differential equations and non-linear ordinary differential equations.

Chan, W. L.↗

Mission analysis flow sequencing optimization

This investigation is an extension of a project dealing with the problem of optimal use of ground resources for future space missions. This problem was formulated as a linear programming problem using an indirect approach. Instead of minimizing the inventory level of needed ground resources, the overlapping periods during which the same types of resources are used by various flights are minimized. The model was built upon the assumption that during the time interval under consideration, the costs of various needed resources remain constant. Under other assumptions concerning costs of resources, the objective function, in general, assumes a non-linear form. In this study, one case where the form of the objective function turns out to be quadratic is considered. Also, disadvantages and limitations of the approach used are briefly discussed.

Scott, M.↗

Position Control for Non-linear, Multiple-Link Robots

Several approaches to complicated control problem examined. Report surveys methods for controlling motion of robot manipulators. Applies to coupled, highly nonlinear multiple-link robots. Presents adaptive control technique based on discrete linear model of robot arm, obtained by linearizing nonlinear dynamical equations about nominal trajectory. Parameters of model calculated with adaptive parameter-identification algorithm. Based on system model, optimal control input computed using position and velocity feedback.

Seraji, H.↗

A study of the use of linear programming techniques to improve the performance in design optimization problems

This project has two objectives. The first is to determine whether linear programming techniques can improve performance when handling design optimization problems with a large number of design variables and constraints relative to the feasible directions algorithm. The second purpose is to determine whether using the Kreisselmeier-Steinhauser (KS) function to replace the constraints with one constraint will reduce the cost of total optimization. Comparisons are made using solutions obtained with linear and non-linear methods. The results indicate that there is no cost saving using the linear method or in using the KS function to replace constraints.

Young, Katherine C.↗

A study of equation solvers for linear and non-linear finite element analysis on parallel processing computers

Concurrent computing environments provide the means to achieve very high performance for finite element analysis of systems, provided the algorithms take advantage of multiple processors. The authors have examined several algorithms for both linear and nonlinear finite element analysis. The performance of these algorithms on an Alliant FX/80 parallel supercomputer has been studied. For single load case linear analysis, the optimal solution algorithm is strongly problem dependent. For multiple load cases or nonlinear analysis through a modified Newton-Raphson method, decomposition algorithms are shown to have a decided advantage over element-by-element preconditioned conjugate gradient algorithms.

Watson, Brian C.↗

Mechanically fastened composite laminates subjected to combined bearing-bypass and shear loading

Bolts and rivets provide a means of load transfer in the construction of aircraft. However, they give rise to stress concentrations and are often the source and location of static and fatigue failures. Furthermore, fastener holes are prone to cracks during take-off and landing. These cracks present the most common origin of structural failures in aircraft. Therefore, accurate determination of the contact stresses associated with such loaded holes in mechanically fastened joints is essential to reliable strength evaluation and failure prediction. As the laminate is subjected to loading, the contact region, whose extent is not known, develops between the fastener and the hole boundary through this contact region, which consists of slip and no-slip zones due to friction. The presence of the unknown contact stress distribution over the contact region between the pin and the composite laminate, material anisotropy, friction between the pin and the laminate, pin-hole clearance, combined bearing-bypass and shear loading, and finite geometry of the laminate result in a complex non-linear problem. In the case of bearing-bypass loading in compression, this non-linear problem is further complicated by the presence of dual contact regions. Previous research concerning the analysis of mechanical joints subjected to combined bearing-bypass and shear loading is non-existent. In the case of bearing-bypass loading only, except for the study conducted by Naik and Crews (1991), others employed the concept of superposition which is not valid for this non-linear problem. Naik and Crews applied a linear finite element analysis with conditions along the pin-hole contact region specified as displacement constraint equations. The major shortcoming of this method is that the variation of the contract region as a function of the applied load should be known a priori. Also, their analysis is limited to symmetric geometry and material systems, and frictionless boundary conditions. Since the contact stress distribution and the contact region are not known a priori, they did not directly impose the boundary conditions appropriate for modelling the contact and on-contact regions between the fastener and the hole. Furthermore, finite element analysis is not suitable for iterative design calculations for optimizing laminate construction in the presence of fasteners under complex loading conditions. In this study, the solution method developed by Madenci and Ileri (1992a,b) has been extended to determine the contact stresses in mechanical joints under combined bearing-bypass and shear loading, and bearing-bypass loading in compression resulting in dual contact regions.

Madenci, Erdogan↗

Fluid flow volume measurements using a capacitance/conductance probe system

A probe system has been developed to measure the flow volume of a single fluid passing through an orifice or flow line. The system employs both capacitance and a conductance probe at the orifice, together with phase detection and data acquisition circuitry to measure flow volume and salinity under low or zero gravity conditions. A wide variety of frequencies can be used for the radio frequency (RF) signal source which is chosen primarily by the capacitance of the orifice probe and the fluid passing through the orifice. Rapid measurements are made using the reflected signal from the orifice probe to determine the 'instantaneous' permittivity of the fluid/gas mixture passing through. The 'instantaneous' measurements are integrated over time to determine flow volume. Analysis reveals that a narrow orifice helps to reduce non-linearities caused by differing flow rates. The geometry of 'deflectors' and 'directors' for the flowing fluid are important in obtaining linearity. Measured data shows that a volume measurement accuracy of approximately four percent can be consistently achieved. The prototype hardware system and associated software have been optimized and are available for further applications. The system has immediate application in low or zero gravity environments for measurements of urine or other liquid volumes.

Nguyen, T. X.↗

Development of Curved-Plate Elements for the Exact Buckling Analysis of Composite Plate Assemblies Including Transverse-Shear Effects

The analytical formulation of curved-plate non-linear equilibrium equations including transverse-shear-deformation effects is presented. The formulation uses the principle of virtual work. A unified set of non-linear strains that contains terms from both physical and tensorial strain measures is used. Linearized, perturbed equilibrium equations (stability equations) that describe the response of the plate just after buckling occurs are then derived after the application of several simplifying assumptions. These equations are then modified to allow the reference surface of the plate to be located at a distance z(sub c) from the centroidal surface. The implementation of the new theory into the VICONOPT exact buckling and vibration analysis and optimum design computer program is described as well. The terms of the plate stiffness matrix using both Classical Plate Theory (CPT) and first-order Shear-Deformation Plate Theory (SDPT) are presented. The necessary steps to include the effects of in-plane transverse and in-plane shear loads in the in-plane stability equations are also outlined. Numerical results are presented using the newly implemented capability. Comparisons of results for several example problems with different loading states are made. Comparisons of analyses using both physical and tensorial strain measures as well as CPT and SDPF are also made. Results comparing the computational effort required by the new analysis to that of the analysis currently in the VICONOPT program are presented. The effects of including terms related to in-plane transverse and in-plane shear loadings in the in-plane stability equations are also examined. Finally, results of a design-optimization study of two different cylindrical shells subject to uniform axial compression are presented.

McGowan, David Michael↗