Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “extrapolation”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Chromium-coated cladding analysis under simulated LOCA burst conditions

Of the near-term accident-tolerant fuel concepts identified, chromium coatings have been shown to slow cladding oxidation without altering fuel system geometry or greatly affecting neutronic performance. To identify how coated cladding tubes perform under high-temperature accident conditions, pressurized-tube burst tests were conducted in the Severe Accident Test Station at Oak Ridge National Laboratory. To analyze how the coating affects cladding behavior, coated and uncoated cladding burst tests were simulated by using the BISON fuel performance code. Thermocouple data from these tests were fit into axial and azimuthal profiles and combined to generate 3D cladding surface temperatures, and pressure transducer data were compared until cladding failure. Additionally, the cladding temperatures at failure and the internal gas pressure evolution show relatively good agreement between the simulation and experiment results. Simulations were then extrapolated to demonstrate the coatings effectiveness to increasing the cladding burst temperature by using a parametric evaluation of the initial tube gas pressure. This work demonstrates the possibility of an increased cladding failure margin under transient conditions due to the addition of chromium coating, and more pertinently, an increased cladding failure margin from more accurate experimental characterization.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Effect of neutron beamline irradiation on hydrogen equilibrium vapor pressure over hydrided Zircaloy-4

Here, the effect of neutron irradiation on the hydrogen equilibrium vapor pressure over hydrided samples of nickel-plated Zircaloy-4 was evaluated at the Washington State University Research Reactor’s epithermal beamline facility. To enable this study, a measurement capability was designed and assembled to obtain low pressure (<0.1 Pa) measurements over the temperature range of interest (i.e., ~300° - 400°C). Measurement results were in good agreement with tends extrapolated from previous data obtained at higher temperatures as well as recent thermodynamic data. No irradiation enhancement of vapor pressure was observed under the incident flux of epithermal neutrons over this temperature range.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Modeling occupancy-driven building loads for large and diversified building stocks through the use of parametric schedules

Building energy modeling provides a fundamental tool to assess the potential for energy efficiency to contribute to reducing world energy consumption and global emissions. Occupancy-related operations are a key source of uncertainty for building energy analysis, particularly for aggregated building stocks. At a district or city level, it is critical to estimate aggregated power load profiles for sizing power grid infrastructure, power plant capacity allocation, and energy efficiency measures. The stochastic nature of behavior-related operations complicates the creation of models that accurately capture building load profiles for entire building stocks. This research introduces a new methodology called parametric schedules to model occupancy-driven schedules for large and diverse building stocks. In contrast to computationally expensive methodologies proposed in the literature, our work does not use a recursive time-consuming step. Occupancy is estimated by the extrapolation of operation times directly from metered electric consumption data; occupancy-related schedules are stochastically assigned to each building model, guaranteeing diversity of operation times in the stock. Our procedure has been tested on a large, diversified data-set of 25,000 commercial buildings in Los Angeles, California. It proved to be able to adequately represent the stochastic schedules diversity of the stock and to refine the stock calibration process by 1%. This innovative approach represents a useful asset for utility companies, grid operators, urban planners, and balancing authorities, seeking to improve building stock modeling and better estimate the impact of energy conservation measures. – This work is part of a larger stock modeling tool called ComStock, which is under development by NREL.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Effects of land-based wind turbine upsizing on community sound levels and power and energy density

Multiple technological, social, and market factors of wind power are evolving rapidly. Most notably, significant wind turbine scaling is occurring and is forecasted to continue. While the larger turbines expected to be deployed in the future are more powerful and efficient, they are also expected to operate at higher sound levels and require larger setbacks than those installed in the last decade. These sometimes-competing deployment trends and impacts cannot be understood via simple extrapolations of past trends. This study analyzes the effect of these future larger turbines on wind turbine micro-siting, project-level power and energy density, and community noise impacts. Due to their taller heights, larger rotors, and higher sound power levels, future wind turbines will require larger setbacks from homes and greater inter-turbine spacing, resulting in fewer turbines deployed for a given land area. This research finds these changes more than offset the effect of the higher turbine sound emissions, significantly decreasing the average sound levels that wind plant hosting communities experience. Yet, simultaneously, plant layouts using future turbine designs also result in projects with higher installed capacities and annual energy output for a given land area. These increases will likely lead to increased tax benefits and local income in the community. The deployment of fewer turbines on a smaller number of parcels could have implications on siting flexibility and landowner payments.

