SEARCH · Engineering Papers

Results for “Variational approach”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 919 records · Page 51

All sources NASA NTRS DOE OSTI

NASA NTRS2013

Uncertainty Determination for Aeroheating in Uranus and Saturn Probe Entries by the Monte Carlo Method

The 2013-2022 Decaedal survey for planetary exploration has identified probe missions to Uranus and Saturn as high priorities. This work endeavors to examine the uncertainty for determining aeroheating in such entry environments. Representative entry trajectories are constructed using the TRAJ software. Flowfields at selected points on the trajectories are then computed using the Data Parallel Line Relaxation (DPLR) Computational Fluid Dynamics Code. A Monte Carlo study is performed on the DPLR input parameters to determine the uncertainty in the predicted aeroheating, and correlation coefficients are examined to identify which input parameters show the most influence on the uncertainty. A review of the present best practices for input parameters (e.g. transport coefficient and vibrational relaxation time) is also conducted. It is found that the 2(sigma) - uncertainty for heating on Uranus entry is no more than 2.1%, assuming an equilibrium catalytic wall, with the uncertainty being determined primarily by diffusion and H(sub 2) recombination rate within the boundary layer. However, if the wall is assumed to be partially or non-catalytic, this uncertainty may increase to as large as 18%. The catalytic wall model can contribute over 3x change in heat flux and a 20% variation in film coefficient. Therefore, coupled material response/fluid dynamic models are recommended for this problem. It was also found that much of this variability is artificially suppressed when a constant Schmidt number approach is implemented. Because the boundary layer is reacting, it is necessary to employ self-consistent effective binary diffusion to obtain a correct thermal transport solution. For Saturn entries, the 2(sigma) - uncertainty for convective heating was less than 3.7%. The major uncertainty driver was dependent on shock temperature/velocity, changing from boundary layer thermal conductivity to diffusivity and then to shock layer ionization rate as velocity increases. While radiative heating for Uranus entry was negligible, the nominal solution for Saturn computed up to 20% radiative heating at the highest velocity examined. The radiative heating followed a non-normal distribution, with up to a 3x variation in magnitude. This uncertainty is driven by the H(sub 2) dissociation rate, as H(sub 2) that persists in the hot non-equilibrium zone contributes significantly to radiation.

Palmer, Grant↗

NASA NTRS1987

Photometry of six radar target asteroids

Photoelectric photometry of six earth-approaching asteroids is presented. The selection criterion was that they were close enough in 1986 to be observed by radar. Rotation periods were obtained for 1986 DA, 3199, 3103, and 1983 RD. 1986 JK and 1986 RA showed no detectable brightness variations during the monitoring time on several nights each, and therefore were either seen pole-on or have long rotation periods. Asteroids 1986 JK and 1986 RA are of taxonomic class C, 1986 DA and 3103 of class X, 1983 RD of class Q, and only 3199 of the class S that was previously believed to be predominant among earth-approaching asteroids.

Wisniewski, W. Z.↗

NASA NTRS1990

Distributed memory approaches for robotic neural controllers

The suitability is explored of two varieties of distributed memory neutral networks as trainable controllers for a simulated robotics task. The task requires that two cameras observe an arbitrary target point in space. Coordinates of the target on the camera image planes are passed to a neural controller which must learn to solve the inverse kinematics of a manipulator with one revolute and two prismatic joints. Two new network designs are evaluated. The first, radial basis sparse distributed memory (RBSDM), approximates functional mappings as sums of multivariate gaussians centered around previously learned patterns. The second network types involved variations of Adaptive Vector Quantizers or Self Organizing Maps. In these networks, random N dimensional points are given local connectivities. They are then exposed to training patterns and readjust their locations based on a nearest neighbor rule. Both approaches are tested based on their ability to interpolate manipulator joint coordinates for simulated arm movement while simultaneously performing stereo fusion of the camera data. Comparisons are made with classical k-nearest neighbor pattern recognition techniques.

Jorgensen, Charles C.↗

NASA NTRS1973

