Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “observational 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,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Nimbus 7 CZCS - Reduction of its radiometric sensitivity with time

Preliminary results are described for an effort to quantify the sensitivity decay of a radiometry sensor (the Coastal Zone Color Scanner or CZCS aboard Nimbus 7). The method used in the study is to (1) compute the water-leaving radiance for imagery acquired in regions where this radiance is known or can be independently estimated, and (2) adjust the sensor calibration to force agreement between the two radiances. Decay factors for orbit numbers from 0 to 20,000 are plotted, and surface and space measurements are compared for the Gulf Stream and the Northern Sargasso Sea at different seasons. The fact that a seasonal variability in the chlorophyll a concentration in the Sargasso Sea was found in the sensor analysis (apparently the first such satellite observation) increases confidence in the method.

Gordon, H. R.↗
NASA NTRS1993

Modeling and data assimilation for Mars Observer

The Mars Observer mission will return a substantial set of meteorological data for the atmosphere of Mars, principally in the form of thermal infrared soundings from the Pressure Modulator Infrared Radiometer (PMIRR). The dataset will enable global studies of a range of dynamical and physical phenomena in the martian atmosphere. The solar orbital configuration of Mars Observer, however, means that the coverage of the martian atmosphere and surface by remote sounding instruments will be asynchronous. This leads to difficulties of interpretation when studying synoptic-scale phenomena with timescales of a few days. In an attempt to overcome such problems we propose to use data assimilation techniques, as currently employed for operational weather forecasting on the Earth, in conjunction with a Martian General Circulation Model (MGCM) that is under development at Oxford and Reading Universities. An overview of the model and assimilation scheme will be given, with reference to problems encountered so far in applying such methods to remote-sensing observations of the martian atmosphere.

Lewis, Stephen R.↗
DOE OSTI2026

Data for Clumping Index Estimation With 30°-tilted Cameras in Row Crops: Evaluation of Methods and Segment Size Effects

The clumping index (CI) quantifies the spatial distribution of foliage elements and is essential for accurately estimating the plant area index (PAI), canopy radiative transfer, and photosynthesis. Traditionally, the finite-length averaging method (LX), the gap size distribution method (CC), and a combined approach of CC and LX (CLX) have been applied to instruments like TRAC and digital hemispherical photography to estimate CI. However, a comprehensive evaluation of these methods in row crops remains limited, especially regarding the influence of segment size on CI. Meanwhile, digital cameras offer a cost-effective and user-friendly solution for canopy measurements in row crops, yet their application in this context remains underexplored. In this study, we employed a new approach using a 30°-tilted digital camera to estimate CI in corn and soybean fields, applying the LX, CC, and CLX methods. We systematically assessed the performance of these three methods by combining field measurements in real-world fields with simulations using the LESS 3D radiative transfer model. Our results showed that CLX applied to the whole image and 45° segment offered accurate estimation of CI (bias within ±0.1, RMSE < 0.2) and PAI (bias within ±0.4, RMSE < 1) in real-world fields and LESS simulations. The accuracy of the LX method was highly sensitive to segment size, with the best performance observed at the 15° segment (PAI bias within ±0.4). In contrast, the CC method remained stable across different segment sizes, and its performance was generally comparable to that of LX, except at the 15° segment. Across view zenith angles, CI derived from CC generally showed a continuous increase, while those from LX and CLX followed a rising trend at small zenith angles but began to decline at 68°, likely due to an increasing proportion of no-gap segments. Seasonally, LX tended to show decreasing CI during early growth stages but increased as the canopy matured, whereas CC and CLX showed gradually increasing CI before plateauing at peak PAI. The 30°-tilted camera effectively captured CI variations across different angles and growth stages, making it a practical and robust instrument for row crop canopy structure analysis. Applying these CI methods to digital cameras offers a low-cost and accessible CI estimation alternative, improving canopy structure monitoring accuracy in row crops.

Modeling↗
NASA NTRS1993

Maximum a posteriori classification of multifrequency, multilook, synthetic aperture radar intensity data

