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 1,189 records · Page 66

The multi-zone calculation of turbomachinery flows. II - The multi-zone calculation of the turbulent, two-specie flow through the SSME HPFTP first and second stage cavities

A multi-zone Navier-Stokes methodology to calculate the two-specie flow through the first and second stage cooling cavities of the Space Shuttle Main Engine (SSME) high pressure fuel turbopump (HPFTP) is developed. A simplified two-component fluid formulation is used to model the interaction of coolant and hot gas. Johnston's secant approximation is used to define an appropriate near wall velocity for use in a three-dimensional law of the wall. The basic Navier-Stokes algorithm used is a finite-volume, predictor-corrector algorithm which uses a pressure correction technique. A multi-zone method is used to partition each cavity into easily handled subdomains. The results show that coolant flow is pumped up the turbine wheel for both cavities, creating a region of large temperature gradients on the turbine shank.

Williams, M.↗

Vapor concentration measurement with photothermal deflectometry

Theoretical and experimental results for using the photothermal deflection technique to measure vapor species concentration, while minimizing the disturbance of the transport (material) parameters due to vapor heating, are developed and described. In contrast to common practice, the above constraints require using a pump-beam duty cycle of less than 50 percent. The theoretical description of the shortened heating time is based on a step-function formulation of the pumping cycle. The results are obtained as closed-form solutions of the energy equation for many chopping cycles until steady state is reached, by use of a Green's-function method. The Euler formulation of the Fermat principle is used to calculate the deflection angle. The equations are expanded to include the effects of vapor velocity on both the temperature and temperature gradient profiles. The effects of finite (unfocused) pump and probe beams and thermal (Soret) diffusion are also accounted for. Excellent agreement between theory and experiment is obtained.

Banish, R. Michael↗

Low-speed flowfield characterization by infrared measurements of surface temperatures

An experimental program was aimed at identifying areas in low speed aerodynamic research where infrared imaging systems can make significant contributions. Implementing a new technique, a long electrically heated wire was placed across a laminar jet. By measuring the temperature distribution along the wire with the IR imaging camera, the flow behavior was identified. Furthermore, using Nusselt number correlations, the velocity distribution could be deduced. The same approach was used to survey wakes behind cylinders in a wind-tunnel. This method is suited to investigate flows with position dependent velocities, e.g., boundary layers, confined flows, jets, wakes, and shear layers. It was found that the IR imaging camera cannot accurately track high gradient temperature fields. A correlation procedure was devised to account for this limitation. Other wind-tunnel experiments included tracking the development of the laminar boundary layer over a warmed flat plate by measuring the chordwise temperature distribution. This technique was applied also to the flow downstream from a rearward facing step. Finally, the IR imaging system was used to study boundary layer behavior over an airfoil at angles of attack from zero up to separation. The results were confirmed with tufts observable both visually and with the IR imaging camera.

Gartenberg, E.↗

Thermal calibration and analysis of a Bridgman flight furnace

A combination of computer modeling and experimental measurements is proposed which permits the calibration of the temperature profile for a spacecraft flight furnace. The Fluent code based on the finite difference method uses the furnace set-point temperatures as boundary conditions to model a Bridgman directional solidification furnace. Experimental measurements are employed as a comparison for the numerical analysis, and the radial gradients and wall temperatures are selected when the calculated values are sufficiently close to those from measurements. The calculations are found to agree with experimental data, and some results of the modeling are given which relate to specific furnace parameters. The CFD approach facilitates the transposition of furnace calibration to microgravitational operation in addition to providing precise initial radial temperature gradients and furnace-wall temperatures.

Knuteson, D. J.↗

Microstructural development during directional solidification of peritectic alloys

Despite the widespread commercial use of peritectic alloys (e.g., steels, brass, bronze, intermetallic compounds, Co based superalloys and A3B type superconductors), the characterization of the microstructural development during directional solidification of peritectics has historically lagged behind similar efforts directed towards other types of binary invariant reactions such as eutectic or monotectic. A wide variety of possible microstructures has been shown to form in peritectics depending upon the imposed temperature gradient, G, the solidification velocity, V, as well as the presence or absence of convection in the melt. This has important technological implications since many commercially important alloys exhibit peritectics and processing methods such as casting and welding often involve widely changing conditions. It has been the aim of this project to examine, in a systematic fashion, both experimentally and theoretically, the influence of gravitationally driven convection on segregation and microstructural development during solidification in peritectic systems under terrestrial conditions. The scientific results of the project will be used to establish ground based data in support of a meaningful microgravity flight experiment.