A numerical procedure for the parametric optimization of three dimensional scramjet nozzles

A numerical procedure permitting the rapid determination of the internal performance of a class of scramjet nozzle configurations is presented. The approach developed is based on the construction of quasi two dimensional simple wave networks, where lateral expansion effects are incorporated by one dimensional approximations. A numerical procedure following this approach has has been developed and results obtained are highly comparable to those obtained employing a characteristic procedure. The numerical program developed permits the parametric variation of cowl length, turning angles on the cowl and vehicle undersurface and lateral expansion and is subject to fixed constraints such as the vehicle length and nozzle exit height. The program requires uniform initial conditions at the burner exit station and yields the location of all predominant wave zones, accounting for lateral expansion effects. In addition, the program yields the detailed pressure distribution on the cowl and vehicle undersurface and calculates the nozzle thrust, lift and pitching moment.

Dash, S.↗

NASA NTRS1980

Dynamics and control of a continuum model for a solar power system

An approach for modeling dynamic equations of motion of a plate attached with rigid bodies is presented. The equations of motion are developed using the principle of virtual work. Lagrange multipliers are used as interaction forces and/or moments to maintain prescribed constraints which is the basis of the interconnection between the plate and rigid bodies. The overall approach is unique in the sense that a continuous model described by a family of partial differential equations is established. An approximate formulation by using variational method is established yielding a solution compatible with the assumed degree of approximation. The formulation is useful particularly when parametric study of dynamic response for a satellite power system is desired. As an example, an approximate governing equation of algebraic eigenvalue problem is given for a dual microwave power transmission system. Controller design is discussed.

Juang, J. N.↗

NASA NTRS1981

Canopy reflectance as influenced by solar illumination angle

The interaction of the solar illumination angle and row azimuth angle of the measured reflectance factor (RF) of soybean canopies was investigated. Diurnal changes of nearly 140% were observed in the red wavelength region when canopies covered 64% of the soil. The amount of shadow observed was a function of the plant geometry and row width. As soil cover approached 100%, the diurnal changes diminished. A function tht describes the solar illumination angle with respect to the row azimuth explained most of the diurnal variation in the measured RF. Variation in near infrared response was much less and did not appear to be as strongly related to Sun-row angle interactions. The near infrared/red ratio was highly sensitive to Sun angle-row direction interactions, whereas the greenness function, utilizing all four spectral bands, was not.

Bauer, M. E.↗

NASA NTRS1985

Cyanobacterial mats: Microanalysis of community metabolism

The microbial communities in two sites were studied using several approaches: (1) light microscopy; (2) the measurement of microprofiles of oxygen and sulfide at the surface of the microbial mat; (3) the study of diurnal variation of oxygen and sulfides; (4) in situ measurement of photosynthesis and sulfate reduction and study of the coupling of these two processes; (5) measurement of glutathione in the upper layers of the microbial mat as a possible oxygen quencher; (6) measurement of reduced iron as a possible intermediate electron donor along the established redoxcline in the mats; (7) measurement of dissolved phosphate as an indicator of processes of break down of organic matter in these systems; and (8) measurement of carbon dioxide in the interstitial water and its delta C-13 in an attempt to understand the flow of CO2 through the systems. Microbial processes of primary production and initial degradation at the most active zone of the microbial mat were analyzed.

Cohen, Y.↗

NASA NTRS1985

Application of TVD schemes for the Euler equations of gas dynamics

Highly accurate and yet stable shock-capturing finite difference schemes have been designed for the computation of the Euler equations of gas dynamics. Four different principles for the construction of high resolution total variation diminishing (TVD) schemes are available, including hybrid schemes, a second-order extension of Godunov's scheme by van Leer (1979), the modified flux approach of Harten (1983, 1984), and the numerical fluctuation approach of Roe (1985). The present paper has the objective to review the class of second-order TVD schemes via the modified flux approach. Attention is given to first-order TVD schemes, a second-order accurate explicit TVD scheme, the global order of accuracy of the second-order TVD scheme, extensions to systems and two-dimensional conservation laws, numerical experiments with a second-order explicit TVD scheme, implicit TVD schemes, and second-order implicit TVD schemes.