We present a maximum a posteriori (MAP) classifier for classifying multifrequency, multilook, single polarization SAR intensity data into regions or ensembles of pixels of homogeneous and similar radar backscatter characteristics. A model for the prior joint distribution of the multifrequency SAR intensity data is combined with a Markov random field for representing the interactions between region labels to obtain an expression for the posterior distribution of the region labels given the multifrequency SAR observations. The maximization of the posterior distribution yields Bayes's optimum region labeling or classification of the SAR data or its MAP estimate. The performance of the MAP classifier is evaluated by using computer-simulated multilook SAR intensity data as a function of the parameters in the classification process. Multilook SAR intensity data are shown to yield higher classification accuracies than one-look SAR complex amplitude data. The MAP classifier is extended to the case in which the radar backscatter from the remotely sensed surface varies within the SAR image because of incidence angle effects. The results obtained illustrate the practicality of the method for combining SAR intensity observations acquired at two different frequencies and for improving classification accuracy of SAR data.

Rignot, E.↗
NASA NTRS1994

On the determination of the number of O stars in H II regions and starburst galaxies

The hot star population in H II regions, H II galaxies, and starburst galaxies is often described in terms of the number of 'equivalent' O stars of a single representative subtype and luminosity class needed to produce the ionizing luminosity deduced from the nebular recombination lines in the optical spectra. In this paper we define conversion factors eta(sub 0), eta(sub 1), and zeta(sub 5000) with which the total number of O V stars and their flux contribution at 5000 A can be derived from the number of these 'equivalent' stars. These quantities depend primarily on three parameters: the slope and upper mass limit of the stellar mass function and the metallicity of the region. Using the latest stellar atmosphere and evolution models, we calculate eta(sub 0), eta(sub 1), and zeta(sub 5000) for a large number of values of these parameters. The results are presented in tabular as well as graphical form. We apply our results to two H II regions for which the hot star population are known and find that the predicted numbers of O stars agree well the observed counts. In addition, we describe a method by which the values of eta(sub 0) and eta(sub 1) and the observed emission-line fluxes can be used to place constraints on the allowed values of the slope and upper mass limit of the stellar mass function in a region.

Vacca, William D.↗
NASA NTRS1998

Microstructural Development of Directionally Solidified Hg(1-x)Zn(x)Se Alloys

Hg(1-x)Zn(x)Se alloys have been studied as an alternative to HgCdTe for the detection of electromagnetic radiation, because the shorter ZnSe and HgSe bonds were predicted to improve lattice stability. Previous studies showed that optical properties were stable with time and that the microhardness was higher than for HgTe based alloys. However, for this material to be commercially viable, compositional variations must be reduced, and the distribution of defects induced by the growth process must be understood. The basic structural properties of bulk materials were characterized, and dislocation etch pit densities were measured to determine whether the addition of Zn to the HgSe lattice reduced the potential for dislocation formation. Hg(0.9)Zn(0.1)Se alloys were directionally solidified using a modified Bridgman-Stockbarger method and in an applied magnetic field. Radial compositional variations were greatly reduced when solidification occurred in an applied magnetic field. The combination of a convex liquid-solid interface shape and a method for reduced wetting produced boules that were a single crystal after reaching full diameter. Observed surface features indicated ampoule wetting could be eliminated using a graphite getter. Microstructural characteristics were greatly improved over HgCdTe alloys, In six boules, a total of only one twin was observed, and no inclusions were found. A method for polishing and producing dislocation etch pits was developed for these alloys, revealing dislocation etch pit densities one to two orders of magnitude less than HgTe based alloys. A kink in the thermal profile during processing of one boule generated more dislocations than did lattice mismatch due to compositional variations. The results of this investigation indicate this alloy has improved microstructural properties and resistance to dislocation formation compared with similar II-VI alloys, and should be further investigated for the detection of electromagnetic radiation.

Cobb, S. D.↗
NASA NTRS2020

Reducing Errors in Velocity–Azimuth Display (VAD) Wind and Deformation Retrievals from Airborne Doppler Radars in Convective Environments

The present study describes methods to reduce the uncertainty of velocity–azimuth display (VAD) wind and deformation retrievals from downward-pointing, conically scanning, airborne Doppler radars. These retrievals have important applications in data assimilation and real-time data processing. Several error sources for VAD retrievals are considered here, including violations to the underlying wind field assumptions, Doppler velocity noise, data gaps, temporal variability, and the spatial weighting function of the VAD retrieval. Specific to airborne VAD retrievals, we also consider errors produced due to the radar scans occurring while the instrument platform is in motion. While VAD retrievals are typically performed using data from a single antenna revolution, other strategies for selecting data can be used to reduce retrieval errors. Four such data selection strategies for airborne VAD retrievals are evaluated here with respect to their effects on the errors. These methods are evaluated using the second hurricane nature run numerical simulation, analytic wind fields, and observed Doppler radar radial velocities. The proposed methods are shown to reduce the median absolute error of the VAD wind retrievals, especially in the vicinity of deep convection embedded in stratiform precipitation. The median absolute error due to wind field assumption violations for the along-track and for the across-track wind is reduced from 0.36 to 0.08 m s−1 and from 0.35 to 0.24 m s−1, respectively. Although the study focuses on Doppler radars, the results are equally applicable to conically scanning Doppler lidars as well.

