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,729 records · Page 96

Multivariate Analysis of the Influence of Process Parameters on Suspension Plasma-Sprayed Coatings

Suspension plasma spray (SPS) is far more complicated than conventional plasma spray and requires a deep knowledge about the influence of process parameters and their correlations. In this study, YSZ coatings were manufactured by SPS with six different process parameters such as plasma power, suspension mass load, original powder size, substrate surface topology, spray distance, and spray step. Afterwards, the porosity of as-prepared coatings was investigated by image method and X-ray transmission technique. A multivariate analysis on the collected experimental data was carried out by employing mathematical statistics methods. The results showed that: 1) Coating porosity has a negative correlation with plasma power and suspension mass load and a positive correlation with the original powder size, spray distance, spray step, and substrate roughness; 2) Spraying distance is the main factor affecting to coating porosity, followed by suspension mass load and substrate surface roughness, respectively. A linear model for porosity prediction was developed and was verified by experiments. The mechanism by which process parameters influence coating porosity is also discussed.

Zhao, Yongli↗

Strong Lensing Parameter Estimation on Ground-Based Imaging Data Using Simulation-Based Inference

Current ground-based cosmological surveys, such as the Dark Energy Survey (DES), are predicted to discover thousands of galaxy-scale strong lenses, while future surveys, such as the Vera Rubin Observatory Legacy Survey of Space and Time (LSST) will increase that number by 1-2 orders of magnitude. The large number of strong lenses discoverable in future surveys will make strong lensing a highly competitive and complementary cosmic probe. To leverage the increased statistical power of the lenses that will be discovered through upcoming surveys, automated lens analysis techniques are necessary. We present two Simulation-Based Inference (SBI) approaches for lens parameter estimation of galaxy-galaxy lenses. We demonstrate the successful application of Neural Posterior Estimation (NPE) to automate the inference of a 12-parameter lens mass model for DES-like ground-based imaging data. We compare our NPE constraints to a Bayesian Neural Network (BNN) and find that it outperforms the BNN, producing posterior distributions that are for the most part both more accurate and more precise; in particular, several source-light model parameters are systematically biased in the BNN implementation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Search For Excited Cascade Hypersons (¿*-) Using the CLAS12 Spectrometer at Jefferson Laboratory

The number of experimentally observed doubly strange Cascade states to date is far fewer than has been predicted theoretically. The CLAS12 Very Strange physics program (E12- 11-005A) at the Thomas Jefferson National Accelerator Facility aims to study the electroproduction of these states using the newly upgraded CEBAF Large Acceptance Spectrometer for 12 GeV (CLAS12) in experimental Hall B. In this project, the reaction ep ! e0K+K+??? ! e0K+K+K?(?=?0) is studied using CLAS12 Run Group A (RG-A) data sets taken by impinging electron beams of 10:2 and 10:6 GeV energies on an LH2 target. Scattered electrons are detected with either the Forward Detector (FD), covering a polar angle range of 5? to 35? to study electroproduction, or with the Forward Tagger (FT), covering a polar angle range of 2:5? to 4:5?, to study quasi-real photoproduction. The CLAS12 detector with nearly a 4? solid angle coverage is used to detect scattered electrons and charged kaons in the ?final state. ?=?0 hyperons are reconstructed using the missing mass technique to explore intermediate doubly-strange hyperons (???) that decay to K? and ?=?0. No statistically significant ??? states other than the ???(1530) were found in the missing mass spectra based on the currently available statistics in the CLAS12-FD acceptance only. A maximum log-likelihood method is implemented to determine the statistical significance and the upper limits on the ???(1820) electroproduction cross section by constructing 95% confidence level boundaries on the ???(1820) yields. Additional research was conducted to investigate the upper limit electroproduction cross section of the reaction ep ! e0K+K+??? ! e0K+K+K?(?=?0) as a function of the electroproduced ??? mass and the differential cross section ( d?dMM(e0K+K+) ) as a function of missing mass MM(e0K+K+) for electroproduction and quasi-real photoproduction processes.

Khanal, Achyut↗

State densities of heavy nuclei in the static-path plus random-phase approximation

