Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistically designed experiments”

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 127 records · Page 7

Comparing different generations of idealized solar geoengineering simulations in the Geoengineering Model Intercomparison Project (GeoMIP)

Abstract. Solar geoengineering has been receiving increased attention in recent years as a potential temporary solution to offset global warming. One method of approximating global-scale solar geoengineering in climate models is via solar reduction experiments. Two generations of models in the Geoengineering Model Intercomparison Project (GeoMIP) have now simulated offsetting a quadrupling of the CO2 concentration with solar reduction. This simulation is idealized and designed to elicit large responses in the models. Here, we show that energetics, temperature, and hydrological cycle changes in this experiment are statistically indistinguishable between the two ensembles. Of the variables analyzed here, the only major differences involve highly parameterized and uncertain processes, such as cloud forcing or terrestrial net primary productivity. We conclude that despite numerous structural differences and uncertainties in models over the past two generations of models, including an increase in climate sensitivity in the latest generation of models, the models are consistent in their aggregate climate response to global solar dimming.

54 ENVIRONMENTAL SCIENCES↗

High throughput battery failure experimental platform

Addressing the need to increase the sample set to understand the causes of lithium battery thermal runaway, we conceived of an experimental platform with capability to increase the number of runaway experiments (currently 3-5 per week), while also collecting detailed electrochemical impedance spectroscopy measurements (EIS). Once expanded, the platform would enable data collection on 10s to 100s of cells that all experience runaway, thereby creating a statistical database necessary to identify early indication of risk. A primary containment unit to house cylindrical cells of variety NMC811 and of size 21700 (21 mm diameter by 70 mm length) was designed with features such as debris containment, preloaded cells in an exchangeable port, nitrogen ventilation, and exhaust containment. We performed the first overcharge abuse experiments of several 21700 cells, handpicked because of different initial EIS, and demonstrated that EIS changes dramatically during early stages of overcharge, but in a different manner than previous pouch cell experiments. The abuse experiments also revealed the discharge pattern and debris field created during runaway, as well as the cell temperature control and overheat, that must be considered in the primary containment apparatus. We designed and built a switching relay system to permit measurement of EIS without an active charging circuit, and created instrument control software for charging, EIS, and cell temperature control. The late-start funding was insufficient to fully construct the primary containment unit, but the foundational design and knowhow is available for any future work.

25 ENERGY STORAGE↗

Continuous and fast calibration of the CMS experiment: design of the automated workflows and operational experience

The exploitation of the full physics potential of the LHC experiments requires fast and efficient processing of the largest possible dataset with the most refined understanding of the detector conditions. To face this challenge, the CMS collaboration has setup an infrastructure for the continuous unattended computation of the alignment and calibration constants, allowing for a refined knowledge of the most time-critical parameters already a few hours after the data have been saved to disk. This is the prompt calibration framework which, since the beginning of the LHC Run-I, enables the analysis and the High Level Trigger of the experiment to consume the most up-to-date conditions optimizing the performance of the physics objects. In the Run-II this setup has been further expanded to include even more complex calibration algorithms requiring higher statistics to reach the needed precision. This imposed the introduction of a new paradigm in the creation of the calibration datasets for unattended workflows and opened the door to a further step in performance. The paper reviews the design of these automated calibration workflows, the operational experience in the Run-II and the monitoring infrastructure developed to ensure the reliability of the service.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Preliminary Design of Critical Experiments Involving Commercially Available B 4 C Neutron Absorber Plates with Low-Enriched UO 2 Fuel

