Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “gradient methods”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS1978

Inviscid, nonadiabatic flow fields over blunt, sonic corner bodies for outer planet entry conditions by a method of integral relations

An investigation has been made into the ability of a method of integral relations to calculate inviscid zero degree angle of attack, radiative heating distributions over blunt, sonic corner bodies for some representative outer planet entry conditions is investigated. Comparisons have been made with a more detailed numerical method, a time asymptotic technique, using the same equilibrium chemistry and radiation transport subroutines. An effort to produce a second order approximation (two-strip) method of integral relations code to aid in this investigation is also described and a modified two-strip routine is presented. Results indicate that the one-strip method of integral relations cannot be used to obtain accurate estimates of the radiative heating distribution because of its inability to resolve thermal gradients near the wall. The two-strip method can sometimes be used to improve these estimates; however, the two-strip method has only a small range of conditions over which it will yield significant improvement over the one-strip method.

Gnoffo, P. A.↗
NASA NTRS2018

Solution of Ordinary Differential Equations in Gradient-Based Multidisciplinary Design Optimization

A gradient-based approach to multidisciplinary design optimization enables efficient scalability to large numbers of design variables. However, the need for derivatives causes difficulties when integrating ordinary differential equations (ODEs) in models. To simplify this, we propose the use of the general linear methods framework, which unifies all Runge-Kutta and linear multistep methods. This approach enables rapid implementation of integration methods without the need to differentiate each one, even in a gradient-based optimization context. We also develop a new parallel time integration algorithm that enables vectorization across time steps. We present a set of benchmarking results using a stiff ODE, a non-stiff nonlinear ODE, and an orbital dynamics ODE, and compare integration methods. In a modular gradient-based multidisciplinary design optimization context, we find that the new parallel time integration algorithm with high-order implicit methods, especially Gauss-Legendre collocation, is the best choice for a broad range of problems.

Hwang, John T.↗
NASA NTRS2006

Single- and Multiple-Objective Optimization with Differential Evolution and Neural Networks

Genetic and evolutionary algorithms have been applied to solve numerous problems in engineering design where they have been used primarily as optimization procedures. These methods have an advantage over conventional gradient-based search procedures became they are capable of finding global optima of multi-modal functions and searching design spaces with disjoint feasible regions. They are also robust in the presence of noisy data. Another desirable feature of these methods is that they can efficiently use distributed and parallel computing resources since multiple function evaluations (flow simulations in aerodynamics design) can be performed simultaneously and independently on ultiple processors. For these reasons genetic and evolutionary algorithms are being used more frequently in design optimization. Examples include airfoil and wing design and compressor and turbine airfoil design. They are also finding increasing use in multiple-objective and multidisciplinary optimization. This lecture will focus on an evolutionary method that is a relatively new member to the general class of evolutionary methods called differential evolution (DE). This method is easy to use and program and it requires relatively few user-specified constants. These constants are easily determined for a wide class of problems. Fine-tuning the constants will off course yield the solution to the optimization problem at hand more rapidly. DE can be efficiently implemented on parallel computers and can be used for continuous, discrete and mixed discrete/continuous optimization problems. It does not require the objective function to be continuous and is noise tolerant. DE and applications to single and multiple-objective optimization will be included in the presentation and lecture notes. A method for aerodynamic design optimization that is based on neural networks will also be included as a part of this lecture. The method offers advantages over traditional optimization methods. It is more flexible than other methods in dealing with design in the context of both steady and unsteady flows, partial and complete data sets, combined experimental and numerical data, inclusion of various constraints and rules of thumb, and other issues that characterize the aerodynamic design process. Neural networks provide a natural framework within which a succession of numerical solutions of increasing fidelity, incorporating more realistic flow physics, can be represented and utilized for optimization. Neural networks also offer an excellent framework for multiple-objective and multi-disciplinary design optimization. Simulation tools from various disciplines can be integrated within this framework and rapid trade-off studies involving one or many disciplines can be performed. The prospect of combining neural network based optimization methods and evolutionary algorithms to obtain a hybrid method with the best properties of both methods will be included in this presentation. Achieving solution diversity and accurate convergence to the exact Pareto front in multiple objective optimization usually requires a significant computational effort with evolutionary algorithms. In this lecture we will also explore the possibility of using neural networks to obtain estimates of the Pareto optimal front using non-dominated solutions generated by DE as training data. Neural network estimators have the potential advantage of reducing the number of function evaluations required to obtain solution accuracy and diversity, thus reducing cost to design.

