Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “applied 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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

Global Intercomparison of hyper-resolution ECOSTRESS coastal sea surface temperature measurements from the Space Station with VIIRS-N20

The ECOSTRESS multi-channel thermal radiometer on the Space Station has an unprecedented spatial resolution of 70 m and a return time of hours to 5 days. It resolves details of oceanographic features not detectable in imagery from MODIS or VIIRS, and has open-ocean coverage, unlike Landsat. We calibrated two years of ECOSTRESS sea surface temperature observations with L2 data from VIIRS-N20 (2019–2020) worldwide but especially focused on important upwelling systems currently undergoing climate change forcing. Unlike operational SST products from VIIRS-N20, the ECOSTRESS surface temperature algorithm does not use a regression approach to determine temperature, but solves a set of simultaneous equations based on first principles for both surface temperature and emissivity. We compared ECOSTRESS ocean temperatures to well-calibrated clear sky satellite measurements from VIIRS-N20. Data comparisons were constrained to those within 90 min of one another using co-located clear sky VIIRS and ECOSTRESS pixels. ECOSTRESS ocean temperatures have a consistent 1.01 ◦C negative bias relative to VIIRS-N20, although deviation in brightness temperatures within the 10.49 and 12.01 µm bands were much smaller. As an alternative, we compared the performance of NOAA, NASA, and U.S. Navy operational split-window SST regression algorithms taking into consideration the statistical limitations imposed by intrinsic SST spatial autocorrelation and applying corrections on brightness temperatures. We conclude that standard bias-correction methods using already validated and well-known algorithms can be applied to ECOSTRESS SST data, yielding highly accurate products of ultra-high spatial resolution for studies of biological and physical oceanography in a time when these are needed to properly evaluate regional and even local impacts of climate change.

Nicolas Weidberg↗
DOE OSTI2022

Large-scale application of free energy perturbation calculations for antibody design

Abstract Alchemical free energy perturbation (FEP) is a rigorous and powerful technique to calculate the free energy difference between distinct chemical systems. Here we report our implementation of automated large-scale FEP calculations, using the Amber software package, to facilitate antibody design and evaluation. In combination with Hamiltonian replica exchange, our FEP simulations aim to predict the effect of mutations on both the binding affinity and the structural stability. Importantly, we incorporate multiple strategies to faithfully estimate the statistical uncertainties in the FEP results. As a case study, we apply our protocols to systematically evaluate variants of the m396 antibody for their conformational stability and their binding affinity to the spike proteins of SARS-CoV-1 and SARS-CoV-2. By properly adjusting relevant parameters, the particle collapse problems in the FEP simulations are avoided. Furthermore, large statistical errors in a small fraction of the FEP calculations are effectively reduced by extending the sampling, such that acceptable statistical uncertainties are achieved for the vast majority of the cases with a modest total computational cost. Finally, our predicted conformational stability for the m396 variants is qualitatively consistent with the experimentally measured melting temperatures. Our work thus demonstrates the applicability of FEP in computational antibody design.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Statistical Mechanics on Lattices (smol) v0.0.1

smol is a Python implementation of the Cluster Expansion Method and extensions of this methodology; which is used to fit applied lattice models from first principle calculations for subsequent use in Monte Carlo simulations for computing thermodynamic and statistical properties of atomic configuration. smol has a several notable advantages over other similar software packages. The first being its Python implementation which makes it easy to use even for users with little software and coding skills. Despite being implemented in Python, essential Monte Carlo routines are implemented in Cython such that performance is not compromised and is competitive even to similar software implemented in C/C++. Furthermore, the package has specific functionality for handling complex ionic materials (such as cathodes and electrolytes) that is not readily available in other packages. Lastly, the software has a flexible and modular design with the intention of making it fast and efficient to develop and extend the methodology.

Barroso-Luque, Luis↗
NASA NTRS1980

Some adaptive filtering techniques applied to the passive remote sensing problem

