Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,459 records · Page 81

Statistical chronometry of Meteorites: II. Initial abundances and homogeneity of short-lived radionuclides

Astrophysical models of planet formation require accurate radiometric dating of meteoritic components by short-lived (Al-Mg, Mn-Cr, Hf-W) and long-lived (U-Pb) chronometers, to develop a timeline of such events in the solar nebula as formation of Ca-rich, Al-rich Inclusions (CAIs), chondrules, planetesimals, etc. CAIs formed mostly around a time (“t=0”) when the short-lived radionuclide 26 Al (t 1/2 = 0.72 Myr) was present and presumably homogeneously distributed at a known level we define as ( 26 Al/ 27 Al) SS = 5.23 × 10 -5 . The time of formation after t=0 of another object can be found by determining its initial ( 26 Al/ 27 Al) 0 ratio and comparing it to ( 26 Al/ 27 Al) SS . Dating of meteoritic objects using the Mn-Cr or Hf-W systems is hindered because the abundances ( 53 Mn/ 55 Mn) SS and ( 182 Hf/ 180 Hf) SS at t=0 are not known precisely. To constrain these quantities, we compile literature Al-Mg, Mn-Cr, Hf-W and Pb-Pb data for 13 achondrites and use novel statistical techniques to minimize the discrepancies between their times of formation across these systems. We find that for ( 53 Mn/ 55 Mn) SS = (7.80±0.36)×10 -6 , ( 182 Hf/ 180 Hf) SS = (10.41 ± 0.12) × 10-5, t SS = 4568.65 ± 0.10 Myr, and a 53 Mn half-life of 3.98 ± 0.22 Myr, these four free parameters make concordant 18 formation times recorded by the different systems in all six known volcanic achondrites (the quenched angrites D’Orbigny, Sahara 99555, Northwest Africa 1670; the psuedo-eucrites Asuka 881394 and Ibitira; and the basaltic achondrite Northwest Africa 7325). These parameters also make concordant the ages derived for chondrules from CB/CH achondrites, formed simultaneously in an impact. The other seven achondrites are not quite concordant, but are plutonic angrites or ‘carbonaceous achondrites’ for which simultaneous closure of the isotopic systems might not be expected. Our findings provide very strong support for homogeneity of 26 Al, 53 Mn, and 182 Hf in the solar nebula, and our approach offers a path for more precise chronometry

79 ASTRONOMY AND ASTROPHYSICS↗

Quality Control of High-Dose-Rate Brachytherapy: Treatment Delivery Analysis Using Statistical Process Control

Purpose: Statistical process control (SPC) is a quality control method used to ensure that a process is well controlled and operates with little variation. This study determined whether SPC was a viable technique for evaluating the proper operation of a high-dose-rate (HDR) brachytherapy treatment delivery system. Methods and Materials: A surrogate prostate patient was developed using Vyse ordnance gelatin. A total of 10 metal oxide semiconductor field-effect transistors (MOSFETs) were placed from prostate base to apex. Computed tomography guidance was used to accurately position the first detector in each train at the base. The plan consisted of 12 needles with 129 dwell positions delivering a prescribed peripheral dose of 200 cGy. Sixteen accurate treatment trials were delivered as planned. Subsequently, a number of treatments were delivered with errors introduced, including wrong patient, wrong source calibration, wrong connection sequence, single needle displaced inferiorly 5 mm, and entire implant displaced 2 mm and 4 mm inferiorly. Two process behavior charts (PBC), an individual and a moving range chart, were developed for each dosimeter location. Results: There were 4 false positives resulting from 160 measurements from 16 accurately delivered treatments. For the inaccurately delivered treatments, the PBC indicated that measurements made at the periphery and apex (regions of high-dose gradient) were much more sensitive to treatment delivery errors. All errors introduced were correctly identified by either the individual or the moving range PBC in the apex region. Measurements at the urethra and base were less sensitive to errors. Conclusions: SPC is a viable method for assessing the quality of HDR treatment delivery. Further development is necessary to determine the most effective dose sampling, to ensure reproducible evaluation of treatment delivery accuracy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Statistical Modeling Approach to Quantitative Analysis of Interobserver Variability in Breast Contouring

