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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Alkali-silica reaction expansion model for confined concrete with stress-dependency and casting direction anisotropy

Alkali-silica reaction (ASR) is a deleterious chemical reaction between alkali hydroxyl ions and types of silica found in some aggregates of concrete. Owners and regulators of nuclear power plants aim to ensure the safety of the concrete structures with optimal maintenance strategies. A new model was developed for predicting the expansion of concrete structures affected by alkali-silica reaction. The model includes a novel combination of existing models as an alkali-silica reaction advancement model, a casting direction anisotropic expansion model, a stress-dependent anisotropic expansion model, and a material property evolution model dependent on the degree of ASR expansion. The model parameters were calibrated based on existing literature data and data generated by previous efforts of this study. The calibrated model was then validated with the experiments carried out in previous efforts of this study. The model was shown to accurately predict the ASR-expansion of large-scale reinforced concrete specimens with confinement.

36 MATERIALS SCIENCE↗
DOE OSTI2022

3D simulations of spinlike flames in Co/Al multilayers with enhanced conduction losses

Reactive Co/Al multilayers are uniformly structured materials that may be ignited to produce rapid and localized heating. Prior studies varying the bilayer thickness (i.e., sum of two individual layers of Co and Al) have revealed different types of flame morphologies, including: (a) steady/planar, (b) wavy/periodic, and (c) transverse bands, originating in the flame front. These instabilities resemble the “spin waves” first observed in the early studies of solid combustion (i.e., Ti cylinder in a N 2 atmosphere), and are likewise thought to be due to the balance of heat released by reaction and heat conduction forward into the unreacted multilayer. However, the multilayer geometry and three-dimensional (3D) edge effects are relatively unexplored. In this work, a new diffusion-limited reaction model for Co/Al multilayers was implemented in large, novel 3D finite element analysis (FEA) simulations, in order to study the origins of these spinlike flames. This reaction model builds upon previous work by introducing three new phase-dependent property models for: (1) the diffusion coefficient, (2) anisotropic thermal conductivity tensor, and (3) bulk heat capacity, as well as one additional model for the bilayer-dependent heat of reaction. These novel 3D simulations are the first to predict both steady and unsteady flames in Co/Al multilayers. Moreover, two unsteady modes of flame propagation are identified, which depend on the enhanced conduction losses with slower flames, as well as flame propagation around notched edges. Finally, future work will consider the generality of the current modeling approach and also seek to define a more generalized set of stability criteria for additional multilayer systems.

42 ENGINEERING↗
DOE OSTI2023

Forward sensitivity analysis and mode dependent control for closure modeling of Galerkin systems

Model reduction by projection-based approaches is often associated with losing some of the important features that contribute towards the dynamics of the retained scales. As a result, a mismatch occurs between the predicted trajectories of the original system and the truncated one. We put forth a framework to apply a continuous time control signal in the latent space of the reduced order model (ROM) to account for the effect of truncation. We set the control input using parameterized models by following energy transfer principles. Our methodology relies on observing the system behavior in the physical space and using the projection operator to restrict the feedback signal into the latent space. Then, we leverage the forward sensitivity method (FSM) to derive relationships between the feedback and the desired mode-dependent control. We test the performance of the proposed approach using two test cases, corresponding to viscous Burgers and vortex merger problems at high Reynolds number. Results show that the ROM trajectory with the applied FSM control closely matches its target values in both the data-dense and data-sparse regimes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Modeling Kilonova Light Curves: Dependence on Nuclear Inputs