A review is made of recursive statistical regression techniques incorporating past or past and future observations through smoothing and Kalman filtering, respectively; with results for the cases of the Tiros-N/MSU and Nimbus-6/Scams remote sensing satellite experiments. In response to the lack of a satisfactory model for the medium sounded, which is presently a major limitation on retrieval technique performance, a novel, global approach is proposed which casts the retrieval problem into the framework of adaptive filtering. A numerical implementation of such an adaptive system is presented, with a multilayer, semi-spectral general circulation model for the atmosphere being used to fine-tune the sensor as well as the dynamical equations of a Kalman filter. It is shown that the assimilation of radiometric data becomes a straightforward subproblem.

Toldalagi, P. M.↗
NASA NTRS1990

A method of classification for multisource data in remote sensing based on interval-valued probabilities

An axiomatic approach to intervalued (IV) probabilities is presented, where the IV probability is defined by a pair of set-theoretic functions which satisfy some pre-specified axioms. On the basis of this approach representation of statistical evidence and combination of multiple bodies of evidence are emphasized. Although IV probabilities provide an innovative means for the representation and combination of evidential information, they make the decision process rather complicated. It entails more intelligent strategies for making decisions. The development of decision rules over IV probabilities is discussed from the viewpoint of statistical pattern recognition. The proposed method, so called evidential reasoning method, is applied to the ground-cover classification of a multisource data set consisting of Multispectral Scanner (MSS) data, Synthetic Aperture Radar (SAR) data, and digital terrain data such as elevation, slope, and aspect. By treating the data sources separately, the method is able to capture both parametric and nonparametric information and to combine them. Then the method is applied to two separate cases of classifying multiband data obtained by a single sensor. In each case a set of multiple sources is obtained by dividing the dimensionally huge data into smaller and more manageable pieces based on the global statistical correlation information. By a divide-and-combine process, the method is able to utilize more features than the conventional maximum likelihood method.

Kim, Hakil↗
DOE OSTI2022

A phase-shift-periodic parallel boundary condition for low-magnetic-shear scenarios

Abstract We formulate a generalized periodic boundary condition as a limit of the standard twist-and-shift parallel boundary condition that is suitable for simulations of plasmas with low magnetic shear. This is done by applying a phase shift in the binormal direction when crossing the parallel boundary. While this phase shift can be set to zero without loss of generality in the local flux-tube limit when employing the twist-and-shift boundary condition, we show that this is not the most general case when employing periodic parallel boundaries, and may not even be the most desirable. A non-zero phase shift can be used to avoid the convective cells that plague simulations of the three-dimensional Hasegawa–Wakatani system, and is shown to have measurable effects in periodic low-magnetic-shear gyrokinetic simulations. We propose a numerical program where a sampling of periodic simulations at random pseudo-irrational flux surfaces are used to determine physical observables in a statistical sense. This approach can serve as an alternative to applying the twist-and-shift boundary condition to low-magnetic-shear scenarios, which, while more straightforward, can be computationally demanding.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2023

Spacesuit Fit and Mobility Assessments by Digital Human Modeling

Spacesuit Fit and Mobility Assessments by Digital Human Modeling K. Han Kim (Leidos, Inc.) Elizabeth A. Benson (KBR, Inc.) Sudhakar L. Rajulu (NASA Johnson Space Center) Spacesuits are required to accommodate safe operations for astronauts across gender and a wide variety of body shapes and sizes. This goal has been of particular importance given the increasing diversity of NASA crewmembers for upcoming Missions. While testing with design prototypes is a critical step for spacesuit development, iterative mockup design, fabrication, and human subject tests can be extremely costly and time consuming. Moreover, testing with a limited subject pool has often raised questions for validity, as test subjects need to represent the entire astronaut population, not only of the past or current, but also the future. This study is aimed at demonstrating how digital human modeling (DHM) tools have been built and directly supported NASA spacesuit developments. With DHM, the computer aided design (CAD) model of a spacesuit was integrated with human body models. Two use scenarios are presented in this paper, namely fit and mobility. For fit assessments, the suit-to-body contact and compression patterns were estimated using 3D human body scan models virtually wearing a spacesuit model. A statistical fit classifier was made from the contact patterns and applied against a large database of body scans (N=2,500). With this technique, the NASA reference design spacesuit Exploration EVA Mobility Unit (xEMU)was verified to accommodate 90% of the astronaut-like population, with the critical dimensions covering 1stto 99thpercentiles of the target body measurements. A similar technique assessed XEMU mobility. The maximum reach envelopes were considered, within which the work objects and critical hardware interfaces should be located for safe and ergonomic operations. While the reach envelope geometry varies significantly with the suit wearer’s body size and strength, the existing test data did not include the subjects critically required to define suit mobility requirements, such as very small females or large males. Using the xEMU virtual model kinematically simulated and permuted for a hypothetical wearer, however, the existing data were statistically transformed and scaled. This method enabled for a parametric estimation of the reach envelopes from the 1stpercentile female or 99thpercentile male. The outcome was successfully incorporated for requirement developments. With the new DHM tools, human integration of the spacesuit was structurally simulated and predictively assessed, which would have been otherwise impossible. Also, the needs for iterative mockup and subject tests were significantly reduced, which resulted in time and cost saving. Additional work is in progress to integrate additional vehicles, tools, and hardware with DHM framework.

