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

Data-driven methods for diffusivity prediction in nuclear fuels

The growth rate of structural defects in nuclear fuels under irradiation is intrinsically related to the diffusion rates of the defects in the fuel lattice. The generation and growth of atomistic structural defects can significantly alter the performance characteristics of the fuel. This alteration of functionality must be accurately captured to qualify a nuclear fuel for use in reactors. Predicting the diffusion coefficients of defects and how they impact macroscale properties such as swelling, gas release, and creep is therefore of significant importance in both the design of new nuclear fuels and the assessment of current fuel types. In this article, we apply data-driven methods focusing on machine learning (ML) to determine various diffusion properties of two nuclear fuels—uranium oxide and uranium nitride. We show that using ML can increase, often significantly, the accuracy of predicting diffusivity in nuclear fuels in comparison to current analytical models. We also illustrate how ML can be used to quickly develop fuel models with parameter dependencies that are more complex and robust than what is currently available in the literature. In conclusion, these results suggest there is potential for ML to accelerate the design, qualification, and implementation of nuclear fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Bio-inspired alula-based winglet design for enhanced heat transfer in high temperature fin-and-tube heat exchangers

Fin-and-tube heat exchangers (FTHEs) are widely used for high-temperature flue-gas heat recovery, but their performance is often limited by wake regions and non-uniform fin-surface temperatures. This study proposes and numerically evaluates four bio-inspired longitudinal vortex generator (VG) configurations in a high-temperature FTHE with flue-gas inlet temperature ∼1230 K: double-delta, curved double-delta, alula, and a new curved-alula geometry. The reference fin is not hydraulically plain; it already incorporates leading-edge separation columns and convex protrusions, so the alula-type winglets are assessed as downstream add-ons acting on a strongly disturbed flow. In a second step, perforations (one, two and three circular holes) are introduced into the curved-alula VGs to further tailor the flow field. Three-dimensional simulations with the Shear Stress Transpor (SST) $k - ω$ model, temperature-dependent flue-gas properties and conjugate conduction are carried out for gas-side Reynolds numbers $Re_g ≈ 8.0$ x $10^2 - 3.6$ x $10^3$ (mass flow rates 0.5 – 2.5 g/s), and the designs are compared in terms of surface heat flux, Nusselt number, friction factor and hydrothermal performance factor (HTPF). For this already-promoted fin, the additional downstream winglets provide moderate, incremental hydrothermal gains. At the highest Reynolds number, the best non-perforated design (curved-alula) increases surface heat flux from 1630.9 to 1794.7 kW/m² (∼ 10 % gain) and the Nusselt number from 227.6 to 242.6 (∼ 7 % gain), while the friction factor rises from 0.26 to about 0.30, yielding HTPF values close to unity (∼ 0.9 – 1.0). Introducing circular perforations into the curved-alula winglets acts mainly as a wake-bleeding refinement: the three-hole configuration provides a heat flux of 1824.7 kW/m² and a pressure drop of 127.9 Pa, with HTPF in the range ∼ 1.03 – 1.14 and a small (∼ 1 – 3 %) improvement over the solid curved-alula design. Flow-field analysis shows that the perforated curved-alula VGs shrink tube-wake regions, thin the thermal boundary layer and homogenize the fin-surface temperature (outlet-gas temperature ∼ 510 – 520 K and fin-surface temperature ∼ 420 – 421 K for the three-hole case). An optimal flue-gas mass flow rate of ∼ 1 g/s ($Re_g ≈ 1.5$ x $10^3$) is identified, beyond which additional heat-transfer gains are offset by rapidly increasing pressure losses. Overall, the results highlight that initial fin geometry and VG placement are as important as VG shape: alula-based winglets are expected to yield larger relative gains on simpler flat-fin layouts or when positioned closer to the fin leading edge and tube

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

An experimental system and procedure of unsteady-state relative permeability test for gas hydrate-bearing sediments