We report that nuclear state densities are important inputs to statistical models of compound-nucleus reactions. State densities are often calculated with self-consistent mean-field approximations that do not include important correlations and must be augmented with empirical collective enhancement factors. Here, we benchmark the static-path plus random-phase approximation (SPA + RPA) to the state density in a chain of samarium isotopes 148–155 Sm against exact results (up to statistical errors) obtained with the shell-model Monte Carlo (SMMC) method. The SPA + RPA method incorporates all static fluctuations beyond the mean field together with small-amplitude quantal fluctuations around each static fluctuation. Using a pairing plus quadrupole interaction, we show that the SPA + RPA state densities agree well with the exact SMMC densities for both the even- and odd-mass isotopes. For the even-mass isotopes, we also compare our results with mean-field state densities calculated with the finite-temperature Hartree-Fock-Bogoliubov (HFB) approximation. We find that the SPA + RPA repairs the deficiencies of the mean-field approximation associated with broken rotational symmetry in deformed nuclei and with the violation of particle-number conservation in the pairing condensate. In particular, in deformed nuclei the SPA + RPA reproduces the rotational enhancement of the state density relative to the mean-field state density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

High–Yield Deterministic Focused Ion Beam Implantation of Quantum Defects Enabled by In Situ Photoluminescence Feedback

Focused ion beam implantation is ideally suited for placing defect centers in wide bandgap semiconductors with nanometer spatial resolution. However, the fact that only a few percent of implanted defects can be activated to become efficient single photon emitters prevents this powerful capability to reach its full potential in photonic/electronic integration of quantum defects. Here an industry adaptive scalable technique is demonstrated to deterministically create single defects in commercial grade silicon carbide by performing repeated low ion number implantation and in situ photoluminescence evaluation after each round of implantation. An array of 9 single defects in 13 targeted locations is successfully created—a ≈70% yield which is more than an order of magnitude higher than achieved in a typical single pass ion implantation. The remaining emitters exhibit non-classical photon emission statistics corresponding to the existence of at most two emitters. This approach can be further integrated with other advanced techniques such as in situ annealing and cryogenic operations to extend to other material platforms for various quantum information technologies.

36 MATERIALS SCIENCE↗

Development of a Time-Height Profile Assimilation Technique for Large-Eddy Simulation

Mesoscale-to-microscale coupling (MMC) aims to address the limited scope of traditional large-eddy simulations by driving the microscale flow with information concerning large-scale weather patterns provided by mesoscale models. We present a new offline MMC technique for horizontally homogeneous microscale flow conditions, in which internal forcing terms are computed based on mesoscale time–height profiles of mean-flow quantities. The advantage of such an approach is that it can be used to drive a microscale simulation with either mesoscale or observational data, and that it does not rely on specific terms in the mesoscale budget equations, which are typically not part of the default output of a mesoscale solver. The performance of the proposed profile assimilation technique is assessed based on the simulation of a typical diurnal cycle over the Scaled Wind Farm Technology site in west Texas. Results indicate that simple data assimilation techniques lead to unphysically high levels of shear and turbulence caused by the algorithm’s inability to cope with inaccuracies in the mesoscale time–height profiles. Modifying the algorithm to account for vertical coherence in the mesoscale source terms gives the microscale solver a greater ability to correct the provided mesoscale time–height profiles, leading to improved predictions of shear and turbulence statistics. The resulting turbulence statistics are in good agreement with meteorological tower observations and simulation results obtained with state-of-the-art coupling techniques using mesoscale budget components.

data assimilation↗

Statistical approach to nuclear level density

We discuss the level density in a finite many-body system with strong interaction between the constituents. Our primary object of applications is the atomic nucleus but the same techniques can be applied to other mesoscopic systems. We calculate and compare nuclear level densities for given quantum numbers obtained by different methods, such as nuclear shell model (the most successful microscopic approach), our main instrument - moments method (statistical approach), and Fermi-gas model; the calculation with the moments method can use any shell-model Hamiltonian excluding the spurious states of the center-of-mass motion. Our goal is to investigate statistical properties of nuclear level density, define its phenomenological parameters, and offer an affordable and reliable way of calculation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

UNSUPERVISED TRANSIENT LIGHT CURVE ANALYSIS VIA HIERARCHICAL BAYESIAN INFERENCE

Historically, light curve studies of supernovae (SNe) and other transient classes have focused on individual objects with copious and high signal-to-noise observations. In the nascent era of wide field transient searches, objects with detailed observations are decreasing as a fraction of the overall known SN population, and this strategy sacrifices the majority of the information contained in the data about the underlying population of transients. A population level modeling approach, simultaneously fitting all available observations of objects in a transient sub-class of interest, fully mines the data to infer the properties of the population and avoids certain systematic biases. We present a novel hierarchical Bayesian statistical model for population level modeling of transient light curves, and discuss its implementation using an efficient Hamiltonian Monte Carlo technique. As a test case, we apply this model to the Type IIP SN sample from the Pan-STARRS1 Medium Deep Survey, consisting of 18,837 photometric observations of 76 SNe, corresponding to a joint posterior distribution with 9176 parameters under our model. Our hierarchical model fits provide improved constraints on light curve parameters relevant to the physical properties of their progenitor stars relative to modeling individual light curves alone. Moreover, we directly evaluate the probability for occurrence rates of unseen light curve characteristics from the model hyperparameters, addressing observational biases in survey methodology. We view this modeling framework as an unsupervised machine learning technique with the ability to maximize scientific returns from data to be collected by future wide field transient searches like LSST.

