Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical techniques”

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,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Fluctuating Hydrodynamic Methods for Fluid-Structure Interactions in Confined Channel Geometries

Computational methods are developed for the study of fluid-structure interactions subject to thermal fluctuations when conned within channels having slit-like geometry. The methods take into account the hydrodynamic coupling and diusivity of microstructures when influenced by their proximity to no-slip walls. We develop stochastic numerical methods subject to no-slip boundary conditions using a staggered nite volume discretization. We show that by imposing an exact fluctuation-dissipation condition to discretize the stochastic driving elds combined with using an exact projection method to enforce incompressibility is sucient to ensure results consistent with statistical mechanics. We demonstrate our methods by investigating how the proximity of ellipsoidal colloids to the channel wall eects their active hydrodynamic responses and passive diusivity. We also study the collective dynamics of a large number of particles by considering the intermediate scattering functions for the relaxation of density fluctuations. We expect our introduced stochastic computational methods to be applicable broadly to applications in which confinement effects play an important role in the dynamics of microstructures subject to hydrodynamic coupling and thermal fluctuations. uctuations.

Wang, Yaohong↗
DOE OSTI2025

Ensemble Simulation Techniques and Fast Randomized Algorithms

The major goals of the project were to develop and analyze new ensemble simulation techniques, including trajectory stratification and preconditioned MCMC techniques, as well as develop fast numerical linear algebra techniques closely related to ensemble simulation ideas. The trajectory stratification techniques involve simulating in parallel short trajectory fragments of a Markov process confined to a specific region of space‐time and then patching together the statistics gathered to assemble estimates of very general dynamical properties. We have also developed this approach for rare event simulation and extended the techniques to applications requiring a more general framework (such as electronic structure calculations). The preconditioned MCMC techniques involve simulating multiple Markov chains in parallel and then using information from the ensemble to speed the mixing of each individual chain. The fast randomized linear algebra methods are motivated by the diffusion Monte Carlo technique, but are applicable to finding the dominant eigenvalue of (almost) general matrices. For most non‐negative matrices, the schemes result in an error (compared to the power method) that is constant in the dimension of the problem. For more general matrices, we see a very clear sublinear cost trend in computational tests.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1976

A demonstration of a transportable radio interferometric surveying system with 3-cm accuracy on a 307-m base line

A precision geodetic measurement system (Aries, for Astronomical Radio Interferometric Earth Surveying) based on the technique of very long base line interferometry has been designed and implemented through the use of a 9-m transportable antenna and the NASA 64-m antenna of the Deep Space Communications Complex at Goldstone, California. A series of experiments designed to demonstrate the inherent accuracy of a transportable interferometer was performed on a 307-m base line during the period from December 1973 to June 1974. This short base line was chosen in order to obtain a comparison with a conventional survey with a few-centimeter accuracy and to minimize Aries errors due to transmission media effects, source locations, and earth orientation parameters. The base-line vector derived from a weighted average of the measurements, representing approximately 24 h of data, possessed a formal uncertainty of about 3 cm in all components. This average interferometry base-line vector was in good agreement with the conventional survey vector within the statistical range allowed by the combined uncertainties (3-4 cm) of the two techniques.

Ong, K. M.↗
DOE OSTI2014

Comparative quantification and statistical analysis of η′ and η precipitates in aluminum alloy AA7075-T651 by TEM and AFM

Quantification of nanometric precipitates in metallic alloys has been traditionally performed using transmission electron microscopy, which is nominally a low throughput technique. This work presents a comparative study of quantification of η′ and η precipitates in aluminum alloy AA7075-T651 using transmission electron microscopy (TEM) and non-contact atomic force microscopy (AFM). AFM quantification was compared with 2-D stereological results reported elsewhere. Also, a method was developed, using specialized software, to characterize nanometric size precipitates observed in dark-field TEM micrographs. Statistical analysis of the quantification results from both measurement techniques supports the use of AFM for precipitate characterization. Once the precipitate stoichiometry has been determined by appropriate analytical techniques like TEM, as it is the case for η′ and η in AA7075-T651, the relative ease with which specimens are prepared for AFM analysis could be advantageous in product and process development, and quality control, where a large number of samples are expected for analysis on a regular basis. - Highlights: • Nanometric MgZn{sub 2} precipitates in AA7075-T651 were characterized using AFM and TEM. • Phase-contrast AFM was used to differentiate metal matrix from MgZn{sub 2} precipitates. • TEM and AFM micrographs were analyzed using commercially available software. • AFM image analysis and TEM 2-D stereology render statistically equivalent results.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Closing the scale gap between land surface parameterizations and GCMs with a new scheme, SiB3-Bins: SOIL MOISTURE SCALE GAP

