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 685 records · Page 38

Optimal trajectories for aeroassisted orbital transfer

Consideration is given to classical and minimax problems involved in aeroassisted transfer from high earth orbit (HEO) to low earth orbit (LEO). The transfer is restricted to coplanar operation, with trajectory control effected by means of lift modulation. The performance of the maneuver is indexed to the energy expenditure or, alternatively, the time integral of the heating rate. Firist-order optimality conditions are defined for the classical approach, as are a sequential gradient-restoration algorithm and a combined gradient-restoration algorithm. Minimization techniques are presented for the aeroassisted transfer energy consumption and time-delay integral of the heating rate, as well as minimization of the pressure. It is shown that the eigenvalues of the Jacobian matrix of the differential system is both stiff and unstable, implying that the sequential gradient restoration algorithm in its present version is unsuitable. A new method, involving a multipoint approach to the two-poing boundary value problem, is recommended.

Miele, A.↗

Aerothermal environment of a blunted three-dimensional nonaxisymmetric body at Mach 6.8

Measurements of surface pressure and heat flux at Mach 6.8 for a nonaxisymmetric body, typical of the forward fuselage of an entry vehicle, have been obtained in the Langley 8-ft high-temperature tunnel. Comparisons of these data with predictions obtained from two numerical methods, which treat the test medium as an ideal gas with a laminar boundary layer, are presented. In general, the numerical results are consistent with each other and follow the data in regions of laminar flow; however, the data for angle of attack 14.9 deg indicate transitional flow on most of the windward surface. Also, the two numerical methods disagree slightly in regions of zero or nearly zero pressure gradients.

Martin, F. W., Jr.↗

Spatial dependence of the local diffusion coefficient measured upstream of the November 12, 1978 interplanetary traveling shock

Characteristics of wuprathermal particles accelerated by quasi-parallel interplanetary traveling shocks have been generally explained in terms of a first order Fermi mechanism. Such models require diffusive scattering of particles upstream of the shock. This scattering is characterized by a local diffusion coefficient, kappa, which is determined by the local power density of waves in the upstream region. The dependence of the diffusion coefficient of suprathermal upstream protons on distance from the November 12, 1978 interplanetary traveling shock using a different approach is studied. Unlike previous studies this method, which is based on measurements of particle streaming and intensity gradients, does not rely on predictions. The local spatial variations of Kappa upstream of the November 12, 1978 shock have been chosen for study because the characteristics of this quasi-parallel shock have been extensively studied, and also because of its favorable geometry (i.e. B field nearly radial).

Gloeckler, G.↗

Prediction of high-resolution flowfields for rotorcraft aeroacoustics

A new approach to the vortex wake of a helicopter is described. This approach makes it possible to obtain highly accurate resolution of the flowfield experienced by the tail rotor with modest computational effort relative to alternative models. A full-span free wake model of the main rotor in a scheme is presented which reconstructs high resolution flow solutions from preliminary simulations with coarse resolution. A new reconstruction method uses local velocity correction terms to capture the steep velocity gradients characteristic of the the vortex-dominated incident flow. The approach is considered to be successful in capturing the high resolution flows necessary for analysis of rotor-wake/rotor interactions with substantially reduced computational cost.

Quackenbush, Todd R.↗

Multidisciplinary Design Techniques Applied to Conceptual Aerospace Vehicle Design

