Engineering Papers⌕ Search

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.

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

At least 325 records · Page 18

Dependence of Ru2O3 Activity on Composition of Silicate Melts: Using Statistical Correlations to Infer Thermodynamic Behavior in the Melt

Understanding variations in activity with composition is an essential step in improving prediction of partition coefficients during magma evolution. Variations in activity with composition are complex and do not generally exhibit ideal behavior relative to a traditional melt-component set. Although deviations from component ideality can be modeled numerically by simply fitting to compositional variables (such as in a regular or subregular solution model), such models have not been particularly successful for describing variations in trace component activities. A better approach might be to try to identify components that do a better job of describing the behavior of the species in the melt. Electrochemical Measurement of Ru2O3 activities: Electrodes were inserted into silicate melt beads of various compositions (Table 1) suspended on Ptwire loops in a 1-atm gas mixing furnace. An electrical potential was imposed between the electrodes, the imposed potential increasing along a step ramp with a pulse imposed on each step (Fig. 1). Current flows between electrodes when electroactive species in the melt are oxidized or reduced at the electrodes. The resulting current was measured at the top and bottom of the voltage pulse, and the difference (the differential current) was plotted against potential. The peak of the resulting curve is related to the activity coefficient for the particular electroactive species (Ru2O3) in the melt [1, 2, 3]. A significant part of the nonideal contribution to activity is due not to intrinsic properties of the component in the melt, but to our ignorance about the state and mixing properties of the component in the melt.

Colson, R. O.↗

Parametric Geometry, Structured Grid Generation, and Initial Design Study for REST-Class Hypersonic Inlets

Computational design and analysis of three-dimensional hypersonic inlets with shape transition has been a significant challenge due to the complex geometry and grid required for three-dimensional viscous flow calculations. Currently, the design process utilizes an inviscid design tool to produce initial inlet shapes by streamline tracing through an axisymmetric compression field. However, the shape is defined by a large number of points rather than a continuous surface and lacks important features such as blunt leading edges. Therefore, a design system has been developed to parametrically construct true CAD geometry and link the topology of a structured grid to the geometry. The Adaptive Modeling Language (AML) constitutes the underlying framework that is used to build the geometry and grid topology. Parameterization of the CAD geometry allows the inlet shapes produced by the inviscid design tool to be generated, but also allows a great deal of flexibility to modify the shape to account for three-dimensional viscous effects. By linking the grid topology to the parametric geometry, the GridPro grid generation software can be used efficiently to produce a smooth hexahedral multiblock grid. To demonstrate the new capability, a matrix of inlets were designed by varying four geometry parameters in the inviscid design tool. The goals of the initial design study were to explore inviscid design tool geometry variations with a three-dimensional analysis approach, demonstrate a solution rate which would enable the use of high-fidelity viscous three-dimensional CFD in future design efforts, process the results for important performance parameters, and perform a sample optimization.

Ferlemann, Paul G.↗

Soil Moisture Limitations on Monitoring Boreal Forest Regrowth Using Spaceborne L-Band SAR Data

A study was carried out to investigate the utility of L-band SAR data for estimating aboveground biomass in sites with low levels of vegetation regrowth. Data to estimate biomass were collected from 59 sites located in fire-disturbed black spruce forests in interior Alaska. PALSAR L-band data (HH and HV polarizations) collected on two dates in the summer/fall of 2007 and one date in the summer of 2009 were used. Significant linear correlations were found between the log of aboveground biomass (range of 0.02 to 22.2 t ha-1) and (L-HH) and (L-HV) for the data collected on each of the three dates, with the highest correlation found using the LHV data collected when soil moisture was highest. Soil moisture, however, did change the correlations between L-band and aboveground biomass, and the analyses suggest that the influence of soil moisture is biomass dependent. The results indicate that to use L-band SAR data for mapping aboveground biomass and monitoring forest regrowth will require development of approaches to account for the influence that variations in soil moisture have on L-band microwave backscatter, which can be particularly strong when low levels of aboveground biomass occur

Kasischke, Eric S.↗

Variable-Speed Power-Turbine Research at Glenn Research Center