Reliable estimations of the relative permeability of gas and water in hydrate-bearing sediments (HBS) and the dependency of the relative permeability on hydrate saturation are critical to predict the productivity of a hydrate reservoir. Yet, this remains poorly estimated owing to lack of experimental data associated with difficulties in conducting multiphase flow experiments in HBS. Recognizing the experimental challenges, this study intends to develop and validate a new experimental system and procedure of unsteady-state relative permeability test that can generate reliable and reproducible flow measurements in HBS. Furthermore, gas hydrate is considered as a part of solid matrix in the sediment, so one of the challenges is to maintain a constant hydrate saturation, which is achieved in this experimental study using tight pressure-temperature (P-T) control near the hydrate stability boundary. The measured differential pressure across the specimen, methane injection flow rate, and volume of displaced brine are used to calculate the relative permeability by adopting a conventional Buckley-Leverett theory-based interpretation method. Residual brine saturation calculated for the hydrate-bearing specimen is higher than that of hydrate-free specimen, presumably due to decrease in pore size, increase in heterogeneity of solid matrix, and increase in size distribution of solid matrix and pore in the presence of hydrates. Further studies are necessary to represent the results of the unsteady-state flow experiment in HBS with a gas hydrate-dependent relative permeability model.

03 NATURAL GAS↗
DOE OSTI2025

Combining geometric-optical and spectral invariants theories for modeling canopy fluorescence anisotropy

The spectral invariants theory ( p -theory) has received much attention in the field of quantitative remote sensing over the past few decades and has been adopted for modeling of canopy solar-induced chlorophyll fluorescence (SIF). However, the spectral invariant properties (SIP) in simple analytical formulae have not been applied for modeling canopy fluorescence anisotropy primarily because they are parameterized in terms of leaf total scattering, which precludes the differentiation between forward and backward leaf SIF emissions. In this study, we have developed the canopy-SIP SIF model by combining geometric-optical (GO) theory to account for asymmetric leaf SIF forward and backward emissions at the first-order scattering and by modeling multiple scattering based on the p-theory, thus avoiding the dependence on radiative transfer models. The applicability of the model simulations especially over 3D heterogeneous canopies was improved by incorporating canopy structure through multi-angular clumping index, and by modeling single scattering from the four components of the scene in view according to the GO approach. The results show good consistency with both the state-of-the-art SIF models and multi-angular field SIF observations over grass and chickpea canopies. Further, the coefficient of determination (R²) between the simulated SIF and field measurements was 0.75 (red) and 0.74 (far-red) for chickpea, and 0.65 (both red and far-red) for grass. The average relative error was approximately 3% for 1D homogeneous scenes when comparing the canopy-SIP SIF model simulations to the SCOPE model simulations, and around 4% for the 3D heterogeneous scene when comparing to the LESS model simulations. The results indicate that the proposed approach for separating asymmetric leaf SIF emissions is a robust way to keep a balance between satisfactory simulation accuracy and efficiency. Model simulations suggest that neglecting the leaf SIF asymmetry can lead to an underestimation of canopy red SIF by 6.3% to 42.6% for various leaf biochemical and canopy structural parameters. This study presents a simple but efficient analytical approach for canopy fluorescence modeling, with potential for large-scale canopy fluorescence simulations.

3D heterogeneous structure↗
DOE OSTI2023

Relationships between maternally-transferred mercury and hatchling development, behavior, and survival in the American alligator ( Alligator mississippiensis )