Multidisciplinary design optimization (MDO) is an emerging discipline within aerospace engineering. Its goal is to bring structure and efficiency to the complex design process associated with advanced aerospace launch vehicles. Aerospace vehicles generally require input from a variety of traditional aerospace disciplines - aerodynamics, structures, performance, etc. As such, traditional optimization methods cannot always be applied. Several multidisciplinary techniques and methods were proposed as potentially applicable to this class of design problem. Among the candidate options are calculus-based (or gradient-based) optimization schemes and parametric schemes based on design of experiments theory. A brief overview of several applicable multidisciplinary design optimization methods is included. Methods from the calculus-based class and the parametric class are reviewed, but the research application reported focuses on methods from the parametric class. A vehicle of current interest was chosen as a test application for this research. The rocket-based combined-cycle (RBCC) single-stage-to-orbit (SSTO) launch vehicle combines elements of rocket and airbreathing propulsion in an attempt to produce an attractive option for launching medium sized payloads into low earth orbit. The RBCC SSTO presents a particularly difficult problem for traditional one-variable-at-a-time optimization methods because of the lack of an adequate experience base and the highly coupled nature of the design variables. MDO, however, with it's structured approach to design, is well suited to this problem. The result of the application of Taguchi methods, central composite designs, and response surface methods to the design optimization of the RBCC SSTO are presented. Attention is given to the aspect of Taguchi methods that attempts to locate a 'robust' design - that is, a design that is least sensitive to uncontrollable influences on the design. Near-optimum minimum dry weight solutions are determined for the vehicle. A summary and evaluation of the various parametric MDO methods employed in the research are included. Recommendations for additional research are provided.

Olds, John Robert↗

Refraction of high frequency noise in an arbitrary jet flow

Refraction of high frequency noise by mean flow gradients in a jet is studied using the ray-tracing methods of geometrical acoustics. Both the two-dimensional (2D) and three-dimensional (3D) formulations are considered. In the former case, the mean flow is assumed parallel and the governing propagation equations are described by a system of four first order ordinary differential equations. The 3D formulation, on the other hand, accounts for the jet spreading as well as the axial flow development. In this case, a system of six first order differential equations are solved to trace a ray from its source location to an observer in the far field. For subsonic jets with a small spreading angle both methods lead to similar results outside the zone of silence. However, with increasing jet speed the two prediction models diverge to the point where the parallel flow assumption is no longer justified. The Doppler factor of supersonic jets as influenced by the refraction effects is discussed and compared with the conventional modified Doppler factor.

Khavaran, Abbas↗

Formulation for Simultaneous Aerodynamic Analysis and Design Optimization

An efficient approach for simultaneous aerodynamic analysis and design optimization is presented. This approach does not require the performance of many flow analyses at each design optimization step, which can be an expensive procedure. Thus, this approach brings us one step closer to meeting the challenge of incorporating computational fluid dynamic codes into gradient-based optimization techniques for aerodynamic design. An adjoint-variable method is introduced to nullify the effect of the increased number of design variables in the problem formulation. The method has been successfully tested on one-dimensional nozzle flow problems, including a sample problem with a normal shock. Implementations of the above algorithm are also presented that incorporate Newton iterations to secure a high-quality flow solution at the end of the design process. Implementations with iterative flow solvers are possible and will be required for large, multidimensional flow problems.

Hou, G. W.↗

Aerodynamic Shape Sensitivity Analysis and Design Optimization of Complex Configurations Using Unstructured Grids

A three-dimensional unstructured grid approach to aerodynamic shape sensitivity analysis and design optimization has been developed and is extended to model geometrically complex configurations. The advantage of unstructured grids (when compared with a structured-grid approach) is their inherent ability to discretize irregularly shaped domains with greater efficiency and less effort. Hence, this approach is ideally suited for geometrically complex configurations of practical interest. In this work the nonlinear Euler equations are solved using an upwind, cell-centered, finite-volume scheme. The discrete, linearized systems which result from this scheme are solved iteratively by a preconditioned conjugate-gradient-like algorithm known as GMRES for the two-dimensional geometry and a Gauss-Seidel algorithm for the three-dimensional; similar procedures are used to solve the accompanying linear aerodynamic sensitivity equations in incremental iterative form. As shown, this particular form of the sensitivity equation makes large-scale gradient-based aerodynamic optimization possible by taking advantage of memory efficient methods to construct exact Jacobian matrix-vector products. Simple parameterization techniques are utilized for demonstrative purposes. Once the surface has been deformed, the unstructured grid is adapted by considering the mesh as a system of interconnected springs. Grid sensitivities are obtained by differentiating the surface parameterization and the grid adaptation algorithms with ADIFOR (which is an advanced automatic-differentiation software tool). To demonstrate the ability of this procedure to analyze and design complex configurations of practical interest, the sensitivity analysis and shape optimization has been performed for a two-dimensional high-lift multielement airfoil and for a three-dimensional Boeing 747-200 aircraft.