The main rotors of the NASA Large Civil Tilt-Rotor (LCTR) notional vehicle operate over a wide speed-range, from 100 percent at takeoff to 54 percent at cruise. The variable-speed power turbine (VSPT) offers one approach by which to effect this speed variation. VSPT aerodynamics challenges include high work factors at cruise, wide (40 to 60 ) incidence-angle variations in blade and vane rows over the speed range, and operation at low Reynolds numbers. Rotordynamics challenges include potential responsiveness to shaft modes within the 50 percent VSPT speed-range. A research effort underway at NASA Glenn Research Center, intended to address these key aerodynamic and rotordynamic challenges, is described. Conceptual design and 3-D multistage RANS and URANS analyses, conducted internally and under contract, provide expected VSPT sizing, stage-count, performance and operability information, and maps for system studies. Initial steps toward experimental testing of incidence-tolerant blading in a transonic linear cascade are described, and progress toward development/improvement of a simulation capability for multistage turbines with low Reynolds number transitional flow is summarized. Preliminary rotordynamics analyses indicate that viable concept engines with 50 percent VSPT shaft-speed range. Assessments of potential paths toward VSPT component-level testing are summarized.

Welch, Gerard E.↗

100 kHz High-Spectral-Resolution NO-PLIF Measurements for Compressible Flows

In the present work, we use a burst-mode laser and optical parametric oscillator system to perform high-spectral resolution NO-PLIF measurements of an underexpanded jet at a repetition-rate of 100 kHz, with the motivation of multi-parameter measurements of temperature, pressure, and velocity. The laser frequency of the 1064 nm seed laser for the burst-mode laser was scanned during the burst to cover two neighboring absorption line pairs near 226 nm. The peak PLIF signal intensity varies along the axial (z) direction of the underexpanded jet as the laser frequency is scanned, which we attribute to the collisional shift induced by the flow and only revealed due to the narrow linewidth of the laser. A pseudo-Voight fit is applied to the LIF excitation spectra on a pixel-by-pixel basis to measure the spectral position of the peak intensity for the two transition pairs and their amplitude. The spectral position of the peak intensity is used to derive a frequency shift, which is separated into its collisional and Doppler components using the axisymmetric nature of the flow field. The amplitude of the measured peaks is used for two-line rotational thermometry. Challenges for quantitative measurements using such an approach are discussed, including measuring the spatial variations in the energy distribution of the laser sheet at a 100 kHz repetition rate and uncertainty/variability in the step size during the fast frequency scan.

laser induced fluorescence↗

100 kHz High-Spectral-Resolution NO-PLIF Measurements for Compressible Flows

In the present work, we use a burst-mode laser and optical parametric oscillator system to perform high-spectral resolution NO-PLIF measurements of an underexpanded jet at a repetition-rate of 100 kHz, with the motivation of multi-parameter measurements of temperature, pressure, and velocity. The laser frequency of the 1064 nm seed laser for the burst-mode laser was scanned during the burst to cover two neighboring absorption line pairs near 226 nm. The peak PLIF signal intensity varies along the axial (z) direction of the underexpanded jet as the laser frequency is scanned, which we attribute to the collisional shift induced by the flow and only revealed due to the narrow linewidth of the laser. A pseudo-Voight fit is applied to the LIF excitation spectra on a pixel-by-pixel basis to measure the spectral position of the peak intensity for the two transition pairs and their amplitude. The spectral position of the peak intensity is used to derive a frequency shift, which is separated into its collisional and Doppler components using the axisymmetric nature of the flow field. The amplitude of the measured peaks is used for two-line rotational thermometry. Challenges for quantitative measurements using such an approach are discussed, including measuring the spatial variations in the energy distribution of the laser sheet at a 100 kHz repetition rate and uncertainty/variability in the step size during the fast frequency scan.

laser induced fluorescence↗

Reference Mission Version 3.0 Addendum to the Human Exploration of Mars: The Reference Mission of the NASA Mars Exploration Study Team