The interaction of land with the atmosphere is sensitive to soil moisture (W). Evapotranspiration (ET) reacts to soil moisture in a nonlinear way, f(W), as soils dry from saturation to wilt point. This nonlinear behavior and the fact that soil moisture varies on scales as small as 1–10 m in nature, while numerical general circulation models (GCMs) have grid cell sizes on the order of 1 to 100s of kilometers, makes the calculation of grid cell-average ET problematic. It is impractical to simulate the land in GCMs on the small scales seen in nature, so techniques have been developed to represent subgrid scale heterogeneity, including: (1) statistical-dynamical representations of grid subelements of varying wetness, (2) relaxation of f(W), (3) moderating f(W) with approximations of catchment hydrology, (4) “tiling” the landscape into vegetation types, and (5) hyperresolution. Here we present an alternative method for representing subgrid variability in W, one proven in a conceptual framework where landscape-scale W is represented as a series of “Bins” of increasing wetness from dry to saturated. The grid cell-level f(W) is defined by the integral of the fractional area of the wetness bins and the value of f(W) associated with each. This approach accounts for the spatiotemporal dynamics of W. We implemented this approach in the SiB3 land surface parameterization and then evaluated its performance against a control, which assumes a horizontally uniform field of W. We demonstrate that the Bins method, with a physical basis, attenuates unrealistic jumps in model state and ET seen in the control runs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

On Bayes Equality and Related Issues

To make calculations in the Bayesian analysis, the formalism of which is based on the layering of a probability measure defined on the product of measurable spaces, it is useful to have a summary of the properties of this layering. In this paper we formulate and prove those of them that are used in calculations more often than others. Particularly, we prove Hoeffding-type inequalities using direct elementary techniques.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Pilot Point Optimization of Mining Boundaries for Lateritic Metal Deposits: Finding the Trade-off Between Dilution and Ore Loss

Geological contacts in lateritic metal deposits (footwall topographies) often delineate the orebody boundaries. Spatial variations seen in such contacts are frequently higher than those for the metal grades of the deposit. Therefore, borehole spacing chosen based on the grade variations cannot adequately capture the geological contact variability. Consequently, models created using such boreholes cause high volumetric uncertainties in the actual and targeted ore extraction volumes, which, in turn, lead to high unplanned dilution and ore losses. In this paper, a method to design optimum ore/mining boundaries for lateritic metal deposits is presented. The proposed approach minimizes the dilution/ore losses and comprises two main steps. First, the uncertainty on the orebody boundary is represented using a set of stochastic realizations generated with a multiple-point statistics algorithm. Then, the optimal orebody boundary is determined using an optimization technique inspired by a model calibration method called Pilot Points. The pilot points represent synthetic elevation values, and they are used to construct smooth mining boundaries using the multilevel B-spline technique. The performance of a generated surface is evaluated using the expected sum of losses in each of the stochastic realizations. The simulated annealing algorithm is used to iteratively determine the pilot point values which minimize the expected losses. The results show a significant reduction in the dilution volume as compared to those obtained from the actual mining operation.

58 GEOSCIENCES↗
DOE OSTI2017

Influences of granular constraints and surface effects on the heterogeneity of elastic, superelastic, and plastic responses of polycrystalline shape memory alloys

Deformation heterogeneities at the microstructural length-scale developed in polycrystalline shape memory alloys (SMAs) during superelastic loading are studied using both experiments and simulations. In situ X-ray diffraction, specifically the far-field high energy diffraction microscopy (ff-HEDM) technique, was used to non-destructively measure the grain-averaged statistics of position, crystal orientation, elastic strain tensor, and volume for hundreds of austenite grains in a superelastically loaded nickel-titanium (NiTi) SMA. These experimental data were also used to create a synthetic microstructure within a finite element model. The development of intragranular stresses were then simulated during tensile loading of the model using anisotropic elasticity. Driving forces for phase transformation and slip were calculated from these stresses. The grain-average responses of individual austenite crystals examined before and after multiple stress-induced transformation events showed that grains in the specimen interior carry more axial stress than the surface grains as the superelastic response "shakes down". Examination of the heterogeneity within individual grains showed that regions near grain boundaries exhibit larger stress variation compared to the grain interiors. As a result, this intragranular heterogeneity is more strongly driven by the constraints of neighboring grains than the initial stress state and orientation of the individual grains.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Non-targeted discovery-based analysis for gas chromatography with mass spectrometry: A comparison of peak table, tile, and pixel-based Fisher ratio analysis