Taylor, Arthur C., III↗

Wall relaxation in growing stems: comparison of four species and assessment of measurement techniques

This study was carried out to develop improved methods for measuring in-vivo stress relaxation of growing tissues and to compare relaxation in the stems of four different species. When water uptake by growing tissue is prevented, in-vivo stress relaxation occurs because continued wall loosening reduces wall stress and cell turgor pressure. With this procedure one may measure the yield threshold for growth (Y), the turgor pressure in excess of the yield threshold (P-Y), and the physiological wall extensibility (phi). Three relaxation techniques proved useful: "turgor-relaxation", "balance-pressure" and "pressure-block". In the turgor-relaxation method, water is withheld from growing tissue and the reduction in turgor is measured directly with the pressure probe. This technique gives absolute values for P and Y, but requires tissue excision. In the balance-pressure technique, the excised growing region is sealed in a pressure chamber, and the subsequent reduction in water potential is measured as the applied pressure needed to return xylem sap to the cut surface. This method is simple, but only measures (P-Y), not the individual values of P and Y. In the pressure-block technique, the growing tissue is sealed into a pressure chamber, growth is monitored continuously, and just sufficient pressure is applied to the chamber to block growth. The method gives high-resolution kinetics of relaxation and does not require tissue excision, but only measures (P-Y). The three methods gave similar results when applied to the growing stems of pea (Pisum sativum L.), cucumber (Cucumis sativus L.), soybean (Glycine max (L.) Merr.) and zucchini (Curcubita pepo L.) seedlings. Values for (P-Y) averaged between 1.4 and 2.7 bar, depending on species. Yield thresholds averaged between 1.3 and 3.0 bar. Compared with the other methods, relaxation by pressure-block was faster and exhibited dynamic changes in wall-yielding properties. The two pressure-chamber methods were also used to measure the internal water-potential gradient (between the xylem and the epidermis) which drives water uptake for growth. For the four species it was small, between 0.3 and 0.6 bar, and so did not limit growth substantially.

NASA Program Space Biology↗

Automation of POST Cases via External Optimizer and "Artificial p2" Calculation

During early conceptual design of complex systems, speed and accuracy are often at odds with one another. While many characteristics of the design are fluctuating rapidly during this phase there is nonetheless a need to acquire accurate data from which to down-select designs as these decisions will have a large impact upon program life-cycle cost. Therefore enabling the conceptual designer to produce accurate data in a timely manner is tantamount to program viability. For conceptual design of launch vehicles, trajectory analysis and optimization is a large hurdle. Tools such as the industry standard Program to Optimize Simulated Trajectories (POST) have traditionally required an expert in the loop for setting up inputs, running the program, and analyzing the output. The solution space for trajectory analysis is in general non-linear and multi-modal requiring an experienced analyst to weed out sub-optimal designs in pursuit of the global optimum. While an experienced analyst presented with a vehicle similar to one which they have already worked on can likely produce optimal performance figures in a timely manner, as soon as the "experienced" or "similar" adjectives are invalid the process can become lengthy. In addition, an experienced analyst working on a similar vehicle may go into the analysis with preconceived ideas about what the vehicle's trajectory should look like which can result in sub-optimal performance being recorded. Thus, in any case but the ideal either time or accuracy can be sacrificed. In the authors' previous work a tool called multiPOST was created which captures the heuristics of a human analyst over the process of executing trajectory analysis with POST. However without the instincts of a human in the loop, this method relied upon Monte Carlo simulation to find successful trajectories. Overall the method has mixed results, and in the context of optimizing multiple vehicles it is inefficient in comparison to the method presented POST's internal optimizer functions like any other gradient-based optimizer. It has a specified variable to optimize whose value is represented as optval, a set of dependent constraints to meet with associated forms and tolerances whose value is represented as p2, and a set of independent variables known as the u-vector to modify in pursuit of optimality. Each of these quantities are calculated or manipulated at a certain phase within the trajectory. The optimizer is further constrained by the requirement that the input u-vector must result in a trajectory which proceeds through each of the prescribed events in the input file. For example, if the input u-vector causes the vehicle to crash before it can achieve the orbital parameters required for a parking orbit, then the run will fail without engaging the optimizer, and a p2 value of exactly zero is returned. This poses a problem, as this "non-connecting" region of the u-vector space is far larger than the "connecting" region which returns a non-zero value of p2 and can be worked on by the internal optimizer. Finding this connecting region and more specifically the global optimum within this region has traditionally required the use of an expert analyst.

