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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

An Efficient Reactive Force Field without Explicit Coordination Dependence for Studying Caustic Aluminum Chemistry

Reactive force fields (RFFs) are an expedient approach to sample chemical reaction paths in complex systems, relative to density functional theory. However, there is continued need to improve efficiencies, specifically in systems that have slow transverse degrees of freedom, as in highly viscous and superconcentrated solutions. Here, we present an RFF that is differentiated from current models (e.g., ReaxFF) by omitting explicit dependence on the atom coordination and employing a small parameter set based on Lennard-Jones, Gaussian, and Stillinger–Weber potentials. In this study, the model was parametrized from AIMD simulation data and is used to model aluminate reactivity in sodium hydroxide solutions with extensive validation against experimental radial distribution functions, computed free energy profiles for oligomerization, and formation energies. In conclusion, the model enables simulation of early stage Al(OH) 3 nucleation which has significant relevance to industrial processing of aluminum and has a computational cost that is reduced by 1 order of magnitude relative to ReaxFF.

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

The Role of Phytoplankton Biomacromolecules in Controlling Ocean Surface Roughness

Satellite altimetric data routinely map sea surface topography by measuring the ocean return signal. One source of altimeter measurement contamination occurs when the radar ocean backscatter becomes unusually large, a situation termed a Sigma-0 bloom. Past research suggests Sigma-0 blooms are associated with weak wind and natural surface slick conditions where capillary waves at the air–sea interface are suppressed. To date, no explicit connection between these conditions and Sigma-0 bloom presence has been provided. Using a series of simplified equations, our reduced model determines capillary wave heights from estimates of planktonic carbon concentrations and regional wind speed. Our results suggest that the radar signal reflection increases as capillary wave height decreases. This relationship depends on surfactant concentration, surfactant composition, and wind speed. Model sensitivity analysis indicates that the interface reflectivity depends on biological activity and wind speed. Our proposed simplified model provides a method to identify potential Sigma-0 bloom regions. We conclude that because of the demonstrated impact of biological surfactants on ocean roughness, it is necessary to consider the biological activity, i.e., phytoplankton bloom events, when interpreting signals from radar altimetry and when developing ocean hydrology models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Phylogenetic distribution and experimental characterization of corrinoid production and dependence in soil bacterial isolates

Abstract Soil microbial communities impact carbon sequestration and release, biogeochemical cycling, and agricultural yields. These global effects rely on metabolic interactions that modulate community composition and function. However, the physicochemical and taxonomic complexity of soil and the scarcity of available isolates for phenotypic testing are significant barriers to studying soil microbial interactions. Corrinoids—the vitamin B12 family of cofactors—are critical for microbial metabolism, yet they are synthesized by only a subset of microbiome members. Here, we evaluated corrinoid production and dependence in soil bacteria as a model to investigate the ecological roles of microorganisms involved in metabolic interactions. We isolated and characterized a taxonomically diverse collection of 161 soil bacteria from a single study site. Most corrinoid-dependent bacteria in the collection prefer B12 over other corrinoids, while all tested producers synthesize B12, indicating metabolic compatibility between producers and dependents in the collection. Furthermore, a subset of producers release B12 at levels sufficient to support dependent isolates in laboratory culture at estimated ratios of up to 1000 dependents per producer. Within our isolate collection, we did not find strong phylogenetic patterns in corrinoid production or dependence. Upon investigating trends in the phylogenetic dispersion of corrinoid metabolism categories across sequenced bacteria from various environments, we found that these traits are conserved in 47 out of 85 genera. Together, these phenotypic and genomic results provide evidence for corrinoid-based metabolic interactions among bacteria and provide a framework for the study of nutrient-sharing ecological interactions in microbial communities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Analysis of SPND Material Candidates for Real-Time Neutron Spectrum Monitoring Poster