17 WIND ENERGY↗
DOE OSTI2020

On the prediction of critical heat flux using a physics-informed machine learning-aided framework

The critical heat flux (CHF) corresponding to the departure from nucleate boiling (DNB) crisis is essential to the design and safety of a two-phase flow boiling system. Despite the abundance of predictive tools available to the thermal engineering community, the path for an accurate, robust CHF model remains elusive due to lack of consensus on the DNB triggering mechanism. This work aims to apply a physics-informed machine learning (ML)-aided hybrid framework to achieve superior predictive capabilities. Such a hybrid approach takes advantage of existing understanding in the field of interest (i.e., domain knowledge) and uses ML to capture undiscovered information from the mismatch between the actual and domain knowledge-predicted target. A detailed case study is carried out with an extensive DNB-specific CHF database to demonstrate (1) the improved performance of the hybrid approach as compared to traditional domain knowledge-based models, and (2) the hybrid model’s superior generalization capabilities over standalone ML methods across a wide range of flow conditions. The hybrid framework could also readily extend its applicability domain and complexity on the fly, showing an elevated level of flexibility and robustness. Based on the case study conclusions, the window-type extrapolation mapping methodology is further proposed to better inform high-cost experimental work.

42 ENGINEERING↗
DOE OSTI2020

Improving air leakage prediction of buildings using the fan pressurization method with the Weighted Line of Organic Correlation

In many countries, the fan pressurization method is the most frequently chosen approach for measuring the air leakage of houses. The measurements are usually performed at pressures that far exceed pressures to which buildings are exposed to under normal conditions. A fit of these tests to the power-law formulation allows an extrapolation to data points outside the measured pressure range. With the Ordinary Least Square (OLS) fitting method, the pressure exponent and flow coefficient can be determined. However, the measurement results are highly sensitive to uncertainties induced by external factors like changing wind conditions during the tests, which is neglected by OLS. This may lead to errors in the prediction of flows at lower pressures. The Weighted Line of Organic Correlation (WLOC) is an alternative approach and takes measurement uncertainty into account. In this paper, a statistical analysis of an extensive data set of pressurization measurements has been performed. Both regression techniques have been compared for almost 7500 fan pressurization measurements of six houses in 109 different house leak configurations. The variability in predicting pressure exponent and flow coefficient for both WLOC and OLS regression was analyzed using probability density functions. It was found that the Weighted Line of Organic Correlation significantly decreases the uncertainty in predicting pressure exponent, flow coefficient, and other low-pressure air leakage metrics compared to the Ordinary Least Square fitting. Here, the authors highly recommend the implementation of WLOC in current measurement standards and test equipment.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Thermodynamic modeling of KCl-PrCl 3 and KCl-LiCl-PrCl 3 systems

Molten salt electrolysis can recover the actinides from spent nuclear fuels, and it involves a eutectic LiCl-KCl in molten form as an electrolyte. During reprocessing, the concentration of fission products such as La, Nd, Pr, etc., increases and affects the recovery efficiency of the electrolyte. In this work, thermodynamic modeling of KCl-PrCl 3 and KCl-LiCl-PrCl 3 systems was carried out using the CALPHAD (Calculation of Phase Diagrams) approach for the first time. The thermodynamic functions for the pure salts were taken from the SGTE (Scientific Group Thermodata Europe) Substances (SSUB) database. The experimental thermochemical and phase equilibria data available in the literature were used as input for the assessment of KCl-PrCl 3 and KCl-LiCl-PrCl 3 systems. The model parameters for the KCl-LiCl system were adjusted to include the new Gibbs energy descriptions for the pure salts. In addition, the sublattice model for the liquid phase in the LiCl-PrCl 3 system was modified and reassessed to ensure the model compatibility for higher-order extrapolation. There is a good agreement between the experimental and calculated thermochemical and phase diagram data for all the optimized constituent binaries and the ternary system. Furthermore, this work will be beneficial for determining the solubility limit of PrCl 3 in molten LiCl-KCl electrolytes and their thermodynamic properties for improving the efficiency of the pyrochemical process.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Revisiting the SGTE lattice stability of bcc aluminum