This Addendum to the Mars Reference Mission was developed as a companion document to NASA Special Publication 6107, “Human Exploration of Mars: The Reference Mission of the NASA Exploration Study Team.” The Addendum summarizes changes and updates to the Mars Reference Mission the Exploration Office developed since the final draft of SP 6107 was printed in early 1998. The Reference Mission is a tool used by the Exploration Team and the exploration community to compare and evaluate approaches to mission and system concepts that could be used for human exploration missions. It is intended to identify and clarify system “drivers,” or significant sources of cost, performance, risk, and schedule variation. It does not represent a final or recommended approach to human Mars missions. Several alternative scenarios, including human exploration missions to the Moon, asteroids, or other targets beyond Earth orbit as well as employing different technical approaches to solving mission and technology challenges, are currently under study by the Exploration Team. Comparing alternative approaches provides the basis for continual improvement to technology investment plans and general understanding of future human exploration missions. The Addendum represents a “snapshot” of work in progress in support of planning for future human exploration missions through May 1998. Annual publications of revisions to the Reference Mission are planned beginning in late 1998.

Mars↗

From minerals to rocks: Toward modeling lithologies with remote sensing

High spectral resolution imaging spectroscopy will play an important role in future planetary missions. Sophisticated approaches will be needed to unravel subtle, super-imposed spectral features typically of natural systems, and to maximize the science return of these instruments. Carefully controlled laboratory investigations using homogeneous mineral separates have demonstrated that variations due to solid solution, changes in modal abundances, and the effects of particle size are well understood from a physical basis. In many cases, these variations can be modeled quantitatively using photometric models, mixing approaches, and deconvolution procedures. However, relative to the spectra of individual mineral components, reflectance spectra of rocks and natural surfaces exhibit a reduced spectral contrast. In addition, soils or regolith, which are likely to dominate any natural planetary surface, exhibit spectral properties that have some similarities to the parent materials, but due to weathering and alteration, differences remain that cannot yet be fully recreated in the laboratory or through mixture modeling. A significant challenge is therefore to integrate modeling approaches to derive both lithologic determinations and include the effects of alteration. We are currently conducting laboratory investigations in lithologic modeling to expand upon the basic results of previous analyses with our initial goal to more closely match physical state of natural systems. The effects of alteration are to be considered separately.

Mustard, John F.↗

Neptune's cloud structure in 1989 - Photometric variations and correlation with ground-based images

Ground-based photoelectric photometry in b, y, and the 6190 and 7250 A methane-bands, as well as spectrum scans of the methane 6190 A band and CCD images at 6190 and 8900 A, were obtained for Neptune during Voyager 2's approach of that planet on August 24, 1989. Photometric variations are presently correlated with the disk transit of bright planetary features, and the changes in feature distribution and brightness noted in the results are evaluated for implications bearing on long-term variability. It is suggested that the long-term secular variation is related to a slow change in a size of location of both the bright companion and the Great Dark Spot.

Lockwood, G. W.↗

Assessment of Low Cycle Fatigue Behavior of Powder Metallurgy Alloy U720

The fatigue lives of modem powder metallurgy disk alloys are influenced by variabilities in alloy microstructure and mechanical properties. These properties can vary as functions of variables the different steps of materials/component processing: powder atomization, consolidation, extrusion, forging, heat treating, and machining. It is important to understand the relationship between the statistical variations in life and these variables, as well as the change in life distribution due to changes in fatigue loading conditions. The objective of this study was to investigate these relationships in a nickel-base disk superalloy, U720, produced using powder metallurgy processing. Multiple strain-controlled fatigue tests were performed at 538 C (1000 F) at limited sets of test conditions. Analyses were performed to: (1) assess variations of microstructure, mechanical properties, and LCF failure initiation sites as functions of disk processing and loading conditions; and (2) compare mean and minimum fatigue life predictions using different approaches for modeling the data from assorted test conditions. Significant variations in life were observed as functions of the disk processing variables evaluated. However, the lives of all specimens could still be combined and modeled together. The failure initiation sites for tests performed at a strain ratio R(sub epsilon) = epsilon(sub min)/epsilon(sub max) of 0 were different from those in tests at a strain ratio of -1. An approach could still be applied to account for the differences in mean and maximum stresses and strains. This allowed the data in tests of various conditions to be combined for more robust statistical estimates of mean and minimum lives.

Gabb, Tomothy P.↗

SABER (TIMED) and MLS (UARS) Temperature Observations of Mesospheric and Stratospheric QBO and Related Tidal Variations

