Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sufficient statistics”

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 379 records · Page 21

Pressure Effects on the Relaxation of an Excited Ethane Molecule in High-Pressure Bath Gases

Here, we use molecular dynamics to calculate the rotational and vibrational energy relaxation of C 2 H 6 in Ar, Kr, and Xe bath gases over a pressure range of 10 to 400 atm and at temperatures of 300 K and 800 K. The C 2 H 6 is instantaneously excited by 80 kcal/mol randomly distributed into both vibrational and rotational modes. The computed relaxation rates show little sensitivity to the identity of the noble gas in the bath. Vibrational relaxation rates show a non-linear pressure dependence at 300 K. At 800 K the reduced range of bath gas densities covered by the range of pressures do not yet show any non-linearity in the pressure dependence. Rotational relaxation is characterized with two relaxation rates. The slower rate is comparable to the vibrational relaxation rate. The faster rate has a linear pressure dependence at 300 K but an irregular, nonlinear pressure dependence at 800 K. To understand this, a model was developed based on approximating the periodic box used in the molecular dynamics simulations by an equal-volume collection of cubes where each cube is sized to allow only single occupancy by the noble gas or the molecule. Combinatorial statistics then leads to a pressure and temperature dependent analytic distribution of the bath gas species the molecule encounters in a collision. This distribution, the dissociation energy of molecule/bath gas complexes and bath gas clusters, and the computed energy release per collision combine to show that only at 300 K is the energy release sufficient to dissociate likely complexes and clusters. This suggests that persistent and pressure-dependent clusters and complexes at 800 K may be responsible for the non-linear pressure dependence of rotational relaxation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The GLAS physical inversion method for analysis of TIROS N data

The Goddard Laboratory for Atmospheric Sciences approach to temperature sounding from satellite observations is fundamentally different from that of NESS in that heavy reliance is placed on the ability to model accurately the instrumental response to atmospheric and surface conditions, while no use is made of statistical relationships between satellite observations and atmospheric temperature profiles. The method involves starting with a guess set of atmospheric and surface conditions from which expected brightness temperatures for the satellite observations are computed. Then, iterative relaxation of atmospheric and surface conditions is performed according to the difference between observed and computed brightness temperatures until sufficient agreement is reached. In the absence of sufficient agreement, no retrieval is produced for that location.

Susskind, J.↗

Impact of Defect Distributions on the Reliability of Structural Materials

The purpose of this study is to develop a probabilistic model to estimate the fatigue life of gas turbine disks made of powder superalloy metals. Such metals are known to have high strength and ductility under extremely high temperatures and known to have high fatigue lives. However, despite the high quality controls instituted during metal processing, sufficiently high numbers of ceramic inclusions materialize within the metal matrix creating weak links that lead to reduced fatigue lives. To study the effect of such ceramic inclusions on the fatigue life of gas turbine disks, an experimental program has been initiated at the NASA Glenn Research. In this phase of the study, a statistical analysis of the test results is undertaken which served as the basis of the probabilistic simulation model. The study is a part of the Crack Resistant Disk Materials task within the Ultra Safe Propulsion Project.

Ghosn, Michel J.↗

Some statistical problems inherent in measuring precipitation

A formula permitting calculation of the mean-square error of the mean value of a random variable due to periodic sampling is derived and applied to estimating the sampling error for satellite observation of the mean rainfall during the GARP Atlantic Tropical Experiment (GATE). The effects of both spatial resolution and frequency of observation on the sampling error are summarized in graphs. It is found that four observations per day are sufficient to determine the monthly mean rainfall over an area of 2.5 deg square (280km x 280km) to within a standard deviation of 5 percent of its mean value; two samples per day would yield an error with a standard deviation slightly less than 10 percent of the mean. A satellite instrument with less frequent sampling may produce significantly greater error in the estimate of monthly mean rainfall.

Flueck, J. A.↗

Filtered Rayleigh-Ritz is all you need