The ability to discover minute differences between samples or sample classes for gas chromatography coupled to mass spectrometry (GC-MS) can be a challenging endeavor, especially when those differences are not a priori . Fisher ratio (F-ratio) analysis is an apt technique to probe the differences between GC-MS chromatograms. F-ratio analysis is a supervised, non-targeted, discovery-based method that compares two different samples (or sample classes) to reduce the GC-MS dataset into a hit list composed of class distinguishing compounds. Three different F-ratio techniques, peak table, tile, and pixel-based were used to “discover” nine non-native analytes that were spiked into gasoline at four different nominal concentrations of 250, 85, 25, 5 parts-per-million (ppm). For the tile and pixel-based F-ratio calculations, a novel methodology is introduced to improve the sensitivity of the F-ratio calculations while reducing false positives. Furthermore, we use a combinatorial technique using null class comparisons, termed null distribution analysis, to determine a statistical F-ratio cutoff for analysis of the hit lists. The pixel-based algorithm was the most sensitive method and was able to “discover” all nine spiked analytes at a nominal concentration of 250 ppm albeit with one false positive interspersed towards the bottom of the hit list. The pixel-based software was also able to “discover” more of the spiked analytes at the lower concentrations with seven of the spiked analytes “discovered” at 85 ppm, four of the spiked analytes “discovered” at 25 ppm, and one analyte “discovered” at 5 ppm.

47 OTHER INSTRUMENTATION↗
DOE OSTI2017

Improved molecular level identification of organic compounds using comprehensive two-dimensional chromatography, dual ionization energies and high resolution mass spectrometry

A new analytical methodology combining comprehensive two-dimensional gas chromatography (GC×GC), dual ionization energies and high resolution time of flight mass spectrometry has been developed that improves molecular level identification of organic compounds in complex environmental samples. GC×GC maximizes compound separation providing cleaner mass spectra by minimizing erroneous fragments from interferences and co-eluting peaks. Traditional electron ionization (EI, 70 eV) provides MS fragmentation patterns that can be matched to published EI MS libraries while vacuum ultraviolet photoionization (VUV, 10.5 eV) yields MS with reduced fragmentation enhancing the abundance of the molecular ion providing molecular formulas when combined with high resolution mass spectrometry. We demonstrate this new approach by applying it to a sample of organic aerosol. In this sample, 238 peaks were matched to EI MS library data with FM ≥ 800 but a fifth (42 compounds) were determined to be incorrectly identified because the molecular formula was not confirmed by the VUV MS data. This highlights the importance of using a complementary technique to confirm compound identifications even for peaks with very good matching statistics. In total, 171 compounds were identified by EI MS matching to library spectra with confirmation of the molecular formula from the high resolution VUV MS data and were not dependent on the matching statistics being above a threshold value. A large number of unidentified peaks were still observed with FM < 800, which in routine analysis would typically be neglected. Where possible, these peaks were assigned molecular formulas from the VUV MS data (211 in total). In total, the combination of EI and VUV MS data provides more than twice as much molecular level peak information than traditional approaches and improves confidence in the identification of individual organic compounds. The molecular formula data from the VUV MS data was used, in conjunction with GC×GC retention times and the observed EI MS, to generate a new, searchable EI MS library compatible with the standard NIST MS search program. This library is deliberately dynamic and editable so that other end users can add new entries and update existing entries as new information becomes available.A new analytical methodology has been developed to improve molecular level identification of organic compounds in complex samples.

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

SDSS/SEGUE spectral feature analysis for stellar atmospheric parameter estimation

Large-scale and deep sky survey missions are rapidly collecting a large amount of stellar spectra, which necessitate the estimation of atmospheric parameters directly from spectra and make it feasible to statistically investigate latent principles in a large data set. We present a technique for estimating parameters T{sub eff}, log g, and [Fe/H] from stellar spectra. With this technique, we first extract features from stellar spectra using the LASSO algorithm; then, the parameters are estimated from the extracted features using the support vector regression. On a subsample of 20,000 stellar spectra from the Sloan Digital Sky Survey (SDSS) with reference parameters provided by the SDSS/SEGUE Spectroscopic Parameter Pipeline, estimation consistency are 0.007458 dex for log T{sub eff} (101.609921 K for T{sub eff}), 0.189557 dex for log g, and 0.182060 for [Fe/H], where the consistency is evaluated by mean absolute error. Prominent characteristics of the proposed scheme are sparseness, locality, and physical interpretability. In this work, each spectrum consists of 3821 fluxes, and 10, 19, and 14 typical wavelength positions are detected, respectively, for estimating T{sub eff}, log g, and [Fe/H]. It is shown that the positions are related to typical lines of stellar spectra. This characteristic is important in investigating physical indications from analysis results. Then, stellar spectra can be described by the individual fluxes on the detected positions (PD) or local integration of fluxes near them (LI). The aforementioned consistency is the result based on features described by LI. If features are described by PD, consistency is 0.009092 dex for log T{sub eff} (124.545075 K for T{sub eff}), 0.198928 dex for log g, and 0.206814 dex for [Fe/H].

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Reconstructing cosmic polarization rotation with ResUNet-CMB