The International Criticality Safety Benchmark Evaluation Project (ICSBEP) is an initiative to provide high-quality benchmark data in a standardized format for criticality safety analysts to use to validate calculational tools and nuclear data. In the last released 2022 version, the ICSBEP handbook contains over 5,000 critical experiment descriptions, results, and associated models. This paper describes the critical configurations proposed and calculations performed as part of the conceptual design phase of Integral Experiment Request (IER) 554 (IER-554). The goal of IER-554 is to add to the benchmark experiments available to the criticality safety community by designing critical experiments that can analyze how adding commercially available Boralcan plates to an assembly of low-enriched UO2 fuel rods affects the effective neutron multiplication factor (k eff ). Boralcan, a neutron absorber product developed by Rio Tinto, is made of B4C in an Al 1100 matrix and is commonly used for criticality suppression in fuel storage pools. In the current version of the ICSBEP handbook, only a few dated critical experiments involve B4C materials, and none involve Boralcan in a thin plate shape. The experiments designed as part of IER-554 are intended to be performed at the Sandia Pulsed Reactor Facility/Critical Experiments (SPRF/CX) apparatus at Sandia National Laboratories. SPRF/CX is a well-characterized assembly considered trustworthy by the benchmarking community because of the numerous high-quality evaluations with very low experimental uncertainties included in the ICSBEP handbook. Figure 1 illustrates the assembly for a particular configuration from LCT-078, one of the published benchmarks in the ICSBEP handbook. The experiments will be moderated and reflected by light water at ambient atmospheric pressure. Before the plates are inserted in the critical assembly, they will be characterized with x-ray computed tomography (XCT) to identify the sizes and distribution of the B 4 C powder particles inside the plates. The preliminary design calculations were performed with SCALE 6.3.0 using the KENO V.a sequence for the criticality calculations, the CE-TSUNAMI-3D sequence for the sensitivity and uncertainty studies, and the ENDF/B-VII.1 Continuous Energy cross section library. All the calculation results presented have 10 pcm statistical uncertainties. A set of critical experiments with Boralcan plates could increase the confidence of the criticality safety community in its modeling methods when using this type of neutron absorber material. Experiments would also help validate k eff calculations, and the industry could use these validations to change the B loading credit limits from the US Nuclear Regulatory Commission standard review plan for dry cask storage of spent nuclear fuel. This paper summarizes only part of the analysis documented in the preliminary design report.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Quinoa, Version 0.2

Quinoa is a set of computational tools that enables research and numerical analysis in fluid dynamics. At this time it remains a test-bed to experiment with various algorithms using fully asynchronous runtime systems. Currently, Quinoa consists of the following tools: (1) Walker, a numerical integrator for systems of stochastic differential equations in time. It is a mathematical tool to analyze and design the behavior of stochastic differential equations. It allows the estimation of arbitrary coupled statistics and probability density functions and is currently used for the design of statistical moment approximations for multiple mixing materials in variable-density turbulence. (2) Inciter, an overdecomposition-aware finite element field solver for partial differential equations using 3D unstructured grids. Inciter is used to research asynchronous mesh-based algorithms and to experiment with coupling asynchronous to bulk-synchronous parallel code. Two planned new features of Inciter, compared to the previous release (LA-CC-16-015), to be implemented in 2017, are (a) a simple Navier-Stokes solver for ideal single-material compressible gases, and (b) solution-adaptive mesh refinement (AMR), which enables dynamically concentrating compute resources to regions with interesting physics. Using the NS-AMR problem we plan to explore how to scale such high-load-imbalance simulations, representative of large production multiphysics codes, to very large problems on very large computers using an asynchronous runtime system. (3) RNGTest, a test harness to subject random number generators to stringent statistical tests enabling quantitative ranking with respect to their quality and computational cost. (4) UnitTest, a unit test harness, running hundreds of tests per second, capable of testing serial, synchronous, and asynchronous functions. (5) MeshConv, a mesh file converter that can be used to convert 3D tetrahedron meshes from and to either of the following formats: Gmsh, http://www.geuz.org/gmsh, Netgen, http://sourceforge.net/apps/mediawiki/netgen-mesher, ExodusII, http://sourceforge.net/projects/exodusii, HyperMesh, http://www.altairhyperworks.com/product/HyperMesh.

Bakosi, Jozsef↗

Correlated materials design: prospects and challenges