Yee, H. C.↗

NASA NTRS1986

The influence of free-stream turbulence on separation of turbulent boundary layers in incompressible, two-dimensional flow

Experiments were conducted to determine if free-stream turbulence scale affects separation of turbulent boundary layers. In consideration of possible interrelation between scale and intensity of turbulence, the latter characteristic also was varied and its role was evaluated. Flow over a 2-dimensional airfoil in a subsonic wind tunnel was studied with the aid of hot-wire anemometry, liquid-film flow visualization, a Preston tube, and static pressure measurements. Profiles of velocity, relative turbulence intensity, and integral scale in the boundary layer were measured. Detachment boundary was determined for various angles of attack and free-stream turbulence. The free-stream turbulence intensity and scale were found to spread into the entire turbulent boundary layer, but the effect decreased as the airfoil surface was approached. When the changes in stream turbulence were such that the boundary layer velocity profiles were unchanged, detachment location was not significantly affected by the variations of intensity and scale. Pressure distribution remained the key factor in determining detachment location.

Potter, J. Leith↗

NASA NTRS1992

A unified approach for numerical simulation of viscous compressible and incompressible flows over adiabatic and isothermal walls

A new formulation (including the choice of variables, their non-dimensionalization, and the form of the artificial viscosity) is proposed for the numerical solution of the full Navier-Stokes equations for compressible and incompressible flows with heat transfer. With the present approach, the same code can be used for constant as well as variable density flows. The changes of the density due to pressure and temperature variations are identified and it is shown that the low Mach number approximation is a special case. At zero Mach number, the density changes due to the temperature variation are accounted for, mainly through a body force term in the momentum equation. It is also shown that the Boussinesq approximation of the buoyancy effects in an incompressible flow is a special case. To demonstrate the new capability, three examples are tested. Flows in driven cavities with adiabatic and isothermal walls are simulated with the same code as well as incompressible and supersonic flows over a wall with and without a groove. Finally, viscous flow simulations of an oblique shock reflection from a flat plate are shown to be in good agreement with the solutions available in literature.

Hafez, M.↗

NASA NTRS1956

Theoretical and Experimental Investigation of the Subsonic-Flow Fields Beneath Swept and Unswept Wings with Tables of Vortex-Induced Velocities

The flow-field characteristics beneath swept and unswept wings as determined by potential-flow theory are compared with the experimentally determined flow fields beneath swept and unswept wing-fuselage combinations. The potential-flow theory utilized considered both spanwise and chordwise distributions of vorticity as well as the wing-thickness effects. The perturbation velocities induced by a unit horseshoe vortex are included in tabular form. The results indicated that significant chordwise flow gradients existed beneath both swept and unswept wings at zero lift and throughout the lift range. The theoretical predictions of the flow-field characteristics were qualitatively correct in all cases considered, although there were indications that the magnitudes of the downwash angles tended to be overpredicted as the tip of the swept wing was approached and that the sidewash angles ahead of the unswept wing were underpredicted. The calculated effects of compressibility indicated that significant increases in the chordwise variation of flow angles and dynamic-pressure ratios should be expected in going from low to high subsonic speeds.

Alford, William J., Jr.↗

NASA NTRS1957

Theoretical and Experimental Investigation of the Subsonic-Flow Fields Beneath Swept and Unswept Wings with Tables or Vortex-induced Velocities

The flow-field characteristics beneath swept and unswept wings as determined by potential-flow theory are compared with the experimentally determined flow fields beneath swept and unswept wing-fuselage combinations. The potential-flow theory utilized considered both spanwise and chordwise distributions of vorticity as well as the wing-thickness effects. The perturbation velocities induced by a unit horseshoe vortex are included in tabular form. The theoretical predictions of the flow-field characteristics were qualitatively correct in all cases considered, although there were indications that the magnitudes of the downwash angles tended to be overpredicted as the tip of the swept wing was approached and that the sidewash angles ahead of the unswept wing were underpredicted. The calculated effects of compressibility indicated that significant increases in the chordwise variation of flow angles and dynamic-pressure ratios should be expected in going from low to high subsonic speeds.

