Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dependability modeling”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Restoration Strategy for Active Distribution Systems Considering Endogenous Uncertainty in Cold Load Pickup

Cold load pickup (CLPU) phenomenon is identified as the persistent power inrush upon a sudden load pickup after an outage. Under the active distribution system (ADS) paradigm, where distributed energy resources (DERs) are extensively installed, the decreased outage duration can induce a strong interdependence between CLPU pattern and load pickup decisions. In this paper, we propose a novel modelling technique to tractably capture the decision-dependent uncertainty (DDU) inherent in the CLPU process. Subsequently, a two-stage stochastic decision-dependent service restoration (SDDSR) model is constructed, where first stage searches for the optimal switching sequences to decide step-wise network topology, and the second stage optimizes the detailed generation schedule of DERs as well as the energization of switchable loads. Further, to tackle the computational burdens introduced by mixed-integer recourse, the progressive hedging algorithm (PHA) is utilized to decompose the original model into scenario-wise subproblems that can be solved in parallel. The numerical test on modified IEEE 123-node test feeders has verified the efficiency of our proposed SDDSR model and provided fresh insights into the monetary and secure values of DDU quantification.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Investigation of grain size and geometrically necessary dislocation density dependence of flow stress in Mg-4Al by using nanoindentation

Models for grain size dependence of flow stress typically ignore contributions from geometrically necessary dislocations (GNDs) in the grain boundary regions. Ashby had proposed the following equation for flow stress (τ) as a function of grain size (d) and effective plastic strain ($\bar{ε}$): $τ = τ_{0} + C'G\sqrt{^{b\bar{ε}}_{d}}$ using an unknown empirical coefficient, C’, In this study, the coefficient, C’, and equation for the grain size dependence of flow stress were investigated in Mg-4Al using nanoindentation. C’ ranges from 1.4 to 7.11 that is related to the hardness difference (∆H0) between the adjacent grains as $C' + 0.62\sqrt{ΔH_{0}}$. Ashby equation was used to predict uniaxial stress-strain curves in grain sizes of 55 µm, 187 µm, and 333 µm. The results showed lower flow stress and lower hardening rate compared to experimentally measured tensile stress-strain response. The discrepancy seems to be a result of hardening effect of SSDs. Here, we propose a modified Ashby equation that takes into account the hardening effect of SSDs and GNDs separately, resulting in a better agreement between calculated and experimentally measured values. This analysis demonstrates estimation of uniaxial stress-strain response using nanoindentation measurements that captures effects of grain size, and strain hardening from SSDs and GNDs. The combination of the revised equation and nanoindentation method will enable rapid exploration of stress-strain response of a diverse range of single-phase polycrystalline alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2023

ENSO-Induced Teleconnection: Process-Oriented Diagnostics to Assess Rossby Wave Sources and Ambient Flow Properties in Climate Models

