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

Adaptive Online Model Update Algorithm for Predictive Control in Networked Systems

In this article, we introduce an adaptive on-line model update algorithm designed for predictive control applications in networked systems, particularly focusing on power distribution systems. Unlike traditional methods that depend on historical data for offline model identification, our approach utilizes real-time data for continuous model updates. This method integrates seamlessly with existing online control and optimization algorithms and provides timely updates in response to real-time changes. This methodology offers significant advantages, including a reduction in the communication network bandwidth requirements by minimizing the data exchanged at each iteration and enabling the model to adapt after disturbances. Furthermore, our algorithm is tailored for non-linear convex models, enhancing its applicability to practical scenarios. The efficacy of the proposed method is validated through a numerical study, demonstrating improved control performance using a synthetic IEEE test case.

data-driven model predictive control↗
DOE OSTI2022

Nramp: Deprive and conquer?

Solute carriers 11 (Slc11) evolved from bacterial permease (MntH) to eukaryotic antibacterial defense (Nramp) while continuously mediating proton (H+)-dependent manganese (Mn2+) import. Also, Nramp horizontal gene transfer (HGT) toward bacteria led to mntH polyphyly. Prior demonstration that evolutionary rate-shifts distinguishing Slc11 from outgroup carriers dictate catalytic specificity suggested that resolving Slc11 family tree may provide a function-aware phylogenetic framework. Hence, MntH C (MC) subgroups resulted from HGTs of prototype Nramp (pNs) parologs while archetype Nramp (aNs) correlated with phagocytosis. PHI-Blast based taxonomic profiling confirmed MntH B phylogroup is confined to anaerobic bacteria vs. MntH A (MA)’s broad distribution; suggested niche-related spread of MC subgroups; established that MA-variant MH, which carries ‘eukaryotic signature’ marks, predominates in archaea. Slc11 phylogeny shows MH is sister to Nramp. Site-specific analysis of Slc11 charge network known to interact with the protonmotive force demonstrates sequential rate-shifts that recapitulate Slc11 evolution. 3D mapping of similarly coevolved sites across Slc11 hydrophobic core revealed successive targeting of discrete areas. The data imply that pN HGT could advantage recipient bacteria for H+-dependent Mn2+ acquisition and Alphafold 3D models suggest conformational divergence among MC subgroups. It is proposed that Slc11 originated as a bacterial stress resistance function allowing Mn2+-dependent persistence in conditions adverse for growth, and that archaeal MH could contribute to eukaryogenesis as a Mn2+ sequestering defense perhaps favoring intracellular growth-competent bacteria.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

SABATH: Surrogate AI Benchmarking Applications' Testing Harness

SABATH provides benchmarking infrastructure for evaluating scientific ML/AI models. It contains support for scientific machine learning surrogates from external repositories. The software dependences are explicitly exposed in the surrogate model definition, which allows the use of advanced optimization, communication, and hardware features.