Dees, Patrick D.↗

Exergy-based Sensitivity Analysis of the Generic Hypersonic Vehicle using FUN3D

In this paper, the implementation of an exergy-based objective function and its adjoint gradient into NASA’s FUN3D solver is discussed and verified. In order to verify that the exergy-based functional is properly implemented, it is used to predict the drag of the Generic Hypersonic Vehicle (GHV), which is then compared to more traditional force-based drag predictions. In addition to the functional implementation, FUN3D’s adjoint capability was extended to obtain sensitivities. Results were verified using FUN3D’s native complex step method for di↵erentiation using a generic wing configuration. The complex and adjoint gradients yielded discrete agreement demonstrating correct implementation and that the functional can be used for gradient-based multidisciplinary analysis and optimization. Next, various trade studies are conducted on the GHV to understand the design space of the vehicle. Finally, an inverse design problem is solved to verify the utilized design optimization framework which is ready to be deployed for exergy-based optimizations in future work.

Neal L. Novotny↗

The Control of Compressible Dynamic Stall Using Adaptive Airfoils

Ongoing research has shown that compressible dynamic stall can be caused by vastly different mechanisms even for small changes in flow conditions. For example, at low Mach numbers (M less than 0.3) and Reynolds numbers (less than 1 x 10(exp 6)), the bursting of the laminar separation bubble induces dynamic stall. At a slightly higher Mach number (M = 0.45) and around the same Reynolds number, shock induced separation can cause dynamic stall. Also, the mechanism changes from that of laminar separation bubble bursting to that due to excessive adverse pressure gradient, with increase of Reynolds number. Complex interactions occurring between the local supersonic flow and the bubble lead to another possible mechanism of dynamic stall. Since all these mechanisms may be encountered by a rotor-blade during a single cycle, there is a strong need to devise a reliable flow control method for use under such dramatically varying conditions. Whereas it is a Ali 'Beneficent challenge, identification of the fact that these mechanisms are all in some way related to the airfoil leading edge adverse pressure gradient, has resulted in an effective method of flow control to be described.

Chandasekhara, M. S.↗

Simple computer method provides contours for radiological images

Computer is provided with information concerning boundaries in total image. Gradient of each point in digitized image is calculated with aid of threshold technique; then there is invoked set of algorithms designed to reduce number of gradient elements and to retain only major ones for definition of contour.

Newell, J. D.↗

Shadowgraph Study of Gradient Driven Fluctuations