The design of correlated materials challenges researchers to combine the maturing, high throughput framework of DFT-based materials design with the rapidly-developing first-principles theory for correlated electron systems. We review the field of correlated materials, distinguishing two broad classes of correlation effects, static and dynamics, and describe methodologies to take them into account. We introduce a material design workflow, and illustrate it via examples in several materials classes, including superconductors, charge ordering materials and systems near an electronically driven metal to insulator transition, highlighting the interplay between theory and experiment with a view towards finding new materials. We review the statistical formulation of the errors of currently available methods to estimate formation energies. We formulate an approach for estimating a lower-bound for the probability of a new compound to form. Correlation effects have to be considered in all the material design steps. Furthermore, these include bridging between structure and property, obtaining the correct structure and predicting material stability. We introduce a post-processing strategy to take them into account.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

On the importance of phase saturation heterogeneity in the analysis of laboratory studies of hydrate dissociation

Methane hydrates (MHs) have been considered as the future source of energy because of its vast resource volume and high energy density. Energy recovery from MH-bearing sediments has attracted intensifying research activities. Fundamentally, heat transfer, fluid flow through porous media, and the kinetics of hydrate reaction are the three key processes controlling the behavior of MH dissociation and the associated fluid production. Earlier studies have suggested that heterogeneous spatial distribution of SH is inevitable in MH-bearing samples synthesized in laboratory. In this paper, we extend our study to analyze numerically the simulation results from the two realizations of the samples (homogeneous and heterogeneous) to identify differences in the fluid production and to determine if they are sufficiently different. Additionally, we conduct a sensitivity analysis and a statistical analysis on the key transport and kinetic rate parameters that could affect hydrate dissociation and fluid production in the context of a heterogeneous hydrate-bearing sample, in an effort to provide insights that could lead to improved designs for laboratory experiments and (possibly) field applications. Our results suggest that the approximation of an artificial hydrate-bearing core with heterogeneous phase saturations by an assumption of uniform phase saturation distributions results in practically similar fluid production profile except for the very early stage with maximum 20.0% deviation in the water production. From the sensitivity and statistical analysis, we determine that gas production depends strongly on the kinetic rate constant, K d 0 and the composite thermal conductivity of the hydrate-bearing sediments, λ θ ; while, water production is very sensitive to K d 0 and the absolute permeability of the sandy medium, k. Understanding the effect of phase heterogeneity and the relative importance of key parameters on the production behavior of hydrate-bearing sediments could provide basis for novel production technologies that lead to enhanced gas production and energy efficiency in the energy recovery process.

42 ENGINEERING↗

An Efficient Pipeline for Abdomen Segmentation in CT Images

Computed tomography (CT) scans usually include some disadvantages due to the nature of the imaging procedure, and these handicaps prevent accurate abdomen segmentation. Discontinuous abdomen edges, bed section of CT, patient information, closeness between the edges of the abdomen and CT, poor contrast, and a narrow histogram can be regarded as the most important handicaps that occur in abdominal CT scans. Currently, one or more handicaps can arise and prevent technicians obtaining abdomen images through simple segmentation techniques. In other words, CT scans can include the bed section of CT, a patient’s diagnostic information, low-quality abdomen edges, low-level contrast, and narrow histogram, all in one scan. These phenomena constitute a challenge, and an efficient pipeline that is unaffected by handicaps is required. In addition, analysis such as segmentation, feature selection, and classification has meaning for a real-time diagnosis system in cases where the abdomen section is directly used with a specific size. A statistical pipeline is designed in this study that is unaffected by the handicaps mentioned above. Intensity-based approaches, morphological processes, and histogram-based procedures are utilized to design an efficient structure. Performance evaluation is realized in experiments on 58 CT images (16 training, 16 test, and 26 validation) that include the abdomen and one or more disadvantage(s). The first part of the data (16 training images) is used to detect the pipeline’s optimum parameters, while the second and third parts are utilized to evaluate and to confirm the segmentation performance. The segmentation results are presented as the means of six performance metrics. Thus, the proposed method achieves remarkable average rates for training/test/validation of 98.95/99.36/99.57% (jaccard), 99.47/99.67/99.79% (dice), 100/99.91/99.91% (sensitivity), 98.47/99.23/99.85% (specificity), 99.38/99.63/99.87% (classification accuracy), and 98.98/99.45/99.66% (precision). In summary, a statistical pipeline performing the task of abdomen segmentation is achieved that is not affected by the disadvantages, and the most detailed abdomen segmentation study is performed for the use before organ and tumor segmentation, feature extraction, and classification.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