79 ASTRONOMY AND ASTROPHYSICS↗

The effect of sensor sheltering and averaging techniques on wind measurements at the Shuttle Landing Facility

This document presents results of a field study of the effect of sheltering of wind sensors by nearby foliage on the validity of wind measurements at the Space Shuttle Landing Facility (SLF). Standard measurements are made at one second intervals from 30-feet (9.1-m) towers located 500 feet (152 m) from the SLF centerline. The centerline winds are not exactly the same as those measured by the towers. A companion study, Merceret (1995), quantifies the differences as a function of statistics of the observed winds and distance between the measurements and points of interest. This work examines the effect of nearby foliage on the accuracy of the measurements made by any one sensor, and the effects of averaging on interpretation of the measurements. The field program used logarithmically spaced portable wind towers to measure wind speed and direction over a range of conditions as a function of distance from the obstructing foliage. Appropriate statistics were computed. The results suggest that accurate measurements require foliage be cut back to OFCM standards. Analysis of averaging techniques showed that there is no significant difference between vector and scalar averages. Longer averaging periods reduce measurement error but do not otherwise change the measurement in reasonably steady flow regimes. In rapidly changing conditions, shorter averaging periods may be required to capture trends.

Merceret, Francis J.↗

Radiation dose rate distributions of spent fuel dry casks estimated with MAVRIC based on detailed geometry and continuous-energy models

This study presents a detailed comparison of the dose rate distributions of a dry TN–32 fuel cask with two geometry models and two cross-sectional datasets. The accuracies of radiation dose rate estimation and computational efficiencies of each geometry model with two cross-sectional data-sets are compared. The use of automated variance reduction techniques can significantly improve the computational efficiency of such a realistic, deep penetration problem that involves radiation transport from different volumetric sources, thereby eliciting only a small statistical error. Monaco with Automated Variance Reduction using Importance Calculations (MAVRIC) is a computational sequence within the SCALE 6.2 code package based on consistent adjoint driven importance sampling (CADIS), a type of automated variance reduction technique. Homogenous and full fuel assembly models are built herein, and two nuclear cross-section libraries (V7–200N47G and continuous energy) are applied in this work. Based on the detailed comparisons, we found that neutron dose rate estimation is more dependent on geometry modeling than on cross-section data. For neutron-induced gamma rays, the dose rate distribution depends on both the spatial self-shielding effect and the cross-section library. The primary gamma rays respectively contribute to the total dose rate by ~91% and ~ 99% on side and top surfaces, and the dose rate accuracy is more dependent on the cross-section library than on geometry modeling. In terms of the computation efficiency and efforts spent on geometry modeling, the homogenous fuel assembly model with the MG library can produce an acceptable dose rate distribution, but the detailed fuel assembly model with the continuous-energy library is required for more precise dose rate estimation.

61 RADIATION PROTECTION AND DOSIMETRY↗

Evaluating the Benefits of Bayesian Hierarchical Methods for Analyzing Heterogeneous Environmental Datasets: A Case Study of Marine Organic Carbon Fluxes

Large compilations of heterogeneous environmental observations are increasingly available as public databases, allowing researchers to test hypotheses across datasets. Statistical complexities arise when analyzing compiled data due to unbalanced spatial sampling, variable environmental context, mixed measurement techniques, and other reasons. Hierarchical Bayesian modeling is increasingly used in environmental science to describe these complexities, however few studies explicitly compare the utility of hierarchical Bayesian models to simpler and more commonly applied methods. Here we demonstrate the utility of the hierarchical Bayesian approach with application to a large compiled environmental dataset consisting of 5,741 marine vertical organic carbon flux observations from 407 sampling locations spanning eight biomes across the global ocean. We fit a global scale Bayesian hierarchical model that describes the vertical profile of organic carbon flux with depth. Profile parameters within a particular biome are assumed to share a common deviation from the global mean profile. Individual station-level parameters are then modeled as deviations from the common biome-level profile. The hierarchical approach is shown to have several benefits over simpler and more common data aggregation methods. First, the hierarchical approach avoids statistical complexities introduced due to unbalanced sampling and allows for flexible incorporation of spatial heterogeneitites in model parameters. Second, the hierarchical approach uses the whole dataset simultaneously to fit the model parameters which shares information across datasets and reduces the uncertainty up to 95% in individual profiles. Third, the Bayesian approach incorporates prior scientific information about model parameters; for example, the non-negativity of chemical concentrations or mass-balance, which we apply here. We explicitly quantify each of these properties in turn. We emphasize the generality of the hierarchical Bayesian approach for diverse environmental applications and its increasing feasibility for large datasets due to recent developments in Markov Chain Monte Carlo algorithms and easy-to-use high-level software implementations.