K Han Kim↗
DOE OSTI2024

Simulation-trained machine learning models for Lorentz transmission electron microscopy

Understanding the collective behavior of complex spin textures, such as lattices of magnetic skyrmions, is of fundamental importance for exploring and controlling the emergent ordering of these spin textures and inducing phase transitions. It is also critical to understand the skyrmion–skyrmion interactions for applications such as magnetic skyrmion-enabled reservoir or neuromorphic computing. Magnetic skyrmion lattices can be studied using in situ Lorentz transmission electron microscopy (LTEM), but quantitative and statistically robust analysis of the skyrmion lattices from LTEM images can be difficult. In this work, we show that a convolutional neural network, trained on simulated data, can be applied to perform segmentation of spin textures and to extract quantitative data, such as spin texture size and location, from experimental LTEM images, which cannot be obtained manually. This includes quantitative information about skyrmion size, position, and shape, which can, in turn, be used to calculate skyrmion–skyrmion interactions and lattice ordering. We apply this approach to segmenting images of Néel skyrmion lattices so that we can accurately identify skyrmion size and deformation in both dense and sparse lattices. The model is trained using a large set of micromagnetic simulations as well as simulated LTEM images. This entirely open-source training pipeline can be applied to a wide variety of magnetic features and materials, enabling large-scale statistical studies of spin textures using LTEM.

McCray, Arthur R. C. (ORCID:0000000160774698)↗
NASA NTRS1975

Statistical analysis of trim maneuvers in low thrust interplanetary navigation

An analytical technique for the statistical analysis of a low thrust trim maneuver in interplanetary navigation has been developed. The maneuver is treated rigorously as a nonlinear function of the trajectory errors to be removed, and the corresponding nongaussian statistics are developed. The method determines the statistical properties of the maneuver direction and duration, as well as the final state dispersion covariance. The analytical technique is applied to the terminal maneuver in a 1980 Encke slow flyby mission and is found to yield accurate statistical results with much less effort than Monte Carlo simulation. The method permits consideration of a shorter terminal maneuver than does linear analysis, with a consequent improvement in delivery accuracy.

Rinker, G. C.↗
NASA NTRS1995

Analysis of the progressive failure of brittle matrix composites

This report investigates two of the most common modes of localized failures, namely, periodic fiber-bridged matrix cracks and transverse matrix cracks. A modification of Daniels' bundle theory is combined with Weibull's weakest link theory to model the statistical distribution of the periodic matrix cracking strength for an individual layer. Results of the model predictions are compared with experimental data from the open literature. Extensions to the model are made to account for possible imperfections within the layer (i.e., nonuniform fiber lengths, irregular crack spacing, and degraded in-situ fiber properties), and the results of these studies are presented. A generalized shear-lag analysis is derived which is capable of modeling the development of transverse matrix cracks in material systems having a general multilayer configuration and under states of full in-plane load. A method for computing the effective elastic properties for the damaged layer at the global level is detailed based upon the solution for the effects of the damage at the local level. This methodology is general in nature and is therefore also applicable to (0(sub m)/90(sub n))(sub s) systems. The characteristic stress-strain response for more general cases is shown to be qualitatively correct (experimental data is not available for a quantitative evaluation), and the damage evolution is recorded in terms of the matrix crack density as a function of the applied strain. Probabilistic effects are introduced to account for the statistical nature of the material strengths, thus allowing cumulative distribution curves for the probability of failure to be generated for each of the example laminates. Additionally, Oh and Finney's classic work on fracture location in brittle materials is extended and combined with the shear-lag analysis. The result is an analytical form for predicting the probability density function for the location of the next transverse crack occurrence within a crack bounded region. The results of this study verified qualitatively the validity of assuming a uniform crack spacing (as was done in the shear-lag model).