Recent work has shown that the (block) Lanczos algorithm can be used to extract approximate energy spectra and matrix elements from (matrices of) correlation functions in quantum field theory, and identified exact coincidences between Lanczos analysis methods and others. In this work, we note another coincidence: the Lanczos algorithm is equivalent to the well-known Rayleigh-Ritz method applied to Krylov subspaces. Rayleigh-Ritz provides optimal eigenvalue approximations within subspaces; we find that spurious-state filtering allows these optimality guarantees to be retained in the presence of statistical noise. We explore the relation between Lanczos and Prony's method, their block generalizations, generalized pencil of functions (GPOF), and methods based on the generalized eigenvalue problem (GEVP), and find they all fall into a larger "Prony-Ritz equivalence class", identified as all methods which solve a finite-dimensional spectrum exactly given sufficient correlation function (matrix) data. This equivalence allows simpler and more numerically stable implementations of (block) Lanczos analyses.

97 MATHEMATICS AND COMPUTING↗

Detection of Digital Elevation Model Errors Using X-band Weather Radar

Flight in Instrument Meteorological Conditions requires pilots to manipulate flight controls while referring to a Primary Flight Display. The Primary Flight Display indicates aircraft attitude along with, in some cases, many other state variables such as altitude, speed, and guidance cues. Synthetic Vision Systems have been proposed that overlay the traditional information provided on Primary Flight Displays onto a scene depicting the location of terrain and other geo-spatial features.Terrain models used by these displays must have sufficient quality to avoid providing misleading information. This paper describes how X-band radar measurements can be used as part of a monitor, and/or maintenance system, to quantify the integrity of terrain models that are used by systems such as Synthetic Vision. Terrain shadowing effects, as seen by the radar, are compared in a statistical manner against estimated shadow feature elements extracted from the stored terrain model from the perspective of the airborne observer. A test statistic is defined that enables detection of errors as small as the range resolution of the radar. Experimental results obtained from two aircraft platforms hosting certified commercial-off-the-shelf X-band radars test the premise and illustrate its potential.

Young, Steven D.↗

The Mira–Titan Universe – IV. High-precision power spectrum emulation

Modern cosmological surveys are delivering data sets characterized by unprecedented quality and statistical completeness; this trend is expected to continue in the future as new ground- and space-based surveys come online. In order to maximally extract cosmological information from these observations, matching theoretical predictions are needed. At low redshifts, the surveys probe the non-linear regime of structure formation where cosmological simulations are the primary means of obtaining the required information. The computational cost of sufficiently resolved large-volume simulations makes it prohibitive to run very large ensembles. Nevertheless, precision emulators built on a tractable number of high-quality simulations can be used to build very fast prediction schemes to enable a variety of cosmological inference studies. We have recently introduced the Mira–Titan Universe simulation suite designed to construct emulators for a range of cosmological probes. This gravity-only set of simulations covers the standard six cosmological parameters {ω m , ω b , σ 8 , $h, n_s, w_0$} and, in addition, includes massive neutrinos and a dynamical dark energy equation of state {ω ν , $w_a$}. In this paper, we present the final emulator for the matter power spectrum based on 111 cosmological simulations, each covering a (2.1 Gpc) 3 volume and evolving 3200 3 particles. In this work, an additional set of 1776 lower resolution simulations and TimeRG perturbation theory results for the power spectrum are used to cover scales straddling the linear to mildly non-linear regimes (maximum wavenumber $\textit{k}$ = 5 Mpc –1 ). The emulator provides predictions at the 2–3 percent level of accuracy over a wide range of cosmological parameters and is publicly released as part of this paper.

79 ASTRONOMY AND ASTROPHYSICS↗

Using Temporal Deep Learning Models to Estimate Daily Snow Water Equivalent Over the Rocky Mountains

Abstract In this study we construct and compare three different deep learning (DL) models for estimating daily snow water equivalent (SWE) from high‐resolution gridded meteorological fields over the Rocky Mountain region. To train the DL models, Snow Telemetry (SNOTEL) station‐based SWE observations are used as the prediction target. All DL models produce higher median Nash‐Sutcliffe Efficiency (NSE) values than a conceptual SWE model and interpolated gridded data sets, although mean squared errors also tend to be higher. Sensitivity of the SWE prediction to the model's input variables is analyzed using an explainable artificial intelligence (XAI) method, yielding insight into the physical relationships learned by the models. This method reveals the dominant role precipitation and temperature play in snowpack dynamics. In applying our models to estimate SWE throughout the Rocky Mountains, an extrapolation problem arises since the statistical properties of SWE (e.g., annual maximum) and geographical properties of individual grid points (e.g., elevation) differ from the training data. This problem is solved by normalizing the SWE with its historical maximum value to alleviate extrapolation for all tested DL models. Our work shows that the DL models are promising tools for estimating SWE, and sufficiently capture relevant physical relationships to make them useful for spatial and temporal extrapolation of SWE values.