Charles N Helms↗
NASA NTRS1984

Evaluation of glass model precipitation forecasts for North America during SOP-1

A study has been made of the accuracy of the GLAS fourth order forecast model's precipitation forecasts over North America during the first Special Observing Period (SOP-1) of FGGE and of the impact of the FGGE special observing systems on these forecasts. Fourteen 120 h predictions were generated using the coarse 4 deg latitude by 5 deg longitude version of the GLAS model from both the FGGE and NOSAT assimilation cycles. These forecasts were then verified against a detailed set of precipitation observations. Separate verifications were performed for precipitation accumulations or = .01 in., or = .1 in., and or = 1 in. occuring in 6, 12 and 24 time periods. In addition, three different methods of verification were applied. In the first method, model precipitation forecasts are verified against the greatest precipitation amount observed within a 4 deg latitude by 5 deg longitude gridbox centered on each gridpoint. For the prediction at a gridpoint to be verified only one observation within the gridbox is required. In the second method, a weighted mean of the closest observations to the gridpoint is used, provided that at least two observations on opposite sides of a gridpoint are available. For the third method, the model precipitation forecasts were interpolated to the observation locations.

Atlas, R.↗
DOE OSTI2022

Understanding Recent Global Hydroclimate Change Using Multivariate Detection and Attribution Techniques and GCM Experiments

The observed slowdown in the rate of global average surface temperature (T) warming during the post-1998 decade (referred to at the time as the “global warming hiatus”) posed a challenge for climate researchers because it was not captured in global climate models projection. Observational and modeling evidence hinted that that leveling of the T trend was linked with the simultaneous appearance of cold surface temperature in the eastern equatorial Pacific that formed a protracted La Niña-like state, a phenomenon known as the Inter-decadal Pacific Oscillation (IPO). It was however not clear whether the concurrent equatorial Pacific cooling was a regional response to the rise in atmospheric greenhouse gas levels or anthropogenic aerosols that is not captured by the climate models or whether it was an unrelated, serendipitous manifestation of a natural phenomenon, internal to the climate system. To resolve this issue, we proposed in this project to develop and apply fingerprinting methods to differentiate between the observed climate response to various external forcing agents and the manifestations of internal climate variability. The project also sought to find out if climate models are able to reproduce the change signals detected in key observational phenomena and if they can be used to investigate the underlying processes. As such we examined changes in global hydroclimate variability and in an outstanding example of regional hydroclimate change, namely the multi-decadal 20 th century change in Sahel rainfall.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1995

Fast Multipole Methods for Three-Dimensional N-body Problems

We are developing computational tools for the simulations of three-dimensional flows past bodies undergoing arbitrary motions. High resolution viscous vortex methods have been developed that allow for extended simulations of two-dimensional configurations such as vortex generators. Our objective is to extend this methodology to three dimensions and develop a robust computational scheme for the simulation of such flows. A fundamental issue in the use of vortex methods is the ability of employing efficiently large numbers of computational elements to resolve the large range of scales that exist in complex flows. The traditional cost of the method scales as Omicron (N(sup 2)) as the N computational elements/particles induce velocities at each other, making the method unacceptable for simulations involving more than a few tens of thousands of particles. In the last decade fast methods have been developed that have operation counts of Omicron (N log N) or Omicron (N) (referred to as BH and GR respectively) depending on the details of the algorithm. These methods are based on the observation that the effect of a cluster of particles at a certain distance may be approximated by a finite series expansion. In order to exploit this observation we need to decompose the element population spatially into clusters of particles and build a hierarchy of clusters (a tree data structure) - smaller neighboring clusters combine to form a cluster of the next size up in the hierarchy and so on. This hierarchy of clusters allows one to determine efficiently when the approximation is valid. This algorithm is an N-body solver that appears in many fields of engineering and science. Some examples of its diverse use are in astrophysics, molecular dynamics, micro-magnetics, boundary element simulations of electromagnetic problems, and computer animation. More recently these N-body solvers have been implemented and applied in simulations involving vortex methods. Koumoutsakos and Leonard (1995) implemented the GR scheme in two dimensions for vector computer architectures allowing for simulations of bluff body flows using millions of particles. Winckelmans presented three-dimensional, viscous simulations of interacting vortex rings, using vortons and an implementation of a BH scheme for parallel computer architectures. Bhatt presented a vortex filament method to perform inviscid vortex ring interactions, with an alternative implementation of a BH scheme for a Connection Machine parallel computer architecture.