Thomas, David J.↗
DOE OSTI2021

Towards testing the theory of gravity with DESI: summary statistics, model predictions and future simulation requirements

Shortly after its discovery, General Relativity (GR) was applied to predict the behavior of our Universe on the largest scales, and later became the foundation of modern cosmology. Its validity has been verified on a range of scales and environments from the Solar system to merging black holes. However, experimental confirmations of GR on cosmological scales have so far lacked the accuracy one would hope for — its applications on those scales being largely based on extrapolation and its validity there sometimes questioned in the shadow of the discovery of the unexpected cosmic acceleration. Future astronomical instruments surveying the distribution and evolution of galaxies over substantial portions of the observable Universe, such as the Dark Energy Spectroscopic Instrument (DESI), will be able to measure the fingerprints of gravity and their statistical power will allow strong constraints on alternatives to GR. In this paper, based on a set of N-body simulations and mock galaxy catalogs, we study the predictions of a number of traditional and novel summary statistics beyond linear redshift distortions in two well-studied modified gravity models — chameleon f(R) gravity and a braneworld model — and the potential of testing these deviations from GR using DESI. These summary statistics employ a wide array of statistical properties of the galaxy and the underlying dark matter field, including two-point and higher-order statistics, environmental dependence, redshift space distortions and weak lensing. We find that they hold promising power for testing GR to unprecedented precision. The major future challenge is to make realistic, simulation-based mock galaxy catalogs for both GR and alternative models to fully exploit the statistic power of the DESI survey (by matching the volumes and galaxy number densities of the mocks to those in the real survey) and to better understand the impact of key systematic effects. Using these, we identify future simulation and analysis needs for gravity tests using DESI.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1996

NASA Rat Acoustic Tolerance Test 1994-1995: 8 kHz, 16 kHz, 32 kHz Experiments

Adult male Sprague-Dawley rats were exposed to chronic applied sound (74 to 79 dB, SPL) with octave band center frequencies of either 8, 16 or 32 kHz for up to 60 days. Control cages had ambient sound levels of about 62 dB (SPL). Groups of rats (test vs. control; N=9 per group) were euthanized after 0. 5. 14, 30, and 60 days. On each euthanasia day, objective evaluation of their physiology and behavior was performed using a Stress Assessment Battery (SAB) of measures. In addition, rat hearing was assessed using the brain stem auditory evoked potential (BAER) method after 60 days of exposure. No statistically significant differences in mean daily food use could be attributed to the presence of the applied test sound. Test rats used 5% more water than control rats. In the 8 kHz and 32 kHz tests this amount was statistically significant(P less than .05). This is a minor difference of questionable physiological significance. However, it may be an indication of a small reaction to the constant applied sound. Across all test frequencies, day 5 test rats had 6% larger spleens than control rats. No other body or organ weight differences were found to be statistically significant with respect to the application of sound. This spleen effect may be a transient adaptive process related to adaptation to the constant applied noise. No significant test effect on differential white blood cell counts could be demonstrated. One group demonstrated a low eosinophil count (16 kHz experiment, day 14 test group). However this was highly suspect. Across all test frequencies studied, day 5 test rats had 17% fewer total leukocytes than day 5 control rats. Sound exposed test rats exhibited 44% lower plasma corticosterone concentrations than did control rats. Note that the plasma corticosterone concentration was lower in the sound exposed test animals than the control animals in every instance (frequency exposure and number of days exposed).