Mercury is a toxic and pervasive environmental contaminant that can be transferred from mother to offspring during development. Consequences of maternally-transferred mercury have been observed in vertebrate taxa, including reduced clutch viability, reduced offspring size, and behavioral alterations. These sublethal effects have been assumed to decrease survivorship, though this is seldom assessed. Here, we examined how maternally-transferred mercury interacts with incubation temperature to influence reproductive success, offspring behavior, and subsequent survival in the American alligator (Alligator mississippiensis). We collected nine clutches of eggs from a mercury contaminated reservoir on the Savannah River Site, South Carolina, and incubated eggs at either female- or male-promoting temperatures. Clutch-averaged mercury in egg yolk was high relative to other studies in crocodilians and ranged from 0.248 to 0.554 ppm compared to 0.018–0.052 ppm at a site with low levels of mercury contamination; mercury levels in hatchling blood ranged from 0.090 to 0.490 ppm ($\overline{x}$= 0.240 ppm, n = 158). We found few, mostly negligible correlations between life history traits and mercury but noted a positive relationship with egg mass, possibly mediated by correlated maternal effects such as resource provisioning. Incubation temperature exerted strong effects on hatchling phenotypes, with warmer, male-promoting temperatures producing larger and bolder hatchlings. Presumptive females, produced from cooler incubation temperatures, spent more time in warm areas during behavior trials. Hatchlings were released 10–15 days post-hatch and surveyed over eight months to assess survival. Survivorship was positively correlated with hatchling size and negatively correlated with proportional time spent in warm areas. Presumptive females had much lower survival, and overall survivorship for the eight-month period was 0.185–0.208, depending on the modelling approach. Finally, our study suggests that, within the range of concentrations we observed, incubation temperature has a stronger effect on offspring behavior and survival than maternally-transferred mercury pollution in American alligators.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

MCTrans++: a 0-D model for centrifugal mirrors

The centrifugal mirror confinement scheme incorporates supersonic rotation of a plasma into a magnetic mirror device. This concept has been shown experimentally to drastically decrease parallel losses and increase plasma stability as compared with prior axisymmetric mirrors. MCTrans++ is a dimensionless (0-D) scoping tool which rapidly models experimental operating points in the Centrifugal Mirror Fusion Experiment (CMFX) at the University of Maryland. In the low-collisionality regime, parallel losses can be modelled analytically. A confining potential is set up that is partially ambipolar and partially centrifugal. Due to the stabilizing effects of flow shear, the perpendicular losses can be modelled as classical. Radiation losses such as bremsstrahlung and cyclotron emission are taken into account. A neutrals model is included, and, in some circumstances, charge-exchange losses are found to exceed all other loss mechanisms. We use the SUNDIALS ARKODE library to solve the underlying equations of this model; the resulting software is suitable for scanning large parameter spaces, and can also be used to model time-dependent phenomena such as a capacitive discharge. MCTrans++ has been used to verify results from prior centrifugal mirrors, create an experimental plan for CMFX and find configurations for future reactor-scale fusion devices.

Schwartz, Nick R. (ORCID:0000000296846435)↗
DOE OSTI2024

Trust Not Verify? The Critical Need for Data Curation Standards in Materials Informatics

The importance of data curation has been recognized in multiple areas of research; however, the discussion of this important issue is only beginning to emerge in materials science. In this Perspective, we highlight the benefits of using the standardized data curation protocols in materials science and discuss current gaps in accurate and reproducible data reporting using case studies drawn from high-impact materials science papers and well-known databases such as the Crystallography Open Database (COD) and the Cambridge Structural Database (CSD). We argue that both experimental and computational materials scientists need to embrace a culture of rigorous data curation as part of modern research data management. We propose a sample data curation pipeline for materials chemistry and illustrate its use by creating two new materials chemistry databases. Here, we hope that this perspective will serve to catalyze further discussion and promote the continuous development of rigorous data curation practices within the materials science research community. We posit that adherence to best practices of data curation will promote and enhance the reliability, reproducibility, and integrity of materials research and enable the development of reliable AI and machine learning models that critically depend on the use of quality data.

Chemical structure↗
DOE OSTI2026

Deciphering the Solvation Structure of Aqueous ZnCl 2 Solutions from X-ray Absorption Spectra Using the Interpretable Graph Neural Network

Machine learning (ML) provides powerful pathways for predicting spectroscopic observables from atomic structures, but its broader impact depends on making model predictions interpretable in terms of physical and chemical principles. Here, we introduce a physics-guided graph neural network (GNN) model that predicts Zn K-edge X-ray spectroscopy (XAS) spectra of aqueous ZnCl 2 solutions. Training data are generated from ab initio XAS calculations on molecular dynamics snapshots obtained using a machine learning interatomic potential. The GNN reproduces experimental spectra across concentrations from dilute (<0.1 m) to highly concentrated (30 m, “water-in-salt”) regimes and scales efficiently to large, disordered liquid systems beyond the reach of conventional ab initio approaches. Gradient-based attribution analysis reveals that the model learns physically meaningful structure-spectrum relationships. Ligand-specific attributions reflect orbital hybridization patterns and the origin of the excitations derived from the density functional theory. Bond-length attributions recover spectral shifts consistent with multiple-scattering theory. Finally, this work bridges data-driven prediction with electronic-structure theory, establishing a general paradigm for interpretable ML that links atomic structure, electronic structure, and spectroscopic observables.