Lograsso, Thomas A.↗

Effects of spatial order of accuracy on the computation of vortical flowfields

The objective of the present investigation is to assess the effect of the spatial order of accuracy used for the evaluation of the inviscid fluxes on the resolution of higher order quantitites, such as velocity gradients. The viscous terms are computed as second-order accurate with central difference formulas, even though for the explicit part of the algorithm higher order approximations may be used. A viscous/inviscid method is used, and the outer part of the flowfield is computed with the inviscid flow equations. The viscous boundary-layer type flow region close to the body surface is computed with an algebraic eddy viscosity model. Results obtained with the conservative and nonconservative formulations and the viscous/inviscid approach are compared with available experimental data. The effect of grid refinement on the accuracy of the solution is also presented.

Ekaterinaris, J. A.↗

Temperature Measurements Reveal Accretion Of Ice On Airfoils

Method of detecting accretion of ice on airfoil based on measurement of difference between temperatures at icing and nonicing locations on airfoil. Release of latent heat of fusion gives rise to characteristic gradient of temperature. Temperatures and difference between them measured by use of contact sensors or remote sensor. Measurements used to trigger alarms and/or deicing equipment.

Hansman, Robert J., Jr.↗

Pioneer 10 and 11 Data Analysis

This report finishes the work of NASA Grant NAS2-153, which supported data analysis for the UCSD instruments on Pioneers 10 and 11. The data analyzed under this grant span 22 years of interplanetary measurements in the inner and outer heliosphere. The UCSD instruments made their mark in cosmic ray research based upon their high energy thresholds, directional responses, and reliable data streams. one of their primary scientific objectives concerns the size, configuration, and time behavior of the heliosphere. The size scale is inferred from the radial intensity gradient, which is measured between the two spacecraft and extrapolated to interstellar intensity levels at the cosmic ray modulation boundary. This boundary still eludes us, and its position, motion, and the best method of extrapolation are ongoing problems. Current projections place the boundary beyond 100 AU, which may be beyond the termination shock, and raises the question of possible modulation in the heliosheath. Probably our only hope of seeing this region in the immediate future rides on the possibility that the boundary will move inward. Our instruments have recorded many Forbush, or transient, decreases in the outer heliosphere. These observations led us to a model that attributes many of the decreases to solar wind stream-stream interactions, and relates the cosmic ray variations to the locally observed magnetic field magnitude. As the cosmic ray variations in this model result only from topological changes in the modulation integral, the model is a tool for studying the possibility that the 11 year cosmic ray modulation cycle can be accounted for by a superposition of Forbush decreases. The cosmic ray angular distribution function is measurable, given a good telemetry rate, by the UCSD Cerenkov detector which counts particles of energy greater than 500 MeV/n. We obtained statistically significant samples from 1 to 9 AU, at 13 AU, and at 34 AU. The anisotropy tends to be a few tenths of a per cent at all radial distances. A quasiperiodic variation in the east-west anisotropy with period of about 50 days remains unexplained.

Fillius, Walker↗

Noninvasive estimation of transmitral pressure drop across the normal mitral valve in humans: importance of convective and inertial forces during left ventricular filling