A fluid or fluid mixture, subjected to a vertical temperature and/or concentration gradient in a gravitational field, exhibits greatly enhanced light scattering at small angles. This effect is caused by coupling between the vertical velocity fluctuations due to thermal energy and the vertically varying refractive index. Physically, small upward or downward moving regions will be displaced into fluid having a refractive index different from that of the moving region, thus giving rise to the enhanced scattering. The scattered intensity is predicted to vary with scattering wave vector q, as q(sup -4), for sufficiently large q, but the divergence is quenched by gravity at small q. In the absence of gravity, the long wavelength fluctuations responsible for the enhanced scattering are predicted to grow until limited by the sample dimensions. It is thus of interest to measure the mean-squared amplitude of such fluctuations in the microgravity environment for comparison with existing theory and ground based measurements. The relevant wave vectors are extremely small, making traditional low-angle light scattering difficult or impossible because of stray elastically scattered light generated by optical surfaces. An alternative technique is offered by the shadowgraph method, which is normally used to visualize fluid flows, but which can also serve as a quantitative tool to measure fluctuations. A somewhat novel shadowgraph apparatus and the necessary data analysis methods will be described. The apparatus uses a spatially coherent, but temporally incoherent, light source consisting of a super-luminescent diode coupled to a single-mode optical fiber in order to achieve extremely high spatial resolution, while avoiding effects caused by interference of light reflected from the various optical surfaces that are present when using laser sources. Results obtained for a critical mixture of aniline and cyclohexane subjected to a vertical temperature gradient will be presented. The sample was confined between two horizontal parallel sapphire plates with a vertical spacing of 1 mm. The temperatures of the sapphire plates were controlled by independent circulating water loops that used Peltier devices to add or remove heat from the room air as required. For a mixture with a temperature gradient, two effects are involved in generating the vertical refractive index gradient, namely thermal expansion and the Soret effect, which generates a concentration gradient in response to the applied temperature gradient. For the aniline/cyclohexane system, the denser component (aniline) migrates toward the colder surface. Consequently, when heating from above, both effects result in the sample density decreasing with altitude and are stabilizing in the sense that no convective motion occurs regardless of the magnitude of the applied temperature gradient. The Soret effect is strong near a binary liquid critical point, and thus the dominant effect is due to the induced concentration gradient. The results clearly show the divergence at low q and the predicted gravitational quenching. Results obtained for different applied temperature gradients at varying temperature differences from the critical temperature, clearly demonstrate the predicted divergence of the thermal diffusion ratio. Thus, the more closely the critical point is approached, the smaller becomes the temperature gradient required to generate the same signal. Two different methods have been used to generate pure concentration gradients. In the first, a sample cell was filled with a single fluid, ethylene glycol, and a denser miscible fluid, water, was added from below thus establishing a sharp interface to begin the experiment. As time went on the two fluids diffused into each other, and large amplitude fluctuations were clearly observed at low q. The effects of gravitational quenching were also evident. In the second method, the aniline/cyclohexane sample was used, and after applying a vertical temperature gradient for several hours, the top and bottom temperatures were set equal and the thermal gradient died on a time scale of seconds, leaving the Soret induced concentration gradient in place. Again, large-scale fluctuations were observed and died away slowly in amplitude as diffusion destroyed the initial concentration gradient.

Cannell, David↗

Skin friction drag measurements by LDV

A laser Doppler velocimeter employing a microscope objective as the receiving lens has been developed for measuring fluid velocity inside the boundary layer flow field with a spatial resolution of 40 microns. The method was applied for direct measurement of aerodynamic skin friction drag from the measured velocity gradient at the wall. Experimental results obtained on skin friction and on velocity components in a turbulent boundary layer on a low speed wind tunnel showed good agreement with previously reported data using conventional instruments such as hot-wire anemometers and Preston tubes. The method thus provides a tool for measurement and control of skin friction on aerodynamic bodies without perturbing the flow field.

Mazumder, M. K.↗

Compression Stockings May Ameliorate Orthostatic Intolerance in Astronauts After Short-Duration Space Flight