The CALPHAD framework rests on the ability to assign well-defined free energies to all phases, including mechanically unstable ones. Given the experimental difficulties associated with the latter, researchers have turned to computational methods to obtain a theoretically justified definition of the free energy of mechanically unstable phases. Recently, a method called “inflection detection” has been shown to yield free energies that are broadly consistent with the current SGTE (Scientific Group Thermodata Europe) standard for pure elements, which were obtained through an extrapolation of stable solid solutions free energies towards unstable end members. Despite this general agreement, apparent failures for a few elements stand out, notably for the bcc–fcc enthalpy difference for the common element Al. Furthermore we revisit the definition of the bcc Al enthalpy and show that a careful CALPHAD-based analysis of enthalpy trends in five different mechanically stable binary solid solutions in fact does agree very well with the inflection-detection approach.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Evolution of PTCDA-derived seeds prior to graphene nanoribbon growth on Ge(001)

The synthesis of graphene nanoribbons (GNRs) can be realized via CH4 chemical vapor deposition (CVD) on substrates such as Ge(001) that promote highly anisotropic growth kinetics. Small polycyclic aromatic hydrocarbons (PAHs) are first sublimed onto Ge near room temperature to form graphene-like seeds that subsequently initiate GNR growth with CH4 exposure at 1173 K. The behaviors of PAHs in the time in between their sublimation onto Ge and GNR growth are unclear. Here, we study an archetypical PAH, perylene-3,4,9,10-tetracarboxyl acid dianhydride (PTCDA), on Ge(001) using both scanning tunneling microscopy (STM) and density functional theory (DFT) to characterize PAH conformation, surface diffusivity, and clustering, versus temperature. PTCDA becomes mobile above 673 K (consistent with a DFT calculated diffusion barrier of 2.01 eV). The mobile PTCDA molecules meet, cluster, and fuse at Ge step edges. These clusters have an apparent height of 0.4 nm, similar to small graphene islands, and grow laterally with increasing temperature – reaching 1.2 – 2.1 nm in width at 1173 K, consistent with the seed size extrapolated from CVD experiments. These results provide a plausible picture for how PTCDA forms seeds for anisotropic GNR CVD and show that PAHs with reduced surface diffusivity and inter-molecular reactivity are needed to enable more monodisperse PAH-seeded GNR synthesis.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Accurate hydrodynamic force and torque correlations for prolate spheroids from Stokes regime to high Reynolds numbers

Detailed simulations of flow around various prolate spheroids are performed using multi-relaxation time lattice Boltzmann method (LBM). The simulations are performed in the Reynolds numbers range 0.1 ≤ Re ≤ 2000 at different incident angles 0°≤ φ ≤ 90° for various prolate spheroids of aspect ratios 1 ≤ λ ≤ 8. The Re is based on the volume equivalent sphere diameter. From the simulations, accurate correlations for average drag, lift and torque coefficients (C D , C L , and C T respectively) are proposed. The mean deviations between the correlations and the simulation data for C D , C L , and C T are 2%, 6.5%, and 4.7% respectively. Furthermore, the correlations are fitted with the physics in mind such that they can be extrapolated beyond the regimes simulated. The fits are designed such that correlations reduce to analytical (C D and C L ) solutions available for prolate spheroids in the limit Re ≈ 0. At high Re, the correlations mimic the asymptotic flattening of C D as for spheres and therefore, we expect the correlations to be valid until the critical Re (≈ 10 5 ). Furthermore, extensive comparison of our correlations with other high-fidelity simulations from literature demonstrate the accuracy of our work compared to other correlations from literature. This work enables researchers to perform accurate unresolved Euler-Lagrangian simulations for a wide range of elongated particles.

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

Theoretical and kinetic modeling study of hydrazine oxidation