OBJECTIVES: We hypothesized that color M-mode (CMM) images could be used to solve the Euler equation, yielding regional pressure gradients along the scanline, which could then be integrated to yield the unsteady Bernoulli equation and estimate noninvasively both the convective and inertial components of the transmitral pressure difference. BACKGROUND: Pulsed and continuous wave Doppler velocity measurements are routinely used clinically to assess severity of stenotic and regurgitant valves. However, only the convective component of the pressure gradient is measured, thereby neglecting the contribution of inertial forces, which may be significant, particularly for nonstenotic valves. Color M-mode provides a spatiotemporal representation of flow across the mitral valve. METHODS: In eight patients undergoing coronary artery bypass grafting, high-fidelity left atrial and ventricular pressure measurements were obtained synchronously with transmitral CMM digital recordings. The instantaneous diastolic transmitral pressure difference was computed from the M-mode spatiotemporal velocity distribution using the unsteady flow form of the Bernoulli equation and was compared to the catheter measurements. RESULTS: From 56 beats in 16 hemodynamic stages, inclusion of the inertial term ([deltapI]max = 1.78+/-1.30 mm Hg) in the noninvasive pressure difference calculation significantly increased the temporal correlation with catheter-based measurement (r = 0.35+/-0.24 vs. 0.81+/-0.15, p< 0.0001). It also allowed an accurate approximation of the peak pressure difference ([deltapc+I]max = 0.95 [delta(p)cathh]max + 0.24, r = 0.96, p<0.001, error = 0.08+/-0.54 mm Hg). CONCLUSIONS: Inertial forces are significant components of the maximal pressure drop across the normal mitral valve. These can be accurately estimated noninvasively using CMM recordings of transmitral flow, which should improve the understanding of diastolic filling and function of the heart.

Non-NASA Center↗

A Comparative Study of Probability Collectives Based Multi-agent Systems and Genetic Algorithms