25 ENERGY STORAGE↗
DOE OSTI2020

One-Dimensional Glassy Behavior of Ultraconfined Water Strings

Water is renowned for its anomalous behaviors, which can be linked to a distributed H-bond network in bulk water. Ultraconfinement of the water molecule can remove H-bonding, leaving only molecular water. In natural cordierite crystals, water is trapped in an orthorhombic channel with an average diameter of 5.7 Å, running through the center of the unit cell parallel to the c-axis. Calorimetric measurements reveal the existence of a one-dimensional (1D) glass linked to this water. In these channels, water molecules in truncated, sparse 1D strings interact only via dipole–dipole correlations. Here, a physical 1D glass is formed from these strings. This unusual state can be explained by a modified Ising model. This model predicts a dependence of the glass transition temperature, Tg, on the size of these domains. This is confirmed experimentally.

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

Hydrate is a Nonwetting Phase in Porous Media

In porous media containing gas hydrate, the hydrate morphology impacts the flow behavior of the fluid phases. We hypothesize that hydrate emplaces itself as a nonwetting phase and use this idea to describe relative permeability of water in a hydrate/water system. We perform steady-state relative permeability measurements in hydrate-bearing samples with a range of hydrate saturation. We measure and compare water relative permeability in the presence of gas and in the presence of hydrate and find that the water relative permeability is the same in both cases. This suggests that (1) relative permeability for hydrate systems can be obtained without performing difficult experiments on hydrate bearing sediments, (2) flow properties are porous-media dependent, and (3) models that assume a fixed pore or tube geometry are inadequate.

58 GEOSCIENCES↗
DOE OSTI2024

Computed entropy spectra for grain boundary segregation in polycrystals

Abstract Grain boundary solute segregation influences most bulk material properties, and understanding solute thermodynamics at grain boundaries is critical for engineering them. However, the vast grain boundary space in polycrystals is challenging to evaluate due to its size, especially for the intrinsically hard-to-compute segregation excess entropy. Here data science methods are used to generate a database of site-wise grain boundary segregation entropy spectra for 155 dilute binary alloys within the harmonic approximation. The spectral framework allows scale bridging between the calculated atomistic site-wise energy-entropy spectra and macroscopic segregation entropy estimates. The results affirm that macroscopic averaging is not sufficient: a spectral treatment of grain boundary segregation is needed to accurately model bulk temperature dependence of grain boundary solute segregation. The calculated spectral entropy database and thermodynamic framework can be applied for both understanding segregation experiments and alloy design exercises, paving the way to a finite-temperature grain boundary genome.

Chemistry↗
DOE OSTI2021

Gaussian processes for autonomous data acquisition at large-scale synchrotron and neutron facilities

The execution and analysis of complex experiments are challenged by the vast dimensionality of the underlying parameter spaces. Although an increase in data-acquisition rates should allow broader querying of the parameter space, the complexity of experiments and the subtle dependence of the model function on input parameters remains daunting owing to the sheer number of variables. New strategies for autonomous data acquisition are being developed, with one promising direction being the use of Gaussian process regression (GPR). GPR is a quick, non-parametric and robust approximation and uncertainty quantification method that can be applied directly to autonomous data acquisition. We review GPR-driven autonomous experimentation and illustrate its functionality using real-world examples from large experimental facilities in the USA and France. We introduce the basics of a GPR-driven autonomous loop with a focus on Gaussian processes, and then shift the focus to the infrastructure that needs to be built around GPR to create a closed loop. Finally, the case studies we discuss show that Gaussian-process-based autonomous data acquisition is a widely applicable method that can facilitate the optimal use of instruments and facilities by enabling the efficient acquisition of high-value datasets.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Raman spectrum and polarizability of liquid water from deep neural networks