Orthostatic intolerance following spaceflight has been observed since the early days of manned spaceflight, and no countermeasure has been 100% effective. During re-entry NASA astronauts currently wear an inflatable anti-gravity suit (AGS) which compresses the legs and abdomen, but this device is uncomfortable and loses effectiveness upon egress from the Space Shuttle. We previously reported that foot-to-thigh, gradient compression stockings were comfortable and effective during standing after Shuttle missions. More recently we showed in a ground-based model of spaceflight that the addition of splanchnic compression to the foot-to-thigh compression stockings, creating foot-to-breast high compression, improved orthostatic tolerance in hypovolemic subjects to a level similar to the AGS. Purpose: To evaluate a new three-piece, foot-to-breast high gradient compression garment as a countermeasure to post-spaceflight orthostatic intolerance. Methods: Fourteen astronauts completed this experiment (7 control, 7 treatment) following Space Shuttle missions lasting 12-16 days. Treatment subjects were custom-fitted for a three-piece, foot-to-breast high compression garment consisting of shorts and foot-to-thigh stockings. The garments were constructed to provide 55 mmHg compression at the ankle and decreased gradually to 15 mmHg over the abdomen. Orthostatic testing occurred ~30 days before flight (without garments) and ~2 hours after flight (with garments for treatment group only) on landing day. Blood pressure (BP) and heart rate (HR) were acquired for 2 minutes while the subject lay prone and then for 3.5 minutes after the subject stood. Data are reported as mean +/- SE. Results: The compression garment successfully prevented the tachycardia and hypotension typically seen post-spaceflight. On landing day, treatment subjects had a smaller change in HR (11+/-1 vs. 21+/-4 beats/min, p< or =0.05) and no decrease in systolic BP (2+/-4 vs. -9+/-2 mmHg, p< or =0.05). Garments also received good comfort ratings and were relatively easy to don. Conclusion: In this small group of astronauts, foot-to-breast high gradient compression garments seem to have prevented these negative effects of spaceflight on the cardiovascular responses to standing.

Platts, Steven H.↗

Design optimization of natural laminar flow bodies in compressible flow

An optimization method has been developed to design axisymmetric body shapes such as fuselages, nacelles, and external fuel tanks with increased transition Reynolds numbers in subsonic compressible flow. The new design method involves a constraint minimization procedure coupled with analysis of the inviscid and viscous flow regions and linear stability analysis of the compressible boundary-layer. In order to reduce the computer time, Granville's transition criterion is used to predict boundary-layer transition and to calculate the gradients of the objective function, and linear stability theory coupled with the e(exp n)-method is used to calculate the objective function at the end of each design iteration. Use of a method to design an axisymmetric body with extensive natural laminar flow is illustrated through the design of a tiptank of a business jet. For the original tiptank, boundary layer transition is predicted to occur at a transition Reynolds number of 6.04 x 10(exp 6). For the designed body shape, a transition Reynolds number of 7.22 x 10(exp 6) is predicted using compressible linear stability theory coupled with the e(exp n)-method.

Dodbele, Simha S.↗

The NTF Inlet Guide Vanes Thermal Gradient Problem and Its Mitigation

The National Transonic Facility (NTF) utilizes Inlet Guide Vanes (IGV) to provide precise, quick response Mach number control for the tunnel. During cryogenic operations, the massive IGV structure can experience large thermal gradients, measured as "Delta T or (Delta)T", between the IGV ring and its support structure called the transfer case. If these temperature gradients are too large, the IGV structure can be stressed beyond its safety limit and cease operation. In recent years, (Delta)T readings exceeding the prescribed safety limits were observed frequently during cryogenic operations, particularly during model access. The tactical operation methods of the tunnel to minimize (Delta)T did not always succeed. One obvious option to remedy this condition is to warm up the IGV structure by disabling the main drive operation, but this "natural" warm up method can takes days in some cases, resulting in productivity loss. This paper documents the thermal gradient problem associated with the IGV structure during cryogenic operation and how the facility has recently achieved an acceptable mitigation which has resulted in improved efficiency of operations.

Venkat, Venki S.↗