Climate model fidelity in representing ENSO-induced teleconnection is assessed with process-oriented diagnostics that examine a chain of processes, from equatorial Pacific precipitation to the midlatitude circulation pattern over the Pacific–North American regions. Such processes are rarely addressed during model development. Using an upper-tropospheric divergent level, local vorticity gradient of the ambient zonal flow (∂ 2 U¯/∂y 2 ) and a restoring force for Rossby waves (β * ) are estimated, the equivalent barotropic vorticity equation is solved, and an anomalous Rossby wave source (RWS') quantified. The analysis is applied to AMIP5 and AMIP6 simulations. For a realistic circulation response representation, the hypothesis that models accurately represent the strength and location of RWS', and spatial variations in β * is tested. Compared to AMIP5, in AMIP6 there are clear improvements in representing RWS' and β * . To validate the hypothesis, the analysis identifies two metrics: spatially coherent RWS' in the subtropical North Pacific, and longitudes of negative β * over the western-central North Pacific. By projecting these metrics in two and three-dimensional views, improvements or degradations in model versions are apparent. If a model’s fidelity in representing ∂ 2 U¯/∂y 2 and RWS' are compromised, then radiated Rossby waves are reflected more equatorward, resulting in zonally elongated circulation anomalies over the central North Pacific. Thus, during climate model development, applying this analysis frequently will keep a regular check on the fidelity of the modeled response to anomalous El Niño convection in conjunction with changing model ambient flow dependencies. As a result, this analysis is intended to form a process-oriented diagnostics package, a community contribution to the NOAA Model Diagnostics Task Force.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

An Improved Fracture Mechanics-Informed Multiscale Thermomechanical Damage Model for Ceramic Matrix Composites

This paper extends recent work done by the authors in modeling length scale-dependent damage behavior of ceramic matrix composites (CMCs) to include effects of local anisotropy introduced by matrix cracking. This model captures scale-dependent damage initiation and propagation behavior of the brittle matrix by employing internal state variable (ISV) theory within a multiscale modeling framework to obtain damaged matrix stress/strain constitutive relationships at each length scale. The damage ISV captures the effects of matrix cracking and growth by using fracture mechanics and the self-consistent scheme to determine the reduced stiffness of the cracked matrix. Matrix cracks, which activate when stress intensity factors near manufacturing induced cavities exceed the fracture toughness of the material, are assumed to be transversely isotropic in the plane of the crack, and matrix anisotropy occurs when the damaged stiffness tensor is rotated from the crack plane to the global axes. The crack progression and temporal evolution of the damage ISV are governed by fracture mechanics and crack growth kinetics. The model effectively captures first matrix cracking, which is the first significant deviation from linear elasticity. The nonlinear predictive capabilities of the material model are demonstrated for monolithic silicon carbide (SiC) and a 2D woven five-harness satin (5HS) carbon fiber SiC matrix (C/SiC) CMC.

Skinner, Travis↗
DOE OSTI2020

A new method to construct a horizontal resolution-dependent wind speed adjustment factor for tropical cyclones in global climate model simulations

A new method to construct a horizontal resolution-dependent wind speed adjustment factor for evaluating tropical cyclones (TCs) in global climate models (GCMs) is presented here. In contrast to the previous studies that used idealized axisymmetric wind fields, this study analyzes 48 hours of 10-second surface wind fields from 1-km TC simulations. The adjustment factor is derived from filtering the simulated TC wind fields onto various horizontal grid spacings typical of those used in GCMs. The new adjustment factor leads to TCs with greater intensity than the existing adjustment factors for horizontal grid spacings smaller than 30 km. This difference is attributed to more realistic wind fields in the TC simulations that contain highly asymmetric, localized patches of higher wind speeds instead of axisymmetric wind fields. Applying the new adjustment factor to select GCM simulations suggests the common interpretation of low-intensity bias in GCM-simulated TCs might be slightly exaggerated.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Determining the gluonic gravitational form factors of the proton

The proton is one of the main building blocks of all visible matter in the Universe. Among its intrinsic properties are its electric charge, mass and spin. These properties emerge from the complex dynamics of its fundamental constituents—quarks and gluons—described by the theory of quantum chromodynamics. The electric charge and spin of protons, which are shared among the quarks, have been investigated previously using electron scattering. An example is the highly precise measurement of the electric charge radius of the proton. By contrast, little is known about the inner mass density of the proton, which is dominated by the energy carried by gluons. Gluons are hard to access using electron scattering because they do not carry an electromagnetic charge. Here we investigated the gravitational density of gluons using a small colour dipole, through the threshold photoproduction of the J/ψ particle. We determined the gluonic gravitational form factors of the proton from our measurement. Here we used a variety of models and determined, in all cases, a mass radius that is notably smaller than the electric charge radius. In some, but not all cases, depending on the model, the determined radius agrees well with first-principle predictions from lattice quantum chromodynamics. This work paves the way for a deeper understanding of the salient role of gluons in providing gravitational mass to visible matter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

CMIP6 skill at predicting interannual to multi-decadal summer monsoon precipitation variability

Monsoons affect the economy, agriculture, and human health of two thirds of the world’s population. Therefore, predicting variations in monsoon precipitation is societally important. We explore the ability of climate models from the sixth phase of the Climate Model Intercomparison Project to predict summer monsoon precipitation variability by using hindcasts from the Decadal Climate Prediction Project (Component A). The multi-model ensemble-mean shows significant skill at predicting summer monsoon precipitation from one year to 6–9 years ahead. However, this skill is dependent on the model, monsoon domain, and lead-time. In general, the skill of the multi-model ensemble-mean prediction is low in year 1 but increases for longer-lead times and is largely consistent with externally forced changes. The best captured region is northern Africa for the 2–5 and 6–9 year forecast lead times. In contrast, there is no significant skill using the ensemble-mean over East and South Asia and, furthermore, there is significant spread in skill among models for these domains. By sub-sampling the ensemble we show that the difference in skill between models is tied to the simulation of the externally forced response over East and South Asia, with models with a more skilful forced response capable of better predictions. A further contribution is from skilful prediction of Pacific Ocean temperatures for the South Asian summer monsoon at longer lead-times. Therefore, these results indicate that predictions of the East and South Asian monsoons could be significantly improved.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Quantum Interference in Two-Atom Resonant X-ray Scattering of an Intense Attosecond Pulse

We theoretically investigate resonant x-ray scattering from two noninteracting Ne+ ions driven by an intense attosecond pulse using a nonrelativistic, QED-based time-dependent framework. Our model includes Rabi oscillations, photoionization, Auger decay, and quantum interference among elastic scattering and resonance fluorescence pathways. We analyze how the total scattering signal depends on pulse intensity, atomic configuration, and initial electronic state. We find that the total resonant scattering yield exceeds its nonresonant counterpart; the angular dependence of the signal qualitatively resembles a two-atom structure factor; and the visibility of interference fringes is sensitive to pulse area and the initial electronic state. Only a subset of final states reached via resonance fluorescence exhibits interference, determined by the indistinguishability of photon emission pathways. Fringe visibility is maximized in the linear scattering regime, where ionization is minimal and resonance fluorescence pathways can be largely indistinguishable. These results highlight optimal conditions for applying ultrafast resonant x-ray scattering to single-particle imaging.

Venkatesh, Akilesh↗
DOE OSTI2025

Relaxation times for disoriented isospin condensates in high energy heavy ion collisions

Fluctuations between charged and neutral kaons measured by the ALICE Collaboration in Pb-Pb collisions at the CERN Large Hadron Collider (LHC) exceed conventional explanations. Previously it was shown that if the scalar condensate is accompanied by an electrically neutral isospin-1 field then the combination can produce large equilibrium fluctuations where ⟨$\overline{𝑢}$ ⁢𝑢⟩≠⟨ $\overline{𝑑}$𝑑⟩. Hadronizing strange and antistrange quarks might then strongly fluctuate between charged ($\overline{𝑢}$⁢𝑠 or $\overline{𝑠}$⁢𝑢) and neutral (𝑑⁢$\overline{𝑠}$⁢ or 𝑠⁢$\overline{𝑑}$) kaons. Here, we estimate the times for the condensates to achieve their equilibrium probability distributions within causal volumes in high energy heavy ion collisions. This is achieved by modeling the temperature dependence of the condensates, mesonic collective excitations, decay rates of the associated fields, and employing the Langevin and Fokker-Planck equations. Within this model, we find that the equilibration times are short compared with the expansion time, making disoriented isospin condensates a viable explanation for the anomalous fluctuations observed at the LHC.

Quark-gluon plasma↗
DOE OSTI2022

Design Methods, Tools, and Data for Ceramic Solar Receivers Year 1 Continuation Report

This report describes the first year of work on a project to develop the methods, tools, and data required to analyze high temperature ceramic Concentrating Solar Power (CSP) components. This first year focused on developing the methods and data required for a time-independent assessment of potential components, focusing in particular on ceramic solar receivers. The report describes both model development and testing work focused on accomplishing this goal. Our overall conclusion is that high temperature ceramic CSP components are viable and could provide a means to overcome the expected low reliability and short service life for equivalent components constructed from Ni-based superalloys. Based on the results reported here, we recommend the project continue to Phase II which will develop more sophisticated, realistic models for time-dependent failure of ceramics operating in expected CSP component conditions and develop the time-dependent ceramic test data needed to parameterize these models, using commercial SiC as a reference material.

14 SOLAR ENERGY↗
DOE OSTI2024

DECOVALEX-2023: Task C Final Report

The Full-scale Emplacement (FE) heater experiment at the Mont Terri Underground Rock Laboratory (URL) was designed and conducted by Nagra to replicate an emplacement tunnel of Nagra’s reference repository design at 1:1 scale. Alongside testing the technical feasibility of constructing disposal tunnels, emplacing waste containers in the tunnels and then backfilling them, the main goals of the FE experiment are (1) to obtain a better understanding of the coupled effects of induced thermo-hydro-mechanical (THM) processes that may occur and (2) to validate existing coupled THM models (Müller et al., 2017). A key aspect of ensuring safety for repositories located in low-permeability rock involves minimizing any damage to the rock itself, thereby preserving its integrity and promoting a stable environment Amongst a number of processes that could damage the rock is the increase in pore pressure due to thermal loading caused by heat emitted from the waste. To reduce the potential damage of the rock, it is important to analyse the evolution of heat over time due to the heat load of the containers and assess possible consequences by coupled THM models. The aim of Task C of DECOVALEX-2023 was to build 3D numerical models of the FE experiment, focussing in particular on the heating induced pore pressure change in the Opalinus Clay. Data from a large number of sensors were available from the FE experiment for model comparison. These sensors measured temperature and relative humidity in the bentonite around the heaters, and temperature, pressure and displacement/strain in the surrounding Opalinus clay. Data were available from the start of excavation (April 2012) up to August 2020 for most sensors (more than 5 years from the start of heating in December 2014). To fulfil the overall aim of the task, the work was broken down into a number of steps, starting with simpler models to build confidence in each team’s approach and then moving to more complex models that better represent the FE experiment. Step 0 consisted of 2D benchmark models, gradually increasing the number of processes that are represented from thermal (T) only models in Step 0a, to coupled thermal hydraulic (TH) models in Step 0b with a representation of changing porosity, to coupled thermo-hydro-mechanical (THM) models in Step 0c, where porosity changes are calculated by the mechanical model. A detailed specification of processes, parameters, initial and boundary conditions was provided for this step, with the ambition that all teams would work towards close agreement in their model results, thus building confidence in the model implementations. vi It was not straightforward to achieve agreement between the teams, so additional steps (Step 0b2, 0b3, 0c2, 0c3) were added along with derivation of some analytical solutions against which the models could be compared. The reasons for the differences between teams were investigated and found to be caused primarily by different conceptual model assumptions (including temperature dependence of the thermal expansion of water), different model formulations (including porosity evolution) and differences in modelled domain sizes, boundary conditions and grid discretisation. This demonstrates that comparisons between multiple modelling teams and/or comparison with analytical results and experimental data are highly beneficial in providing an indication of uncertainty in model predictions. At the conclusion of Step 0, almost all teams had achieved a close agreement in model results and those that had not achieved an agreement knew the reason for this. Step 1 moved from 2D models to 3D models of the FE experiment without adding technical features like shotcrete or EDZ, and only considering the heating phase. Initially the 3D model was tightly specified to continue to build confidence in the model implementations (Step 1a). The results of Step 1a were compared to the data from the FE-experiment without the teams seeing the data. The teams were then provided with a sub-set of the data from the FE-experiment and invited to consider how best to use the large dataset for model comparison (Step 1b). Teams were then asked to use the data provided to calibrate their models, only changing material property values rather than adding features or processes to their models (Step 1c). In Step 1, teams were asked to only model the heating phase of the experiment, so pressure in the Opalinus Clay was reported as change in pressure since the initial conditions were specified rather than modelled. The change from 2D to 3D models was accompanied by an increase in the dispersion of results between the teams. Some of this was resolved during the task, but some remained and is potentially due to model discretisation. Calibration of parameters was useful in improving the fit of the models to the data but the remaining differences indicated that the models were missing features or processes. In Step 2, the teams were asked to update their models with additional features and processes as well as calibrating parameters to try and improve the fit of the models to the data. Teams were encouraged to represent ventilation of the open FE tunnel prior to backfilling with heaters and bentonite and in Step 2, the absolute pressure in the Opalinus Clay was compared between the teams. Teams took different approaches, but there was consideration of adding shotcrete and an EDZ into the model, representing stress change during excavation and different approaches to modelling ventilation of the FE tunnel. Overall, the documented results showed a very good agreement for temperature. The results for porewater pressure evolution showed a significant improvement for most teams compared to Step 1c with a good agreement to the measurements for several teams whereas some teams overpredicted the pressure increase and others overpredicted the drainage effect especially for the sensors close to the heater. Step 3 was an opportunity for teams to use the models developed in Step 1 and Step 2 to make predictions about the temperature and pressure changes that will be expected at the FE experiment over the next few years in light of the planned changes in thermal output of the heaters.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

A Model to Account for the Effects of Load Ratio and Hydrogen Pressure on the Fatigue Crack Growth Behavior of Pressure Vessel Steels

A phenomenological model for estimating the effects of load ratio R and hydrogen pressure P H2 on the hydrogen-assisted fatigue crack growth rate (HA–FCGR) behavior in the transient and steady-state regimes of pressure vessel steels is described. The “transient regime” is identified with crack growth within a severely embrittled zone of intense plasticity at the crack tip. The “steady-state” behavior is associated with the crack growing into a region of comparatively lower hydrogen concentration located further away from the crack tip. The model treats the effects of R and P H2 as being functionally separable. In the transient regime, the effects of the hydrogen pressure on the HA–FCGR behavior were negligible but were significant in the steady-state regime. The hydrogen concentration in the steady-state region is modeled as being dependent on the kinetics of lattice diffusion, which is sensitive to pressure. Experimental HA–FCGR data from the literature were used to validate the model. The new model was shown to be valid over a wide range of conditions that ranged between –1 ≤ R ≤ 0.8 and 0.02 ≤ P H2 ≤ 103 MPa for pressure vessel steels.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Constraining multi-field inflation using the SPHEREx all-sky survey power spectra

Abstract We investigate how well the SPHEREx all-sky survey can constrain local primordial non-Gaussianity beyond the parameterf NL using galaxy power spectra. We forecast joint constraints on the parametersf NL ,g NL andτ NL obtained assuming a simple two-field curvaton model of inflation. The parametersf NL andg NL characterise the squeezed limits of the primordial bispectrum and trispectrum respectively, and lead to a characteristic scale-dependence of the galaxy bias that increases out to arbitrarily large scales. Values of the parameterτ NL > (6/5f NL ) 2 cause the galaxy power spectrum to have a stochastic component which also increases out to arbitrarily large scales. Our MCMC forecasts indicate that SPHEREx can provide joint constraints on any two of the three parametersf NL ,g NL andτ NL . Due to strong degeneracies among these parameters, measurements of the galaxy power spectra alone may not be sufficient to jointly constrain all three. Constraints onf NL ,g NL andτ NL obtained from galaxy power spectrum observations depend on the modelling of underlying nuisance parameters. We study the robustness of our forecast constraints to modelling choices and note that even with relatively conservative modelling assumptions, SPHEREx galaxy power spectra can provide strong evidence of local non-Gaussianity, even if the particular values off NL andg NL cannot be measured precisely.

Astronomy & Astrophysics↗
DOE OSTI2021

Defense Waste Processing Facility Nitric-Glycolic Flowsheet Chemical Process Cell Chemistry: Part 2

The conversions of nitrite to nitrate, the destruction of glycolate, and the conversion of glycolate to formate and oxalate were modeled for the Nitric-Glycolic flowsheet using data from Chemical Process Cell (CPC) simulant runs conducted by Savannah River National Laboratory (SRNL) from 2011 to 2016. The goal of this work was to develop empirical correlation models to predict these values from measurable variables from the chemical process so that these quantities could be predicted a-priori from the sludge or simulant composition and measurable processing variables. The need for these predictions arises from the need to predict the REDuction/OXidation (REDOX) state of the glass from the Defense Waste Processing Facility (DWPF) melter. This report summarizes the work on these correlations based on the aforementioned data. Previous work on these correlations was documented in a technical report covering data from 2011-2015. This current report supersedes this previous report. Further refinement of the models as additional data are collected is recommended. The glass REDOX depends on the concentrations of nitrate and manganese (oxidants), and of glycolate, formate, oxalate, carbon, and antifoam (reductants) in the melter feed. The waste sludge contains nitrite, nitrate, manganese (Mn), and oxalate. Virtually all of the nitrite is converted to nitrate or NO+NO 2 +N 2 O gases in the CPC. The portion of the nitrite converted to nitrate increases the amount of nitrate in the sludge. The amount of glycolate in the final melter feed depends on the amount of the glycolic acid feed that is destroyed. Similarly, the amounts of formate and oxalate formed during the decomposition of glycolic acid are required. The material balance on carbon was found to not close in most cases. Generally, there was less carbon at the end of testing compared to the inputs. The most uncertain product variable was glycolate, so material balances were performed where the glycolate concentration was adjusted, usually upward, to close the balance. Correlation versus the original, as-measured, data was generally poor, but correlation against the material balance adjusted values was greatly improved. It was also shown that the correlation of the measured REDOX versus the predicted REDOX was much better when the material balance adjusted glycolate values were used. Three data series were primarily used during the regressions of the data; these series were 1) Sludge Batch 9 NG flowsheet simulant runs NG51-62 (SB9-NG); 2) Scaled Runs + Bounding Hydrogen Runs (SR+BH); and 3) Runs GN43-50 and 57 (43-50,57). The glycolate destruction was found to correlate with acid stoichiometry (AS), percent reducing acid (PRA), and for some data series, headspace to simulant volume ratio (HSV), mercury (Hg), and nitrate. Although glycolate destruction for pairs of data series (e.g., [SB9-NG] and [SR+BH]) were found to depend on HSV, the combination of all three data series was not found to have significant dependence on this variable. The best model for glycolate destruction depended on AS, nitrate, and Hg. This model predicted the product glycolate compositions of the data to within 92-106%. The conversion of glycolate to formate was high when noble metals and Hg were not present, with values up to 100%. When noble metals and Hg were present, this conversion ranged from zero to 7%, and was dependent on AS. Lower AS gave higher conversions to formate. The conversion to oxalate was found to depend on the AS and the initial concentration of nitrite. An alternative fit versus AS and the form of ruthenium (Ru) used is a possible alternative. This fit was somewhat less statistically significant. This second model predicts that more oxalate is formed when Ru-nitrosyl nitrate is used rather than Ru chloride. The conversion of glycolate to oxalate ranged from zero to 6%. The conversion of nitrite to nitrate depended primarily on AS and PRA, with HSV and Hg being significant when these variables were varied. For multiple series of data, nitrite was also needed to SRNL-STI-2017-00172 5HYLVLRQ viL distinguish between data series, and the effect of HSV became insignificant. The best model for nitrite to nitrate conversion depended on AS, PRA, nitrite, and Hg. The 95% confidence intervals on the predicted values of glycolate destruction, glycolate to oxalate conversion, and nitrite to nitrate conversion were used to determine the uncertainty in the predicted REDOX when starting with only the composition of the sludge, AS, and PRA. Using the 95% confidences on an individual value (that is the confidence in getting a particular value for one single test as opposed to what the mean would be for multiple tests), the uncertainty in the predicted REDOX was calculated. The uncertainty in the actual product composition glycolate, oxalate, formate, and nitrate concentrations translated to an uncertainty in the REDOX value of ±0.1,which is approximately the uncertainty claimed in the REDOX model itself.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Impact of light output on the timing resolution of organic glass scintillator bars in a dual-ended readout configuration