More than three years of temperature observations from the SABER (TIMED) and MLS WARS) instruments are analyzed to study the annual and inter-annual variations extending from the stratosphere into the upper mesosphere. The SABER measurements provide data from a wide altitude range (15 to 95 km) for the years 2002 to 2004, while the MLS data were taken in the 16 to 55 km altitude range a decade earlier. Because of the sampling properties of SABER and MLS, the variations with local solar time must be accounted for when estimating the zonal mean variations. An algorithm is thus applied that delineates with Fourier analysis the year-long variations of the migrating tides and zonal mean component. The amplitude of the diurnal tide near the equator shows a strong semiannual periodicity with maxima near equinox, which vary from year to year to indicate the influence from the Quasi-biennial Oscillation (QBO) in the zonal circulation. The zonal mean QBO temperature variations are analyzed over a range of latitudes and altitudes, and the results are presented for latitudes from 48"s to 48"N. New results are obtained for the QBO, especially in the upper stratosphere and mesosphere, and at mid-latitudes. At Equatorial latitudes, the QBO amplitudes show local peaks, albeit small, that occur at different altitudes. From about 20 to 40 km, and within about 15" of the Equator, the amplitudes can approach 3S K for the stratospheric QBO or SQBO. For the mesospheric QBO or MQBO, we find peaks near 70 km, with temperature amplitudes reaching 3.5"K, and near 85 km, the amplitudes approach 2.5OK. Morphologically, the amplitude and phase variations derived from the SABER and MLS measurements are in qualitative agreement. The QBO amplitudes tend to peak at the Equator but then increase again pole-ward of about 15" to 20'. The phase progression with altitude varies more gradually at the Equator than at mid-latitudes. A comparison of the observations with results from the Numerical Spectral Model (NSM) reveals that there is qualitative agreement. The NSM generates the QBO extending from the stratosphere into the upper mesosphere, with temperature variations extending to mid latitudes, but the predicted amplitudes are smaller than those observed.

Huang, Frank T.↗

Spectral Ratio Imaging with Hyperion Satellite Data for Geological Mapping

Since the advent of LANDSAT I in 1972, many different multispectral satellites have been orbited by the U.S. and other countries. These satellites have varied from 4 spectral bands in LANDSAT I to 14 spectral bands in the ASTER sensor aboard the TERRA space platform. Hyperion is a relatively new hyperspectral sensor with over 220 spectral bands. The huge increase in the number of spectral bands offers a substantial challenge to computers and analysts alike when it comes to the task of mapping features on the basis of chemical composition, especially if little or no ground truth is available beforehand from the area being mapped. One approach is the theoretical approach of the modeler, where all extraneous information (atmospheric attenuation, sensor electronic gain and offset, etc.) is subtracted off and divided out, and laboratory (or field) spectra of materials are used as training sets to map features in the scene of similar composition. This approach is very difficult to keep accurate because of variations in the atmosphere, solar illumination, and sensor electronic gain and offset that are not always perfectly recorded or accounted for. For instance, to apply laboratory or field spectra of materials as data sets from the theoretical approach, the header information of the files must reflect the correct, up-to-date sensor electronic gain and offset and the analyst must pick the exact atmospheric model that is appropriate for the day of data collection in order for classification procedures to accurately match pixels in the scene with the laboratory or field spectrum of a desired target on the basis of the hyperspectral data. The modeling process is so complex that it is difficult to tell when it is operating well or determine how to fix it when it is incorrect. Recently RSI has announced that the latest version of their ENVI software package is not performing atmospheric corrections correctly with the FLAASH atmospheric model. It took a long time to determine that it was wrong, and may take an equally long time (or longer) to fix.

Vincent, Robert K.↗

Document analysis with neural net circuits

Document analysis is one of the main applications of machine vision today and offers great opportunities for neural net circuits. Despite more and more data processing with computers, the number of paper documents is still increasing rapidly. A fast translation of data from paper into electronic format is needed almost everywhere, and when done manually, this is a time consuming process. Markets range from small scanners for personal use to high-volume document analysis systems, such as address readers for the postal service or check processing systems for banks. A major concern with present systems is the accuracy of the automatic interpretation. Today's algorithms fail miserably when noise is present, when print quality is poor, or when the layout is complex. A common approach to circumvent these problems is to restrict the variations of the documents handled by a system. In our laboratory, we had the best luck with circuits implementing basic functions, such as convolutions, that can be used in many different algorithms. To illustrate the flexibility of this approach, three applications of the NET32K circuit are described in this short viewgraph presentation: locating address blocks, cleaning document images by removing noise, and locating areas of interest in personal checks to improve image compression. Several of the ideas realized in this circuit that were inspired by neural nets, such as analog computation with a low resolution, resulted in a chip that is well suited for real-world document analysis applications and that compares favorably with alternative, 'conventional' circuits.