The present work constitutes the first theoretical and kinetic modeling study of hydrazine oxidation, which may be important for burnout in ammonia-fueled combustion. The kinetics of the oxidation of N 2 H 4 , N 2 H 3 and tHNNH by molecular oxygen were investigated via a quantum chemistry/canonical transition state theory approach. Geometries and anharmonic frequencies were obtained with density functional theory, and energies from coupled cluster calculations (CCSD(T)) extrapolated to the infinite basis set limit, with corrections for core-valence electron correlation, scalar relativistic effects, and higher level correlation up to lambda coupled cluster, CCSDT(Q) Λ . The key reactions occurred on the N 2 H 4 O 2 potential energy surface, where the results indicated a fast reaction of N 2 H 3 with HO 2 via singlet adducts to yield tHNNH + H 2 O 2 and HNN(H)O + H 2 O, while reaction on the triplet surface proceeds via a bound complex followed by a tight, submerged barrier to yield N 2 H 4 + O 2 . The results were incorporated in a detailed reaction mechanism, which was used to interpret the shock tube results from Michel and Wagner (1965) on the effect of O 2 on hydrazine conversion at 1100–1400 K. The kinetic model captured qualitatively the observed behavior, but underestimated the reaction rate under oxidizing conditions. The hydrazine pyrolysis chemistry dominated conversion at reducing conditions and/or high temperature. At oxidizing conditions and intermediate temperatures (≲ 1400 K), reactions of N 2 -amines with HO 2 and O 2 were important for the oxidation rate.

Ab initio calculations↗
DOE OSTI2021

Rheological characterization of 3D printable geopolymers

This study demonstrates a two-step approach that enables quantification of concrete printability through dynamics mode rheological measurements; I) modeling shearing history during extrusion; II) monitoring the hardening evolution of deposited material by applying a strain smaller than the critical strain. It is shown that the shearing history of the material is removed by imposing a pre-shearing above the critical strain, and zero value yield stress is measured for all specimens. At step II, a linear extrapolation of the green strength development can quantify the static yield stress at the origin, which we used to quantify the material printability. As far as this yield stress surpasses the stress level that exists in the printed structure, the material retains its shape stability. We show the performance of the model on a series of geopolymer mortars with a wide range of rheological properties as a function of material composition, aging, and pre-shearing consequences.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Hydrothermal solubility of Dy hydroxide as a function of pH and stability of Dy hydroxyl aqueous complexes from 25 to 250 °C

The rare earth elements (REE) have important applications in green energy technologies. The formation of mineral deposits in geologic systems commonly involves hydrothermal fluids which can mobilize the REE. However, the REE speciation is not well known as a function of pH. The thermodynamic properties of REE hydroxyl complexes used in geochemical models are based on the Helgeson-Kirkham-Flowers (HKF) equation of state parameters which were derived by extrapolation of low temperature experimental and estimated data. In this study, Dy hydroxide solubility experiments are combined with available literature data to improve these models from 25 to 250 °C and optimize the thermodynamic properties of Dy 3+ and Dy hydroxyl complexes using GEMSFITS. Batch-type solubility experiments were conducted from 150 to 250 °C and at saturated water vapor pressure in perchloric acid solutions with initial pH values of 2 to 5 in 0.5 pH unit increments. The measured solubility of Dy hydroxide is retrograde with temperature and decreases with pH. The logarithm of total dissolved Dy molality ranges from –2.3 to –5.3 at 150 °C (pH 4.7–5.5), from –2.4 to –5.6 at 200 °C (pH 3.9–5.1), and from –3.7 to –6.9 at 250 °C (pH of 3.4 and 5.0). The optimized standard partial molal Gibbs energies of formation (Δ f G° T ) derived for Dy 3+ and DyOH 2+ display a close to linear relationship with temperature, fitting with previous optimizations based on DyPO 4 solubility data in the literature. A comparison of the optimized ΔfG°T values for aqueous Dy species with predictions from available HKF parameters indicates significant differences ranging from +11 to –26 kJ/mol between 25 and 250 °C. The experimental fits are used to derive the Dy hydroxide solubility products (K s0 ) and formation constants for the hydrolysis of Dy (β n with n = 1 to 3; Dy 3+ + nOH – = DyOH n 3-n ) as a function of temperature. The optimization method presented yields accurate thermodynamic properties for the Dy 3+ aqua ions and the DyOH 2+ species at the acidic to mildly acidic pH studied whereas more experimental work is needed at near-neutral and alkaline conditions to better constrain the other hydroxyl complexes. Furthermore, the optimized thermodynamic data have a significant impact on geochemical modeling of the mobility and solubility of REE minerals in acidic hydrothermal fluids.

58 GEOSCIENCES↗
DOE OSTI2023