The effectiveness of detector system modeling and validation depends on the quality of characterization performed on the system. For nuclear nonproliferation applications, which require detectors to address complex and variable field conditions, access to accurate system models is essential. This study presents the timing characterization of organic glass scintillator bars used in an imaging system developed at the University of Michigan. 137 Cs and 60 Co gamma-ray sources were used to determine the coincidence time resolution as a function of measured light output for two Compton scatter events between two organic glass scintillator bars. Timing resolution for two events, each with light output ranging from 100 to 1100 keVee were characterized. The resulting data were fit to a parametric function that was validated to agree with experimental results within ±25 ps. Simulations were performed to establish appropriate calibration points for each organic scintillator, thereby ensuring accurate calibration of light output values during analysis. This work successfully characterized the light output dependence of the coincidence timing resolution of a pair of organic glass detectors for use in an imaging system. Depending on the amount of scintillation light output from the events, the full width at half maximum of the measured coincidence time resolution ranged from 653.2 ± 16.1 ps for low light yields, 100 keVee, to 121.0 ± 8.4 ps for higher yields at 1100 keVee. The insights obtained from the timing resolution behavior will allow for accurate simulation capabilities in future security and verification efforts using scatter-based imaging systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Region and cloud regime dependence of parametric sensitivity in E3SM atmosphere model