Graf, Hans Peter↗

Novel Estimation of Pilot Performance Characteristics

Two mechanisms internal to the pilot that affect performance during a tracking task are: 1) Pilot equalization (i.e. lead/lag); and 2) Pilot gain (i.e. sensitivity to the error signal). For some applications McRuer's Crossover Model can be used to anticipate what equalization will be employed to control a vehicle's dynamics. McRuer also established approximate time delays associated with different types of equalization - the more cognitive processing that is required due to equalization difficulty, the larger the time delay. However, the Crossover Model does not predict what the pilot gain will be. A nonlinear pilot control technique, observed and coined by the authors as 'amplitude clipping', is shown to improve stability, performance, and reduce workload when employed with vehicle dynamics that require high lead compensation by the pilot. Combining linear and nonlinear methods a novel approach is used to measure the pilot control parameters when amplitude clipping is present, allowing precise measurement in real time of key pilot control parameters. Based on the results of an experiment which was designed to probe workload primary drivers, a method is developed that estimates pilot spare capacity from readily observable measures and is tested for generality using multi-axis flight data. This paper documents the initial steps to developing a novel, simple objective metric for assessing pilot workload and its variation over time across a wide variety of tasks. Additionally, it offers a tangible, easily implementable methodology for anticipating a pilot's operating parameters and workload, and an effective design tool. The model shows promise in being able to precisely predict the actual pilot settings and workload, and observed tolerance of pilot parameter variation over the course of operation. Finally, an approach is proposed for generating Cooper-Harper ratings based on the workload and parameter estimation methodology.

Human pilot models↗

Kinematic Modeling of Separation Compression for Paired Approaches to Closely-Spaced Parallel Runways

In a simultaneous paired approach to closely-spaced parallel runways, a pair of aircraft flies in close proximity on parallel approach paths. The longitudinal separation between the aircraft must be maintained within a range that avoids wake encounters and, if one of the aircraft blunders, avoids collision. To increase operational availability, the approach procedure must accommodate a mixture of aircraft sizes and, consequently, approach speeds. In these procedures, the slower aircraft is placed in the lead position. The faster aircraft maintains separation from the slow aircraft in a dependent operation until final approach and flies independently afterward. Due to the higher approach speed of the fast aircraft, longitudinal separation will decrease during final approach. Therefore, the fast aircraft must position itself before the final approach so that it will remain within the safe range of separation as separation decreases. Given the approach geometry and speed schedule for each aircraft, one can use kinematics to estimate the separation loss between a pair of aircraft. A kinematic model can complement fast-time Monte-Carlo simulations of the approach by enabling a tailored reduction in the variation of starting position for the fast aircraft. One could also implement the kinematic model in ground-based or on-board decision support tools to compute the optimal initial separation for a given pair of aircraft. To better match the auto-coupled flight of real aircraft, the paper derives a kinematic model where the speed schedule is flown using equivalent airspeed. The predicted time of flight using the equivalent airspeed kinematic model compares well against a high-fidelity aircraft simulation performing the same approach. This model also demonstrates a modest increase in the predicted loss of separation when contrasted against a kinematic model that assumes the scheduled speed is true airspeed.

Madden, Michael M.↗

Efficient Gradient-Based Shape Optimization Methodology Using Inviscid/Viscous CFD

