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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Coupled Multiphysics Modeling of Lithium-Ion Batteries for Automotive Crashworthiness Applications

Considerable advances have been made in battery safety models, but achieving predictive accuracy across a wide range of conditions continues to be challenging. Interactions between dynamically evolving mechanical, electrical, and thermal state variables make model prediction difficult during mechanical abuse scenarios. In this study, we develop a physics-based modeling approach that allows for choosing between different mechanical and electrochemical models depending on the required level of analysis. We demonstrate the use of this approach to connect cell-level abuse response to electrode-level and particle-level transport phenomena. A pseudo-two-dimensional model and simplified single-particle models are calibrated to electrical-thermal cycling data and applied to mechanically induced short-circuit scenarios to understand how the choice of electrochemical model affects the model prediction under abuse scenarios. These models are implemented using user-defined subroutines on ls-dyna finite element software and can be coupled with existing automotive crash safety models.

analysis and design of components↗
DOE OSTI2022

An Electron-scattering Time Delay in Black Hole Accretion Disks

Abstract Universal to black hole X-ray binaries, the high-frequency soft lag gets longer during the hard-to-intermediate state transition, evolving from ≲1 to ∼10 ms. The soft lag production mechanism is thermal disk reprocessing of nonthermal coronal irradiation. X-ray reverberation models account for the light-travel time delay external to the disk, but assume instantaneous reprocessing of the irradiation inside the electron-scattering-dominated disk atmosphere. We model this neglected scattering time delay as a random walk within an α -disk atmosphere, with approximate opacities. To explain soft lag trends, we consider a limiting case of the scattering time delay that we dub the thermalization time delay, t th ; this is the time for irradiation to scatter its way down to the effective photosphere, where it gets thermalized, and then scatter its way back out. We demonstrate that t th plausibly evolves from being inconsequential for low mass accretion rates m ̇ characteristic of the hard state, to rivaling or exceeding the light-travel time delay for m ̇ characteristic of the intermediate state. However, our crude model confines t th to a narrow annulus near peak accretion power dissipation, so cannot yet explain in detail the anomalously long-duration soft lags associated with larger disk radii. We call for time-dependent models with accurate opacities to assess the potential relevance of a scattering delay.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Multifidelity uncertainty quantification with models based on dissimilar parameters

Multifidelity uncertainty quantification (MF UQ) sampling approaches have been shown to significantly reduce the variance of statistical estimators while preserving the bias of the highest-fidelity model, provided that the low-fidelity models are well correlated. However, maintaining a high level of correlation can be challenging, especially when models depend on different input uncertain parameters, which drastically reduces the correlation. Existing MF UQ approaches do not adequately address this issue. In this work, we propose a new sampling strategy that exploits a shared space to improve the correlation among models with dissimilar parameterization. We achieve this by transforming the original coordinates onto an auxiliary manifold using the adaptive basis (AB) method (Tipireddy and Ghanem, 2014). The AB method has two main benefits: (1) it provides an effective tool to identify the low-dimensional manifold on which each model can be represented, and (2) it enables easy transformation of polynomial chaos representations from high- to low-dimensional spaces. This latter feature is used to identify a shared manifold among models without requiring additional evaluations. Here we present two algorithmic flavors of the new estimator to cover different analysis scenarios, including those with legacy and non-legacy high-fidelity (HF) data. We provide numerical results for analytical examples, a direct field acoustic test, and a finite element model of a nuclear fuel assembly. For all examples, we compare the proposed strategy against both single-fidelity and MF estimators based on the original model parameterization.

42 ENGINEERING↗
DOE OSTI2024

Low- ν method with LHC neutrinos

The Forward Physics Facility (FPF) plans to use neutrinos produced at the Large Hadron Collider to make a variety of measurements at previously unexplored TeV energies. Its primary goals include precision measurements of the neutrino cross section and using the measured neutrino flux both to uncover information about far-forward hadron production and to search for various beyond standard model scenarios. However, these goals have the potential to conflict: Extracting information about the flux or cross section relies upon an assumption about the other. In this paper, we demonstrate that the FPF can use the low- ν method—a technique for constraining the flux shape by isolating neutrino interactions with low energy transfer to the nucleus—to break this degeneracy. We show that the low- ν method is effective for extracting the ν μ flux shape, in a model-independent way. We discuss its application for extracting the ν ¯ μ flux shape but find that this is significantly more model dependent. Finally, we explore the precision to which the ν μ flux shape could be constrained at the FPF for a variety of proposed detector options. We find that the precision would be sufficient to discriminate between various realistic flux models. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