Abstract The Department of Energy (DOE)’s Energy Exascale Earth System Model (E3SM), including its atmosphere model (EAM), has many relatively new features. In a previous study we conducted a systematic parametric sensitivity analysis for EAM based on short, perturbed parameter ensemble (PPE) simulations, mainly focusing on global mean climate features and metrics. While parameter values in global climate models are generally invariant in space and time, model response to parameters perturbation may vary by regions and climate regimes, which motivates the need to better understand the EAM model behaviors and physics at regional scale and process level. In this study, using the same set of PPE simulations and a similar sensitivity analysis framework, we identify parameters that cause largest sensitivities over different regions and compare model responses in fast atmospheric processes to the parameters across different cloud regimes for several important cloud-related fidelity metrics. We find that cloud forcing has opposite response to some parameters over mid-latitude vs. tropical land. We also analyze how the parametric sensitivity varies as stratocumulus transitions to shallow convection and to deep convection over ocean. Low cloud forcing and shortwave cloud forcing in the subtropical eastern Pacific are most sensitive to the parameters controlling the width of the probability density function (PDF) of the subgrid vertical velocity ( w’ ) ( gamma ) and the damping of the w’ skewness ( c8 ) near the coast but become more sensitive to the parameter affecting the damping of the w’ variance ( c1 ) further offshore. Detailed interpretation of the spatial dependence of parametric sensitivity is provided. We also investigate how the parametric sensitivity evolves with prediction duration. This study improves our process-level understanding of cloud physics and parameterization and provides insights for developing more advanced regime-aware parameterization schemes in global climate model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