Rai, Man Mohan↗
NASA NTRS1974

Heat flow and geothermal potential of the East Mesa KGRA, Imperial Valley, California

The East Mesa KGRA (Known Geothermal Resource Area) is located in the southeast part of the Imperial Valley, California, and is roughly 150 kilometers square in areal extent. A new heat flow technique which utilizes temperature gradient measurements across best clays is presented and shown to be as accurate as conventional methods for the present study area. Utilizing the best clay gradient technique, over 70 heat flow determinations have been completed within and around the East Mesa KGRA. Background heat flow values range from 1.4 to 2.4 hfu (1 hfu = .000001 cal. per square centimeter-second) and are typical of those throughout the Basin and Range province. Heat flow values for the northwest lobe of the KGRA (Mesa anomaly) are as high as 7.9 hfu, with the highest values located near gravity and seismic noise maxima and electrical resistivity minima. An excellent correlation exists between heat flow contours and faults defined by remote sensing and microearthquake monitoring.

Swanberg, C. A.↗
NASA NTRS1975

Numerical computation of viscous flows on the lee side of blunt shapes flying at supersonic speeds

A numerical method for solving the parabolic approximation to the steady-state compressible Navier-Stokes equations is examined. The approximation neglects only the streamwise gradients of shear stress. An implicit finite difference method is used which advances the solution downstream from an initial data surface and determines the complete viscous-inviscid flow between the body and bow shock wave. It is necessary that the inviscid portion of the flow field be supersonic. Crossflow separation is determined as part of the solution. The method is applied to a 15 deg sphere-cone at 15 deg angle of attack, and the results are compared with an inviscid method-of-characteristics calculation.

Rakich, J. V.↗
NASA NTRS1983

Mixed solar wind originating from coronal regions of different temperatures

Ionization states of elements in the solar wind are often used to determine thermal gradients in the lower corona. This method is based on the assumption, that in the beginning, solar wind material has a homogeneous temperature determining the original charge state of elements. Features in M/Q-spectra which might appear if the above assumption is violated are investigated and compared with observational evidence.

Bochsler, P.↗
NASA NTRS1984

Simulated flight through JAWS wind shear - In-depth analysis results

The Joint Airport Weather Studies (JAWS) field experiment was carried out in 1982 near Denver. An analysis is presented of aircraft performance in the three-dimensional wind fields. The fourth dimension, time, is not considered. The analysis seeks to prepare computer models of microburst wind shear from the JAWS data sets for input to flight simulators and for research and development of aircraft control systems and operational procedures. A description is given of the data set and the method of interpolating velocities and velocity gradients for input to the six-degrees-of-freedom equations governing the motion of the aircraft. The results of the aircraft performance analysis are then presented, and the interpretation classifies the regions of shear as severe, moderate, or weak. Paths through the severe microburst of August 5, 1982, are then recommended for training and operational applications. Selected subregions of the flow field defined in terms of planar sections through the wind field are presented for application to simulators with limited computer storage capacity, that is, for computers incapable of storing the entire array of variables needed if the complete wind field is programmed.

Frost, W.↗
NASA NTRS1988

Transient natural convection heat and mass transfer in a rectangular enclosure - A numerical analysis

A numerical analysis of transient heat and solute transport across a rectangular cavity with combined horizontal temperature and concentration gradients is performed by a numerical method based on the SIMPLE. Numerical results show that the average Nusselt and Sherwood numbers both decrease markedly when the solutal and thermal buoyancy forces act in the opposite directions. When the solutal and thermal buoyancy forces act in the same directions, however, the average Sherwood number increases significantly and yet the average Nusselt number decreases slightly.

Han, Samuel S.↗
NASA NTRS1989

Identification of atmospheric fronts over the ocean with microwave measurements of water vapor and rain

This paper describes some basic research techniques and algorithms developed to diagnose fronts in cyclonic storms over the ocean with data from satellite-borne microwave radiometers. Methods are developed for flagging strong gradients in integrated atmospheric water vapor and the presence of rain by using data from the SSMR on board the polar orbiting Seasat and Nimbus-7 satellites. Examination of 65 frontal systems showed that the water vapor gradient flag correctly identified 86 percent of the fronts, while the precipitation flagged 91 percent. The two types of flags emphasize different portions of the cyclone and are therefore complementary. Ultimately, these techniques are intended for operational use with data from the Special Sensor Microwave Imager which was launched in June 1987 on a satellite in the Defense Meteorological Satellite Program (DMSP).