Purpose: To develop a new approach for interobserver variability analysis. Methods and Materials: Eight radiation oncologists specializing in breast cancer radiation therapy delineated a patient's left breast “from scratch” and from a template that was generated using deformable image registration. Three of the radiation oncologists had previously received training in Radiation Therapy Oncology Group consensus contouring for breast cancer atlas. The simultaneous truth and performance level estimation algorithm was applied to the 8 contours delineated “from scratch” to produce a group consensus contour. Individual Jaccard scores were fitted to a beta distribution model. We also applied this analysis to 2 or more patients, which were contoured by 9 breast radiation oncologists from 8 institutions. Results: The beta distribution model had a mean of 86.2%, standard deviation (SD) of ±5.9%, a skewness of −0.7, and excess kurtosis of 0.55, exemplifying broad interobserver variability. The 3 RTOG-trained physicians had higher agreement scores than average, indicating that their contours were close to the group consensus contour. One physician had high sensitivity but lower specificity than the others, which implies that this physician tended to contour a structure larger than those of the others. Two other physicians had low sensitivity but specificity similar to the others, which implies that they tended to contour a structure smaller than the others. With this information, they could adjust their contouring practice to be more consistent with others if desired. When contouring from the template, the beta distribution model had a mean of 92.3%, SD ± 3.4%, skewness of −0.79, and excess kurtosis of 0.83, which indicated a much better consistency among individual contours. Similar results were obtained for the analysis of 2 additional patients. Conclusions: The proposed statistical approach was able to measure interobserver variability quantitatively and to identify individuals who tended to contour differently from the others. The information could be useful as feedback to improve contouring consistency.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Influence of numerical schemes on statistical properties of computed charged particle trajectories in turbulent electromagnetic fields

A class of numerical schemes is developed for the study of charged particle transport in complex stationary electromagnetic fields and tested for fields obtained from a numerical solution of the magneto-hydrodynamic equation. The performances of these schemes are evaluated by analyzing the conservation of energy and the statistical properties of the trajectories. Energy conservation is affected by the interpolation technique used to estimate the field value at the particle position. However, the particle transport properties are more robust, except in the limit of low energy when a significant fraction of the particles are trapped.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Enhancing statistical moment calculations for stochastic Galerkin solutions with Monte Carlo techniques

Here, we provide a method for enhancing stochastic Galerkin moment calculations to the linear elliptic equation with random diffusivity using an ensemble of Monte Carlo solutions. This hybrid approach combines the accuracy of low-order stochastic Galerkin and the computational efficiency of Monte Carlo methods to provide statistical moment estimates which are significantly more accurate than performing each method individually. The hybrid approach involves computing a low-order stochastic Galerkin solution, after which Monte Carlo techniques are used to estimate the residual. We show that the combined stochastic Galerkin solution and residual is superior in both time and accuracy for a one-dimensional test problem and a more computational intensive two-dimensional linear elliptic problem for both the mean and variance quantities.

97 MATHEMATICS AND COMPUTING↗

A probabilistic drought forecasting framework: A combined dynamical and statistical approach

In order to improve drought forecasting skill, this study develops a probabilistic drought forecasting framework comprised of dynamical and statistical modeling components. The novelty of this study is to seek the use of data assimilation to quantify initial condition uncertainty with the Monte Carlo ensemble members, rather than relying entirely on the hydrologic model or land surface model to generate a single deterministic initial condition, as currently implemented in the operational drought forecasting systems. Next, the initial condition uncertainty is quantified through data assimilation and coupled with a newly developed probabilistic drought forecasting model using a copula function. The initial condition at each forecast start date are sampled from the data assimilation ensembles for forecast initialization. Finally, seasonal drought forecasting products are generated with the updated initial conditions. This study introduces the theory behind the proposed drought forecasting system, with an application in Columbia River Basin, Pacific Northwest, United States. Results from both synthetic and real case studies suggest that the proposed drought forecasting system significantly improves the seasonal drought forecasting skills and can facilitate the state drought preparation and declaration, at least three months before the official state drought declaration.

data assimilation↗

Capturing the fractocohesive length scale in elastomers through a statistical mechanics-based gradient enhanced damage model