The mergers of binary neutron stars, as well as black hole–neutron star systems, are expected to produce an electromagnetic counterpart that can be analyzed to infer the element synthesis that occurred in these events. We investigate one source of uncertainties pertinent to lanthanide-rich outflows: the nuclear inputs to rapid neutron capture nucleosynthesis calculations. We begin by examining 32 different combinations of nuclear inputs: eight mass models, two types of spontaneous fission rates, and two types of fission daughter product distributions. We find that such nuclear physics uncertainties typically generate at least one order of magnitude uncertainty in key quantities such as the nuclear heating (one and a half orders of magnitude at 1 day post-merger), the bolometric luminosity (one order of magnitude at 5 days post-merger), and the inferred mass of material from the bolometric luminosity (factor of 8 when considering the 8–10 day region). Since particular nuclear processes are critical for determining the electromagnetic signal, we provide tables of key nuclei undergoing β-decay, α-decay, and spontaneous fission important for heating at different times, identifying decays that are common among the many nuclear input combinations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Collision-Energy Dependence in Heavy-Ion Collisions from Nonlinear QCD Evolution

We explore the effects of including the energy dependence determined from evolution equations within the color glass condensate framework on observables in ultrarelativistic heavy-ion collisions. This amounts to integrating the JIMWLK evolution equations into the impact parameter-dependent glasma model, which is then coupled to viscous relativistic hydrodynamics. This methodology allows for a systematic representation of nuclei at specific Bjorken-x values, which are probed at different center-of-mass energies of the collision and rapidities of final state particles. Comparing the methodology to the conventional impact parameter-dependent glasma model, we find significant effects on multiplicity distributions and particle spectra, especially in smaller collision systems at the highest center-of-mass energies. Our results highlight the importance of incorporating nonlinear QCD evolution in the description of heavy-ion collisions at varying center-of-mass energies, as the precise extraction of transport coefficients will be affected. This Letter establishes a robust framework for understanding the quark gluon plasma and nuclear structure at high energy, integrating small-x physics into the initial conditions of heavy-ion collisions.

Mäntysaari, Heikki↗
DOE OSTI2026

Multiphysics Analysis of Li Cooled Divertor Substrate During Loss of Coolant Accident (LOCA)

In the ongoing study of potential designs for liquid metal (LM) plasma-facing components (PFCs), so-called “slow” and “fast” Li flow divertor concepts are under investigation. In the previous studies on design and analysis of the slow Li flow divertor and comparison with the fast Li flow divertor, the magnetohydrodynamics (MHD)/heat transfer effects of the Li flowing inside the substrate as a second coolant were comprehensively investigated under the normal steady-state operation conditions. Here, in the present study, the multiphysics analysis is extended to the unsteady abnormal divertor scenario where the Li layer on top of the substrate does not provide full coverage or even totally disappears for a certain period of time, so that the substrate becomes directly exposed to the incident high plasma heat flux. Such an unwanted event may happen regardless of the concept of the divertor and is worth detailed investigations, typically referred to as a loss of coolant accident (LOCA). To address this situation, a simplified scoping analysis is conducted first in 2-D, and then an integrated 3-D modeling is performed using a time-dependent multiphysics model in COMSOL Multiphysics that integrates LM MHD, heat transfer, and solid mechanics. The main goal is to evaluate conditions under which the major material limits, such as the maximum allowable temperature, stress, and displacement of the substrate, can still be met. It was shown that the maximum time over which the substrate of RAFM steel can retain structural integrity during the LOCA is around 0.2 ∼ 0.3 s. Any divertor concept that utilizes RAFM steel as a substrate material and liquid Li as a second coolant should take such a permitted time into consideration.

divertor↗
DOE OSTI2025

Dark Energy Survey Year 6 Results: Point-Spread Function Modeling

