Search NASA⌕ Search

SEARCH · Search NASA

Results for “Methods:statistical”

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,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

The influence of temperature on the strain-hardening behavior of Fe-22/25/28Mn-3Al-3Si TRIP/TWIP steels

In this work, the influence of temperature and stacking fault energy (SFE) on the strain-hardening behavior and critical resolved shear stress for twinning was investigated for three Fe–22/25/28Mn–3Al–3Si wt.% transformation- and twinning-induced plasticity (TRIP/TWIP) steels. The SFEs were calculated by two different methods, density functional theory and statistical thermodynamic modeling. The dislocation structure, observed at low levels of plastic deformation, transitions from “planar” to “wavy” dislocation glide with an increase in temperature, Mn content, and/or SFE. The change in dislocation glide mechanisms from planar to wavy reduces the strain hardening rate, in part due to fewer planar obstacles and greater cross slip activity. Additionally, the alloys exhibit a large decrease in strength and ductility with increasing temperature from 25 to 200 °C, attributed to a substantial reduction in the thermally activated component of the flow stress, predominate suppression of TRIP and TWIP, and a significant increase in the critical resolved shear stress for mechanical twinning. Interestingly, the increase in SFE with temperature had a rather minor influence on the critical resolved shear stress for mechanical twinning, and other temperature dependent factors which likely play a more dominant role are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Machine Learning technique for isotopic determination of radioisotopes using HPGe γ-ray spectra

γ-ray spectroscopy is a quantitative, non-destructive technique that may be utilized for the identification and quantitative isotopic estimation of radionuclides. Traditional methods of isotopic determination have various challenges that contribute to statistical and systematic uncertainties in the estimated isotopics. Furthermore, these methods typically require numerous pre-processing steps, and have only been rigorously tested in laboratory settings with limited shielding. Here, in this work, we examine the application of a number of machine learning based regression algorithms as alternatives to conventional approaches for analyzing γ-ray spectroscopy data in the Emergency Response arena. This approach not only eliminates many steps in the analysis procedure, and therefore offers potential to reduce this source of systematic uncertainty, but is also shown to offer comparable performance to conventional approaches in the Emergency Response Application.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

A Bayesian analysis of classical shadows

Abstract The method of classical shadows proposed by Huang, Kueng, and Preskill heralds remarkable opportunities for quantum estimation with limited measurements. Yet its relationship to established quantum tomographic approaches, particularly those based on likelihood models, remains unclear. In this article, we investigate classical shadows through the lens of Bayesian mean estimation (BME). In direct tests on numerical data, BME is found to attain significantly lower error on average, but classical shadows prove remarkably more accurate in specific situations—such as high-fidelity ground truth states—which are improbable in a fully uniform Hilbert space. We then introduce an observable-oriented pseudo-likelihood that successfully emulates the dimension-independence and state-specific optimality of classical shadows, but within a Bayesian framework that ensures only physical states. Our research reveals how classical shadows effect important departures from conventional thinking in quantum state estimation, as well as the utility of Bayesian methods for uncovering and formalizing statistical assumptions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Learning reduced kinetic Monte Carlo models of complex chemistry from molecular dynamics

Here, we propose a novel statistical learning framework for automatically and efficiently building reduced kinetic Monte Carlo (KMC) models of large-scale elementary reaction networks from data generated by a single or few molecular dynamics simulations (MD). Existing approaches for identifying species and reactions from molecular dynamics typically use bond length and duration criteria, where bond duration is a fixed parameter motivated by an understanding of bond vibrational frequencies. Conversely, we show that for highly reactive systems, bond duration should be a model parameter that is chosen to maximize the predictive power of the resulting statistical model. We demonstrate our method on a high temperature, high pressure system of reacting liquid methane, and show that the learned KMC model is able to extrapolate more than an order of magnitude in time for key molecules. Additionally, our KMC model of elementary reactions enables us to isolate the most important set of reactions governing the behavior of key molecules found in the MD simulation. We develop a new data-driven algorithm to reduce the chemical reaction network which can be solved either as an integer program or efficiently using L1 regularization, and compare our results with simple count-based reduction. For our liquid methane system, we discover that rare reactions do not play a significant role in the system, and find that less than 7% of the approximately 2000 reactions observed from molecular dynamics are necessary to reproduce the molecular concentration over time of methane. Furthermore, we describe a framework in this work that paves the way towards a genomic approach to studying complex chemical systems, where expensive MD simulation data can be reused to contribute to an increasingly large and accurate genome of elementary reactions and rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2012

Quasi-equilibrium electron density along a magnetic field line