Luszczek, Piotr [Univ. of Tennessee, Knoxville, TN↗
DOE OSTI2020

Evaluating the Potential and Challenges of an Uncertainty Quantification Method for Long Short–Term Memory Models for Soil Moisture Predictions

Recently, recurrent deep networks have shown promise to harness newly available satellite–sensed data for long–term soil moisture projections. However, to be useful in forecasting, deep networks must also provide uncertainty estimates. Here we evaluated Monte Carlo dropout with an input–dependent data noise term (MCD+N), an efficient uncertainty estimation framework originally developed in computer vision, for hydrologic time series predictions. MCD+N simultaneously estimates a heteroscedastic input–dependent data noise term (a trained error model attributable to observational noise) and a network weight uncertainty term (attributable to insufficiently constrained model parameters). Although MCD+N has appealing features, many heuristic approximations were employed during its derivation, and rigorous evaluations and evidence of its asserted capability to detect dissimilarity were lacking. To address this, we provided an in–depth evaluation of the scheme's potential and limitations. We showed that for reproducing soil moisture dynamics recorded by the Soil Moisture Active Passive (SMAP) mission, MCD+N indeed gave a good estimate of predictive error, provided that we tuned a hyperparameter and used a representative training data set. The input–dependent term responded strongly to observational noise, while the model term clearly acted as a detector for physiographic dissimilarity from the training data, behaving as intended. However, when the training and test data were characteristically different, the input–dependent term could be misled, undermining its reliability. Additionally, due to the data–driven nature of the model, data noise also influences network weight uncertainty, and therefore the two uncertainty terms are correlated. Altogether, this approach has promise, but care is needed to interpret the results.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Atmospheric River Response to Arctic Sea Ice Loss in the Polar Amplification Model Intercomparison Project

The atmospheric river (AR) response to Arctic sea ice loss in the Northern hemisphere winter is investigated using simulations from the Polar Amplification Model Intercomparison Project. Results have shown that the midlatitude responses are dominated by dynamic effects. Poleward of around 60°N, the dynamic and thermodynamic effects cancel each other, resulting in relatively small responses. The response uncertainty can be characterized by leading uncertainty modes, with the responses over the Pacific and Atlantic projecting onto the northeastward extension and equatorward shift mode, respectively. In addition, the responses seem to be mean state-dependent: under the same forcing, models with more poleward-located climatological ARs tend to show stronger equatorward shifts over the Atlantic; over the Pacific, models with more westward-located climatological AR core tend to show stronger northeastward extensions. Finally, these relationships highlight the importance of improving the AR climatology representation on reducing the response uncertainty to Arctic sea ice loss.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

On plasticity-enhanced interfacial toughness in bonded joints

Here, the performance and reliability of many structures and components depend on the integrity of interfaces between dissimilar materials. Interfacial toughness Γ is the key material parameter that characterizes resistance to interfacial crack growth, and Γ is known to depend on many factors including temperature. For example, previous work showed that the toughness of an epoxy/aluminum interface decreased 40 % as the test temperature was increased from -60 °C to room temperature (RT). Interfacial integrity at elevated temperatures is of considerable practical importance. Recent measurements show that instead of continuing to decrease with increasing temperature, Γ increases when test temperature is above RT. Cohesive zone finite element calculations of an adhesively bonded, asymmetric double cantilever beam specimen of the type used to measure Γ suggest that this increase in toughness may be a result of R-curve behavior generated by plasticity-enhanced toughening during stable subcritical crack growth with interfacial toughness defined as the critical steady-state limit value. In these calculations, which used an elastic-perfectly plastic epoxy model with a temperature-dependent yield strength, the plasticity-enhanced increase in Γ above its intrinsic value Γ o depended on the ratio of interfacial strength σ* to the yield strength σyb of the bond material. There is a nonlinear relationship between Γ/Γ o and σ*/σ yb with the value Γ/Γ o increasing rapidly above a threshold value of σ*/σ yb . The predicted increase in toughness can be significant. For example, there is nearly a factor of two predicted increase in Γ/Γ o during micrometer-scale crack-growth when σ*/σ yb = 2 (a reasonable choice for σ*/σ yb ). Furthermore, contrary to other reported results, plasticity-enhanced toughening can occur prior to crack advance as the cohesive zone forms and the peak stress at the tip of the original crack tip translates to the tip of the fully formed cohesive zone. These results suggest that plasticity-enhanced toughening should be considered when modeling interfaces at elevated temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Pair momentum dependence of a tilted source in heavy-ion collisions

In noncentral heavy-ion collisions, the particle-emitting source can be tilted away from the beam direction, an effect that becomes particularly significant at collision energies of a few GeV and lower. This phenomenon, manifest itself in many observables such as directed flow, polarization, and vorticity, is therefore important to investigate. Here, in this paper, we study the consistency between the tilt extracted directly from the freeze-out distribution of pions and the tilt parameter obtained using the azimuthally sensitive femtoscopy (asHBT) method. Using the Ultrarelativistic Quantum Molecular Dynamics model, we demonstrate a strong dependence of the tilt parameter extracted with asHBT on the momentum of the particle pair. Considering the experimental challenges in accessing low particle momenta—where the tilt parameter extracted with asHBT closely matches the tilt of the freeze-out distribution of pions—we propose an exponential extrapolation method to obtain the tilt of the entire freeze-out distribution. This approach aims to enhance the accuracy of experimental measurements of tilt in noncentral heavy-ion collisions.

particle correlations & fluctuations↗
DOE OSTI2023

Tail Dependence as a Measure of Teleconnected Warm and Cold Extremes of North American Wintertime Temperatures

Current models for spatial extremes are concerned with the joint upper (or lower) tail of the distribution at two or more locations. Such models cannot account for teleconnection patterns of 2-m surface air temperature ( T 2m ) in North America, where very low temperatures in the contiguous United States may coincide with very high temperatures in Alaska in the wintertime. This dependence between warm and cold extremes motivates the need for a model with opposite-tail dependence in spatial extremes. This work develops a statistical modeling framework that has flexible behavior in all four pairings of high and low extremes at pairs of locations. In particular, we use a mixture of rotations of common Archimedean copulas to capture various combinations of four-corner tail dependence. We study teleconnected T 2m extremes using ERA5 of daily average 2-m temperature during the boreal winter. Further, the estimated mixture model quantifies the strength of opposite-tail dependence between warm temperatures in Alaska and cold temperatures in the midlatitudes of North America, as well as the reverse pattern. These dependence patterns are shown to correspond to blocked and zonal patterns of midtropospheric flow. This analysis extends the classical notion of correlation-based teleconnections to considering dependence in higher quantiles.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Simulating Sea-Ice Deformation in Viscous-Plastic Sea-Ice Models With CD-Grids

Linear kinematic features (LKFs) are found everywhere in the Arctic sea-ice cover. They are strongly localized deformations often associated with the formation of leads and pressure ridges. In viscous-plastic (VP) sea-ice models, the simulation of LKFs depends on several factors such as the grid resolution, the numerical solver convergence, and the placement of the variables on the mesh. In this study, we compare two recently proposed discretization with a CD-grid placement with respect to their ability to reproduce LKFs. The first (CD1) is based on a nonconforming finite element discretization, whereas the second (CD2) uses a conforming subgrid discretization. To analyze their resolution properties, we evaluate runs from different models (e.g., FESOM, MPAS) on a benchmark problem using quadrilateral, hexagonal and triangular meshes. Our findings show that the CD1 setup simulates more deformation structure than the CD2 setup. This highlights the importance of the type of spatial discretization for the simulation of LKFs. Due to the higher number of degrees of freedom, both CD-grids resolve more LKFs than traditional A, B, and C-grids at fixed mesh level. This is an advantage of the CD-grid approach, as high spatial mesh resolution is needed in VP sea-ice models to simulate LKFs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Impacts of Pasture Conversion to Sugarcane on Water Fluxes and Water Use Efficiency in the Southeastern US

The expansion of sugarcane (cane), a high-yielding perennial crop, will likely reshape the bioenergy landscape in the Southeastern US. However, its ecohydrological implications, particularly following conversion from grazed pastures, a dominant land use in the region, remain highly uncertain. We investigated the impact of cane expansion on evapotranspiration (ET) and its partitioning, and the mechanisms influencing both ET components and water use efficiency (WUE) across multiple scales and growth cycles in subtropical Florida. We combined eddy covariance, biometric measurements, and process-based stomatal conductance (g s ) models. ET was 1.7% lower in cane than in improved pasture (IMP) but exceeded that in semi-native pasture (SN) by 21%. Transpiration (T) followed a similar pattern, consistent with lower g s in cane relative to IMP. Cane had more conservative water use and greater sensitivity of g s to vapor pressure deficit (VPD) compared to IMP pasture, suggesting cane may be more tolerant of increasing atmospheric water demand. In contrast, SN showed lower g s and weaker stomatal sensitivity to VPD compared to cane, resulting in lower T. In cane, stomatal regulation and T varied across growth cycles, with stomata becoming less water conservative as stands matured, highlighting the importance of incorporating stand age-dependent stomatal regulation into hydrological models. Evaporation (E) was higher in cane than pastures (19%–26%), partially offsetting WUE gains. Cane exhibited higher intrinsic WUE (GPP/g s ; Gross Primary Productivity), ecosystem WUE (GPP/ET), and harvest WUE (harvest/ET) than both pasture types. Large-scale pasture-to-cane conversion could produce widely contrasting hydrological outcomes. The net regional impact will depend on the proportion of each pasture type converted and on cane's high g s sensitivity to VPD, which triggers tight stomatal regulation and conservative water use, both of which will become increasingly consequential under intensifying atmospheric water demand.

bioenergy↗
DOE OSTI2023

Pressure effects on reactivity and extinction of n -dodecane diffusion cool flame

Low-temperature combustion at high pressure has gained increasing interest owing to the growing demand for clean and efficient propulsion and energy systems. In particular, cool flame structure, extinction limits, and fuel reactivities at high pressure are critical in affecting the engine performance at near-limit conditions. Here, this paper investigates the effects of pressure on the cool flame extinction limit, structure, radical index, reactivity, and oxygen concentration dependence by experiments, analysis, and modeling. A large n-alkane, namely n-dodecane (n C 12 H 26 ), is selected to study its diffusion cool flame dynamics and reactivity in a high-pressure counterflow burner up to 10 atm. The experimental results show that higher pressure increases the cool flame extinction strain rates, and that the pressure-weighted extinction strain rate (aP) is proportional to the square of pressure, P 2 . A scaling analysis explains the relationship between the dependence of flame structure, heat release rate, and pressure-weighted cool flame strain rate on pressure. Furthermore, radical indexes at different pressures are measured by isolating the thermal and transport effects from the chemical contribution to diffusion cool flame extinction. The radical index clearly shows that the low-temperature reactivity increases with pressure. In addition, due to the critical role of multiple oxygen addition reactions in low-temperature chemistry, the relationship between the cool flame extinction limit and the oxygen concentration is explored. It is found that the cool flame extinction limits are proportional to the n th power of the oxygen concentration, [O 2 ] n , and increasing pressure leads to stronger extinction limit dependence (larger n) on the oxygen concentration. The present experiment and detailed kinetic analysis show clearly that increasing pressure promotes the low-temperature chemistry including the oxygen addition reactions, while the scaling analysis explains well the experimental results.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

The effect of initial conditions on mixing transition of the Richtmyer–Meshkov instability

We investigate the late-time Richtmyer–Meshkov instability (RMI) growth of sinuous perturbations on an air/sulphur hexafluoride interface (Atwood number, A ~ 0.67 ) subjected to a Mach 1.2 planar shock wave at Los Alamos National Laboratory's vertical shock tube facility. Interface perturbations are established using a novel membraneless technique where cross-flowing air and SF 6 separated by an oscillating splitter plate create a perturbed density interface. The interface formed has multi-modal features and residual small perturbations, however, a dominant mode is still noticeable. The late-time perturbation growths scale with $ka_0$ initial conditions (where $k$ is the wavenumber and $a_0$ is the initial amplitude of the dominant mode) as measured at the pre-shock interface. Past nonlinear models based on potential-flow theory, heuristic/interpolation approaches, Padé approximants and numerical simulations are evaluated against present experimental results. Accounting for an explicit $ka_0$ dependence in Sadot et al.'s model, we propose an empirical rational function that captures the asymptotic behaviour of perturbation growth for a broad range of initial conditions ( 0.30 $≤ ka_0 ≤$ 0.86 ). The onset of mixing transition and its initial condition dependence are investigated with respect to the minimum state criterion ($Re =$ 1.6 $\times$ 10 5 ) for unsteady flows by Zhou. Earlier mixing transitions for higher ka0 initial conditions are noted from local and global Reynolds number estimates which are corroborated by the existence of an inertial sub-range and formation of mixing regions indicating the physical significance of the minimum state criterion in RMI flows. The transition is accompanied by the increasing teapot-like appearance of joint probability density functions of $p–q$ (invariants of the reduced velocity gradient tensor), establishing the technique as a useful tool for turbulence detection in two-dimensional diagnostics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Machine Learning for First Principles Calculations of Material Properties for Ferromagnetic Materials

The investigation of finite temperature properties using Monte-Carlo (MC) methods requires a large number of evaluations of the system’s Hamiltonian to sample the phase space needed to obtain physical observables as function of temperature. DFT calculations can provide accurate evaluations of the energies, but they are too computationally expensive for routine simulations. To circumvent this problem, machine-learning (ML) based surrogate models have been developed and implemented on high-performance computing (HPC) architectures. In this paper, we describe two ML methods (linear mixing model and HydraGNN) as surrogates for first principles density functional theory (DFT) calculations with classical MC simulations. These two surrogate models are used to learn the dependence of target physical properties from complex compositions and interactions of their constituents. We present the predictive performance of these two surrogate models with respect to their complexity while avoiding the danger of overfitting the model. An important aspect of our approach is the periodic retraining with newly generated first principles data based on the progressive exploration of the system’s phase space by the MC simulation. The numerical results show that HydraGNN model attains superior predictive performance compared to the linear mixing model for magnetic alloy materials.

Eisenbach, Markus↗
DOE OSTI2022

PyOECP: A flexible open-source software library for estimating and modeling the complex permittivity based on the open-ended coaxial probe (OECP) technique

Here, we present PyOECP, a Python-based flexible open-source software for estimating and modeling the complex permittivity obtained from the open-ended coaxial probe (OECP) technique. The transformation of the measured reflection coefficient to complex permittivity is performed based on three different methods. The software library contains the dielectric spectra of common reference liquids, which can be used to transform the reflection coefficient into the dielectric spectra. Several Python routines that are commonly employed (e.g., SciPy and NumPy) in the field of science and engineering are required only so that the users can alter the software structure depending on their needs. The modeling algorithm exploits the Markov Chain Monte Carlo method for the data regression. The discrete relaxation models can be built by a proper combination of well-known relaxation models. In addition to these models, electrode polarization, a typical measurement artifact for interpreting dielectric spectra, can be incorporated into the modeling algorithm. A continuous relaxation model, which solves the Fredholm integral equation of the first kind (a mathematically ill-posed problem), is also included. This open-source software enables users to freely adjust the physical parameters to obtain physical insight into their materials under test and will be consistently updated for more accurate measurement and interpretation of dielectric spectra in an automated manner. This work describes the theoretical and mathematical background of the software, lays out the workflow, and validates the software functionality based on both synthetic and empirical data included in the software.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Development, characterization, and modeling of a high-performance Ru/B2CA catalyst for ammonia synthesis

This paper documents the development and performance of a nano-phase Ru catalyst on a (BaO) x (CaO) y (Al 2 O 3 ) support. Extensive screening of the support’s ternary composition shows the best stoichiometry is (BaO) 2 (CaO)(Al 2 O 3 ), denoted B2CA. The paper first describes catalyst preparation and characterization. The paper reports a detailed 12-step reaction mechanism that represents ammonia synthesis over wide ranges of temperature, pressure, space velocity, and feed composition. Additionally, the mechanism is developed and validated using results of packed-bed experiments. The elementary reaction pathways consider surface adsorbates, including catalyst-poisoning behaviors. The rate expressions include important coverage-dependent activation barriers. Machine learning models assist interpretation of the catalyst-support interactions. The detailed chemistry is much more predictive than is possible with global representations (N 2 + 3H 2 ⇌ 2NH 3 ). The validated models can be applied to assist optimizing reactor design and operating conditions.

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

Combining synchrotron X-ray diffraction, mechanistic modeling and machine learning for in situ subsurface temperature quantification during laser melting

Laser melting, such as that encountered during additive manufacturing, produces extreme gradients of temperature in both space and time, which in turn influence microstructural development in the material. Qualification and model validation of the process itself and the resulting material necessitate the ability to characterize these temperature fields. However, well established means to directly probe the material temperature below the surface of an alloy while it is being processed are limited. To address this gap in characterization capabilities, a novel means is presented to extract subsurface temperature-distribution metrics, with uncertainty, from in situ synchrotron X-ray diffraction measurements to provide quantitative temperature evolution data during laser melting. Temperature-distribution metrics are determined using Gaussian process regression supervised machine-learning surrogate models trained with a combination of mechanistic modeling (heat transfer and fluid flow) and X-ray diffraction simulation. The trained surrogate model uncertainties are found to range from 5 to 15% depending on the metric and current temperature. The surrogate models are then applied to experimental data to extract temperature metrics from an Inconel 625 nickel superalloy wall specimen during laser melting. The maximum temperatures of the solid phase in the diffraction volume through melting and cooling are found to reach the solidus temperature as expected, with the mean and minimum temperatures found to be several hundred degrees less. The extracted temperature metrics near melting are determined to be more accurate because of the lower relative levels of mechanical elastic strains. However, uncertainties for temperature metrics during cooling are increased due to the effects of thermomechanical stress.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Flyby reaction trajectories: Chemical dynamics under extrinsic force

Dynamic effects are an important determinant of chemical reactivity and selectivity, but the deliberate manipulation of atomic motions during a chemical transformation is not straightforward. In this work, we demonstrate that extrinsic force exerted upon cyclobutanes by stretching pendant polymer chains influences product selectivity through force-imparted nonstatistical dynamic effects on the stepwise ring-opening reaction. The high product stereoselectivity is quantified by carbon-13 labeling and shown to depend on external force, reactant stereochemistry, and intermediate stability. Computational modeling and simulations show that, besides altering energy barriers, the mechanical force activates reactive intramolecular motions nonstatistically, setting up “flyby trajectories” that advance directly to product without isomerization excursions. A mechanistic model incorporating nonstatistical dynamic effects accounts for isomer-dependent mechanochemical stereoselectivity.

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

Gravitational waves from global cosmic strings and cosmic archaeology

Global cosmic strings are predicted in many motivated extensions to the Standard Model of particle physics, with close connections to axion dark matter physics. Recent studies suggest that, although subdominant relative to Goldstone emission, gravitational wave (GW) signals from global strings can be detectable with current and planned GW detectors such as LIGO, LISA, DECIGO/BBO, ET/CE and AEDGE/AION, as well as pulsar timing arrays such as PPTA, NANOGrav and SKA. This work is an extensive, updated study on GWs from a global cosmic string network, taking into account of the most recent developments related to the subject. The main analysis is based on the analytical Velocity-dependent One-Scale (VOS) model calibrated with recent simulation results, which provides a generic protocol for such calculations with details given. We also demonstrate how the GW signal can be influenced with variations to the baseline model: this includes considering the uncertainties of model parameters and the potential deviation from the conventional VOS model prediction (i.e. the scaling behavior) as suggested by some of the recent simulation results. Furthermore, we investigated in detail the effect of a non-standard cosmology (e.g. early matter domination or kination) or new particle species on the GW signals from global strings. We demonstrate that the frequency spectrum of GW background from global cosmic strings can be used to probe the cosmic history prior to the Big Bang nucleosynthesis (BBN) (i.e. the primordial dark age) up to a temperature of T ~ 10 8 GeV.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