Katsaros, Kristina B.↗
NASA NTRS1990

Constrained trajectory optimization for kinematically redundant arms

Two velocity optimization schemes for resolving redundant joint configurations are compared. The Extended Moore-Penrose Technique minimizes the joint velocities and avoids obstacles indirectly by adjoining a cost gradient to the solution. A new method can incorporate inequality constraints directly to avoid obstacles and singularities in the workspace. A four-link arm example is used to illustrate singularity avoidance while tracking desired end-effector paths.

Carignan, Craig R.↗
NASA NTRS1990

Validation of parameterization scheme for eddy diffusion from satellite data

The eddy diffusion coefficient K(yy) has been calculated usign LIMS for the months of December 1978 and January and February 1979. Two methods have been used. The first implements the suggestion made by Tung (1987) to parameterize the eddy transport as a diffusive process along isentropes. The second method integrates the equation relating the parcel displacements to the eddy velocity fields. The latter method uses a filtering on both space and time domains to isolate transients and is referred to as the 'spectral method'. Results from the first method are shown to be reliable only for quiescent periods, breaking down when the meridional gradient of potential vorticity is negligible. Results from the two methods are in agreement only for very disturbed conditions, when transience is readily isolated. It is concluded that the parameterizations suggested for eddy transport and calculated in this paper may be meaningful for quiet periods, but are not reliable for unsteady and very large amplitude disturbances.

Sassi, F.↗
NASA NTRS1995

A deterministic particle method for one-dimensional reaction-diffusion equations

We derive a deterministic particle method for the solution of nonlinear reaction-diffusion equations in one spatial dimension. This deterministic method is an analog of a Monte Carlo method for the solution of these problems that has been previously investigated by the author. The deterministic method leads to the consideration of a system of ordinary differential equations for the positions of suitably defined particles. We then consider the time explicit and implicit methods for this system of ordinary differential equations and we study a Picard and Newton iteration for the solution of the implicit system. Next we solve numerically this system and study the discretization error both analytically and numerically. Numerical computation shows that this deterministic method is automatically adaptive to large gradients in the solution.

Mascagni, Michael↗
NASA NTRS2000

Global Design Optimization for Aerodynamics and Rocket Propulsion Components

Modern computational and experimental tools for aerodynamics and propulsion applications have matured to a stage where they can provide substantial insight into engineering processes involving fluid flows, and can be fruitfully utilized to help improve the design of practical devices. In particular, rapid and continuous development in aerospace engineering demands that new design concepts be regularly proposed to meet goals for increased performance, robustness and safety while concurrently decreasing cost. To date, the majority of the effort in design optimization of fluid dynamics has relied on gradient-based search algorithms. Global optimization methods can utilize the information collected from various sources and by different tools. These methods offer multi-criterion optimization, handle the existence of multiple design points and trade-offs via insight into the entire design space, can easily perform tasks in parallel, and are often effective in filtering the noise intrinsic to numerical and experimental data. However, a successful application of the global optimization method needs to address issues related to data requirements with an increase in the number of design variables, and methods for predicting the model performance. In this article, we review recent progress made in establishing suitable global optimization techniques employing neural network and polynomial-based response surface methodologies. Issues addressed include techniques for construction of the response surface, design of experiment techniques for supplying information in an economical manner, optimization procedures and multi-level techniques, and assessment of relative performance between polynomials and neural networks. Examples drawn from wing aerodynamics, turbulent diffuser flows, gas-gas injectors, and supersonic turbines are employed to help demonstrate the issues involved in an engineering design context. Both the usefulness of the existing knowledge to aid current design practices and the need for future research are identified.

Shyy, Wei↗
NASA NTRS2003

Optimal Design and Damage Tolerance Verification of an Isogrid Structure for Helicopter Application

A composite isogrid panel design for application to a rotorcraft fuselage is presented. An optimum panel design for the lower fuselage of the rotorcraft that is subjected to combined in-plane compression and shear loads was generated using a design tool that utilizes a smeared-stiffener theory in conjunction with a genetic algorithm. A design feature was introduced along the edges of the panel that facilitates introduction of loads into the isogrid panel without producing undesirable local bending gradients. A low-cost manufacturing method for the isogrid panel that incorporates these design details is also presented. Axial compression tests were conducted on the undamaged and low-speed impact damaged panels to demonstrate the damage tolerance of this isogrid panel. A combined loading test fixture was designed and utilized that allowed simultaneous application of compression and shear loads to the test specimen. Results from finite element analyses are presented for the isogrid panel designs and these results are compared with experimental results. This study illustrates the isogrid concept to be a viable candidate for application to the helicopter lower fuselage structure.