54 ENVIRONMENTAL SCIENCES↗

Approximate Bayesian computation for forward modeling in cosmology

Bayesian inference is often used in cosmology and astrophysics to derive constraints on model parameters from observations. This approach relies on the ability to compute the likelihood of the data given a choice of model parameters. In many practical situations, the likelihood function may however be unavailable or intractable due to non-gaussian errors, non-linear measurements processes, or complex data formats such as catalogs and maps. In these cases, the simulation of mock data sets can often be made through forward modeling. We discuss how Approximate Bayesian Computation (ABC) can be used in these cases to derive an approximation to the posterior constraints using simulated data sets. This technique relies on the sampling of the parameter set, a distance metric to quantify the difference between the observation and the simulations and summary statistics to compress the information in the data. We first review the principles of ABC and discuss its implementation using a Population Monte-Carlo (PMC) algorithm and the Mahalanobis distance metric. We test the performance of the implementation using a Gaussian toy model. We then apply the ABC technique to the practical case of the calibration of image simulations for wide field cosmological surveys. We find that the ABC analysis is able to provide reliable parameter constraints for this problem and is therefore a promising technique for other applications in cosmology and astrophysics. Our implementation of the ABC PMC method is made available via a public code release.

79 ASTRONOMY AND ASTROPHYSICS↗

Advanced Models for Nondestructive Evaluation of Aging Nuclear Power Plant Cables

Background: Nuclear power plant (NPP) cables must function correctly for proper power supply and control of the nuclear reactor. This is especially important for public and environmental safety at critical times such as during reactor shut-down following loss of coolant. Over extended periods of service, the cable insulation and jacket polymer materials suffer from degradation due to exposure to environmental stressors such as heat and ionizing radiation. The focus of this project has been the development of detailed descriptions of the way in which insulation polymers degrade, from a materials science perspective, and the development of new electrical techniques to monitor cable polymer condition nondestructively. Objectives: The objectives of this project have been to i) perform detailed characterization experiments upon cable insulation polymers that have experienced thermal and radiation exposure mimicking that received during NPP operation, including the exploration of novel cable nondestructive evaluation methods; capacitive, THz and infrared methods; ii) identify the most sensitive indicators of aging in these materials, and iii) develop advanced, validated models describing microstructural and chemical changes that are observed in cable insulation polymers due to thermal and radiation exposure. Methods: Accelerated aging was performed on the two most common United States NPP cable insulation polymers (cross-linked polyethylene (XLPE) and ethylene propylene rubber (EPR)) under various conditions of elevated temperature and gamma radiation dose. Aging mechanisms were elucidated through extensive thermal, mechanical and microscopic characterization of the aged samples, and materials modeling. Dielectric and IR spectra were measured on pristine and aged samples. Data was analyzed statistically, via multiple linear regression, partial least squares and probabilistic neural network approaches. Impact: Extant cable monitoring techniques do not deal successfully with all polymer types. This research has laid a foundation upon which new electrical nondestructive test methods for cable polymer condition monitoring may be developed, hence launching the development of new test methods for cable insulation monitoring based on measuring electrical rather than mechanical parameters. Collaboration: This project has benefitted from collaboration between experts in materials science and engineering, materials modeling, electrical engineering, and nondestructive evaluation. The majority of the characterization measurements were carried out at Iowa State University. Materials modeling was carried out primarily at Washington State University. Irradiation of samples took place at Pacific Northwest National Laboratory (PNNL), where a student intern from Iowa State was hosted for one summer, setting up the aging experiments. A research student from the University of Bologna, Italy, joined the project for their thesis research at Iowa State University with samples provided by PNNL. The project benefitted from regular involvement of senior engineers in US-based industry partners. Project investigators presented findings at various national and international conferences and workshops.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Investigation of energy management strategies for photovoltaic systems - An analysis technique