A search for microscopic black holes, string balls, and sphalerons in proton-proton collisions at $\sqrt{s}=13$ TeV

A search for microscopic black holes, string balls, and electroweak sphalerons using proton-proton collisions at $\sqrt{s}=13$ TeV recorded with the CMS detector at the CERN LHC during the 2016–2018 data taking, and corresponding to an integrated luminosity of 138 fb −1 , is presented. Two search strategies based on control samples in data are used. Model-independent limits on the cross section of physics phenomena with multiple energetic jets, leptons, and photons are set using a method that relies on the shape invariance of the scalar sum of the transverse momenta of all objects in the event. Model-dependent limits on black hole and sphaleron production are set using a newly introduced method that has been developed for the identification of collider events with distinct kinematic features by separating them into classes based on phase space proximity. In the context of models with large extra dimensions, semiclassical black holes and string balls with masses below 8.4–11.4 TeV and 9.0–10.7 TeV, respectively, are excluded at 95% confidence level, significantly extending the reach beyond previous searches. Results of a dedicated search for electroweak sphalerons are used to derive an upper limit of 0.0034 at 95% confidence level on the fraction of quark-quark interactions, with a center-of-mass energy above the nominal sphaleron transition energy threshold of 9 TeV, that result in a sphaleron transition.

Beyond Standard Model↗
DOE OSTI2020

Effects of Void Morphology on Detonation Initiation of PBX 9502

Accurate predictive modeling of high explosives in abnormal conditions, such as a fire, is critical to personnel safety. Modeling depends heavily on precise physical characterization of the high explosive in question. The aim of this thesis was to characterize the detonation response to shock and the microstructure at high-temperature (250°C) of PBX 9502 and use these data to attempt to model these results using a new reactive burn model called SURF (Scaled Uniform Reactive Front). In order to characterize the shock response of PBX 9502 at high-temperature a 1D gas gun experiment was designed and executed that provided shock to detonation transition data. These data were then used to do a preliminary calibration of SURF for PBX 9502 at 250°C. Small-angle neutron scattering (SANS) was employed to characterize the void morphology of PBX 9502 as a function of temperature. This was done in order to feed a modified version of SURF called physically informed - SURF (π-SURF) which uses the void size and quantity distribution data as measured by SANS along with the deflagration and detonation characteristics of a high explosive to inform the SURF model. The π-SURF model was then used to model both ambient and high-temperature PBX 9502 SDT response and the results were compared to the physically measured results. It was found that π-SURF works in ambient cases. At high-temperature, however, it was determined that not enough information is known about hot spot formation to use the π-SURF model for high-temperature scenarios.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2023

Evaluating skill in predicting the Interdecadal Pacific Oscillation in initialized decadal climate prediction hindcasts in E3SMv1 and CESM1 using two different initialization methods and a small set of start years

Abstract It is a daunting challenge to conduct initialized hindcasts with enough ensemble members and associated start years to form a drifted climatology from which to compute the anomalies necessary to quantify the skill of the hindcasts when compared to observations. This limits the ability to experiment with case studies and other applications where only a few initial years are needed. Here we run a set of hindcasts with CESM1 and E3SMv1 using two different initialization methods for a limited set of start years and use the respective uninitialized free-running historical simulations to form the model climatologies. Since the drifts from the observed initial states in the hindcasts toward the uninitialized model state are large and rapid, after a few years the drifted initialized models approach the uninitialized model climatological errors. Therefore, hindcasts from the limited start years can use the uninitialized climatology to represent the drifted model states after about lead year 3, providing a means to compute forecast anomalies in the absence of a large hindcast sample. There is comparable skill for predicting spatial patterns of multi-year Pacific sea surface temperature anomalies in the domain of the Interdecadal Pacific Oscillation using this method compared to the conventional methodology with a large hindcast data set, though there is a model dependence to the drifts in the two initialization methods.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Search for a CP-odd Higgs boson decaying into a heavy CP-even Higgs boson and a Z boson in the ${\ell}^{+}{\ell}^{-}t\overline{t}$ and $\nu\overline{\nu}b\overline{b}$ final states using 140 fb –1 of data collected with the ATLAS detector