Alford, William J , Jr↗

NASA NTRS1995

Observations of the CO bulge on Venus and implications for mesospheric winds

Observations of CO at 2.6 mm (115.27 GHz) were made with the Owens Valley Radio Observatory (OVRO) millimeter interferometer in 1986 and 1988, yielding high-quality disk-resolved spectra which were inverted to determine the CO mixing ratio profile from distinct regions on the disk, allowing us to map the distribution of CO in the upper mesosphere of Venus both horizontally and vertically. The 1986 observations were of the morning terminator and were particularly useful in searching for a suspected CO maximum ('bulge') on the nightside. The resulting CO mixing ratio profiles were mapped for various altitudes as functions of latitude and local time, and we report that we have resolved the previously inferred CO bulge. The bulge increases in magnitude from a small day-night variation at 90 km to an extensive nightside peak at 100 km, the upper limit of our observations. The peak bulge-to-dayside ratio approached 20-30 at 100 km in 1986 and may have been as large as 50-100, assuming late-afternoon CO abundances found in 1988 were similar to those in 1986. Three-dimensional mapping shows that in the upper mesosphere the bulge was displaced from local midnight toward the morning equator, centered at 3:30 AM local time. Using the qualitative model of mesospheric circulation on Venus proposed by Clancy and Muhleman (1985, 1991), we explain this shift in terms of strong retrograde zonal winds throughout the mesosphere, matching the directly detected mesospheric circulation (Shah, K., D. O. Muhleman, and G. L. Berge 1991) observed with the same dataset in 1988.

Gurwell, Mark A.↗

NASA NTRS1997

Observed Changes in Atmospheric Boundary Layer Properties at Memphis International Airport During August 1995

As part of the NASA Terminal Area Productivity Program, Langley Research Center embarked on a series of field measurements of wake vortex characteristics and associated atmospheric boundary layer properties. One measurement period was at the Memphis International Airport in August 1995. Atmospheric temperature, humidity, winds, turbulence, radiation, and soil properties were measured from a variety of sensor systems and platforms including sodars, profilers, aircraft and towers. This research focused on: (1) changes that occurred in tower data during sunrise and sunset transitions, (2) vertical variation of temperature and cross-head winds at selected times utilizing combinations of sensors, and (3) changes measured by an OV-10 aircraft during approaches and level flights. Significant but not unusual changes are documented and discussed in terms of expected boundary layer behavior. Questions on measurement and prediction of these changes from existing and near-term capabilities are discussed in the context of a future Aircraft Vortex Spacing System.

Zak, J. Allen↗

NASA NTRS1995

Geometric adjustments to account for eye eccentricity in processing horizontal and vertical eye and head movement data

Neglecting the eccentric position of the eyes in the head can lead to erroneous interpretation of ocular motor data, particularly for near targets. We discuss the geometric effects that eye eccentricity has on the processing of target-directed eye and head movement data, and we highlight two approaches to processing and interpreting such data. The first approach involves determining the true position of the target with respect to the location of the eyes in space for evaluating the efficacy of gaze, and it allows calculation of retinal error directly from measured eye, head, and target data. The second approach effectively eliminates eye eccentricity effects by adjusting measured eye movement data to yield equivalent responses relative to a specified reference location (such as the center of head rotation). This latter technique can be used to standardize measured eye movement signals, enabling waveforms collected under different experimental conditions to be directly compared, both with the measured target signals and with each other. Mathematical relationships describing these approaches are presented for horizontal and vertical rotations, for both tangential and circumferential display screens, and efforts are made to describe the sensitivity of parameter variations on the calculated results.

NASA Discipline Neuroscience↗

NASA NTRS2002