54 ENVIRONMENTAL SCIENCES↗

Multidimensional analysis of data obtained in experiments with X-ray emulsion chambers and extensive air showers

Nonparametric statistical methods are used to carry out the quantitative comparison of the model and the experimental data. The same methods enable one to select the events initiated by the heavy nuclei and to calculate the portion of the corresponding events. For this purpose it is necessary to have the data on artificial events describing the experiment sufficiently well established. At present, the model with the small scaling violation in the fragmentation region is the closest to the experiments. Therefore, the treatment of gamma families obtained in the Pamir' experiment is being carried out at present with the application of these models.

Chilingaryan, A. A.↗

The depths of the mare basalt source region

Chemical data for the parental magmas of the nine known VLT Array I pyroclastic glasses show statistically significant trends due to olivine (approximately Fo83) control during their formation by partial melting. The compositional scatter is largely due to compositional variations in the source regions on the + or - 1 percent level. This compositional scatter is small when one considers that the scale of the source region is up to 1000 km, but is sufficient to make positive identification of the residual phase(s) in the source regions difficult. Nevertheless, when the effects of the scatter are properly modeled, it is relatively clear that olivine is the residual phase in the source region. Hence these data and additional constraints indicate that the source regions are at shallow depths in the moon.

Binder, A. B.↗

Theoretical analysis of isotope effects on ozone formation in oxygen photochemistry

In situ measurements of stratospheric ozone and laboratory studies of ozone production in electric discharge through oxygen have shown previously that ozone containing heavy isotopes of oxygen (O-17, O-18) may be formed preferentially. In order to assess the relevance of thee latter experiment to the stratospheric measurements, detailed understanding of the effect of isotopic substitution on the O3 formation reaction O + O2 + M yields O3 + M and on the O atom exchange reaction O + O2 + O yields O2 + O is necessary. In this work, an estimate of the effect of isotopic substitution on the recombination rate is made by us of approximate dynamical theories and statistical mechanics. The results indicate the possibility of isotope effects on the O + O2 recombination rate of the order of several percent at stratospheric temperatures. In general, recombination reactions involving heavy (mass 49, 50) O3 formation are found to be slower than the reaction leading to normal (mass 48) O3 formation. The calculated isotope effects are sufficiently small that the uncertainties in the model input and the approximations in the dynamical theories will probably make the quantitative nature of these results subject to considerable uncertainty. This isotope effect should not be observable in the atmosphere given the precision of the current measurements but may be crucial in the understanding of the laboratory experiments, where observed enhancements are only of the order of several percent. Possible reasons for this discrepancy between the observed enhancement and predicted depletion are presented.

Kaye, J. A.↗

Visualizing turbulent mixing of gases and particles

A physical model and interactive computer graphics techniques have been developed for the visualization of the basic physical process of stochastic dispersion and mixing from steady-state CFD calculations. The mixing of massless particles and inertial particles is visualized by transforming the vector field from a traditionally Eulerian reference frame into a Lagrangian reference frame. Groups of particles are traced through the vector field for the mean path as well as their statistical dispersion about the mean position by using added scalar information about the root mean square value of the vector field and its Lagrangian time scale. In this way, clouds of particles in a turbulent environment are traced, not just mean paths. In combustion simulations of many industrial processes, good mixing is required to achieve a sufficient degree of combustion efficiency. The ability to visualize this multiphase mixing can not only help identify poor mixing but also explain the mechanism for poor mixing. The information gained from the visualization can be used to improve the overall combustion efficiency in utility boilers or propulsion devices. We have used this technique to visualize steady-state simulations of the combustion performance in several furnace designs.

Ma, Kwan-Liu↗

DSN Array Simulator