A methodology is developed to determine the density of high-energy electrons along a magnetic field line for a low-{beta} plasma. This method avoids the expense and statistical noise of traditional particle tracking techniques commonly used for high-energy electrons in bombardment plasma generators. By preserving the magnetic mirror and assuming a mixing timescale, typically the elastic collision frequency with neutrals, a quasi-equilibrium electron distribution can be calculated. Following the transient decay, the analysis shows that both the normalized density and the reduction fraction due to collision converge to a single quasi-equilibrium solution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Robustness analysis of an air heating plant and control law by using polynomial chaos

This paper presents a robustness analysis of an air heating plant with a multivariable closed-loop control law by using the polynomial chaos methodology (MPC). The plant consists of a PVC tube with a fan in the air input (that forces the air through the tube) and a mass flux sensor in the output. A heating resistance warms the air as it flows inside the tube, and a thermo-couple sensor measures the air temperature. The plant has thus two inputs (the fan's rotation intensity and heat generated by the resistance, both measured in percent of the maximum value) and two outputs (air temperature and air mass flux, also in percent of the maximal value). The mathematical model is obtained by System Identification techniques. The mass flux sensor, which is nonlinear, is linearized and the delays in the transfer functions are properly approximated by non-minimum phase transfer functions. The resulting model is transformed to a state-space model, which is used for control design purposes. The multivariable robust control design techniques used is the LQG/LTR, and the controllers are validated in simulation software and in the real plant. Finally, the MPC is applied by considering some of the system's parameters as random variables (one at a time, and the system's stochastic differential equations are solved by expanding the solution (a stochastic process) in an orthogonal basis of polynomial functions of the basic random variables. This method transforms the stochastic equations in a set of deterministic differential equations, which can be solved by traditional numerical methods (That is the MPC). Statistical data for the system (like expected values and variances) are then calculated. The effects of randomness in the parameters are evaluated in the open-loop and closed-loop pole's positions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Elucidating plasma dynamics in Hasegawa–Wakatani turbulence by information geometry

In this work, the impact of adiabatic electrons on drift-wave turbulence, modeled by the Hasegawa–Wakatani equations, is studied using information length. Information length is a novel theoretical method for measuring distances between statistical states represented by different probability distribution functions (PDFs) along the path of a system and represents the total number of statistically different states that a system evolves through in time. Specifically, the time-dependent PDFs of turbulent fluctuations for a given adiabatic index $A$ are computed. The changes in fluctuation statistics are then quantified in time by using information length. The numerical results provide time traces exhibiting intermittent plasma dynamics, and such behavior is identified by a rapid change in the information length. The effects of $A$ are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Diffusion and stochastic island generation in the magnetic field line random walk

The cross-field diffusion of field lines in stochastic magnetic fields described by the 2D+slab model is studied using a semi-analytic statistical approach, the decorrelation trajectory method. We show that field line trapping and the associated stochastic magnetic islands strongly influence the diffusion coefficients, leading to dependences on the parameters that are different from the quasilinear and Bohm regimes. A strong amplification of the diffusion is produced by a small slab field in the presence of trapping. The diffusion regimes are determined and the corresponding physical processes are identified.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

Crossing statistic: reconstructing the expansion history of the universe

We present that by combining Crossing Statistic [1,2] and Smoothing method [3-5] one can reconstruct the expansion history of the universe with a very high precision without considering any prior on the cosmological quantities such as the equation of state of dark energy. We show that the presented method performs very well in reconstruction of the expansion history of the universe independent of the underlying models and it works well even for non-trivial dark energy models with fast or slow changes in the equation of state of dark energy. Accuracy of the reconstructed quantities along with independence of the method to any prior or assumption gives the proposed method advantages to the other non-parametric methods proposed before in the literature. Applying on the Union 2.1 supernovae combined with WiggleZ BAO data we present the reconstructed results and test the consistency of the two data sets in a model independent manner. Results show that latest available supernovae and BAO data are in good agreement with each other and spatially flat ΛCDM model is in concordance with the current data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Estimating covariance matrices for two- and three-point correlation function moments in Arbitrary Survey Geometries

ABSTRACT We present configuration-space estimators for the auto- and cross-covariance of two- and three-point correlation functions (2PCF and 3PCF) in general survey geometries. These are derived in the Gaussian limit (setting higher order correlation functions to zero), but for arbitrary non-linear 2PCFs (which may be estimated from the survey itself), with a shot-noise rescaling parameter included to capture non-Gaussianity. We generalize previous approaches to include Legendre moments via a geometry-correction function calibrated from measured pair and triple counts. Making use of importance sampling and random particle catalogues, we can estimate model covariances in fractions of the time required to do so with mocks, obtaining estimates with negligible sampling noise in ∼10 (∼100) CPU-hours for the 2PCF (3PCF) autocovariance. We compare results to sample covariances from a suite of BOSS DR12 mocks and find the matrices to be in good agreement, assuming a shot-noise rescaling parameter of 1.03 (1.20) for the 2PCF (3PCF). To obtain strongest constraints on cosmological parameters, we must use multiple statistics in concert; having robust methods to measure their covariances at low computational cost is thus of great relevance to upcoming surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

Use of knowledge-based restraints in phenix.refine to improve macromolecular refinement at low resolution

Traditional methods for macromolecular refinement often have limited success at low resolution (3.0–3.5 Å or worse), producing models that score poorly on crystallographic and geometric validation criteria. To improve low-resolution refinement, knowledge from macromolecular chemistry and homology was used to add three new coordinate-restraint functions to the refinement program phenix.refine . Firstly, a `reference-model' method uses an identical or homologous higher resolution model to add restraints on torsion angles to the geometric target function. Secondly, automatic restraints for common secondary-structure elements in proteins and nucleic acids were implemented that can help to preserve the secondary-structure geometry, which is often distorted at low resolution. Lastly, we have implemented Ramachandran-based restraints on the backbone torsion angles. In this method, a φ,ψ term is added to the geometric target function to minimize a modified Ramachandran landscape that smoothly combines favorable peaks identified from nonredundant high-quality data with unfavorable peaks calculated using a clash-based pseudo-energy function. All three methods show improved MolProbity validation statistics, typically complemented by a lowered R free and a decreased gap between R work and R free .

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Representational analysis of extended disorder in atomistic ensembles derived from total scattering data

With the increased availability of high-intensity time-of-flight neutron and synchrotron X-ray scattering sources that can access wide ranges of momentum transfer, the pair distribution function method has become a standard analysis technique for studying disorder of local coordination spheres and at intermediate atomic separations. In some cases, rational modeling of the total scattering data (Bragg and diffuse) becomes intractable with least-squares approaches, necessitating reverse Monte Carlo simulations using large atomistic ensembles. However, the extraction of meaningful information from the resulting atomistic ensembles is challenging, especially at intermediate length scales. Representational analysis is used here to describe the displacements of atoms in reverse Monte Carlo ensembles from an ideal crystallographic structure in an approach analogous to tight-binding methods. Rewriting the displacements in terms of a local basis that is descriptive of the ideal crystallographic symmetry provides a robust approach to characterizing medium-range order (and disorder) and symmetry breaking in complex and disordered crystalline materials. Lastly, this method enables the extraction of statistically relevant displacement modes (orientation, amplitude and distribution) of the crystalline disorder and provides directly meaningful information in a locally symmetry-adapted basis set that is most descriptive of the crystal chemistry and physics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

BraggNet: integrating Bragg peaks using neural networks

Neutron crystallography offers enormous potential to complement structures from X-ray crystallography by clarifying the positions of low-Z elements, namely hydrogen. Macromolecular neutron crystallography, however, remains limited, in part owing to the challenge of integrating peak shapes from pulsed-source experiments. To advance existing software, this article demonstrates the use of machine learning to refine peak locations, predict peak shapes and yield more accurate integrated intensities when applied to whole data sets from a protein crystal. The artificial neural network, based on the U-Net architecture commonly used for image segmentation, is trained using about 100 000 simulated training peaks derived from strong peaks. After 100 training epochs (a round of training over the whole data set broken into smaller batches), training converges and achieves a Dice coefficient of around 65%, in contrast to just 15% for negative control data sets. Integrating whole peak sets using the neural network yields improved intensity statistics compared with other integration methods, including k-nearest neighbours. These results demonstrate, for the first time, that neural networks can learn peak shapes and be used to integrate Bragg peaks. In conclusion, it is expected that integration using neural networks can be further developed to increase the quality of neutron, electron and X-ray crystallography data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

K-space algorithmic reconstruction (KAREN): a robust statistical methodology to separate Bragg and diffuse scattering

Diffuse scattering occurring in the Bragg diffraction pattern of a long-range-ordered structure represents local deviation from the governing regular lattice. However, interpreting the real-space structure from the diffraction pattern presents a significant challenge because of the dramatic difference in intensity between the Bragg and diffuse components of the total scattering function. In contrast to the sharp Bragg diffraction, the diffuse signal has generally been considered to be a weak expansive or continuous background signal. In this paper, using 1D and 2D models, it is demonstrated that diffuse scattering in fact consists of a complex array of high-frequency features that must not be averaged into a low-frequency background signal. To evaluate the actual diffuse scattering effectively, an algorithm has been developed that uses robust statistics and traditional signal processing techniques to identify Bragg peaks as signal outliers which can be removed from the overall scattering data and then replaced by statistically valid fill values. This method, described as a 'K-space algorithmic reconstruction' (KAREN), can identify Bragg reflections independent of prior knowledge of a system's unit cell. KAREN does not alter any data other than that in the immediate vicinity of the Bragg reflections, and reconstructs the diffuse component surrounding the Bragg peaks without introducing discontinuities which induce Fourier ripples or artifacts from underfilling `punched' voids. The KAREN algorithm for reconstructing diffuse scattering provides demonstrably better resolution than can be obtained from previously described punch-and-fill methods. The superior structural resolution obtained using the KAREN method is demonstrated by evaluating the complex ordered diffuse scattering observed from the neutron diffraction of a single plastic crystal of CBr 4 using pair distribution function analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2023

At-wavelength characterization of X-ray wavefronts in Bragg diffraction from crystals

The advent of next-generation synchrotron radiation sources and X-ray free-electron lasers calls for high-quality Bragg-diffraction crystal optics to preserve the X-ray beam coherence and wavefront. This requirement brings new challenges in characterizing crystals in Bragg diffraction in terms of Bragg-plane height errors and wavefront phase distortions. Here, a quantitative methodology to characterize crystal optics using a state-of-the-art at-wavelength wavefront sensing technique and statistical analysis is proposed. The method was tested at the 1-BM-B optics testing beamline at the Advanced Photon Source for measuring silicon and diamond crystals in a self-referencing single-crystal mode and an absolute double-crystal mode. The phase error sensitivity of the technique is demonstrated to be at the λ/100 level required by most applications, such as the characterization of diamond crystals for cavity-based X-ray free-electron lasers.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

PATHA: Performance Analysis Tool for HPC Applications

Large science projects rely on complex workflows to analyze terabytes or petabytes of data. These jobs are often running over thousands of CPU cores and simultaneously performing data accesses, data movements, and computation. It is difficult to identify bottlenecks or to debug the performance issues in these large workflows. In order to address these challenges, we have developed Performance Analysis Tool for HPC Applications (PATHA) using the state-of-art open source big data processing tools. Our framework can ingest system logs to extract key performance measures, and apply the most sophisticated statistical tools and data mining methods on the performance data. Furthermore, it utilizes an efficient data processing engine to allow users to interactively analyze a large amount of different types of logs and measurements. To illustrate the functionality of PATHA, we conduct a case study on the workflows from an astronomy project known as the Palomar Transient Factory (PTF). This study processed 1.6 TB of system logs collected on the NERSC supercomputer Edison. Using PATHA, we were able to identify performance bottlenecks, which reside in three tasks of PTF workflow with the dependency on the density of celestial objects.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2022

Crash Risks Evaluation of Urban Expressways: A Case Study in Shanghai

We report that proactive traffic safety management systems can reduce crashes by identifying crash precursors, evaluating real-time crash risks, and implementing suitable interventions. The basic prerequisite for developing such a system is to propose a reliable crash risk evaluation model that takes real-time traffic flow data as input. Previous studies have primarily focused on real-time crash prediction using some statistical or machine-learning methods. However, further quantitative evaluation and classification of crash risks have been ignored. In this study, we conduct a systematic crash risk evaluation workflow, including crash risk prediction, crash risk quantification, and crash risk classification. Specifically, the crash risk prediction using an extended logit model is proposed, from which CAS, CSD, UAS, DAS, DTV are identified to be contributing factors of crash risks. Then a crash risk quantification model based on the parameter evaluation of the extended logit model is developed. The crash risks of urban expressways and their spatial-temporal evolution trends are quantified. Finally, the crash risks are classified into high crash risk level, moderate crash risk level, and low crash risk level by the k-means cluster algorithm. Then the threshold boundaries of different crash risk levels are determined. The research results provide a proactive guidance for traffic safety management of urban expressways.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Optimization of a Cyclone Using Multiphase Flow Computational Fluid Dynamics

The U.S. Department of Energy National Energy Technology Laboratory's (NETL) 50 kW th chemical looping reactor (CLR) has an underperforming cyclone, which was designed using empirical correlations. To improve the performance of this cyclone using computational fluid dynamics (CFD)-based modeling simulations, four critical design parameters including the vortex tube radius and length, barrel radius, and the inlet width and height were optimized. Here, NETL's open source multiphase flow with interphase exchange (MFiX) CFD code has been used to model a series of cyclones by systematically varying the geometric design parameters. To perform the optimization process, the surrogate modeling and sensitivity analysis followed by the optimization capability in nodeworks was used. The basic methodology for the process is to employ a statistical design of experiments (DOE) method to generate sampling simulations that fill the design space. Corresponding CFD models are then created, executed, and postprocessed. A response surface is created to characterize the relationship between input parameters and the quantities of interest (QoI). Finally, the CFD-surrogate is used by an optimization method to find the optimal design condition based on the objective and constraints prescribed. The resulting optimal cyclone has a larger diameter and longer vortex tube, a larger diameter barrel, and a taller and narrower solids inlet. The improved design has a predicted pressure drop 11 times lower than the original design while reducing the mass loss by a factor of 2.3.

42 ENGINEERING↗