The future of the nuclear industry will be shaped by advanced reactor designs requiring innovative instrumentation for safe operation. A key area of innovation is the ability to accurately and reliably monitor the neutron energy spectrum during operations, which would also enhance the experimental capabilities of existing research reactors. To address this need, our research focuses on the potential use of Self-Powered Neutron Detectors (SPNDs) for real-time neutron spectrum monitoring. We investigated various candidate materials for SPNDs by evaluating their energy-dependent reaction rates. Promising materials were further modeled using MCNP to determine the free electron production per neutron across the neutron energy spectrum, approximating their energy-dependent signal sensitivity. Our findings identified Samarium (Sm) and Gadolinium (Gd) as top candidates for thermal neutrons, Hafnium (Hf) and Thulium (Tm) for epithermal neutrons, and Tantalum (Ta) for fast neutrons. These materials will undergo future experimental validation to confirm the modeled responses and assess their feasibility for real-time neutron spectrum measurements.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗
DOE OSTI2020

Lattice modeling of excavation damage in argillaceous clay formations: Influence of deformation and strength anisotropy

This paper presents modeling of mechanical anisotropy in argillaceous rocks using an irregular lattice modeling approach, namely the rigid-body-spring network. To represent the mechanical anisotropy, new schemes are implemented in the modeling framework. The directionality of elastic deformation is resolved by modifying the element formulation with anisotropic elastic properties. The anisotropy of strength and failure characteristics is facilitated by adopting orientation-dependent failure criteria into the failure model. The verification of the improved modeling procedures is performed against theoretical model predictions for unconfined compression tests with various bedding orientations. Furthermore, excavation damage and fracturing processes in rock formations are simulated for different geomechanical configurations, such as rock anisotropy and tectonic heterogeneity. The simulated excavation damage characteristics are realistic and comparable with the actual field observation at a tunnel located in an argillaceous clay formation. The simulation results provide insights into the excavation damage zone phenomena with an explicit representation of fracturing processes.

58 GEOSCIENCES↗
DOE OSTI2024

Value-added catalog of physical properties for more than 1.3 million galaxies from the DESI survey

We present an extensive catalog of the physical properties of more than a million galaxies investigated with the Dark Energy Spectroscopic Instrument (DESI), one of the largest spectroscopic surveys to date. Spanning a full range of target types, including emission-line galaxies, luminous red galaxies, and quasars, our survey encompasses an unprecedented range of spectroscopic redshifts, all the way from 0 to 6. The physical properties, such as stellar masses and star formation rates, were derived via the CIGALE spectral energy distribution (SED) fitting code accounting for the contribution coming from active galactic nuclei (AGNs). Based on the modeling of the optical-mid-infrared (grz supplemented with WISE photometry) SEDs, we studied the galaxy properties with respect to their location on the main sequence. We have revised the dependence of stellar mass estimates on model choices and on the availability of WISE photometry. Indeed, the WISE data are required to minimize the misclassification of star-forming galaxies as AGNs. The lack of WISE bands in SED fits leads to elevated AGN fractions for 68% of star-forming galaxies identified using emission line diagnostic diagrams, but this does not significantly affect their stellar mass or star formation estimates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Ab initio molecular dynamics (AIMD) simulations of NaCl, UCl 3 molten salts

Ab initio molecular dynamics (AIMD) simulations are used to calculate select thermophysical and thermodynamic properties of NaCl, UCl and NaCl-UCl 3 molten salts. Following established approaches, the AIMD simulations include a model for Van der Waals interactions. The Langreth & Lundqvist (vDW-DF2), DFT-D3, and density-dependent energy correction (DFT-dDsC) dispersion models are first tested for molten NaCl in order to assess their accuracy for density and heat capacity predictions across a range of temperatures. Based on the NaCl results, the vdW-DF2 and DFT-dDsC methods are extended to the UCl 3 system and compared to available experimental data. Next, mixtures of NaCl-UCl 3 are investigated and analyzed with respect to the deviation from ideal solution behavior. For the DFT-dDsC simulations, density deviates by up to 2% from an ideal mixture, with the maximum occurring close to the eutectic composition. The mixing energy also deviates from an ideal solution and exhibits a minimum of about -0.074 to eV per formula unit, again close to the eutectic composition. Finally, the compressibility and species diffusivity of the pure and mixed salt systems are calculated. The diffusivities are slightly reduced in the mixed compared to the pure systems and the compressibilities loosely follow a linear correlation as a function of the UCl 3 composition. The trends observed for mixing properties are rationalized by correlating them to the coordination chemistry, which emphasizes the importance of maintaining the extended network formed by U and Cl ions as the NaCl concentration increases. The concentration at which breakdown of the extended network occurs, roughly coincides with the minimum of both the mixing energy and the deviation from ideal solution behavior for density.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Investigating the Diurnal Radiative, Turbulent, and Biophysical Processes in the Amazonian Canopy‐Atmosphere Interface by Combining LES Simulations and Observations