MC performance study for an early $v_μ$ charged-current inclusive analysis with MicroBooNE

This note describes an analysis performed on Monte Carlo data to evaluate the sensitivity of MicroBooNE for an early ν µ charged-current inclusive cross section measurement. Such an analysis is intended to be done using the first three months of Booster Neutrino Beam data. The event selection is entirely based on an automated event reconstruction. The Monte Carlo prediction for a flux-integrated and single differential cross section measurement with an approximate estimation of statistical and systematic uncertainties for the MicroBooNE detector as designed is presented. This allows the comparison of the sensitivity of MicroBooNE to theory and other experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Formal Experimentation and Analysis of Handheld RFID Readers as a Tool for Nuclear Material Accounting

Commercial, off-the-shelf Radio Frequency Identification (RFID) systems have been successfully deployed for inventory tracking in numerous industries, but their viability in the tracking of complex environments containing nuclear material is less understood. Of primary interest in this setting is a RFID tracking system for nuclear material accounting which can reliably identify as many tags as possible with high accuracy, while also reducing an operator’s exposure to radiation. In this work, we develop a formal statistical approach to identify relevant handheld RFID reader settings which optimize tagging performance. To achieve this goal, we design a full factorial split-plot experiment for a static shelf configuration scene with 50 randomly placed nuclear material containers affixed with RFID tags. We use Bayesian inference to fit a second-order response surface model which expresses the probability of a successful match for each container as a function of the experimental factors. Such effects are allowed to vary across individual containers and the containers’ population in its entirety to estimate overall effects. Uncertainties of estimates and predictions are quantified via their corresponding posterior distributions. Following extensive model checking and validation, the fitted model is used to identify experimental factors which maximize matching probabilities at both the container-level and for the full shelf configuration scene. We also analyze sensitivity of performance to relevant factors.

25 ENERGY STORAGE↗

Plackett-Burman experimental design to facilitate syntactic foam development

The use of an eight-experiment Plackett–Burman method can assess six experimental variables and eight responses in a polysiloxane-glass microsphere syntactic foam. The approach aims to decrease the time required to develop a tunable polymer composite by identifying a reduced set of variables and responses suitable for future predictive modeling. The statistical design assesses the main effects of mixing process parameters, polymer matrix composition, microsphere density and volume loading, and the blending of two grades of microspheres, using a dummy factor in statistical calculations. Responses cover rheological, physical, thermal, and mechanical properties. The cure accelerator content of the polymer matrix and the volume loading of the microspheres have the largest effects on foam properties. These factors are the most suitable for controlling the gel point of the curing foam, and the density of the cured foam. The mixing parameters introduce widespread variability and therefore should be fixed at effective levels during follow-up testing. Some responses may require greater contrast in microsphere-related factors. As a result, compared to other possible statistical approaches, the run economy of the Plackett–Burman method makes it a valuable tool for rapidly characterizing new foams.

36 MATERIALS SCIENCE↗

Making inertial confinement fusion models more predictive

Computer models of inertial confinement fusion (ICF) implosions play an essential role in experimental design and interpretation as well as our understanding of fundamental physics under the most extreme conditions that can be reached in the laboratory. Building truly predictive models is a significant challenge, with the potential to greatly accelerate progress to high yield and ignition. One path to more predictive models is to use experimental data to update the underlying physics in a way that can be extrapolated to new experiments and regimes. We describe a statistical framework for the calibration of ICF simulations using data collected at the National Ignition Facility (NIF). We perform Bayesian inferences for a series of laser shots using an approach that is designed to respect the physics simulation as much as possible and then build a second model that links the individual-shot inferences together. We show that this approach is able to match multiple X-ray and neutron diagnostics for a whole series of NIF “BigFoot” shots. Within the context of 2D radiation hydrodynamic simulations, our inference strongly favors a significant reduction in fuel compression over other known degradation mechanisms (namely, hohlraum issues and engineering perturbations). This analysis is expanded using a multifidelity technique to pick fuel-ablator mix from several candidate causes of the degraded fuel compression (including X-ray preheat and shock timing errors). Finally, we use our globally calibrated model to investigate the extra laser drive energy that would be required to overcome the inferred fuel compression issues in NIF BigFoot implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Discontinuities in effective permeability due to fracture percolation