Cosmic polarization rotation, which may result from parity-violating new physics or the presence of primordial magnetic fields, converts E-mode polarization of the cosmic microwave background (CMB) into B-mode polarization. Anisotropic cosmic polarization rotation leads to statistical anisotropy in CMB polarization and can be reconstructed with quadratic estimator techniques similar to those designed for gravitational lensing of the CMB. At the sensitivity of upcoming CMB surveys, lensing-induced B-mode polarization will act as a limiting factor in the search for anisotropic cosmic polarization rotation, meaning that an analysis which incorporates some form of delensing will be required to improve constraints on the effect with future surveys. In this study we extend the ResUNet-CMB convolutional neural network to reconstruct anisotropic cosmic polarization rotation in the presence of gravitational lensing and patchy reionization, and we show that the network simultaneously reconstructs all three effects with variance that is lower than that from the standard quadratic estimator nearly matching the performance of an iterative reconstruction method.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Decoding the age–chemical structure of the Milky Way disc: an application of copulas and elicitable maps

In the Milky Way, the distribution of stars in the [α/Fe] versus [Fe/H] and [Fe/H] versus age planes holds essential information about the history of star formation, accretion, and dynamical evolution of the Galactic disc. We investigate these planes by applying novel statistical methods called copulas and elicitable maps to the ages and abundances of red giants in the Apache Point Observatory Galactic Evolution Experiment survey. We find that the high- and low-α disc stars have a clean separation in copula space and use this to provide an automated separation of the α sequences using a purely statistical approach. This separation reveals that the high-α disc ends at the same [α/Fe] and age at high [Fe/H] as the low-[Fe/H] start of the low-α disc, thus supporting a sequential formation scenario for the high- and low-α discs. We then combine copulas with elicitable maps to precisely obtain the correlation between stellar age τ and metallicity [Fe/H] conditional on Galactocentric radius R and height z in the range 0 < R < 20 kpc and |z| < 2 kpc. The resulting trends in the age–metallicity correlation with radius, height, and [α/Fe] demonstrate a ≈0 correlation wherever kinematically cold orbits dominate, while the naively expected negative correlation is present where kinematically hot orbits dominate. This is consistent with the effects of spiral-driven radial migration, which must be strong enough to completely flatten the age–metallicity structure of the low-α disc.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Iterative Monte Carlo analysis of spin-dependent parton distributions

We present a comprehensive new global QCD analysis of polarized inclusive deep-inelastic scattering, including the latest high-precision data on longitudinal and transverse polarization asymmetries from Jefferson Lab and elsewhere. The analysis is performed using a new iterative Monte Carlo fitting technique which generates stable fits to polarized parton distribution functions (PDFs) with statistically rigorous uncertainties. Inclusion of the Jefferson Lab data leads to a reduction in the PDF errors for the valence and sea quarks, as well as in the gluon polarization uncertainty at x ≳ 0.1. Furthermore, the study also provides the first determination of the flavor-separated twist-3 PDFs and the d 2 moment of the nucleon within a global PDF analysis.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Evaluating Synthetic Smart Meter Locations For Communication System Modeling

An important part of good communication models for smart grid applications, particularly when wireless protocols are used, is the location information of the nodes. We’ve developed metrics to assist in evaluating whether a given model with such information is representative of feeders found in the world. These metrics are applied to a set of feeder models with known good location information for the smart meter locations and compared to a set of feeder models with similar location information that is not expected to be representative of real-world feeders. The comparison reveals that the suspect models do not pass statistical tests utilizing the developed metrics, providing initial validation of the analysis technique.

smart meter, feeder model, GridLAB-D↗
DOE OSTI2014

SU-E-J-85: Leave-One-Out Perturbation (LOOP) Fitting Algorithm for Absolute Dose Film Calibration