Mele, Gary D.↗
NASA NTRS1999

Real-Time Assimilation of Goes-Derived Products into A Mesoscale Model and It's Impact on Short-Term (06-36hr) Forecasts from 17 October 1998 through the Present

As the parameterizations of surface energy budgets in regional models have become more complete physically, models have the potential to be much more realistic in simulations of coupling between surface radiation, hydrology, and surface energy transfer. Realizing the importance of properly specifying the surface energy budget, many institutions are using land-surface models to represent the lower boundary forcing associated with biophysical processes and soil hydrology. However, the added degrees of freedom due to inclusion of such land-surface schemes require the specification of additional parameters within the model system such as vegetative resistances, green vegetation fraction, leaf area index, soil physical and hydraulic characteristics, stream flow, runoff, and the vertical distribution of soil moisture. Spatial heterogeneity of these parameters makes correct specification problematic since measurements are not routinely available. A technique has been developed for assimilating GOES-IR skin temperature tendencies, solar insolation, and surface albedo into the surface energy budget equation of a mesoscale model so that the simulated rate of temperature change closely agrees with the satellite observations. The technique has been successfully employed in a number of mesoscale models in case-study mode. We have taken the next step and developed a study to determine if assimilating these types of data into mesoscale models in real-time can improve short-term (648h) forecasts of temperature, relative humidity, and QPF on a daily basis over relatively large regions. Therefore, an operational modeling/assimilation system has been developed at the GHCC during the past summer that allows us to produce simulations out to 48 hours in a timely manor. The PSU/NCAR MM5 is used in a nested configuration with a 25 km grid covering the southeastern third of the US. The model has been on-line since 1 July 1998 and forecast products are posted on our web site. The satellite algorithms that generate data to be assimilated came on-line 17 October 1998. Quantitative assessment of the forecast quality is performed via traditional verification statistics. In addition, invaluable qualitative information is obtained through close collaboration with several NWSFO's who are using the MM5 products in real-time on a daily basis. The assimilation technique has been applied in an off-line mode since 17 October. Results based on bulk statistical verification of surface meteorology over the entire Southeastern US show that assimilating the GOES-derived land surface tendencies and solar radiation results in a significant reduction of the shelter air temperature and RH bias on a daily basis. In fact, the assimilation technique has produced improved temperature and RH forecasts for 97% of the 100 simulations performed to date. Work is currently underway to determine the sensitivity of the assimilation procedure to the availability of satellite data, length of assimilation period, model initialization, and synoptic-scale meteorological conditions. In addition, results from a detailed energy budget analysis using the Early Eta, our operational MM5, and the assimilation runs will help us to better understand the satellite assimilation the land-surface energy budge. Research during the spring-summer of 1999 will focus on the impact of the assimilation technique during the warm season where it is hypothesized that it can have a positive impact on QPF during conditions of weak synoptic-scale forcing.

Lapenta, William M.↗
DOE OSTI2021

Adaptively controlled fast production of defect-free beryllium ion crystals using pulsed laser ablation

Trapped atomic ions find wide applications ranging from precision measurement to quantum information science and quantum computing. Beryllium ions are widely used due to the light mass and convenient atomic structure of beryllium; however, conventional ion loading from thermal ovens exerts undesirable gas loads for a prolonged duration. Here, we demonstrate a method to rapidly produce pure linear chains of beryllium ions with pulsed laser ablation, serving as a starting point for large-scale quantum information processing. Our method is fast compared to thermal ovens, reduces the gas load to only 10 -12 Torr (10 -10 Pa) level, yields a short recovery time of a few seconds, and also eliminates the need for a deep ultraviolet laser for photoionization. We also study the loading dynamics, which show non-Poissonian statistics in the presence of sympathetic cooling. In addition, we apply feedback control to obtain defect-free ion chains with desirable lengths.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Energy dependent calculations of fission product, prompt, and delayed neutron yields for neutron induced fission on 235 U, 238 U, and 239 Pu