The formerly developed preconditioned-biconjugate-gradient (PBCG) solvers for the analysis and the sensitivity equations had resulted in very large error reductions per iteration; quadratic convergence was achieved whenever the solution entered the domain of attraction to the root. Its memory requirement was also lower as compared to a direct inversion solver. However, this memory requirement was high enough to preclude the realistic, high grid-density design of a practical 3D geometry. This limitation served as the impetus to the first-year activity (March 9, 1995 to March 8, 1996). Therefore, the major activity for this period was the development of the low-memory methodology for the discrete-sensitivity-based shape optimization. This was accomplished by solving all the resulting sets of equations using an alternating-direction-implicit (ADI) approach. The results indicated that shape optimization problems which required large numbers of grid points could be resolved with a gradient-based approach. Therefore, to better utilize the computational resources, it was recommended that a number of coarse grid cases, using the PBCG method, should initially be conducted to better define the optimization problem and the design space, and obtain an improved initial shape. Subsequently, a fine grid shape optimization, which necessitates using the ADI method, should be conducted to accurately obtain the final optimized shape. The other activity during this period was the interaction with the members of the Aerodynamic and Aeroacoustic Methods Branch of Langley Research Center during one stage of their investigation to develop an adjoint-variable sensitivity method using the viscous flow equations. This method had algorithmic similarities to the variational sensitivity methods and the control-theory approach. However, unlike the prior studies, it was considered for the three-dimensional, viscous flow equations. The major accomplishment in the second period of this project (March 9, 1996 to March 8, 1997) was the extension of the shape optimization methodology for the Thin-Layer Navier-Stokes equations. Both the Euler-based and the TLNS-based analyses compared with the analyses obtained using the CFL3D code. The sensitivities, again from both levels of the flow equations, also compared very well with the finite-differenced sensitivities. A fairly large set of shape optimization cases were conducted to study a number of issues previously not well understood. The testbed for these cases was the shaping of an arrow wing in Mach 2.4 flow. All the final shapes, obtained either from a coarse-grid-based or a fine-grid-based optimization, using either a Euler-based or a TLNS-based analysis, were all re-analyzed using a fine-grid, TLNS solution for their function evaluations. This allowed for a more fair comparison of their relative merits. From the aerodynamic performance standpoint, the fine-grid TLNS-based optimization produced the best shape, and the fine-grid Euler-based optimization produced the lowest cruise efficiency.

Baysal, Oktay↗

Higher-Order Finite Elements for Computing Thermal Radiation

Two variants of the finite-element method have been developed for use in computational simulations of radiative transfers of heat among diffuse gray surfaces. Both variants involve the use of higher-order finite elements, across which temperatures and radiative quantities are assumed to vary according to certain approximations. In this and other applications, higher-order finite elements are used to increase (relative to classical finite elements, which are assumed to be isothermal) the accuracies of final numerical results without having to refine computational meshes excessively and thereby incur excessive computation times. One of the variants is termed the radiation sub-element (RSE) method, which, itself, is subject to a number of variations. This is the simplest and most straightforward approach to representation of spatially variable surface radiation. Any computer code that, heretofore, could model surface-to-surface radiation can incorporate the RSE method without major modifications. In the basic form of the RSE method, each finite element selected for use in computing radiative heat transfer is considered to be a parent element and is divided into sub-elements for the purpose of solving the surface-to-surface radiation-exchange problem. The sub-elements are then treated as classical finite elements; that is, they are assumed to be isothermal, and their view factors and absorbed heat fluxes are calculated accordingly. The heat fluxes absorbed by the sub-elements are then transferred back to the parent element to obtain a radiative heat flux that varies spatially across the parent element. Variants of the RSE method involve the use of polynomials to interpolate and/or extrapolate to approximate spatial variations of physical quantities. The other variant of the finite-element method is termed the integration method (IM). Unlike in the RSE methods, the parent finite elements are not subdivided into smaller elements, and neither isothermality nor other unrealistic physical conditions are assumed. Instead, the equations of radiative heat transfer are integrated numerically over the parent finite elements by use of a computationally efficient Gaussian integration scheme.

Gould, Dana C.↗

Roman-WFI: characterization of widefield chromatic response

Characterization of filter transmission is essential to control systematic chromatic errors in photometric surveys. Widefield surveys employ large filter optics viewed from high angles of incidence. Large optics may have variations in edge wavelength definition across their surface, and their spectral response shifts at high incidence angles. I will describe the characterization of the band definitions for the Wide Field Instrument of the Roman Space Telescope as part of its thermal vacuum test campaigns and through preparatory work. The measurement approach and analysis enable joint retrieval of chromatic response variations across the filter and focal plane. Additional tests probe out-of-band response and thermal straylight paths.

Eric Switzer↗