Koumoutsakos, P.↗
NASA NTRS1994

Calculation of clear-sky outgoing longwave radiation using ECMWF gridded fields and ISCCP C1 cloud data

Clear-Sky OLR (outgoing longwave radiation) has a variety of definitions. Clear-sky OLR can be measured from satellites by examining the data to detect clouds, eliminating those observations that are contaminated with clouds, and averaging the remaining data. The disadvantage of this method is that satellites do not provide time continuous coverage, requiring many observations and careful analysis to remove sampling biases. Alternatively, clear-sky OLR can be found continuously using a radiative flux code, combined with analyzed fields of temperature, humidity and clouds, either from assimilated observations or a GCM simulation. The disadvantage is that the cloud field is typically very low resolution, on the order of 200 km, so they lack important small scale cloud detail. Even if fractional clouds data is available, it is not obvious how GCM and ERBE results should be compared. The purpose of this study is to assess clear-sky OLR calculated from these methods using available data. Cloud data from ISCCP provides the fractional cloud coverage data necessary to test these methods against calculations based on observed data. Clear-sky OLR is calculated using the ECMWF/TOGA archive. Monthly averages are made in the manner of the several methods listed above and compared to each other and to ERBE.

Ramsey, Perry G.↗
DOE OSTI2020

A Dynamic Trajectory Fit to Multisensor Fireball Observations

Meteorites with known orbital origins are key to our understanding of solar system formation and the source of life on Earth. Fireball networks have been developed globally in a unified effort to record and ultimately retrieve these cosmic samples. However, the accuracy of the determined orbit and the likelihood of meteorite recovery depend directly on the accuracy of the chosen meteoroid triangulation method. There are three leading techniques for meteoroid triangulation discussed in the literature: the method of planes, the straight-line least-squares method, and the multiparameter fit method. Here we describe an alternative method to meteoroid triangulation, called the dynamic trajectory fit. This approach uses the meteoroid’s 3D dynamic equations of motion to fit a realistic trajectory directly to multisensor line-of-sight observations. This method has the ability to resolve fragmentation events, fit systematic observatory timing offsets, and determine mass estimates of the meteoroid along its observable trajectory. Through a comprehensive Monte Carlo analysis of over 100,000 trajectory simulations, we find this new method to more accurately estimate meteoroid trajectories of slow entry events (<25 km s{sup −1}) and events observed from low convergence angles (<10°) compared to existing meteoroid triangulation techniques. Additionally, we triangulate an observed fireball event with visible fragmentation using the various triangulation methods to show that the proposed dynamic trajectory fit implementing fragmentation to best match the captured multisensor line-of-sight data.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

Physical and chemical structure of the IC 63 nebula. 1: Millimeter and far-infrared observations

We present results of a (sub)millimeter and far-infrared study of the reflection/emission nebula IC 63, located close to the BO.5p star gamma Cas. The source has been mapped in the (12)CO 2 - 1 and 3 - 2, (13)CO 2 - 1, and CS 2 - 1 lines and shows a small molecular cloud less than 1'x 2' in extent, which coincides with the brightest optical nebulosity and IRAS 100 micrometer emission. IC 63 is therefore an excellent example of a nearby (d approximately = 230 pc), edge-on photon-dominated region (PDR). Various other molecules have been observed at the peak position through their rotational transitions, in order to probe the physical parameters and to derive abundances. The measured CO, HCO(+) HCN, CS and H2CO line ratios suggest that the cloud is warm, T approximately = 50 K, and dense, n (H2) approximately = 5 x 10(exp 4)/cu cm. Excitation of molecules by electrons may play a significant role in this PDR. On the basis of these physical conditions, column densities have been determined from the observed line strengths. Several different methods are discussed to constrain the H2 column density, including the use of measured submillimeter continuum fluxes. The resulting abundances of species such as CN and CS are similar to those found in cold, dark clouds like TMC-1 and L134N. However, the abundances of other simple molecules such as HNC, HCO(+) and possibly C2H are lower by factors of at least three, probably because of the enhanced photodissociation rates at a distance of 1.3 pc from a B star. Surprisingly, only the abundance of the H2S molecule appears enhanced. More complex, volatile molecules such as CH3OH CH3CN and HNCO, and the sulfur-oxides SO and SO2 have not been found in this cloud. Limited observations of molecules in the reflection nebulea NGC 2023 are presented as well, and the resulting molecular abundances are compared with those found for IC 63.