An experimentally validated mesoscale model for the effective thermal conductivity of U-Zr fuels

An experimentally validated mesoscale model was developed to simulate the effective thermal conductivity (ETC) of metallic U and U–Zr nuclear fuels. The effects of microstructure, temperature, composition, and interfacial thermal resistance (ITR) were investigated. Companion experiments were conducted to validate the model. The numerical simulations clearly demonstrate that accounting for the interface (Kapitza) thermal resistance and Zr precipitation is necessary to improve the model predictions. The dependence of the effective Kapitza resistance of depleted-U and U-10Zr on temperature was determined. Finally, for both materials, the largest difference between the model calculations and the experimental data was about 4.5%, which is within the precision of the experimental measurements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Time-resolved site-selective imaging of predissociation and charge transfer dynamics: the CH 3 I B-band

The predissociation dynamics of the 6s (B 2 E) Rydberg state of gas-phase CH 3 I were investigated by time-resolved Coulomb-explosion imaging using extreme ultraviolet (XUV) free-electron laser pulses. Inner-shell ionization at the iodine 4d edge was utilized to provide a site-specific probe of the ensuing dynamics. The combination of a velocity-map imaging (VMI) spectrometer coupled with the pixel imaging mass spectrometry (PImMS) camera permitted three-dimensional ionic fragment momenta to be recorded simultaneously for a wide range of iodine charge states. In accord with previous studies, initial excitation at 201.2 nm results in internal conversion and subsequent dissociation on the lower-lying A-state surface on a picosecond time scale. Examination of the time-dependent yield of low kinetic energy iodine fragments yields mechanistic insights into the predissociation and subsequent charge transfer following multiple ionization of the iodine products. The effect of charge transfer was observed through differing delay-dependencies of the various iodine charge states, from which critical internuclear distances for charge transfer could be inferred and compared to a classical over-the-barrier model. Time-dependent photofragment angular anisotropy parameters were extracted from the central slice of the Newton sphere, without Abel inversion, and highlight the effect of rotation of the parent molecule before dissociation, as observed in previous works. Our results demonstrate the ability to perform three-dimensional ion imaging at high event rates and showcase the potential benefits of this approach, particularly in relation to further time-resolved studies at free-electron laser facilities.

74 ATOMIC AND MOLECULAR PHYSICS↗