Evaluation of dual-weighted residual and machine learning error estimation for projection-based reduced-order models of steady partial differential equations

Projection-based reduced-order models (pROMs) show great promise as a means to accelerate many-query applications such as forward error propagation, solving inverse problems, and design optimization. In order to deploy pROMs in the context of high-consequence decision making, accurate error estimates are required to determine the region(s) of applicability in the parameter space. The following paper considers the dual-weighted residual (DWR) error estimate for pROMs and compares it to another promising pROM error estimate, machine learned error models (MLEM). Here, we show how DWR can be applied to ROMs and then evaluate DWR on two partial differential equations (PDEs): a two-dimensional linear convection–reaction–diffusion equation, and a three-dimensional static hyper-elastic beam. It is shown that DWR is able to estimate errors for pROMs extrapolating outside of their training set while MLEM is best suited for pROMs used to interpolate within the pROM training set.

42 ENGINEERING↗
DOE OSTI2023

A multifidelity deep operator network approach to closure for multiscale systems

Projection-based reduced order models (PROMs) have shown promise in representing the behavior of multiscale systems using a small set of generalized (or latent) variables. Despite their success, PROMs can be susceptible to inaccuracies, even instabilities, due to the improper accounting of the interaction between the resolved and unresolved scales of the multiscale system (known as the closure problem). In the current work, we interpret closure as a multifidelity problem and use a multifidelity deep operator network (DeepONet) framework to address it. In addition, to enhance the stability and/or accuracy of the multifidelity-based closure, we employ the recently developed ``in-the-loop" training approach from the literature on coupling physics and machine learning models. The resulting approach is tested on shock advection for the one-dimensional viscous Burgers equation and vortex merging for the two-dimensional Navier-Stokes equations. Furthermore, the numerical experiments show significant improvement of the predictive ability of the closure-corrected PROM over the un-corrected one both in the interpolative and the extrapolative regimes.

42 ENGINEERING↗
DOE OSTI2024

tLaSDI: Thermodynamics-informed latent space dynamics identification

Here we propose a latent space dynamics identification method, namely tLaSDI, that embeds the first and second principles of thermodynamics. The latent variables are learned through an autoencoder as a nonlinear dimension reduction model. The latent dynamics are constructed by a neural network-based model that precisely preserves certain structures for the thermodynamic laws through the GENERIC formalism. An abstract error estimate is established, which provides a new loss formulation involving the Jacobian computation of autoencoder. The autoencoder and the latent dynamics are simultaneously trained to minimize the new loss. Computational examples demonstrate the effectiveness of tLaSDI, which exhibits robust generalization ability, even in extrapolation. In addition, an intriguing correlation is empirically observed between a quantity from tLaSDI in the latent space and the behaviors of the full-state solution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Mesh-based super-resolution of fluid flows with multiscale graph neural networks

A graph neural network (GNN) approach is introduced in this work which enables mesh-based three-dimensional super-resolution of fluid flows. In this framework, the GNN is designed to operate not on the full mesh-based field at once, but on localized meshes of elements (or cells) directly. To facilitate mesh-based GNN representations in a manner similar to spectral (or finite) element discretizations, a baseline GNN layer (termed a message passing layer, which updates local node properties) is modified to account for synchronization of coincident graph nodes, rendering compatibility with commonly used element-based mesh connectivities. Furthermore, the architecture is multiscale in nature, and is comprised of a combination of coarse-scale and fine-scale message passing layer sequences (termed processors) separated by a graph unpooling layer. The coarse-scale processor embeds a query element (alongside a set number of neighboring coarse elements) into a single latent graph representation using coarse-scale synchronized message passing over the element neighborhood, and the fine-scale processor leverages additional message passing operations on this latent graph to correct for interpolation errors. Demonstration studies are performed using hexahedral mesh-based data from Taylor–Green Vortex and backward-facing step flow simulations at Reynolds numbers of 1600 and 3200. Through analysis of both global and local errors, the results ultimately show how the GNN is able to produce accurate super-resolved fields compared to targets in both coarse-scale and multiscale model configurations. Reconstruction errors for fixed architectures were found to increase in proportion to the Reynolds number. Geometry extrapolation studies on a separate cavity flow configuration show promising cross-mesh capabilities of the super-resolution strategy.

Backward-facing step↗