Abstract We investigate the atmospheric diurnal variability inside and above the Amazonian rainforest for a representative day during the dry season. To this end, we combine high‐resolution large‐eddy simulations that are constrained and evaluated against a comprehensive observation set, including CO 2 concentrations, gathered during GoAmazon2014/15. We design systematic numerical experiments to quantify whether a multilayer approach in solving the explicit canopy improves our canopy‐atmosphere representation. We particularly focus on the relationship between photosynthesis and plant transpiration, and their distribution at leaf and canopy scales. We found the variability of photosynthesis drivers like vapor pressure deficit and leaf temperature to be about 3 times larger for sunlit leaves compared to shaded leaves. This leads to a large spread on leaf stomatal conductance values with minimum and maximum values varying more than 100%. Regarding the turbulent structure, we find wind‐driven stripe‐like shapes at the canopy top and structures resembling convective cells at the canopy. Wind‐related variables provide the best spatiotemporal agreement between model and observations. The potential temperature and heat flux profiles agree with an observed decoupling near the canopy top interface, although with less variability and cold biases of up to 3 K. The increasing complexity on the biophysical processes leads to the largest disagreements for evaporation, CO 2 plant assimilation and soil efflux. The model is able to capture the correct dependences and trends with the magnitudes still differing. We finally discuss the need to revise leaf and soil models and to complete the observations at leaf and canopy levels.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Lifetime renormalization of weakly anharmonic superconducting qubits. I. Role of number nonconserving terms

The dynamics of a weakly anharmonic superconducting qubit in a complex electromagnetic environment is generally well described by an effective multimode Kerr Hamiltonian at sufficiently weak excitation. This Hamiltonian can be embedded in a master equation with losses determined by the details of the electromagnetic environment. Recent experiments indicate, however, that when a superconducting circuit is driven with microwave signals populating the system with sufficiently high excitations, the observed relaxation rates appear to be substantially different from expectations based on the electromagnetic environment of the qubit alone. This issue is a limiting factor in the optimization of superconducting qubit readout schemes. We claim here that an effective master equation with drive-power-dependent parameters is an efficient approach to model such quantum dynamics. In this sequence of papers, we derive effective master equations, whose parameters exhibit nonlinear dependence on the excitation level of the circuit as well as the electromagnetic environment of the qubit. Further, we show that the number nonconserving terms in the qubit nonlinearity generally lead to a renormalization of dissipative parameters of the effective master equation, while the number conserving terms give rise to a renormalization of the system frequencies. Here, in Paper I, we consider the excitation-relaxation dynamics of a transmon qubit that is prepared in a certain initial state, but is not driven otherwise. A unitary transformation technique is introduced to study the renormalization of (i) qubit relaxation due to coupling to a generic bath and (ii) Purcell decay. Analytic expressions are provided for the dependence of the nonlinear dissipative terms on the details of the electromagnetic environment of the qubit. The perturbation technique based on unitary transformations developed here is generalized to the continuously driven case in Paper II.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Microstructural and infiltration properties of woven preforms during chemical vapor infiltration

Interface-resolved direct numerical simulations (DNSs) of chemical vapor infiltration (CVI) have been performed over a range of furnace-operating conditions (Thiele moduli) and for practical woven preform geometries. A level-set method is used to resolve the geometry of the initial preform at tow scale. The interface between the vapor and solid phase is then evolved in time through the entire CVI densification cycle, fully resolving the time-varying topology between the two phases. In contrast to previous level-set methods for CVI simulation, the physical reaction and diffusion processes govern the level-set movement in the current approach. The surface deposition kinetics is described by the usual one-step model. In this paper, the DNS data are used to study the evolving porosity, surface-to-volume ratio, and flow infiltration properties (permeability and effective diffusivities). Additionally, comparisons are made to popularly-assumed structure functions and the standard, Kozeny–Carmen porous media model commonly employed in modeled CFD simulations of CVI. The virtual DNS experiments reveal a Thiele modulus and preform geometry (fabric layup) dependence which the existing microstructural and infiltration models are not able to describe throughout the entire densification process. The DNS-based, woven geometry-specific correlations can be applied directly to mean-field, furnace-scale CFD simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Machine learning predictions for local electronic properties of disordered correlated electron systems