Aeroacoustic Flow Phenomena Accurately Captured by New Computational Fluid Dynamics Method

One of the challenges in the computational fluid dynamics area is the accurate calculation of aeroacoustic phenomena, especially in the presence of shock waves. One such phenomenon is "transonic resonance," where an unsteady shock wave at the throat of a convergent-divergent nozzle results in the emission of acoustic tones. The space-time Conservation-Element and Solution-Element (CE/SE) method developed at the NASA Glenn Research Center can faithfully capture the shock waves, their unsteady motion, and the generated acoustic tones. The CE/SE method is a revolutionary new approach to the numerical modeling of physical phenomena where features with steep gradients (e.g., shock waves, phase transition, etc.) must coexist with those having weaker variations. The CE/SE method does not require the complex interpolation procedures (that allow for the possibility of a shock between grid cells) used by many other methods to transfer information between grid cells. These interpolation procedures can add too much numerical dissipation to the solution process. Thus, while shocks are resolved, weaker waves, such as acoustic waves, are washed out.

Blech, Richard A.↗

NASA NTRS2003

Optimal Feedback Control of Thermal Networks

An improved approach to the mathematical modeling of feedback control of thermal networks has been devised. Heretofore software for feedback control of thermal networks has been developed by time-consuming trial-and-error methods that depend on engineers expertise. In contrast, the present approach is a systematic means of developing algorithms for feedback control that is optimal in the sense that it combines performance with low cost of implementation. An additional advantage of the present approach is that a thermal engineer need not be expert in control theory. Thermal networks are lumped-parameter approximations used to represent complex thermal systems. Thermal networks are closely related to electrical networks commonly represented by lumped-parameter circuit diagrams. Like such electrical circuits, thermal networks are mathematically modeled by systems of differential-algebraic equations (DAEs) that is, ordinary differential equations subject to a set of algebraic constraints. In the present approach, emphasis is placed on applications in which thermal networks are subject to constant disturbances and, therefore, integral control action is necessary to obtain steady-state responses. The mathematical development of the present approach begins with the derivation of optimal integral-control laws via minimization of an appropriate cost functional that involves augmented state vectors. Subsequently, classical variational arguments provide optimality conditions in the form of the Hamiltonian equations for the standard linear-quadratic-regulator (LQR) problem. These equations are reduced to an algebraic Riccati equation (ARE) with respect to the augmented state vector. The solution of the ARE leads to the direct computation of the optimal proportional- and integral-feedback control gains. In cases of very complex networks, large numbers of state variables make it difficult to implement optimal controllers in the manner described in the preceding paragraph.

Papalexandris, Miltiadis↗

NASA NTRS2017

Statistical Characterization of 18650-Format Lithium-Ion Cell Thermal Runaway Energy Distributions

Effective thermal management systems, designed to handle the impacts of thermal runaway (TR) and to prevent cell-to-cell propagation, are key to safe operation of lithium-ion (Li-ion) battery assemblies. Critical factors for optimizing these systems include the total energy released during a single cell TR event and the fraction of the total energy that is released through the cell casing vs. through the ejecta material. A unique calorimeter was utilized to examine the TR behavior of a statistically significant number of 18650-format Li-ion cells with varying manufacturers, chemistries, and capacities. The calorimeter was designed to contain the TR energy in a format conducive to discerning the fractions of energy released through the cell casing vs. through the ejecta material. Other benefits of this calorimeter included the ability to rapidly test of large quantities of cells and the intentional minimization of secondary combustion effects. High energy (270 Wh/kg) and moderate energy (200 Wh/kg) 18650 cells were tested. Some of the cells had an imbedded short circuit (ISC) device installed to aid in the examination of TR mechanisms under more realistic conditions. Other variations included cells with bottom vent (BV) features and cells with thin casings (0.22 1/4m). After combining the data gathered with the calorimeter, a statistical approach was used to examine the probability of certain TR behavior, and the associated energy distributions, as a function of capacity, venting features, cell casing thickness and temperature.

Thermal Runaway Energy Distribution↗