Motivated by a triaxial coreflood experiment with a sample of Utica shale where an abrupt jump in permeability was observed, possibly due to the creation of a percolating fracture network through the sample, we perform numerical simulations based on the experiment to characterize how the effective permeability of otherwise low-permeability porous media depends on fracture formation, connectivity, and the contrast between the fracture and matrix permeabilities. While a change in effective permeability due to fracture formation is expected, the dependence of its magnitude upon the contrast between the matrix permeability and fracture permeability and the fracture network structure is poorly characterized. We use two different high-fidelity fracture network models to characterize how effective permeability changes as percolation occurs. The first is a dynamic two-dimensional fracture propagation model designed to mimic the laboratory settings of the experiment. The second is a static three-dimensional discrete fracture network (DFN) model, whose fracture and network statistics are based on the fractured sample of Utica shale. Once the network connects the inflow and outflow boundaries, the effective permeability increases non-linearly with network density. In most networks considered, a jump in the effective permeability was observed when the embedded fracture network percolated. We characterize how the magnitude of the jump, should it occur, depends on the contrast between the fracture and matrix permeabilities. For small contrasts between the matrix and fracture permeabilities the change is insignificant. However, for larger contrasts, there is a substantial jump whose magnitude depends non-linearly on the difference between matrix and fracture permeabilities. A power-law relationship between the size of the jump and the difference between the matrix and fracture permeabilities is observed. In conclusion, the presented results underscore the importance of fracture network topology on the upscaled properties of the porous medium in which it is embedded.

36 MATERIALS SCIENCE↗

Data mining graphene: Correlative analysis of structure and electronic degrees of freedom in graphenic monolayers with defects

The link between changes in the material crystal structure and its mechanical, electronic, magnetic, and optical functionalities known as the structure-property relationship is the cornerstone of the contemporary materials science research. The recent advances in scanning transmission electron and scanning probe microscopies (STEM and SPM) have opened an unprecedented path towards examining the materials structure property relationships on the single-impurity and atomic-configuration levels. Lacking, however, are the statistics-based approaches for cross-correlation of structure and property variables obtained in different information channels of the STEM and SPM experiments. Here we have designed an approach based on a combination of sliding window Fast Fourier Transform, Pearson correlation matrix, linear and kernel canonical correlation, to study a relationship between lattice distortions and electron scattering from the SPM data on graphene with defects. Our analysis revealed that the strength of coupling to strain is altered between different scattering channels which can explain coexistence of several quasiparticle interference patterns in the nanoscale regions of interest. In addition, the application of the kernel functions allowed us extracting a non-linear component of the relationship between the lattice strain and scattering intensity in graphene. Lastly, the outlined approach can be further utilized to analyzing correlations in various multi-modal imaging techniques where the information of interest is spatially distributed and has usually a complex multidimensional nature.

36 MATERIALS SCIENCE↗

Simulation-based optimal Bayesian experimental design for nonlinear systems

The optimal selection of experimental conditions is essential to maximizing the value of data for inference and prediction, particularly in situations where experiments are time-consuming and expensive to conduct. We propose a general mathematical framework and an algorithmic approach for optimal experimental design with nonlinear simulation-based models; in particular, we focus on finding sets of experiments that provide the most information about targeted sets of parameters. Our framework employs a Bayesian statistical setting, which provides a foundation for inference from noisy, indirect, and incomplete data, and a natural mechanism for incorporating heterogeneous sources of information. An objective function is constructed from information theoretic measures, reflecting expected information gain from proposed combinations of experiments. Polynomial chaos approximations and a two-stage Monte Carlo sampling method are used to evaluate the expected information gain. Stochastic approximation algorithms are then used to make optimization feasible in computationally intensive and high-dimensional settings. These algorithms are demonstrated on model problems and on nonlinear parameter inference problems arising in detailed combustion kinetics.

97 MATHEMATICS AND COMPUTING↗