We present the point-spread function (PSF) modeling for weak lensing shear measurement using the full six years of the Dark Energy Survey (DES Y6) data. We review the PSF estimation procedure using the PIFF (PSFs In the Full FOV) software package and describe the key improvements made to PIFF and modeling diagnostics since the DES year three (Y3) analysis: (i) use of external Gaia and infrared photometry catalogs to ensure higher purity of the stellar sample used for model fitting, (ii) addition of color-dependent PSF modeling, the first for any weak lensing analysis, and (iii) inclusion of model diagnostics inspecting fourth-order moments, which can bias weak lensing measurements to a similar degree as second-order modeling errors. Through a comprehensive set of diagnostic tests, we demonstrate the improved accuracy of the Y6 models evident in significantly smaller systematic errors than those of the Y3 analysis, in which all g band data were excluded due to insufficiently accurate PSF models. For the Y6 weak lensing analysis, we include g band photometry data in addition to the riz bands, providing a fourth band for photometric redshift estimation. Looking forward to the next generation of wide-field surveys, we describe several ongoing improvements to PIFF, which will be the default PSF modeling software for weak lensing analyses for the Vera C. Rubin Observatory’s Legacy Survey of Space and Time.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Time-dependent water vapor desorption isotherm model of hardened cement paste

Highlights: • Short-term and long-term desorption isotherms of hardened cement paste are reviewed. • Water content at RH depends on both thermodynamic and microstructural equilibria. • An extended GAB equation is proposed to fit the cavitation drop during desorption. • A dynamic time-dependent sorption isotherm model is shown and used for simulations. • Numerical results capture anomalous drying in the same RH as seen experimentally. This study aims to enhance moisture transport modeling by elucidating the so-called anomalous water diffusion in cementitious materials. Water desorption isotherms are measured from samples at various drying stages using techniques with different durations to quantify the impact of the drying duration on the water sorption ability. A single water sorption isotherm does not solely give a relationship between the water content and relative humidity but also represents the state of microstructure. The continuous evolution of a desorption isotherm due to drying-induced microstructural rearrangement is demonstrated. For numerical modeling, the microstructural alteration can be explicitly considered through a dynamic desorption isotherm model, which governs the local thermodynamic equilibrium at the capillary meniscus. This approach is implemented into a multiphase transport code, whose ability to predict drying is validated using literature data. Finally, the effect of prolonged drying on the colloidal nature of the calcium-silicate-hydrate gel is discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Large-Eddy Simulations of Idealized Atmospheric Boundary Layers Using Nalu-Wind

Accurate prediction of wind-plant performance relies, in part, on properly characterizing the turbulent atmospheric boundary layer (ABL) flow in which wind turbines operate. Large-eddy simulation (LES) is a powerful tool for simulating ABLs because it resolves the largest, most energetic scales of three-dimensional turbulent motions. Yet LES predictions are well known to depend on modeling choices such as grid resolution, numerical discretization schemes, and closures for unresolved scales of turbulence. Here, we evaluate how these choices influence predictions of ABL winds using Nalu-Wind, a wind-specific fork of the open-source, generalized, unstructured, massively parallel flow solver NaluCFD/Nalu.

17 WIND ENERGY↗
DOE OSTI2020

Distributionally robust facility location problem under decision-dependent stochastic demand

While the traditional facility location problem considers exogenous demand, in some applications, locations of facilities could affect the willingness of customers to use certain types of services, e.g., carsharing, and therefore they also affect realizations of random demand. Moreover, a decision maker may not know the exact distribution of such endogenous demand and how it is affected by location choices. In this paper, we consider a distributionally robust facility location problem, in which we interpret the moments of stochastic demand as functions of facility-location decisions. We reformulate a two-stage decision-dependent distributionally robust optimization model as a monolithic formulation, and then derive exact mixed-integer linear programming reformulation as well as valid inequalities when the means and variances of demand are piecewise linear functions of location solutions. We conduct extensive computational studies, in which we compare our model with a decision-dependent deterministic model, as well as stochastic programming and distributionally robust models without the decision-dependent assumption. Here, the results show superior performance of our approach with remarkable improvement in profit and quality of service under various settings, in addition to computational speed-ups given by formulation enhancements. These results draw attention to the need of considering the impact of location decisions on customer demand within this strategic-level planning problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Reynolds stress tensor measurements using magnetic resonance velocimetry: expansion of the dynamic measurement range and analysis of systematic measurement errors