The DSN Array Simulator (wherein 'DSN' signifies NASA's Deep Space Network) is an updated version of software previously denoted the DSN Receive Array Technology Assessment Simulation. This software (see figure) is used for computational modeling of a proposed DSN facility comprising user-defined arrays of antennas and transmitting and receiving equipment for microwave communication with spacecraft on interplanetary missions. The simulation includes variations in spacecraft tracked and communication demand changes for up to several decades of future operation. Such modeling is performed to estimate facility performance, evaluate requirements that govern facility design, and evaluate proposed improvements in hardware and/or software. The updated version of this software affords enhanced capability for characterizing facility performance against user-defined mission sets. The software includes a Monte Carlo simulation component that enables rapid generation of key mission-set metrics (e.g., numbers of links, data rates, and date volumes), and statistical distributions thereof as functions of time. The updated version also offers expanded capability for mixed-asset network modeling--for example, for running scenarios that involve user-definable mixtures of antennas having different diameters (in contradistinction to a fixed number of antennas having the same fixed diameter). The improved version also affords greater simulation fidelity, sufficient for validation by comparison with actual DSN operations and analytically predictable performance metrics.

Tikidjian, Raffi↗

Dielectric screening and vacancy formation for large neutral and charged Si n H m ( n > 1500 ) nanocrystals using real-space pseudopotentials

A commonly used procedure for computing the properties of defects in crystalline materials is to consider a large supercell that includes the defect of interest. This is a straightforward technique as standard energy band codes can be used for such computations. For neutral defects, the only impediment of such an approach is to avoid defect-defect interactions between adjoining cells. However, this procedure can be complex if the defect of interest is charged as the system at large contains Coulombic divergences. Moreover, some have recently argued that the conventional definition of formation energies for charged defects cannot be reconciled with statistical mechanics. Here, we focus on an alternative approach. We consider large nanocrystals wherein a charged defect can be placed. Since the system is confined, a charged defect within the nanocrystal does not result in a Coulombic divergence. The chief impediment is computational, i.e., while no defect-defect or Coulombic divergences are present, the nanocrystal must be sufficiently large to allow the system to properly replicate a bulklike configuration. With the development of new algorithms and hardware advances, computations for systems of sufficient size to address this issue are feasible. In particular, we solve the Kohn-Sham equation in real space using pseudopotential-density-functional theory for large silicon nanocrystals, which contain thousands of atoms. Further, we focus on (i) the screening of a point charge and (ii) the formation of a charged vacancy in hydrogen-terminated silicon nanocrystals. This approach allows us to examine the role of quantum confinement in addition to exploring the bulk limit. Comparisons to other methods confirm the viability of this approach.

36 MATERIALS SCIENCE↗

Some analyses of the variability of atmospheric parameters at low altitudes significant for aircraft flyover noise measurements.

A number of special radiosondes were examined to study the statistical variabilities of atmospheric parameters on a time and space scale commensurate with aircraft flyover noise measurements. Particular emphasis is placed on those aspects of the data which relate to the temporal variabilities of the temperature and humidity profiles. It is found that during the late morning hours, continuous surface measurements and measurements of airport tower levels are sufficient to monitor the layer of the atmosphere in which the largest variabilities occur.

Chang, D. T.↗

X-ray Binary Luminosity Function Scaling Relations for Local Galaxies Based on Subgalactic Modeling

We present new Chandra constraints on the X-ray luminosity functions(XLFs)of X-ray binary(XRB)populations, as well as their scaling relations, for a sample of 38 nearby galaxies(D = 3.4–29 Mpc). Our galaxy sample is drawn primarily from the Spitzer Infrared Nearby Galaxies Survey(SINGS)and contains a wealth of Chandra(5.8 Ms total)and multiwavelength data, allowing for star formation rates(SFRs)and stellar masses(M⁎)to be measured on subgalactic scales. We divided the 2478 X-ray-detected sources into 21 subsamples in bins of specific SFR(sSFR ≡ SFR/M⁎)and constructed XLFs. To model the XLF dependence on sSFR, we fit a global XLF model, containing contributions from high-mass XRBs(HMXBs), low-mass XRBs(LMXBs), and background sources from the cosmic X-ray background that respectively scale with SFR, M⁎, and sky area. We find an HMXB XLF that is more complex in shape than previously reported and an LMXB XLF that likely varies with sSFR, potentially due to an age dependence. When applying our global model to XLF data for each individual galaxy, we discover a few galaxy XLFs that significantly deviate from our model beyond statistical scatter. Most notably, relatively low-metallicity galaxies have an excess of HMXBs above≈10^(38)erg/s, and elliptical galaxies that have relatively rich populations of globular clusters(GCs)show excesses of LMXBs compared to the global model. Additional modeling of how the XRB XLF depends on stellar age, metallicity, and GC specific frequency is required to sufficiently characterize the XLFs of galaxies.

Bret D. Lehmer↗

Initial Atmospheric Transport of Particles

A new version of the Initial Atmospheric Transport (IAT) model has been developed. This report fully characterizes the new model. IAT is designed to model heat release into the atmosphere that contains contaminated particles. The heat release forms a buoyant puff or plume. Buoyancy is quickly dissipated in 10's to 100's of seconds through drag and entrainment of ambient air. The final location of contaminant particles in the atmosphere after dissipation of abnormal heat is an important input to long term transportation codes such as the National Oceanic and Atmospheric Administration's (NOAA) Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) code. Ultimately the settling of contaminants after long term transport provides a basis for estimating health consequences. In comparison to the old version, the new model includes explicit simulation of particles, newly formulated equations of motion for a rising puff, and a new plume model. The plume model releases many instances of a single puff along a single trajectory in a way that balances energy released over a longer period than would be appropriate for a single puff. The puff model has additions for vorticity, turbulence, virtual mass, drag, and entrainment. The particle model leads to direct interfacing between IAT and the long-term transport code HYSPLIT. Previously, particle locations were given a pattern based on IAT's puff trajectory with no particle simulation. Now, particles are explicitly simulated, and 3-D spatial distributions of particles are output based on final positions. Ballistic particles are thrown outward at high velocity from explosive events and tend to hit the ground before the end of the simulation. Non-ballistic particles tend to get pulled into the rising puff s ring of vorticity with some escaping to ambient conditions as buoyant energy dissipates. A new model verification and validation effort was made for the new plume model. The Wallops Island Solid Rocket Propellant (SRP) fire plume was extracted to a 3-D transient dataset. After calibration, the model performs better than previous single puff attempts. Even so, error levels in the process used to estimate location and ambient weather conditions approach 30%. This leads to the conclusion that unplanned photogrammetry evaluations of the IAT model are useful but not sufficiently accurate to quantify model parameters. Like previous attempts, the vertical motion entrainment parameters come out much too high in comparison to values for other studies in the literature. It is suspected that a new model of shear entrainment would alleviate this problem. Unfortunately, IAT already has more calibration parameters than the data comparisons thus far can resolve with statistical certainty. This underscores the need to verify and validate the IAT model using detailed CFD studies that allow ideal statistical comparison with the IAT model for hundreds of cases. Even so, validation by comparison to carefully planned experiments (i.e. with acceptable error levels) is also necessary because the dynamics of entrainment are at a sub-grid level for CFD. Though better validation and verification are most important, many new enhancements are still envisioned for future IAT development. These include particle rain out due to condensation and agglomeration of water, direct connection to the North American Mesoscale (NAM) weather dataset, dynamic weather conditions, including ambient particle entrainment into the puff, and detrainment at the puff to ambient mixing interface. In general, including detailed aerosol dynamics is the greatest coming challenge which will require an entire new sub-model in IAT. This is very important if contaminant particles change in such a way that they increase or decrease in their potential harm to the environment.

54 ENVIRONMENTAL SCIENCES↗

Remote sensing applications to hydrologic modeling in the southern Sierra Nevada and portions of the San Joaquin Valley, volume 1

The author has identified the following significant results. Characteristics of LANDSAT MSS imagery present problems in using satellite radiation measurements to estimate the shortwave albedo of an alpine snow cover. Every 15 minute USGS quadrangle contains over 100,000 pixels which poses a computation problem if each pixel is to be evaluated individually. The sampling interval may be sufficiently great to mask some effects of terrain and vegetation on reflectance. Three frames of LANDSAT imagery are needed for complete coverage of the study area, yet less than one third of the area coverage from each frame covers an area of interest. Because of distortions inherent in the imagery, information regarding spacecraft altitude, attitude, and position must be statistically derived with respect to ground control points in the image whose geodetic locations are known. An inspection of shade points indicates that up to one third of the most heavily snow covered areas may saturate in bands 4 through 6. LANDSAT's 9 day repeat cycle is not optimum for snow cover reflectance modeling because the most pronounced changes in albedo occur most nearly following a new snowfall. Such a snowfall, occurring between overpasses, is inadequately represented by extrapolation from the previous overpasses.

Dozier, J.↗