A search for a heavy CP-odd Higgs boson, A, decaying into a Z boson and a heavy CP-even Higgs boson, H, is presented. It uses the full LHC Run 2 dataset of pp collisions at $\sqrt{s}$ = 13 TeV collected with the ATLAS detector, corresponding to an integrated luminosity of 140 fb –1 . The search for A → ZH is performed in the $ℓ^+ℓ^–t\bar{t}$ and $ν\bar{ν}b\bar{b}$ final states and surpasses the reach of previous searches in different final states in the region with mH > 350 GeV and mA > 800 GeV. No significant deviation from the Standard Model expectation is found. Upper limits are placed on the production cross-section times the decay branching ratios. Limits with less model dependence are also presented as functions of the reconstructed $m(t\bar{t}$) and $m(b\bar{b}$) distributions in the $ℓ^+ℓ^–t\bar{t}$ and $ν\bar{ν}b\bar{b}$ channels, respectively. In addition, the results are interpreted in the context of two-Higgs-doublet models

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Strong CP and flavor in multi-Higgs theories

We introduce a class of multi-Higgs doublet extensions of the Standard Model that solve the strong problem with profound consequences for the flavor sector. The Yukawa matrices are constrained to have many zero entries by a “Higgs-flavor” symmetry, , that acts on Higgs and quark fields. The violation of both and occurs in the Higgs mass matrix so that, for certain choices of charges, the strong parameter is zero at tree level. Radiative corrections to are computed in this class of theories. They vanish in realistic two-Higgs doublet models with . We also construct realistic three-Higgs models with , where the one-loop results for are model-dependent. Requiring has important implications for the flavor problem by constraining the Yukawa coupling and Higgs mass matrices. Contributions to from higher-dimension operators are computed at one loop and can also be sufficiently small, although the hierarchy problem of this class of theories is worse than in the Standard Model.

Hall, Lawrence↗
DOE OSTI2022

Late day measurement of excised branches results in uncertainty in the estimation of two stomatal parameters derived from response curves in Populus deltoides Bartr. x Populus nigra L.

Many terrestrial biosphere models depend on an understanding of the relationship between stomatal conductance and photosynthesis. However, unlike the measurement of photosynthetic parameters, such as the maximum carboxylation capacity, where standard methods (e.g., CO 2 response or ACi curves) are widely accepted, a consensus method for empirically measuring parameters representing stomatal response has not yet emerged. Most models of stomatal response to environment represent stomatal conductance as being bounded by a lower intercept parameter (g 0 ), and linearly scaled based on a multivariate term described by the stomatal slope parameter (g 1 ). Here we employ the widely used Unified Stomatal Optimization model, to test whether g 1 and g 0 parameters are impacted by the choice of measurement method, either on an intact branch, or a cut branch segment stored in water. We measured paired stomatal response curves on intact and excised branches of a hybrid poplar clone (Populus deltoides Bartr. x Populus nigra L. OP367), measured twice over a diurnal period. We found that predawn branch excision did not significantly affect measured g 0 and g 1 when measured within 4 h of excision. Measurement in the afternoon resulted in significantly higher values of g 1 and lower values of g0, with values changing by 55% and 56% respectively. Excision combined with afternoon measurement resulted in a marked effect on parameter estimates, with g1 increasing 89% from morning to afternoon and a 25% lower g 1 for cut branches than those measured in situ. We also show that in hybrid poplar the differences in parameter estimates obtained from plants measured under different conditions can directly impact models of canopy function, reducing modeled transpiration by 18% over a simulated 12.5-hour period. While these findings are only for a single isohydric woody species, our findings suggest that stomatal optimality parameters may not remain constant throughout the day.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

The Kinetics of the Coherent Stage of Hydroxide Corrosion on Li 2 O-Covered LiH

The capability to model LiOH growth on vacuum-baked LiH is a necessary precursor to making kinetic predictions of undesirable hydrogen outgassing from corroded LiH materials when placed in sealed-system applications. Here, in this work, LiH samples with initial LiOH surface corrosion layers less than 1.1 μm were outgassed by vacuum baking at high temperature to convert LiOH to Li 2 O. Diffuse-reflectance infrared Fourier transform (DRIFT) spectroscopy was used to measure the subsequent LiOH regrowth during exposure to 25–375 ppm of H 2 O vapor at room temperature. Logarithmic kinetics best model the growth of this coherent hydroxide corrosion layer (up to 1.5 μm) on previously vacuum-baked LiH samples. The logarithmic kinetics are attributed to the dissociation of H 2 O on the corrosion layer and the establishment of a H + /OH – electric field across the thickening LiOH corrosion layer. A time-dependent model of LiOH corrosion growth as a function of both initial LiOH thickness before vacuum baking and moisture partial pressures during re-exposure was developed to help researchers better assess the unwanted hydrogen outgassing potential from LiH/LiD materials.