We present a scalable machine learning (ML) model to predict local electronic properties such as on-site electron number and double occupation for disordered correlated electron systems. Our approach is based on the locality principle, or the nearsightedness nature, of many-electron systems, which means local electronic properties depend mainly on the immediate environment. A ML model is developed to encode this complex dependence of local quantities on the neighborhood. We demonstrate our approach using the square-lattice Anderson-Hubbard model, which is a paradigmatic system for studying the interplay between Mott transition and Anderson localization. We develop a lattice descriptor based on the group-theoretical method to represent the on-site random potentials within a finite region. The resultant feature variables are used as input to a multilayer fully connected neural network, which is trained from data sets of variational Monte Carlo (VMC) simulations on small systems. We show that the ML predictions agree reasonably well with the VMC data. Our work underscores the promising potential of ML methods for multiscale modeling of correlated electron systems.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Modeling of gap conductance for LWR fuel rods applied in the BISON code

A new gap conductance model is proposed in this study as a combination of Toptan’s model and the Ross-Stoute model. A variance-based sensitivity analysis is performed to understand how simulation results depend on all input parameters of the proposed model. Additionally, new modeling options (e.g. fill gas thermal conductivity, temperature jump distance, thermal accommodation coefficient, etc.) are added into the nuclear fuel performance code, BISON. The need for further investigation of the gap heat transfer between fuel and cladding in BISON motivated this study to evaluate its impact on the code’s predictions. New gap conductance modeling is proposed. A series of integral-effects validation tests is performed: (1) to demonstrate the impact of the proposed model on the code’s fuel temperature predictions at the beginning of life and through the reactor’s life; (2) to ensure that the proposed model is capable of accurately modeling gap heat transfer characteristics in real-world problems; and (3) to investigate the impact of the estimation of fission gas release on the fuel temperature predictions with the proposed model. Overall, the results indicate that the proposed gap conductance model improves BISON’s predictions.

42 ENGINEERING↗
DOE OSTI2020

Pseudospin versus magnetic dipole moment ordering in the isosceles triangular lattice material K 3 Er(VO 4 ) 2

Spin- 1 / 2 antiferromagnetic triangular lattice models are paradigms of geometrical frustration, revealing very different ground states and quantum effects depending on the nature of anisotropies in the model. Due to strong spin orbit coupling and crystal field effects, rare-earth ions can form pseudospin- 1 / 2 magnetic moments with anisotropic single-ion and exchange properties. Thus, rare-earth-based triangular lattices enable the exploration of this interplay between frustration and anisotropy. Here we study one such case, the rare-earth double vanadate glaserite material K 3 Er ( VO 4 ) 2 , which is a quasi-two-dimensional (2D) isosceles triangular antiferromagnet. Our specific heat and neutron powder diffraction data from K 3 Er ( VO 4 ) 2 reveal a transition to long range magnetic order at T N = 155 ± 5 mK which accounts for all R ln 2 entropy. Furthermore, we observe what appears to be a coexistence of three-dimensional (3D) and quasi-2D order below T N . The quasi-2D order leads to an anisotropic Warren-like peak profile for ( h k 0 ) reflections, while the 3D order is best-described by layers of antiferromagnetic b -aligned moments alternating with layers of zero moment. Our magnetic susceptibility data reveal that Er 3 + takes on a strong X Y single-ion anisotropy in K 3 Er ( VO 4 ) 2 , leading to vanishing moments when pseudospins are oriented along c . Thus, the magnetic structure, when considered from the pseudospin point of view could comprise of alternating layers of b -axis and c -axis aligned antiferromagnetism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Evaluating Recombination Mechanisms in RbF Treated Cu(In${}_\mathrm{x}$Ga$_\mathrm{1-x}$)Se$_{2}$ Solar Cells