Progress is reported in formulating energy management strategies for stand-alone PV systems, developing an analytical tool that can be used to investigate these strategies, applying this tool to determine the proper control algorithms and control variables (controller inputs and outputs) for a range of applications, and quantifying the relative performance and economics when compared to systems that do not apply energy management. The analysis technique developed may be broadly applied to a variety of systems to determine the most appropriate energy management strategies, control variables and algorithms. The only inputs required are statistical distributions for stochastic energy inputs and outputs of the system and the system's device characteristics (efficiency and ratings). Although the formulation was originally driven by stand-alone PV system needs, the techniques are also applicable to hybrid and grid connected systems.

Cull, R. C.↗

Galactic ArchaeoLogIcaL ExcavatiOns (GALILEO) II. t-SNE portrait of local fossil relics and structures

Based on high-quality Apache Point Observatory Galactic Evolution Experiment (APOGEE) DR17 and Gaia DR3 data for 1742 red giants stars within 5 kpc of the Sun and not rotating with the Galactic disk (V φ < 100 km s -1 ), we used the nonlinear technique of unsupervised analysis t-Distributed Stochastic Neighbor Embedding (t-SNE) to detect coherent structures in the space of ten chemical-abundance ratios: [Fe/H], [O/Fe], [Mg/Fe], [Si/Fe], [Ca/Fe], [C/Fe], [N/Fe], [Al/Fe], [Mn/Fe], and [Ni/Fe]. Additionally, we obtained orbital parameters for each star using the nonaxisymmetric gravitational potential GravPot16. Seven structures are detected, including Splash, Gaia-Sausage-Enceladus (GSE), the high-α heated-disk population, N-C-O peculiar stars, and inner disk-like stars, plus two other groups that did not match anything previously reported in the literature, here named Galileo 5 and Galileo 6 (G5 and G6). These two groups overlap with Splash in [Fe/H], with G5 having a lower metallicity than G6, and they are both between GSE and Splash in the [Mg/Mn] versus [Al/Fe] plane, with G5 being in the α-rich in situ locus and G6 on the border of the α-poor in situ one. Nonetheless, their low [Ni/Fe] hints at a possible ex situ origin. Their orbital energy distributions are between Splash and GSE, with G5 being slightly more energetic than G6. We verified the robustness of all the obtained groups by exploring a large range of t-SNE parameters, applying it to various subsets of data, and also measuring the effect of abundance errors through Monte Carlo tests.

79 ASTRONOMY AND ASTROPHYSICS↗

Final Scientific/Technical Report: Integrating Advanced Software and Statistical Methods for Spectroscopic Dark Energy Surveys

“Dark energy” is a phenomenon of unknown origin that pervades the vacuum of space and is causing the present-day expansion of the universe to accelerate. Understanding the nature of dark energy is one of the most fundamental challenges to modern physics, having the potential to revolutionize our theories of both gravity and elementary particle physics. The effects of dark energy are currently only accessible to astronomical survey experiments that probe large volumes of the cosmos. One of the most powerful forms of experiment for studying dark energy is the cosmological “redshift survey,” which maps the large-scale structure of the universe as traced by the positions of galaxies, quasars, and intergalactic hydrogen clouds. The signature of dark energy and other cosmological parameters is imprinted on this large-scale structure as it evolves over cosmic time. This project focused on maximizing the accurate extraction of information from redshift surveys by improving the algorithms and statistics applied across the chain of inference from raw pixels to high-level cosmological parameter estimation. At the low level, this project developed new techniques for calibrating and analyzing high-background, low signal-to-noise data from astronomical survey spectrographs to provide reliable three-dimensional maps of galaxies and quasars throughout the universe. At the high level, this project developed new statistical frameworks for understanding the completeness of redshift surveys and the underlying demographics of the galaxies that they observe, thereby allowing for a more accurate connection between survey results and theoretical models. This work contributes to maximizing the scientific reach of existing and planned surveys, and opens new potential feasibility for more ambitious surveys in the future.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Study to define unsteady flow fields and their statistical characteristics

Preliminary estimates of space shuttle fluctuating pressure environments were made based on analyses of wind tunnel data, and empirical prediction techniques. Particular emphasis was given to the external tank and solid rocket boosters for the transonic speed regime during launch of a parallel-burn space shuttle configuration. Predicted environments are presented as space-averaged zonal profiles with progressive shading from zone to zone to illustrate spatial variations in the magnitude of the fluctuating pressure coefficient over the surfaces of the external tank and solid rocket boosters. Predictions are provided for the transonic Mach number range from 0.8 equal to or less than M sub infinity equal to or less than 1.5, and for supersonic Mach numbers of 2.0 and 3.0.

Robertson, J. E.↗