Here, this study aims to examine modeling flaw sensitivity in elastomers. The direct incorporation of polymer chain statistical mechanics considerations into a continuum stretch-based gradient-enhanced damage formulation, in turn, allows a representation of diffuse chain damage and fracture events that align with known micromechanical mechanisms. Through a series of numerical experiments, we simulate crack propagation and extract the fracture energy as an output of the model, while keeping track of the micromechanical signatures of diffuse chain damage that accommodate fracture propagation and eventually influence flaw sensitivity. Finally, by combining the fracture toughness and the work to rupture, we identify a fractocohesive length of the material, corresponding to the full width of the damage process zone. As the damage-to-fracture cascade in the proposed GED model is influenced by the introduction of a length scale associated with network imperfection and long-range load transfer, the emerging relationship of the two length scales is discussed, providing a potential link between microscopic damage mechanisms and the observed macroscopic fracture response.

36 MATERIALS SCIENCE↗

Application of the immobilized low-activity waste glass corrosion model to the static dissolution of 24 statistically-designed alkali-borosilicate waste glasses

Glass corrosion models that capture the complex mechanisms of the glass-water reaction enable the prediction of nuclear waste glass durability in disposal scenarios. Parameterization of such models is challenging because of the need to capture changes in corrosion behavior with time, reaction conditions, and glass composition. In this report we describe and employ the immobilized low activity waste (ILAW) glass corrosion model (IGCM) in geochemical simulations of static dissolution tests, at two temperatures (40 °C and 90 °C), for a matrix of 24 enhanced low activity waste (eLAW) glasses statistically designed to cover a processable composition space defined by 8 major glass components (Al 2 O 3 , B 2 O 3 , CaO, Na 2 O, SiO 2 , SnO 2 , ZrO 2 , and Others as defined in the text). The IGCM includes a first-order chemical affinity term and an ion-exchange term that represents the net exchange of Na+ ions in the pristine glass with protons in aqueous solution. Constant, time-dependent, and time- and pH-dependent functional forms of the ion-exchange term are evaluated to reproduce the change in corrosion behavior with time in saturated, static dissolution tests. The agreement with measured aqueous concentrations of the main glass components (B, Na, and Si) improved significantly upon addition of a time-dependent term, which therefore constitutes a simple representation of the glass-water reaction progress. Due to the limited changes in pH in the static dissolution tests, past a short initial period of rapid increase, addition of a pH dependent term did not appreciably improve the fits, indicating that comprehensive model parameterization requires more than one type of glass corrosion test to capture a wide range of solution chemistries. IGCM parameters were found to be dependent on glass composition, and the parameter sets generated in this work will enable the development of composition–parameter correlation models that offer the promise of predicting IGCM parameters, and thus glass corrosion behavior, solely based on glass composition.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Failure statistics for commercial lithium ion batteries: A study of 24 pouch cells

There are relatively few publications that assess capacity decline in enough commercial cells to quantify cell-to-cell variation, but those that do show a surprisingly wide variability. Capacity curves cross each other often, a challenge for efforts to measure the state of health and predict the remaining useful life (RUL) of individual cells. We analyze capacity fade statistics for 24 commercial pouch cells, providing an estimate for the time to 5% failure. Our data indicate that RUL predictions based on remaining capacity or internal resistance are accurate only once the cells have already sorted themselves into “better” and “worse” ones. Analysis of our failure data, using maximum likelihood techniques, provide uniformly good fits for a variety of definitions of failure with normal and with 2- and 3-parameter Weibull probability density functions, but we argue against using a 3-parameter Weibull function for our data. pdf fitting parameters appear to converge after about 15 failures, although business objectives should ultimately determine whether data from a given number of batteries provides sufficient confidence to end lifecycle testing. Increased efforts to make batteries with more consistent lifetimes should lead to improvements in battery cost and safety.

25 ENERGY STORAGE↗

A Bayesian statistical approach to inferring particle dynamics from in-situ pump POC and chloropigment data from the Mediterranean Sea