The Implementation of a Super Mirror Polarizer at the SNS Fundamental Neutron Physics Beamline

A new bender supermirror polarizer is used to polarize the cold neutron beam at the Fundamental Neutron Physics Beamline at the Spallation Neutron Source. We present the design of a compensation magnet that was built around the polarizer to minimize the polarizer fringe fields that could compromise the magnetic field requirements of the NPDGamma experiment for the field uniformity in the spin rotator and the field direction in the liquid hydrogen target located downstream from the polarizer. The entire magnetic field environment of the experiment has been analyzed using a finite-element model. Measurements of the magnetic field gradients and field direction have been carried out and the results are less than the upper limits required in the experiment. According to the results the compensated fields meet the stringent magnetic field requirements of the experiment defined by the systematic errors that have to be well below the statistical uncertainty of 10 8 in our main observable, the gamma asymmetry in neutron capture on hydrogen. We describe the design of the magnetic field, the construction of the compensation magnet, and we compare results of the field measurements with the results from the model.

Balascuta, S.↗

Probing Flavor Asymmetry of Anti-quarks in the Proton by Drell-Yan Experiment SeaQuest

A new measurement on the avor asymmetry between d and u in the proton is reported in this thesis. The proton contains a substantial number of antiquarks which arise from dynamical interactions of gluons such as gluon dissociation to a quark-antiquark pair, g ! q + q, and from non-perturbative processes as described by the pion-cloud model, for example. The antiquarks in the proton undertake an important role in determining the dynamic characteristics of the internal structure of the proton, although its distribution in the proton and its origin are not fully understood. Understanding sea quarks in hadron is an important subject for QCD. The SeaQuest experiment at Fermi National Accelerator Laboratory (Fermilab) is a xed target experiment using the 120 GeV proton beam extracted from the Fermilab Main Injector. One of the goals of the experiment is to measure the avor asymmetry between d quark and u quark in the proton as a function of the target Bjorken x using the Drell-Yan process in the p-p or p-d reactions. This process takes place in hadron-hadron collisions when a quark in one hadron in the beam and an antiquark in other hadron in the target annihilate into a virtual photon that decays into a lepton pair. The avor asymmetry between d and u quarks was found by deep-inelastic scattering experiment NMC at CERN. The E866/NuSea experiment at Fermilab obtained the avor asymmetry in the proton for 0:015 < x < 0:35 using the 800 GeV proton beam extracted from the Fermilab Tevatron. The result indicates the dominance of d; it is 70% larger than u at lower x. The SeaQuest experiment was planned to do a new precise measurement at higher x region. The lower energy beam (120 GeV) increases the Drell-Yan cross section and suppresses the background primarily arising from J/ decays. Therefore, SeaQuest will obtain more statistics in a shorter time than the E866 experiment. After detector construction, detector commissioning and accelerator upgrade, physics data taking started in 2013. The SeaQuest spectrometer is designed to detect dimuon from the Drell-Yan process. It consists of targets, two di-pole magnets, and four tracking detector groups. The third tracking detector group has two drift chambers. One was newly fabricated in Japan by the Japanese group in SeaQuest collaboration and was shipped to Fermilab. The other one was constructed by SeaQuest collaborator in Fermilab under the initiative of the Japanese group. I worked on the construction and installation of the detectors, data taking and data analysis in SeaQuest. I extracted the avor asymmetry as a function of Bjorken x using the SeaQuest data for the rst time. This thesis shows the results using a part of data taken in 2014 and 2015. The asymmetry was extracted for much wider Bjorken x region than the previous experiment. The measured Bjorken x range covers up to 0.58. The result shows that the ratio of d=u is always higher than 1 at 0:1 < x < 0:45, in contrast to the E866 result. For 0:45 < x < 0:58, the result shows that the ratio is close to unity. Predictions made by current PDF parameterizations are in agreement with the present result. Also, a prediction obtained by one of the non-perturbative models, pion-cloud model, is closer to the SeaQuest result than the E866 result. This result of d=u asymmetry at the wide Bjorken x region, 0:1 < x < 0:58, is very important information to understand the inner structure of the proton and the origin of the sea quarks in the proton.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