We introduce a scheme based on machine learning and deep neural networks to model the environmental dependence of the electronic polarizability in insulating materials. Application to liquid water shows that training the network with a relatively small number of molecular configurations is sufficient to predict the polarizability of arbitrary liquid configurations in close agreement with ab initio density functional theory calculations. In combination with a neural network representation of the interatomic potential energy surface, the scheme allows us to calculate the Raman spectra along 2-nanosecond classical trajectories at different temperatures for H 2 O and D 2 O. The vast gains in efficiency provided by the machine learning approach enable longer trajectories and larger system sizes relative to ab initio methods, reducing the statistical error and improving the resolution of the low-frequency Raman spectra. Decomposing the spectra into intramolecular and intermolecular contributions elucidates the mechanisms behind the temperature dependence of the low-frequency and stretch modes.

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

nn-PINNs: Non-Newtonian physics-informed neural networks for complex fluid modeling

Time- and rate-dependent material functions in non-Newtonian fluids in response to different deformation fields pose a challenge in integrating different constitutive models into conventional computational fluid dynamic platforms. Considering their relevance in many industrial and natural settings alike, robust data-driven frameworks that enable accurate modeling of these complex fluids are of great interest. The main goal is to solve the coupled Partial Differential Equations (PDEs) consisting of the constitutive equations that relate the shear stress to the deformation and fully capture the behavior of the fluid under various flow protocols with different boundary conditions. Here, in this work, we present non-Newtonian physics-informed neural networks (nn-PINNs) for solving systems of coupled PDEs adopted for complex fluid flow modeling. The proposed nn-PINN method is employed to solve the constitutive models in conjunction with conservation of mass and momentum by benefiting from Automatic Differentiation (AD) in neural networks, hence avoiding the mesh generation step. nn-PINNs are tested for a number of different complex fluids with different constitutive models and for several flow protocols. These include a range of Generalized Newtonian Fluid (GNF) empirical constitutive models, as well as some phenomenological models with memory effects and thixotropic timescales. nn-PINNs are found to obtain the correct solution of complex fluids in spatiotemporal domains with good accuracy compared to the ground truth solution. We also present applications of nn-PINNs for complex fluid modeling problems with unknown boundary conditions on the surface, and show that our approach can successfully recover the velocity and stress fields across the domain, including the boundaries, given some sparse velocity measurements.

42 ENGINEERING↗
DOE OSTI2026

Assessing plasma face component thermal response to rotating 3D magnetic fields for SPARC tokamak

Thermal response simulations of plasma-facing components (PFCs) in the SPARC tokamak, performed with the HEAT code, show that three-dimensional (3D) heat loads resulting from stationary n=1 perturbations require highly radiative scenarios, with up to 95% of the power crossing the separatrix (P SOL ) being radiated, to maintain PFC temperatures within acceptable operational limits, whereas the application of slowly rotating 3D fields substantially reduces the thermal loads. The HEAT module, developed to predict heat loads from non-axisymmetric plasmas, is extended to model time-dependent heat flux patterns generated by rotating 3D fields, and a comprehensive thermal analysis is performed on PFCs subjected to both the maximum and minimum power loads, as well as to rotating heat flux distributions, to evaluate the temperature evolution for varying perturbation amplitudes and rotation frequencies. The extension of this analysis to 3D fields with toroidal mode number n=2 shows that this configuration leads to weaker localized heat flux peaks relative to the n=1 case, enabling safe operation with less than 80% of the power radiated when static 3D fields of low amplitude are applied, while using slowly rotating fields at higher amplitudes. These results indicate that n=2 perturbations are generally less detrimental to divertor power exhaust, emphasizing the strong dependence of divertor power exhaust on the characteristics of the applied 3D fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Secondary electron emission for reticulated carbon foam surfaces using direct measurements and spectroscopic analysis