Jansen, David J.↗
DOE OSTI2021

Measuring Gas Transport and Sorption in Large Intact Geologic Specimens via the Piezometric Method

Uptake of noble gases into heterogeneous geologic core samples was measured using a piezometric methodology. In addition to measuring accessible porosity—as with gas pycnometry—by monitoring the rate of pressure decay, this method can also be used to estimate the gas effective diffusivity in the sample. In contrast to previous applications of this method, where milligram quantities of fractured grains are characterized, here approximately kilogram core samples were left intact when tested. In doing so, a more representative sample of the heterogeneous field geology is provided. Additionally, alteration of the pore structure and connectivity during sample preparation is avoided. Furthermore, to scale the piezometric method from milligrams to kilograms, the system was designed to operate at medium vacuum (1 to 100 Pa) to restrict transport in pores less than approximately 60 µm to large Knudsen numbers. To test the system performance, two samples of interest were selected: a rhyolitic welded tuff from Blue Canyon Dome at the Energetic Materials Research and Testing Center and a zeolitized non-welded rhyolitic tuff from the Nevada National Security Site. Three noble gases were utilized in this test series; Argon and xenon as they are of direct interest to nuclear monitoring efforts and helium as it is a weakly adsorbing reference standard. Additionally, mercury intrusion porosimetry measurements were made on subsamples of the core to compare the observed porosity by the two methods and to discuss gas transport rates in the context of the measured pore distribution.

58 GEOSCIENCES↗
DOE OSTI2020

Matrix Completion Using Alternating Minimization for Distribution System State Estimation

This paper examines the problem of state estimation in power distribution systems under low-observability conditions. The recently proposed constrained matrix completion method which combines the standard matrix completion method and power flow constraints has been shown to be effective in estimating voltage phasors under low-observability conditions using single-snapshot information. However, the method requires solving a semidefinite programming (SDP) problem, which becomes computationally infeasible for large systems and if multiple-snapshot (time-series) information is used. This paper proposes an efficient algorithm to solve the constrained matrix completion problem with time-series data. This algorithm is based on reformulating the matrix completion problem as a bilinear (non-convex) optimization problem, and applying the alternating minimization algorithm to solve this problem. This paper proves the summable convergence of the proposed algorithm, and demonstrates its efficacy and scalability via IEEE 123-bus system and a real utility feeder system. This paper also explores the value of adding more data from the history in terms of computation time and estimation accuracy.

41 EE - Solar Energy Technologies Office (EE-4S)↗
NASA NTRS1984

Observations from varying the lift and drag inputs to a noise prediction method for supersonic helical tip speed propellers

Previous comparisons between calculated and measured supersonic helical tip speed propeller noise show them to have different trends of peak blade passing tone versus helical tip Mach number. It was postulated that improvements in this comparison could be made first by including the drag force terms in the prediction and then by reducing the blade lift terms at the tip to allow the drag forces to dominate the noise prediction. Propeller hub to tip lift distributions were varied, but they did not yield sufficient change in the predicted lift noise to improve the comparison. This result indicates that some basic changes in the theory may be needed. In addition, the noise predicted by the drag forces did not exhibit the same curve shape as the measured data. So even if the drag force terms were to dominate, the trends with helical tip Mach number for theory and experiment would still not be the same. The effect of the blade shock wave pressure rise was approxmated by increasing the drag coefficient at the blade tip. Predictions using this shock wdave approximation did have a curve shape similar to the measured data. This result indicates that the shock pressure rise probably controls the noise at supersonic tip speed and that the linear prediction method can give the proper noise trend with Mach number.

Dittmar, J. H.↗