Baker, Donald J.↗
NASA NTRS2010

Measurements and Simulations of Nadir-Viewing Radar Returns from the Melting Layer at X- and W-Bands

Simulated radar signatures within the melting layer in stratiform rain, namely the radar bright band, are checked by means of comparisons with simultaneous measurements of the bright band made by the EDOP (X-band) and CRS (W-band) airborne Doppler radars during the CRYSTAL-FACE campaign in 2002. A stratified-sphere model, allowing the fractional water content to vary along the radius of the particle, is used to compute the scattering properties of individual melting snowflakes. Using the effective dielectric constants computed by the conjugate gradient-fast Fourier transform (CGFFT) numerical method for X and W bands, and expressing the fractional water content of melting particle as an exponential function in particle radius, it is found that at X band the simulated radar bright-band profiles are in an excellent agreement with the measured profiles. It is also found that the simulated W-band profiles usually resemble the shapes of the measured bright-band profiles even though persistent offsets between them are present. These offsets, however, can be explained by the attenuation caused by cloud water and water vapor at W band. This is confirmed by the comparisons of the radar profiles made in the rain regions where the un-attenuated W-band reflectivity profiles can be estimated through the X- and W band Doppler velocity measurements. The bright-band model described in this paper has the potential to be used effectively for both radar and radiometer algorithms relevant to the TRMM and GPM satellite missions.

Liao, Liang↗
NASA NTRS2012

Validation of Satellite-Based Objective Overshooting Cloud-Top Detection Methods Using CloudSat Cloud Profiling Radar Observations

Two satellite infrared-based overshooting convective cloud-top (OT) detection methods have recently been described in the literature: 1) the 11-mm infrared window channel texture (IRW texture) method, which uses IRW channel brightness temperature (BT) spatial gradients and thresholds, and 2) the water vapor minus IRW BT difference (WV-IRW BTD). While both methods show good performance in published case study examples, it is important to quantitatively validate these methods relative to overshooting top events across the globe. Unfortunately, no overshooting top database currently exists that could be used in such study. This study examines National Aeronautics and Space Administration CloudSat Cloud Profiling Radar data to develop an OT detection validation database that is used to evaluate the IRW-texture and WV-IRW BTD OT detection methods. CloudSat data were manually examined over a 1.5-yr period to identify cases in which the cloud top penetrates above the tropopause height defined by a numerical weather prediction model and the surrounding cirrus anvil cloud top, producing 111 confirmed overshooting top events. When applied to Moderate Resolution Imaging Spectroradiometer (MODIS)-based Geostationary Operational Environmental Satellite-R Series (GOES-R) Advanced Baseline Imager proxy data, the IRW-texture (WV-IRW BTD) method offered a 76% (96%) probability of OT detection (POD) and 16% (81%) false-alarm ratio. Case study examples show that WV-IRW BTD.0 K identifies much of the deep convective cloud top, while the IRW-texture method focuses only on regions with a spatial scale near that of commonly observed OTs. The POD decreases by 20% when IRW-texture is applied to current geostationary imager data, highlighting the importance of imager spatial resolution for observing and detecting OT regions.

Bedka, Kristopher M.↗
NASA NTRS1985

On the electromagnetic scattering from infinite rectangular conducting grids

The study and development of two numerical techniques for the analysis of electromagnetic scattering from a rectangular wire mesh are described. Both techniques follow from one basic formulation and they are both solved in the spectral domain. These techniques were developed as a result of an investigation towards more efficient numerical computation for mesh scattering. These techniques are efficient for the following reasons: (a1) make use of the Fast Fourier Transform; (b2) they avoid any convolution problems by converting integrodifferential equations into algebraic equations; and (c3) they do not require inversions of any matrices. The first method, the SIT or Spectral Iteration Technique, is applied for regions where the spacing between wires is not less than two wavelengths. The second method, the SDCG or Spectral Domain Conjugate Gradient approach, can be used for any spacing between adjacent wires. A study of electromagnetic wave properties, such as reflection coefficient, induced currents and aperture fields, as functions of frequency, angle of incidence, polarization and thickness of wires is presented. Examples and comparisons or results with other methods are also included to support the validity of the new algorithms.

Christodoulou, C.↗