Purpose: To introduce an outliers-recognition fitting routine for film dosimetry. It cannot only be flexible with any linear and non-linear regression but also can provide information for the minimal number of sampling points, critical sampling distributions and evaluating analytical functions for absolute film-dose calibration. Methods: The technique, leave-one-out (LOO) cross validation, is often used for statistical analyses on model performance. We used LOO analyses with perturbed bootstrap fitting called leave-one-out perturbation (LOOP) for film-dose calibration . Given a threshold, the LOO process detects unfit points (“outliers”) compared to other cohorts, and a bootstrap fitting process follows to seek any possibilities of using perturbations for further improvement. After that outliers were reconfirmed by a traditional t-test statistics and eliminated, then another LOOP feedback resulted in the final. An over-sampled film-dose- calibration dataset was collected as a reference (dose range: 0-800cGy), and various simulated conditions for outliers and sampling distributions were derived from the reference. Comparisons over the various conditions were made, and the performance of fitting functions, polynomial and rational functions, were evaluated. Results: (1) LOOP can prove its sensitive outlier-recognition by its statistical correlation to an exceptional better goodness-of-fit as outliers being left-out. (2) With sufficient statistical information, the LOOP can correct outliers under some low-sampling conditions that other “robust fits”, e.g. Least Absolute Residuals, cannot. (3) Complete cross-validated analyses of LOOP indicate that the function of rational type demonstrates a much superior performance compared to the polynomial. Even with 5 data points including one outlier, using LOOP with rational function can restore more than a 95% value back to its reference values, while the polynomial fitting completely failed under the same conditions. Conclusion: LOOP can cooperate with any fitting routine functioning as a “robust fit”. In addition, it can be set as a benchmark for film-dose calibration fitting performance.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

SU-E-T-270: Diffusion Synthetic Acceleration for Linear Boltzmann Transport Equation

Purpose: Linear Boltzmann transport equation (LBTE) is as accurate as the Monte Carlo method (MC) for dose calculation in photon/particle therapy (LBTE is a deterministic and Eulerian formulation and MC is a statistical and Lagrangian description). An advantage of LBTE is that numerous acceleration techniques can be utilized for acceleration. This work is to explore the acceleration of LBTE via diffusion synthetic acceleration (DSA). Methods: For simplicity, two-dimensional, steady-state, and within-group LBTE is considered with two angular dimensions and two spatial dimensions. The discrete ordinate method is developed for solving this integro-differential equation. The angular variables are discretized using a level-symmetric quadrature set on the unit sphere. The spatial variables are discretized on the structured grid based on the diamond scheme. The source-iteration method (SI) is used to solve the discretized system.Since SI is slow in optically thick and highly scattering regime. DSA is developed to accelerate SI. The motivation for DSA is that diffusion equation (DE) is a good approximation of LBTE in the above regime. However, DE is much cheaper than LBTE computationally since DE only involves spatial variables. Thus, in each DSA iteration, DSA adds to the SI step a computationally-negligible DE step, i.e., to first solve DE with the SI residual as source term, and then compensate the SI solution with DE solution. Results: DSA was benchmarked and compared with SI. The difference between two methods was within 0.12% which verifies the accuracy of DSA, while DSA demonstrated the great advantage in speed, e.g., the reduction of iteration number to 6% and 4% respectively for cases with 100 and 1,000 scattering-absorption ratio that commonly occur in clinical dose calculation. Conclusion: DSA has been developed as one of many possible means for accelerating the numerical solver of LBTE for dose calculation. The authors were partially supported by the NSFC (#11405105), the 973 Program (#2015CB856000) and the Shanghai Pujiang Talent Program (#14PJ1404500)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Differentiable Predictions for Large Scale Structure with SHAMNet

In simulation-based models of the galaxy-halo connection, theoretical predictions for galaxy clustering and lensing are typically made based on Monte Carlo realizations of a mock universe. In this paper, we use Subhalo Abundance Matching (SHAM) as a toy model to introduce an alternative to stochastic predictions based on mock population, demonstrating how to make simulation-based predictions for clustering and lensing that are both exact and differentiable with respect to the parameters of the model. Conventional implementations of SHAM are based on iterative algorithms such as Richardson-Lucy deconvolution; here we use the JAX library for automatic differentiation to train SHAMNet, a neural network that accurately approximates the stellar-to-halo mass relation (SMHM) defined by abundance matching. In our approach to making differentiable predictions of large-scale structure, we map parameterized PDFs onto each simulated halo, and calculate gradients of summary statistics of the galaxy distribution by using autodiff to propagate the gradients of the SMHM through the statistical estimators used to measure one- and two-point functions. Our techniques are quite general, and we conclude with an overview of how they can be applied in tandem with more complex, higher-dimensional models, creating the capability to make differentiable predictions for the multi-wavelength universe of galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