Concentrations of chloropigments and particulate organic carbon (POC) in large-volume in-situ pump samples from the Mediterranean Sea were used to estimate rate constants of processes that control the fate of particles, and specifically chloropigments, in the water column. In this paper we introduce a Bayesian statistical inversion method that combines the data with a new box model and has the capacity to infer rate constants for POC respiration/dissolution, chlorophyll and pheopigment degradation, and particle aggregation and disaggregation. We use first-order kinetics to model disaggregation, and use both first-order and second-order kinetics to model aggregation. Using these methods, the estimated small-particle (1–70 μm) POC respiration rate constant was 2.44 -1.00 +1.69 yr -1 (0.41 yr). The estimated disaggregation and second-order aggregation rate constants were 85.6 -36.4 +63.4 yr -1 (1.17×10 -2 yr) and 2.78 -1.17 +2.01 μM -1 yr -1 , respectively. Using the optimal rate constants and the corresponding particle concentrations, disaggregation is ~4.2 times faster than the small-size POC dissolution rate, which indicates that disaggregation is a dominant processes at the time of sampling. More importantly, by comparing our results with those of previous studies, we conclude that sampling methods have less influence than tracers themselves on inferring particle dynamic rate constants. We previously introduced a somewhat similar approach to modeling SV sediment trap data, but large volume pumps are a much more common sample collection method in oceanographic surveys than SV sediment traps, and thus our new model should have a wider applicability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Segmentation, statistical analysis, and modelling of the wall system in ceramic foams

Closed walls in otherwise open foam structures may have a great impact on macroscopic properties of the materials. In this paper, we present two algorithms for the segmentation of such closed walls from micro-computed tomography images of the foam structure. The techniques are compared on simulated data and applied to tomographic images of ceramic filters. This allows for a detailed statistical analysis of the normal directions and sizes of the walls. Finally, we explain how the information derived from the segmented wall system can be included in a stochastic microstructure model for the foam.

36 MATERIALS SCIENCE↗

Determination of grain boundary mobility during recrystallization by statistical evaluation of electron backscatter diffraction measurements

One of the key aspects influencing microstructural design pathways in metallic systems is grain boundary motion. The present work introduces a method by means of which direct measurement of grain boundary mobility vs. misorientation dependence is made possible. The technique utilizes datasets acquired by means of serial electron backscatter diffraction (EBSD) measurements. The experimental EBSD measurements are collectively analyzed, whereby datasets were used to obtain grain boundary mobility and grain aspect ratio with respect to grain boundary misorientation. The proposed method is further validated using cellular automata (CA) simulations. Single crystal aluminium was cold rolled and scratched in order to nucleate random orientations. Subsequent annealing at 300 °C resulted in grains growing, in the direction normal to the scratch, into a single deformed orientation. Growth selection was observed, wherein the boundaries with misorientations close to Σ7 CSL orientation relationship (38° 〈111〉) migrated considerably faster. The obtained boundary mobility distribution exhibited a non-monotonic behavior with a maximum corresponding to misorientation of 38° ± 2° about 〈111〉 axes ± 4°, which was 10–100 times higher than the mobility values of random high angle boundaries. Correlation with the grain aspect ratio values indicated a strong growth anisotropy displayed by the fast growing grains. The observations have been discussed in terms of the influence of grain boundary character on grain boundary motion during recrystallization. - Highlights: • Statistical microstructure method to measure grain boundary mobility during recrystallization • Method implementation independent of material or crystal structure • Mobility of the Σ7 boundaries in 5N Al was calculated as 4.7 × 10{sup –8} m{sup 4}/J ⋅ s. • Pronounced growth selection in the recrystallizing nuclei in Al • Boundary mobility values during recrystallization 2–3 orders of magnitude larger vis-à-vis curvature driven motion.

36 MATERIALS SCIENCE↗

Statistical analysis of HAADF-STEM images to determine the surface coverage and distribution of immobilized molecular complexes

The surface immobilization of molecular catalysts is attractive because it combines the benefits of homogeneous and heterogeneous catalysis. However, determining the surface coverage and distribution of a molecular catalyst on a solid support is often challenging, inhibiting our ability to design improved catalytic systems. Here, in this work, we demonstrate that the combination of scanning transmission electron microscopy (STEM) and image analysis of the individual positions of heavy atoms in transition metal complexes via a convolutional neural network (CNN) allows statistically robust determination of the surface coverage and distribution of immobilized molecular catalysts. These observations provide information about how changes in the functionalization conditions, attachment group, and structure of the molecular catalyst affect the surface coverage and distribution, providing insight into the chemical mechanism of surface immobilization. The method could be generally valuable for correlating the surface coverage and distribution to the activity, selectivity, and stability of a catalytic system.

HAADF-STEM↗

Reliable statistics-based detection and investigation of anomalies in a SMART valve system