This study presents magnetic resonance velocimetry (MRV) Reynolds Stress measurements in a periodic hill channel with a hill Reynolds number of Re = 29,500. The velocity encoding scheme is based on the ICOSA6 method with six icosahedral encoding directions and multiple encoding values are measured to increase the dynamic range. The full Reynolds stress tensor is obtained from a voxel-wise three-dimensional Gaussian fit using the magnitude data of all acquisitions. The MRV results are compared to a wall-resolved large eddy simulation and laser Doppler velocimetry measurements conducted in the same channel. It is shown that the MRV Reynolds stress data have excellent precision and agree qualitatively with the reference data. However, there are apparent systematic deviations. One of the most prominent error contributions is the signal attenuation caused by higher orders of motion, which leads to an overestimation of the turbulence level. Another fundamental error is identified in the assumption that the turbulence is Gaussian distributed. With the presented reconstruction technique, the MRV data are fitted to a statistical model, and depending on the examined flow setup, the Gaussian model can lead to considerable errors. Possible ways of how to reduce all identified errors are presented. In summary, this technique enables Reynolds stress tensor measurements in complex internal flows with high dynamic range and excellent precision. However, several issues need to be resolved to make the turbulence quantification more accurate.

42 ENGINEERING↗
DOE OSTI2022

Modeling the Conformer-Dependent Electronic Absorption Spectra and Photolysis Rates of Methyl Vinyl Ketone Oxide and Methacrolein Oxide

Criegee intermediates are important atmospheric oxidants, formed via the reaction of ozone with volatile alkenes emitted into the troposphere. Small Criegee intermediates (e.g., CH 2 OO and CH 3 CHOO) are highly reactive, and their removal via unimolecular decay or bimolecular chemistry dominates their atmospheric lifetimes. As the molecular complexity of Criegee intermediates increases, their electronic absorption spectra show a bathochromic shift within the solar spectrum relevant to the troposphere. In such cases, solar photolysis may become a competitive contributor to their atmospheric removal. In this article, we report the conformer-dependent simulated electronic absorption spectra of two four-carbon-centered Criegee intermediates, methyl vinyl ketone oxide (MVK-oxide) and methacrolein oxide (MACR-oxide). Both MVK-oxide and MACR-oxide contain four low-energy conformers, which are convoluted in the experimentally measured spectra. Here, we deconvolute each conformer and estimate contributions from each of the four conformers to the experimentally measured spectra. We also estimate the photolysis rates and predict that solar photolysis should be a more competitive removal process for MVK-oxide and MACR-oxide (cf. CH 2 OO and CH 3 CHOO).

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

Atomistically resolved hot exciton relaxation dynamics in CdSe quantum dots: Experiment and theory

Semiconductor quantum dots (QDs) are well known to give rise to a quantum confined structure of excitons. Because of this quantum confinement, new physics of hot exciton relaxation dynamics arises. Decades of work using transient absorption (TA) spectroscopy have yielded initial simple observations, such as estimates of the cooling rate from single pump photon energy experiments. More detailed TA experiments employed variable pump photon energies to measure excitonic state-resolved transition rates. These TA measurements, usually the simplest form, have been employed to characterize QDs and their relaxation dynamics to this day. Yet, these TA measurements are fundamentally lacking in their ability to measure energy-resolved hot exciton cooling, which requires observation of the full cooling history through the real excitonic manifold. Here, we employ coherent multi-dimensional spectroscopy (CMDS) to perform an atomistically directed study of hot exciton cooling in CdSe QDs, revealing energy resolved relaxation dynamics. CMDS experiments are compared with simulations and prior TA measurements and simpler theories. Our findings reveal a hot exciton relaxation dynamics landscape. This relaxation dynamics landscape is a linear or sub-linear function of excess energy for different structures of QDs, with a strong size dependence. Our model simulations parameterized by the empirical pseudopotential model reproduces the experimental functional form and the dependence upon QD diameter and shell.

Atomistic simulations↗