This study investigates secondary electron emission (SEE) characteristics of reticulated foams using direct measurements and analytical modeling. Total SEE was quantified, revealing suppression of up to 44% in carbon foam structures compared to planar graphite surfaces. An optimal geometric configuration was identified and supported by analytical models. SEE angular dependence experiments showed diverse behaviors: fiber-like behavior and directional dependence for pore and ligaments on the mm scale, with fuzz-like characteristics when the foam features are between 10–100 µm. Electron energy analyzer measurements showed that carbon foams preferentially suppress inelastic backscattered electrons (BSEs) more so than true secondary electrons (SEs). The analysis indicated a larger fraction of low-energy SE generation in foams compared to flat surfaces due to increased emission from curved fiber ligaments and tertiary SEs from high-energy BSEs. These findings have implications for design and optimization of materials with tailored electron emission properties for applications like plasma-facing components, spacecraft materials, and accelerator surfaces.

Auger↗
DOE OSTI2020

An accurate perturbative approach to redshift space clustering of biased tracers in modified gravity

We extend the scale-dependent Gaussian Streaming Model (GSM) to produce analytical predictions for the anisotropic redshift-space correlation function for biased tracers in modified gravity models. Employing the Convolution Lagrangian Perturbation Theory (CLPT) re-summation scheme, with a local Lagrangian bias schema provided by the peak-background split formalism, we predict the necessary ingredients that enter the GSM, the real-space halo pairwise velocity and the pairwise velocity dispersion. We further consider effective field theory contributions to the pairwise velocity dispersion in order to model correctly its large scale behavior. We apply our method on two widely-considered modified gravity models, the chameleon-screened ƒ(R) Hu-Sawicki model and the nDGP Vainshtein model and compare our predictions against state-of-the-art N-body simulations for these models. We demonstrate that the GSM approach to predict the monopole and the quadrupole of the redshift-space correlation function for halos, gives very good agreement with the simulation data, for a wide range of screening mechanisms, levels of screening and halo masses at z = 0.5 and z = 1. Our work shows the applicability of the GSM, for cosmologies beyond GR, demonstrating that it can serve as a powerful predictive tool for the next stage of cosmological surveys like DESI, Euclid, LSST and WFIRST.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

3D radiated power analysis of JET SPI discharges using the Emis3D forward modeling tool

Abstract Precise values for radiated energy in tokamak disruption experiments are needed to validate disruption mitigation techniques for burning plasma tokamaks like ITER and SPARC. Control room analysis of radiated power ( P rad ) on JET assumes axisymmetry, since fitting 3D radiation structures with limited bolometry coverage is an under-determined problem. In mitigated disruptions, radiation is toroidally asymmetric and 3D, due to fast-growing 3D MHD modes and localized impurity sources. To address this problem, Emis3D adopts a physics motivated forward modeling (‘guess and check’) approach, comparing experimental bolometry data to synthetic data from user-defined radiation structures. Synthetic structures are observed with the Cherab modeling framework and a best fit chosen using a reduced χ 2 statistic. 2D tomographic inversion models are tested, as well as helical flux tubes and 3D MHD simulated structures from JOREK. Two nominally identical pure neon shattered pellet injection (SPI) mitigated discharges in JET are analyzed. 2D tomographic inversions with added toroidal freedom are the best fits in the thermal quench (TQ) and current quench (CQ). In the pre-TQ, 2D reconstructions are statistically the best fits, but are likely over-optimized and do not capture the 3D radiation structure seen in fast camera images. The next-best pre-TQ fits are helical structures that extend towards the high-field side, consistent with an impurity flow under the magnetic nozzle effect also observed in JOREK simulations. Whole-disruption radiated fractions of 0.98 + 0.03 / − 0.29 and 1.01 + 0.02 / − 0.17 are found, suggesting that the stored energy may have been fully mitigated by each SPI, although mitigation efficiencies well below ITER and SPARC requirements for high energy pulses are still within the large uncertainties. Emis3D is also used to validate JOREK SPI simulations, and confirms improvements in matching experiment from changes to impurity modeling. Time-dependent toroidal peaking factors are calculated and discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