Rubidium fluoride (RbF) postdeposition treatment (PDT) has been shown to improve the performance of Cu(In x Ga 1-x )Se 2 (CIGS) photovoltaic devices. Here, in this study, temperature-dependent current voltage (JVT) and time-resolved photoluminescence (TRPL) experiments were combined with modeling using the solar cell capacitance simulator (SCAPS) computer code to investigate the effect of the RbF PDT. Two devices, one as-deposited and one with RbF PDT, were deposited by a three stage coevaporation process. JVT measurements suggest the dominant recombination mechanism may be tunneling-enhanced recombination via bandtail states, but that defect states in the bandgap can also be important. RbF PDT is shown to decrease the characteristic energy of the bandtails. TRPL data show an increase in the minority carrier lifetime after RbF PDT, leading to an improved open-circuit voltage. SCAPS modeling indicates that the dominant recombination mechanism is dependent on the specific defect makeup of a device, suggesting that small changes in processing conditions can impact device behavior. This explains the observation that, for some devices, defect states in the gap dominate while others, as is the case here, appear to be dominated by bandtails.

14 SOLAR ENERGY↗
DOE OSTI2024

Skeletal Kinetics Reduction for Astrophysical Reaction Networks

A novel methodology is developed to extract accurate skeletal reaction models for nuclear combustion. Local sensitivities of isotope mass fractions with respect to reaction rates are modeled based on the forced optimally time-dependent (f-OTD) scheme. These sensitivities are then analyzed temporally to generate skeletal models. The methodology is demonstrated by conducting skeletal reduction of constant density and temperature burning of carbon and oxygen relevant to Type Ia supernovae (SNe Ia). The 495-isotopes Torch model is chosen as the detailed reaction network. A map of maximum production of 56 Ni in SNe Ia is produced for different temperatures, densities, and proton-to-neutron ratios. The f-OTD simulations and the sensitivity analyses are then performed with initial conditions from this map. A series of skeletal models are derived and their performances are assessed by comparison against currently existing skeletal models. Previous models have been constructed intuitively by assuming the dominance of α-chain reactions. The comparison of the newly generated skeletal models against previous models is based on the predicted energy release and 44 Ti and 56 Ni abundances by each model. The consequences of ye ≠ 0.5 in the initial composition are also explored where ye is the electron fraction. The simulated results show that 56 Ni production decreases by decreasing ye as expected, and that the 43 Sc is a key isotope in proton and neutron channels toward 56 Ni production. It is shown that an f-OTD skeletal model with 150 isotopes can accurately predict the 56Ni abundance in SNe Ia for ye ≲ 0.5 initial conditions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Predicting fracture behavior in single crystal nickel using a coupled crystal plasticity phase field damage approach

Understanding the fracture behavior of single crystal metals is critical for predicting material performance under mechanical loading. Here, in this study, we investigate the fracture characteristics of single crystal nickel tensile bars using a crystal plasticity coupled phase field damage (CP-PFD) model. Experimental tensile tests were conducted on 15 specimens spanning five crystallographic orientations and three thickness variants per orientation. The results revealed two distinct fracture modes: brittle fractures with 45-degree angled surfaces and ductile fractures characterized by significant necking. The CP-PFD model successfully replicated these fracture behaviors, demonstrating strong agreement with experimental observations. The model effectively predicted the strain at which necking and fracture occurred, as well as the orientation-dependent fracture mechanisms. By comparing experimental and simulated fracture surfaces, we establish the CP-PFD model as a robust tool for predicting single crystal behavior and damage evolution. This work provides insight into the microstructural dependence of fracture behavior and establishes a predictive framework for modeling orientation-dependent damage evolution in single-crystal nickel.

Crystal plasticity↗
DOE OSTI2020

Volatilization Behavior of Supported Au Nanoparticle Arrays under H 2 at High Temperature

In this work, the fabrication and experimental study of model Au nanoparticle arrays (nanodots) was used to allow a direct comparison of reactive evaporation at elevated temperatures and nonambient H 2 -containing gas atmospheres with theoretical models. Strong temperature and H 2 -concentration dependencies were confirmed as expected, and the observed dependence in particle diameter with time was effectively described using a reactive evaporation and gas phase diffusion based quantitative model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