In this article, we perform energy-dependent calculations of independent and cumulative fission product yields for 235 U, 238 U, and 239 Pu in the first chance fission region. Starting with the primary fission fragment distributions taken from available experimental data and analytical functions based on assumptions for the excitation energy and spin-parity distributions, the Hauser-Feshbach statistical decay treatment for fission fragment de-excitation is applied to more than 1,000 fission fragments for the incident neutron energies up to 5 MeV. The calculated independent yields are then used as an input of β-decay calculations to produce the cumulative yield, and summation calculations are performed. Model parameters in these procedures are adjusted by applying the Bayesian technique at the thermal energy for 235 U and 239 Pu and in the fast energy range for 238 U. The calculated fission observable quantities, such as the energy-dependent cumulative yields, and prompt and delayed neutron yields, are compared with available experimental data. We also study the possible impact of the second chance fission opening on the energy dependence of the delayed neutron yield by on the opening of second chance fission on the energy dependence calculation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

A Bayesian interpretation of inconsistency measures in cosmology

Abstract Measures of inconsistency and tension between datasets have become an essential part of cosmological analyses. It is important to accurately evaluate the significance of such tensions when present. We propose here a Bayesian interpretation of inconsistency measures that can extract information about physical inconsistencies in the presence of data scatter. This new framework is based on the conditional probability distribution of the level of physical inconsistency given the obtained value of the measure. We use the index of inconsistency as a case study to illustrate the new interpretation framework, but this can be generalized to other metrics. Importantly, there are two aspects in the quantification of inconsistency that behave differently as the number of model parameters increases. The first is the probability for the level of physical inconsistency to reach a threshold which drops with the increase of the number of parameters under consideration. The second is the actual level of physical inconsistency which remains rather insensitive to such an increase in parameters. The difference between these two aspects is often overlooked, which leads to a long-standing ambiguity: when a given inconsistency is found between two constraints, its “significance” seems to be lower when considered in a higher-dimensional parameter space. This ambiguity is resolved by the Bayesian interpretation we introduce in this work because the conditional probability distribution includes all the statistical information of the level of physical inconsistency. Finally, we apply the Bayesian interpretation to examine the (in)consistency between Planck versus the Cepheid-based local measurement, the Dark Energy Survey (DES), the Atacama Cosmology Telescope (ACT) and WMAP. We confirm and revisit the degrees of previous physical inconsistencies and show the stability of the new interpretation with respect to the number of cosmological parameters compared to the commonly used n-σ interpretation when applied to cosmological tensions in multi-parameter spaces.

Astronomy & Astrophysics↗
DOE OSTI2023

Combining cosmic shear data with correlated photo-z uncertainties: constraints from DESY1 and HSC-DR1

An accurate calibration of the source redshift distribution p(z) is a key aspect in the analysis of cosmic shear data. This, one way or another, requires the use of spectroscopic or high-quality photometric samples. However, the difficulty to obtain colour-complete spectroscopic samples matching the depth of weak lensing catalogs means that the analyses of different cosmic shear datasets often use the same samples for redshift calibration. This introduces a source of statistical and systematic uncertainty that is highly correlated across different weak lensing datasets, and which must be accurately characterised and propagated in order to obtain robust cosmological constraints from their combination. In this paper we introduce a method to quantify and propagate the uncertainties on the source redshift distribution in two different surveys sharing the same calibrating sample. The method is based on an approximate analytical marginalisation of the p(z) statistical uncertainties and the correlated marginalisation of residual systematics. We apply this method to the combined analysis of cosmic shear data from the DESY1 data release and the HSC-DR1 data, using the COSMOS 30-band catalog as a common redshift calibration sample. We find that, although there is significant correlation in the uncertainties on the redshift distributions of both samples, this does not change the final constraints on cosmological parameters significantly. The same is true also for the impact of residual systematic uncertainties from the errors in the COSMOS 30-band photometric redshifts. Additionally, we show that these effects will still be negligible in Stage-IV datasets. Finally, the combination of DESY1 and HSC-DR1 allows us to constrain the "clumpiness" parameter to S 8 = ${0.768}_{-0.017}^{+0.021}$. This corresponds to a ~√(2) improvement in uncertainties with respect to either DES or HSC alone.

79 ASTRONOMY AND ASTROPHYSICS↗