Reliable anomaly detection and diagnosis are critical for the safe operation of complex engineered systems. This study presents a unified framework that integrates statistical, model-based, and data-driven techniques for anomaly detection and investigation, demonstrated on SMART valve systems in hybrid energy applications. Four detection methods—mean deviation, seasonal extreme studentized deviate, ARIMA forecasting, and matrix profiling—were implemented and compared. Matrix profiling was particularly effective in revealing subtle deviations and hidden relationships among variables. Anomaly investigation was performed by analyzing variable-level and grouped signal profiles, with system topology incorporated to distinguish primary faults from propagated effects. Grouping signals by type enhanced interpretability, enabling accurate localization of anomalies across multi-dimensional datasets. Experimental results confirmed the framework's capability to consistently detect and isolate anomalies while providing actionable insights into system interdependencies. The proposed methodology offers a robust, interpretable, and scalable solution for condition monitoring, with potential applications in safety-critical domains such as nuclear energy, aerospace, and process industries.

ARIMA models↗

Statistical estimation of the performance of a fast-neutron multiplicity system for nuclear material accountancy

Statistical analyses have been performed to develop bounding estimates of the expected performance of a conceptual fast-neutron multiplicity system (FNMS) for assaying plutonium. The conceptual FNMS design includes 32 cubic liquid scintillator detectors, measuring 7.62 cm per side, configured into 4 stacked rings of 8 detectors each. Expected response characteristics for the individual FNMS detectors, as well as the response characteristics of the entire FNMS, were determined using Monte Carlo simulations based on prior validation experiments. The results from these simulations were then used to estimate the Pu assay capabilities of the FNMS in terms of counting time, assay mass, and assay mass variance, using assay mass variance as a figure of merit. The analysis results are compared against a commonly used thermal-neutron coincidence counter. The advantages of using a fast-neutron counting system versus a thermal-neutron counting system are significant. Most notably, the time required to perform an assay to an equivalent assay mass variance is greatly reduced with a fast-neutron system, by more than an order of magnitude compared with that of the thermal-neutron system, due to the reduced probability of random summing with the fast system. The improved FNMS performance is especially relevant for assays involving Pu masses of 10 g or more.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Classification of events from α -induced reactions in the MUSIC detector via statistical and ML methods

The Multi-Sampling Ionization Chamber (MUSIC) detector is typically used to measure nuclear reaction cross sections relevant for nuclear astrophysics, fusion studies, and other applications. From the MUSIC data produced in one experiment scientists carefully extract an order of 10 3 events of interest from about 10 9 total events, where each event can be represented by an 18-dimensional vector. However, the standard data classification process is based on expert driven, manually intensive data analysis techniques that require several months to identify patterns and classify the relevant events from the collected data. Here, to address this issue, we present a method for the classification of events originating from specific α-induced reactions by combining statistical and machine learning methods that require significantly less input from the domain scientist, relative to the standard technique. Here, we applied the new method to two experimental data sets and compared our results with those obtained using traditional methods. With few exceptions, the number of events classified by our method agrees within ±20% with the results obtained using traditional methods. With the present method, which is the first of its kind for the MUSIC data, we have established the foundation for the automated extraction of physical events of interest from experiments using the MUSIC detector.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

High fidelity multiphysics tightly coupled model for a lead cooled fast reactor concept and application to statistical calculation of hot channel factors

A tightly coupled multiphysics code system is established using the MOOSE framework for hot channel factor (HCF) evaluation on a Lead Fast Reactor (LFR) concept. The coupled system is driven by the Griffin multiphysics coupling capability under which the MOOSE Heat Transfer module and NekRS computational fluid dynamics solver are coupled for conjugate heat transfer using the Cardinal application. The coupled capability is demonstrated on an LFR assembly model based on materials and geometry of a prototypical lead-cooled fast reactor design by Westinghouse Electric Company, LLC. Moreover, the work integrates the Multiphysics Object Oriented Simulation Environment (MOOSE) Stochastic Tools Module (STM) to perform calculations for statistical analysis of HCF. Furthermore, the coupling strategy and workflow demonstrated in this paper is not only useful for predicting accurate hot channel factors for different kinds of advanced reactors but also for other engineering applications such as control rod worth assessment, generation of high-fidelity database for Artificial intelligence (AI)/machine learning (ML) training, design optimization and multi-resolution modeling.

Cardinal↗