Matt, Sarah M. [Lawrence Livermore National Labora↗
DOE OSTI2022

Robust estimation of charge carrier diffusivity using transient photoluminescence microscopy

Transient microscopy has emerged as a powerful tool for imaging the diffusion of excitons and free charge carriers in optoelectronic materials. In many excitonic materials, extraction of diffusion coefficients can be simplified because of the linear relationship between signal intensity and local excited state population. However, in materials where transport is dominated by free charge carriers, extracting diffusivities accurately from multidimensional data is complicated by the nonlinear dependence of the measured signal on the local charge carrier density. To obtain accurate estimates of charge carrier diffusivity from transient microscopy data, statistically robust fitting algorithms coupled to efficient 3D numerical solvers that faithfully relate local carrier dynamics to raw experimental measurables are sometimes needed. Here, we provide a detailed numerical framework for modeling the spatiotemporal dynamics of free charge carriers in bulk semiconductors with significant solving speed reduction and for simulating the corresponding transient photoluminescence microscopy data. To demonstrate the utility of this approach, we apply a fitting algorithm using a Markov chain Monte Carlo sampler to experimental data on bulk CdS and methylammonium lead bromide (MAPbBr 3 ) crystals. Parameter analyses reveal that transient photoluminescence microscopy can be used to obtain robust estimates of charge carrier diffusivities in optoelectronic materials of interest, but that other experimental approaches should be used for obtaining carrier recombination constants. Additionally, simplifications can be made to the fitting model depending on the experimental conditions and material systems studied. Our open-source simulation code and fitting algorithm are made freely available to the scientific community.

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

Towards Anomaly Detection at the CMS High-Level Trigger System

Traditional trigger strategies in CMS typically rely on model-dependent selections or rigid kinematic cuts, risking the omission of unexpected exotic signatures. To address this, we propose a novel anomaly detection (AD) algorithm for the High-Level Trigger (HLT), designed to serve as a complementary second layer of filtering to the Level-1 AXOL1TL AD algorithm. We employ a transformer-based foundation model trained on a diverse ensemble of Standard Model processes. By combining a joint contrastive and classification objective, and using particle kinematics as inputs, the model learns to map events to a physics-informed latent space where anomalous events are isolated from dominant backgrounds. Preliminary results show that this strategy enhances the signal-to-background ratio across a range of rare SM and BSM scenarios. Furthermore, this work constitutes foundational R&D for the potential implementation of an analogous AD algorithm in the Level-1 trigger system for Phase-2.

Cruz, Roy [U. Wisconsin, Madison (main)] (ORCID:00↗
DOE OSTI2026

LeWRON: Agentic Analysis of Electroweak Phase Transitions

The electroweak phase transition (EWPT) is a central topic in particle physics and cosmology, connecting collider phenomenology, baryogenesis, and gravitational-wave observatories. Its analysis requires a technically demanding, convention-sensitive, and model-dependent pipeline, from constructing the finite-temperature effective potential to tracking thermal histories, computing bubble nucleation rates, and predicting gravitational-wave spectra. We present LeWRON (Learning ElectroWeak phase tRansitiON), an agentic framework that orchestrates this pipeline starting from an input Lagrangian. LeWRON combines audited toolbox construction with an Explorer module that uses the generated model-specific code for further analysis, including scans and plots. Intermediate analytic outputs are checked by auditor agents and stored as structured artifacts, enabling reproducible human inspection and downstream use through both a command-line interface and a public Python API. The framework supports a reproduction mode, which infers conventions from the literature and reproduces published results, and a discovery mode, which guides users through structured checkpoints for new models. We demonstrate LeWRON across representative beyond-the-Standard-Model scenarios and release the code on GitHub.

Wang, Isaac R. [Fermilab] (ORCID:000000030789218X)↗
DOE OSTI2022

Development of a three-dimensional APOLLO3 neutrons deterministic scheme for the CABRI reactor

CABRI is an experimental reactor to study the fuel behavior during reactivity injection transients. These transients being highly multiphysics, the development of suitable modeling and simulation tools to simulate them is important for the optimization of the tests and the control of the experimental conditions. This paper focuses on the development of an APOLLO3 deterministic core calculation tool dedicated to the CABRI transient analysis. It represents the first stage of the incremental process for the implementation of a multiphysics time-dependent modeling of the CABRI transient. The neutron calculation scheme is based on a classical two-step approach. The first step consists of a 281-energy group calculation flux with the TDT-MOC (Method Of Characteristics) solver for cross-section space and energy (23 groups) collapsing for the CABRI different assembly clusters. The bias on a 2D core neutron calculation due to the self-shielding calculation and collapsing on a restricted pattern are investigated thanks to a comparison with a direct full 2D calculation on a quarter of core. The second step relies on a pin-resolved transport 3D transport core calculation with the SN solver MINARET. A progressive numerical validation process is followed to quantify the calculation biases on reactivity and reaction rates at each step using reference calculations with the stochastic code TRIPOLI4. The next development stage toward a multiphysics scheme will be the implementation of the 3D-kinetics equation resolution and the coupling with a core thermal-hydraulics model. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Measurement of the $H \rightarrow \gamma \gamma $ and $H \rightarrow ZZ^* \rightarrow 4 \ell $ cross-sections in $pp$ collisions at $\sqrt{s}=13.6$ TeV with the ATLAS detector

The inclusive Higgs boson production crosssection is measured in the di-photon and the Z Z* → 4$\ell$ decay channels using 31.4 and 29.0 fb –1 of pp collision data respectively, collected with the ATLAS detector at a centreof-mass energy of $\sqrt{s}$ = 13.6 TeV. To reduce the model dependence, the measurement in each channel is restricted to a particle-level phase space that closely matches the channel’s detector-level kinematic selection, and it is corrected for detector effects. These measured fiducial cross-sections are σ fid,γγ = $76^{+14}_{–13}$ fb, and $σ_{\text{fid,4}\ell}$ = 2.80 ± 0.74 fb, in agreement with the corresponding Standard Model predictions of 67.6±3.7 fb and 3.67±0.19 fb. Assuming Standard Model acceptances and branching fractions for the two channels, the fiducial measurements are extrapolated to the full phase space yielding total cross-sections of σ (pp → H) = $67^{+12}_{–11}$ pb and 46±12 pb at 13.6 TeV from the di-photon and Z Z* → 4$\ell$ measurements respectively. The two measurements are combined into a total cross-section measurement of σ (pp → H) = 58.2±8.7 pb, to be compared with the Standard Model prediction of σ (pp → H) SM = 59.9 ± 2.6 pb.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Interplay of freeze-in and freeze-out: Lepton-flavored dark matter and muon colliders

We study a lepton-flavored dark matter model and its signatures at a future muon collider. We focus on the less-explored regime of feeble dark matter interactions, which suppresses the dangerous lepton-flavor-violating processes, gives rise to dark matter freeze-in production, and leads to long-lived particle signatures at colliders. We find that the interplay of dark matter freeze-in and its mediator freeze-out gives rise to an upper bound of around TeV scales on the dark matter mass. The signatures of this model depend on the lifetime of the mediator and can range from generic prompt decays to more exotic long-lived particle signals. In the prompt region, we calculate the signal yield, study useful kinematics cuts, and report tolerable systematics that would allow for a 5 σ discovery. In the long-lived region, we calculate the number of charged tracks and displaced lepton signals of our model in different parts of the detector and uncover kinematic features that can be used for background rejection. We show that, unlike in hadron colliders, multiple production channels contribute significantly, which leads to sharply distinct kinematics for electroweakly charged long-lived particle signals. Ultimately, the collider signatures of this lepton-flavored dark matter model are common among models of electroweak-charged new physics, rendering this model a useful and broadly applicable benchmark model for future muon collider studies that can help inform work on detector design and studies of systematics. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2023

A coarse-grained simulation model for colloidal self-assembly via explicit mobile binders

Colloidal particles with mobile binding molecules constitute a powerful platform for probing the physics of self-assembly. Binding molecules are free to diffuse and rearrange on the surface, giving rise to spontaneous control over the number of droplet–droplet bonds, i.e., valence, as a function of the concentration of binders. This type of valence control has been realized experimentally by tuning the interaction strength between DNA-coated emulsion droplets. Optimizing for valence two yields droplet polymer chains, termed ‘colloidomers’, which have recently been used to probe the physics of folding. To understand the underlying self-assembly mechanisms, here we present a coarse-grained molecular dynamics (CGMD) model to study the self-assembly of this class of systems using explicit representations of mobile binding sites. Further, we explore how valence of assembled structures can be tuned through kinetic control in the strong binding limit. More specifically, we optimize experimental control parameters to obtain the highest yield of long linear colloidomer chains. Subsequently tuning the dynamics of binding and unbinding via a temperature-dependent model allows us to observe a heptamer chain collapse into all possible rigid structures, in good agreement with recent folding experiments. Our CGMD platform and dynamic bonding model (implemented as an open-source custom plugin to HOOMD-Blue) reveal the molecular features governing the binding patch size and valence control, and opens the study of pathways in colloidomer folding. This model can therefore guide programmable design in experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