We compare Genetic Algorithms (GA's) with Probability Collectives (PC), a new framework for distributed optimization and control. In contrast to GA's, PC-based methods do not update populations of solutions. Instead they update an explicitly parameterized probability distribution p over the space of solutions. That updating of p arises as the optimization of a functional of p. The functional is chosen so that any p that optimizes it should be p peaked about good solutions. The PC approach works in both continuous and discrete problems. It does not suffer from the resolution limitation of the finite bit length encoding of parameters into GA alleles. It also has deep connections with both game theory and statistical physics. We review the PC approach using its motivation as the information theoretic formulation of bounded rationality for multi-agent systems. It is then compared with GA's on a diverse set of problems. To handle high dimensional surfaces, in the PC method investigated here p is restricted to a product distribution. Each distribution in that product is controlled by a separate agent. The test functions were selected for their difficulty using either traditional gradient descent or genetic algorithms. On those functions the PC-based approach significantly outperforms traditional GA's in both rate of descent, trapping in false minima, and long term optimization.

Huang, Chien-Feng↗

Intelligent Control for Drag Reduction on the X-48B Vehicle

This paper focuses on the development of an intelligent control technology for in-flight drag reduction. The system is integrated with and demonstrated on the full X-48B nonlinear simulation. The intelligent control system utilizes a peak-seeking control method implemented with a time-varying Kalman filter. Performance functional coordinate and magnitude measurements, or independent and dependent parameters respectively, are used by the Kalman filter to provide the system with gradient estimates of the designed performance function which is used to drive the system toward a local minimum in a steepestdescent approach. To ensure ease of integration and algorithm performance, a single-input single-output approach was chosen. The framework, specific implementation considerations, simulation results, and flight feasibility issues related to this platform are discussed.

Griffin, Brian Joseph↗

Heating of the solar corona by the resonant absorption of Alfven waves

An improved method for calculating the resonance absorption heating rate is discussed and the results are compared with observations in the solar corona. To accomplish this, the wave equation for a dissipative, compressible plasma is derived from the linearized magnetohydrodynamic equations for a plasma with transverse Alfven speed gradients. For parameters representative of the solar corona, it is found that a two-scale description of the wave motion is appropriate. The large-scale motion, which can be approximated as nearly ideal, has a scale which is on the order of the width of the loop. The small-scale wave, however, has a transverse scale much smaller than the width of the loop, with a width of about 0.3-250 km, and is highly dissipative. These two wave motions are coupled in a narrow resonance region in the loop where the global wave frequency equals the local Alfven wave frequency. Formally, this coupling comes about from using the method of matched asymptotic expansions to match the inner and outer (small and large scale) solutions. The resultant heating rate can be calculated from either of these solutions. A formula derived using the outer (ideal) solution is presented, and shown to be consistent with observations of heating and line broadening in the solar corona.

Davila, Joseph M.↗

An adaptive remeshing method for finite element thermal analysis

A finite element remeshing approach that makes use of quadrilateral and triangular elements is described. The approach uses the solution on a previous mesh to create a new mesh. Meshes are completely unstructured with highly refined elements in regions of steep gradients and larger elements where gradients are smaller. Studies of convergence rates for heat conduction problems with exact solutions show that for problems with highly localized solution variations, the remeshing approach gives smaller solution errors with fewer unknowns than refinement of uniform, structured meshes.

Thornton, Earl A.↗

A prediction method for velocity and temperature profiles in a two-dimensional nominally steady turbulent boundary layer

This paper describes a recently developed boundary-layer prediction method for a variable property compressible flow, in which heat transfer takes place primarily by forced convection and for which the mainstream Mach number is small. The leading order terms, in asymptotic expansions for large Reynolds numbers, are obtained for the mean velocity and temperature distribution in both the inner and outer layer of the turbulent boundary layer. Closure in the inner layer is achieved using an analytical model for the mean profiles which is based on the observed coherent structure of the time-dependent inner layer flow. For the outer layer, simple eddy viscosity and conductivity models are developed without recourse to the Reynolds analogy. In the prediction method a numerical solution of the outer layer equations is matched to the analytical inner layer profiles as the computation procedes downstream. Calculations are presented for a range of adverse and favorable pressure gradient flows and the predicted results compare well with existing data.

Weigand, G. G.↗

Calculation of skin-friction coefficients for low Reynolds number turbulent boundary layer flows

An analysis is presented of the reliability of various generally accepted empirical expressions for the prediction of the skin-friction coefficient C/sub f/ of turbulent boundary layers at low Reynolds numbers in zero-pressure-gradient flows on a smooth flat plate. The skin-friction coefficients predicted from these expressions were compared to the skin-friction coefficients of experimental profiles that were determined from a graphical method formulated from the law of the wall. These expressions are found to predict values that are consistently different than those obtained from the graphical method over the range 600 Re/sub theta 2000. A curve-fitted empirical relationship was developed from the present data and yields a better estimated value of C/sub f/ in this range. The data, covering the range 200 Re/sub theta 7000, provide insight into the nature of transitional flows. They show that fully developed turbulent boundary layers occur at Reynolds numbers Re/sub theta/ down to 425. Below this level there appears to be a well-ordered evolutionary process from the laminar to the turbulent profiles. These profiles clearly display the development of the turbulent core region and the shrinking of the laminar sublayer with increasing values of Re/sub theta/.

Barr, P. K.↗

Environmentally induced discharges on satellites

A discharge process whose trigger conditions are a negative exposed metallic surface surrounded by a less negative dielectric, and a large voltage gradient at a dielectric/metal interface is proposed. Analysis of SCATHA data for a discharge substantiates the postulation. Surface discharges cause a small transient charge transfer to space which results in voltage transients. A method of computing these transients, based on the charge lost through the capacitance to space and a fraction of charge stored in the dielectric at the discharge source was developed. It gives an estimate of the discharge transients at the discharge site, which is used as input for coupling code analysis of structure/system response. The transient computations were applied to a three-axis stabilized, geosynchronous satellite for both sunlight and eclipse charging. The energy of the transient pulses are about 1 mJ for sunlight discharge and 8 mJ for eclipse. Changing of selected coatings on the satellite would relieve the stress.

Stevens, N. J.↗

Direct simulation of a turbulent boundary layer up to R(sub)(theta)= 1410

The turbulent boundary layer on a flat plate, with zero pressure gradient, is simulated numerically at four stations between R sub theta = 225 and R sub theta = 1410. The three-dimensional time-dependent Navier-Stokes equations are solved using a spectra method with up to about 10 to the 7th power grid points. Periodic spanwise and stream-wise conditions are applied, and a multiple-scale procedure is applied to approximate the slow streamwise growth of the boundary layer. The flow is studied, primarily, from a statistical point of view. The solutions are compared with experimental results. The scaling of the mean and turbulent quantities with Reynolds number is examined and compared with accepted laws, and the significant deviations are documented. The turbulence at the highest Reynolds number is studied in detail. The spectra are compared with various theoretical models. Reynolds-stress budget data are provided for turbulence-model testing.

Spalart, Philippe R.↗

Optimal control concepts in design sensitivity analysis

A close link is established between open loop optimal control theory and optimal design by noting certain similarities in the gradient calculations. The resulting benefits include a unified approach, together with physical insights in design sensitivity analysis, and an efficient approach for simultaneous optimal control and design. Both matrix displacement and matrix force methods are considered, and results are presented for dynamic systems, structures, and elasticity problems.

Belegundu, Ashok D.↗