Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “setup”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

NO formation by N 2 /O 2 plasma catalysis: The impact of surface reactions, gas-phase reactions, and mass transport

Pathways and timescales relevant to facilitate plasma-assisted N 2 -O 2 reactions are assessed by measuring the consumption of plasma-derived N and the formation of NO in the gas phase and over Ag catalytic surfaces. These measurements are enabled by a setup that enables N 2 activation in an atmospheric pressure RF plasma jet, enables O 2 addition in the plasma afterglow, facilitates reactions over an Ag wire catalyst, and allows species density quantification by molecular beam mass spectrometry. Gas-phase reactions consume N but do not form NO with high selectivity. The presence of the non-porous Ag wire catalyst increases the rate of N conversion to NO, though mass transfer processes, not surface reactions, dictate the rate of N consumption. When O 2 concentrations and the ratio of the surface area of the catalyst to the void volume of the reactor are high (3–5 mol% O 2 , 10900 m –1 ), N conversion to NO reaches 100 % selectivity. When both N 2 and O 2 are fed through the plasma jet, gas-phase NO production increases 10×, although plasma and gas-phase processes do not exclusively produce NO. Above a threshold NO density, N cannot diffuse to the catalyst surface faster than it is consumed in the gas phase by reactions with NO. Furthermore, the use of heterogeneous catalysts to enhance plasma-driven N x O y formation and control N x O y product selectivity is limited to cases where diffusive transport of N from the gas phase to the catalyst surface is faster than consumption of N from gas-phase reactions with NO.

Engineering↗
DOE OSTI2025

Evaluation of IrO 2 catalysts doped with Ti and Nb at industrially relevant electrolyzer conditions: A comprehensive study

A series of commercial Oxygen Evolution Reaction (OER) IrO 2 -based materials doped with acid-stable titanium and niobium species were comprehensively characterized by Brunauer-Emmett-Teller (BET), X-ray diffraction analysis (XRD), X-ray photoelectron spectroscopy (XPS), transmission electron microscopy (TEM) with energy dispersive spectroscopy (EDS), and X-ray Scattering. Electrocatalysts were integrated into Membrane Electrode Assembly (MEA) using a fabrication method developed under the US DOE H2NEW consortium. An electrolysis performance in a commercial setup as well as a laboratory screening system was performed at conditions relevant to industrial application. According to the comprehensive characterizations, the studied materials are closer to doped iridium oxides rather than core–shell structures. In an electrolysis cell, the IrO 2 /TiO x catalyst slightly outperforms the IrO 2 /NbO x based on the activity. It was demonstrated that the operation of electrolysis cells at elevated temperatures and the implementation of thinner Nafion-type membranes allows for substantially increased performance, which is consistent with the literature report. Finally, this work provides valuable baselines including characterization and performance for guiding future research in this direction.

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

Capillary flow in UHPFRC with synthetic fibers, under high tensile stresses

A novel Ultra High-Performance Fiber Reinforced Concrete (UHPFRC) mix with synthetic fibers and a low clinker matrix, henceforth referred to as PE-UHPFRC, has been developed for structural applications. It exhibits a high tensile elastic limit above 7 MPa, a tensile strength of more than 10 MPa, and a very high tensile hardening domain of above 3.5%. In order to effectively use this material, its protective properties have been investigated on the basis of the effect of tensile deformation and subsequent cracking on water transport properties, for a wide range of tensile strain. A special setup was developed to measure the capillary absorption of liquids while the specimen is under tension. The results show a considerable reduction in capillary absorption in case of PE-UHPFRC compared with Strain Hardening Cement-based Composites (SHCC). Moreover, the results highlight the considerable effect of the onset of cracking on capillary absorption.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Predicting oxidation damage of ultra high-temperature carbide ceramics in extreme environments using machine learning

Determining the oxidation resistance of UHTC carbides in extreme environments is challenging theoretically and experimentally due to the high dimensional complexity of influencing variables and intricate testing setups. Herein we demonstrate the use of machine learning (ML) models trained with experimental literature data to predict the oxide thickness of UHTC carbides exposed to air based on composition, mean grain size, relative densification, holding time, and temperature. A multi-dimensional database with 76 occurrences is created containing experimental results of Hf, Zr, and Ta carbides plus additives. In this study, the preprocessed database is then used to train ML models to predict their oxidation behavior. The trained model predicts the oxidation damage in the form of an average oxide thickness in UHTC carbides with a Mean Absolute Error (MAE) of ±65.45 μm for samples in the testing set that developed thicknesses up to 1000 μm. The model successfully predicted oxidation damage for a recession rate lower than 60 μm/min. It is noticed that the ensemble method MAE is increased to ±134.34 μm while forecasting the oxidation of samples with a recession rate higher than the threshold. The unprecedented approach is a novel way to predict the damage through the oxidation of carbide compounds before processing for a smarter design with room for improvement.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Dynamic and static corrosion of chromium-alumina refractory in simulated nuclear waste glass

Monofrax® K-3, a chromium–alumina refractory, is widely used in nuclear waste glass melters due to its high durability. However, aggressive vitrification conditions still lead to corrosion that limits melter lifespan. Here, this study investigates two key corrosion modes, subsurface and melt-line corrosion, under static and dynamic conditions, using a custom setup enabling precise control of melt velocity and temperature. Experimental results, interpreted using diffusion-limited dissolution models, show that (i) subsurface corrosion increases by approximately 10% per mm/s increase in melt velocity, while melt-line corrosion remains virtually unaffected, (ii) corrosion rate dependence on temperature is inversely proportional to melt viscosity, and (iii) glass composition affects both corrosion modes similarly, based on Cr 2 O 3 , Al 2 O 3 , and alkali contents. The experimental data also allowed the determination of diffusion coefficients of major dissolving species, which were found to be consistent between the melt-line and subsurface corrosion models, indicating their suitability for future CFD studies.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2025

Quantifying particle movement in a spout-fluidized bed with irregular feedstock morphology

Here, the spout-bed fluidization behavior of nonspherical, 140 μm SiC feedstock was quantified via particle image velocimetry for varying gas distributor geometries. A bench-scale, room-temperature fluidization setup was assembled to model a 50 mm fluidized bed chemical vapor deposition (FB-CVD) system, and fluidized bed motion was captured using a high-speed camera. Modular tips with varying inlet geometries were 3D printed and tested on the bench-scale rig using identical feedstock and gas flow rates in the range of 3.0–9.0 L/min. Fluidization behavior was quantified by extracting parameters of the bed velocity, frequency, dead time, and other measurements, which were ranked for each inlet geometry configuration tested. The results from this work demonstrate that changing the path of inlet gas flow can significantly change the hydrodynamics within a spout-fluidized bed under identical feedstock, loading, and flow rate conditions, potentially enabling experimental control of particle fluidization behavior for a given condition. Moreover, composite rankings of fluidization behavior for unique distributor geometries hold potential to guide the design of FB-CVD experiments for various engineering and scientific applications.

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

Interaction between dissolution and precipitation during olivine carbonation: Implications for CO 2 mineralization

Large-scale carbonation of olivine is considered a promising approach for in situ mineral carbonation, offering a permanent and stable method for CO 2 storage. A critical aspect of this process is understanding how dissolution and precipitation interact, as this could drive fracturing and enhance further reactions. In this study, we conducted carbonation experiments on olivine using CO 2 -saturated aqueous solutions of NaHCO 3 and NaCl. Two experimental setups were used: one representing an open geochemical system and the other a closed system, corresponding to reaction-limited and flow-limited scenarios, respectively. Further, post-reaction textural analysis using scanning electron microscopy (SEM) revealed surface coatings of reaction products in the closed system, while etch pits and etch channels were prevalent in the open system. Although no direct evidence of reaction-driven fracturing was observed, etch pits and etch channels may serve as initiation points for subcritical crack formation and growth, potentially maintaining permeability and exposing new unreacted surfaces. Using linear elastic fracture mechanics (LEFM) model, we estimate that microcracks could propagate under a pressure of 0.1 GPa if reaction products accumulate within the etch pits. Our findings offer new insights into the mechanisms governing olivine carbonation.

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

Melt pool instability detection using coaxial photodiode system validated by in-situ X-ray imaging

A co-axial photodiode monitoring system with high temporal resolution has been integrated into a proven test bench enabling synchronized side-view high speed X-ray imaging of melt pool dynamics and top-view spectral emission characterization of the melt pool. Here, this setup enables direct observation of melt pool phenomena and correlation between the two monitoring systems which can be directly scaled to commercial systems. The work demonstrates a 92% detection rate in keyhole collapse phenomena related to defect generation in SLM. Furthermore, the impact of gas flow on monitoring signals is studied to understand the fundamental importance of gas flow in commercial systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A comprehensive and fair comparison of two neural operators (with practical extensions) based on $\mathrm{FAIR}$ data

Neural operators can learn nonlinear mappings between function spaces and offer a new simulation paradigm for real-time prediction of complex dynamics for realistic diverse applications as well as for system identification in science and engineering. Herein, we investigate the performance of two neural operators, which have shown promising results so far, and we develop new practical extensions that will make them more accurate and robust and importantly more suitable for industrial-complexity applications. The first neural operator, DeepONet, was published in 2019 (Lu et al., 2019), and its original architecture was based on the universal approximation theorem of Chen & Chen (1995). The second one, named Fourier Neural Operator or FNO, was published in 2020, and it is based on parameterizing the integral kernel in the Fourier space. DeepONet is represented by a summation of products of neural networks (NNs), corresponding to the branch NN for the input function and the trunk NN for the output function; both NNs are general architectures, e.g., the branch NN can be replaced with a CNN or a ResNet. According to Kovachki et al. (2021), FNO in its continuous form can be viewed conceptually as a DeepONet with a specific architecture of the branch NN and a trunk NN represented by a trigonometric basis. In order to compare FNO with DeepONet computationally for realistic setups, we develop several extensions of FNO that can deal with complex geometric domains as well as mappings where the input and output function spaces are of different dimensions. We also develop an extended DeepONet with special features that provide inductive bias and accelerate training, and we present a faster implementation of DeepONet with cost comparable to the computational cost of FNO, which is based on the Fast Fourier Transform. Here we consider 16 different benchmarks to demonstrate the relative performance of the two neural operators, including instability wave analysis in hypersonic boundary layers, prediction of the vorticity field of a flapping airfoil, porous media simulations in complex-geometry domains, etc. We follow the guiding principles of FAIR (Findability, Accessibility, Interoperability, and Reusability) for scientific data management and stewardship. The performance of DeepONet and FNO is comparable for relatively simple settings, but for complex geometries the performance of FNO deteriorates greatly. We also compare theoretically the two neural operators and obtain similar error estimates for DeepONet and FNO under the same regularity assumptions.

42 ENGINEERING↗
DOE OSTI2022

Geometric prior of multi-resolution yielding manifolds and the local closest point projection for nearly non-smooth plasticity

Elastoplasticity models often introduce a scalar-valued yield function to implicitly represent the boundary between elastic and plastic material states. Here this paper introduces a new alternative where the yield envelope is represented by a manifold of which the topology and the geometry are learned from a set of data points in a parametric space (e.g. principal stress space, π-plane). Here, deep geometric learning enables us to reconstruct a highly complex yield envelope by breaking it down into multiple coordinate charts. The global atlas that consists of these coordinate charts in return allows us to represent the yield surface via multiple overlapping patches, each with a specific local parametrization. This setup provides several advantages over the classical implicit function representation approach. For instance, the availability of coordinate charts enables us to introduce an alternative stress integration algorithm where the trial stress may project directly on a local patch and hence circumvent the issues related to non-smoothness and the lack of convexity of yield surfaces. Meanwhile, the local parametric approach also enables us to predict hardening/softening locally in the parametric space, even without complete knowledge of the yield surface. Comparisons between the classical yield function approach on the non-smooth plasticity and anisotropic cam-clay plasticity model are provided to demonstrate the capacity of the models for highly precise yield surface and the feasibility of the implementation of the learned model in the local stress integration algorithm.

42 ENGINEERING↗
DOE OSTI2022

A comprehensive study of non-adaptive and residual-based adaptive sampling for physics-informed neural networks

Physics-informed neural networks (PINNs) have shown to be effective tools for solving both forward and inverse problems of partial differential equations (PDEs). PINNs embed the PDEs into the loss of the neural network using automatic differentiation, and this PDE loss is evaluated at a set of scattered spatio-temporal points (called residual points). The location and distribution of these residual points are highly important to the performance of PINNs. However, in the existing studies on PINNs, only a few simple residual point sampling methods have mainly been used. Here, we present a comprehensive study of two categories of sampling for PINNs: non-adaptive uniform sampling and adaptive nonuniform sampling. We consider six uniform sampling methods, including (1) equispaced uniform grid, (2) uniformly random sampling, (3) Latin hypercube sampling, (4) Halton sequence, (5) Hammersley sequence, and (6) Sobol sequence. We also consider a resampling strategy for uniform sampling. To improve the sampling efficiency and the accuracy of PINNs, we propose two new residual-based adaptive sampling methods: residual-based adaptive distribution (RAD) and residual-based adaptive refinement with distribution (RAR-D), which dynamically improve the distribution of residual points based on the PDE residuals during training. Hence, we have considered a total of 10 different sampling methods, including six non-adaptive uniform sampling, uniform sampling with resampling, two proposed adaptive sampling, and an existing adaptive sampling. We extensively tested the performance of these sampling methods for four forward problems and two inverse problems in many setups. Our numerical results presented in this study are summarized from more than 6000 simulations of PINNs. Here, we show that the proposed adaptive sampling methods of RAD and RAR-D significantly improve the accuracy of PINNs with fewer residual points for both forward and inverse problems. Furthermore, the results obtained in this study can also be used as a practical guideline in choosing sampling methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Multi-output multilevel best linear unbiased estimators via semidefinite programming

Multifidelity forward uncertainty quantification (UQ) problems often involve multiple quantities of interest and heterogeneous models (e.g., different grids, equations, dimensions, physics, surrogate and reduced-order models). While computational efficiency is key in this context, multi-output strategies in multilevel/multifidelity methods are either sub-optimal or non-existent. In this paper we extend multilevel best linear unbiased estimators (MLBLUE) to multi-output forward UQ problems and we present new semidefinite programming formulations for their optimal setup. Not only do these formulations yield the optimal number of samples required, but also the optimal selection of low-fidelity models to use. While existing MLBLUE approaches are single-output only and require a non-trivial nonlinear optimization procedure, the new multi-output formulations can be solved reliably and efficiently. Here, we demonstrate the efficacy of the new methods and formulations in practical UQ problems with model heterogeneity.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Hutchinson Trace Estimation for high-dimensional and high-order Physics-Informed Neural Networks

Physics-Informed Neural Networks (PINNs) have proven effective in solving partial differential equations (PDEs), especially when some data are available by seamlessly blending data and physics. However, extending PINNs to high-dimensional and even high-order PDEs encounters significant challenges due to the computational cost associated with automatic differentiation in the residual loss function calculation. Herein, we address the limitations of PINNs in handling high-dimensional and high-order PDEs by introducing the Hutchinson Trace Estimation (HTE) method. Starting with the second-order high-dimensional PDEs, which are ubiquitous in scientific computing, HTE is applied to transform the calculation of the entire Hessian matrix into a Hessian vector product (HVP). This approach not only alleviates the computational bottleneck via Taylor-mode automatic differentiation but also significantly reduces memory consumption from the Hessian matrix to an HVP’s scalar output. We further showcase HTE’s convergence to the original PINN loss and its unbiased behavior under specific conditions. Comparisons with the Stochastic Dimension Gradient Descent (SDGD) highlight the distinct advantages of HTE, particularly in scenarios with significant variability and variance among dimensions. We further extend the application of HTE to higher-order and higher-dimensional PDEs, specifically addressing the biharmonic equation. By employing tensor-vector products (TVP), HTE efficiently computes the colossal tensor associated with the fourth-order high-dimensional biharmonic equation, saving memory and enabling rapid computation. The effectiveness of HTE is illustrated through experimental setups, demonstrating comparable convergence rates with SDGD under memory and speed constraints. Additionally, HTE proves valuable in accelerating the Gradient-Enhanced PINN (gPINN) version as well as the Biharmonic equation. Overall, HTE opens up a new capability in scientific machine learning for tackling high-order and high-dimensional PDEs.

Curse of dimensionality↗
DOE OSTI2022

Diffuse back-illumination temperature imaging (DBI-TI), a novel soot thermometry technique

To meet stringent emissions regulations on soot emissions, it is critical to further advance the fundamental understanding of in-cylinder soot formation and oxidation processes. Among several optical techniques for soot quantification, diffuse back-illumination extinction imaging (DBI-EI) has recently gained traction mainly due to its ability to compensate for beam steering, which if not addressed, can cause unacceptably high measurement uncertainty. Until now, DBI-EI has only been used to measure the amount of soot along the line of sight, and here, we extend the capabilities of a DBI-EI setup to also measure in-cylinder soot temperature. This proof of concept of diffuse back-illumination temperature imaging (DBI-TI) as a soot thermometry technique is presented by implementing DBI-TI in a single cylinder, heavy-duty, optical diesel engine to provide 2-D line-of-sight integrated soot temperature maps. The potential of DBI-TI to be an accurate thermometry technique for use in optical engines is analyzed. The achievable accuracy is due in part to simultaneous measurement of the soot extinction, which circumvents the uncertainty in dispersion coefficients that depend on the optical properties of soot and the wavelength of light utilized. Analysis shows that DBI-TI provides temperature estimates that are closer to the mass-averaged soot temperature when compared to other thermometry techniques that are more sensitive to soot temperature closer to the detector. Furthermore, uncertainty analysis and Monte Carlo (MC) simulations provide estimates of the temperature measurement errors associated with this technique. The MC simulations reveal that for the light intensities and optical densities encountered in these experiments, the accuracy of the DBI-TI technique is comparable or even better than other established optical thermometry techniques. Thus, DBI-TI promises to be an easily implementable extension to the existing DBI-EI technique, thereby extending its ability to provide comprehensive line-of-sight integrated information on soot.

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

A detailed experimental and modeling comparison of molecular radiative heat loss in a spark-ignition engine

Radiative heat transfer has been chiefly considered negligible in internal combustion engines, except for Diesel engines where soot radiation was recognized as a significant radiative transfer source. Only more recently, detailed simulations have shown that molecular radiation can be substantial as well. In extension to this, molecular radiative heat transfer can reach detectable levels of about 5–10% of the total heat transfer in spark-ignited engines. The broadband radiative nature of the significant emitting molecules, carbon dioxide and water, makes it necessary to address whether radiative trapping plays a substantial role in either total heat loss or energy redistribution within the cylinder. An experimental setup that allows measurements of the infrared emissions at a 2-crank-angle-degree resolution was developed to determine the importance of radiative trapping by carbon dioxide and water. Measurements were conducted in the well-characterized and documented TCC-III engine at the University of Michigan. The engine operated on a stoichiometric propane/air mixture at a speed of 1,300 RPM. Large Eddy Simulations with added line-by-line photon Monte-Carlo simulations of the molecular radiation were an integral part of this study to plan and devise the measurements. Post-processing of the simulation data included an accurate representation of the experimental volume from where infrared signals are collected. Additionally, the experimental data were used for validation of the photon Monte-Carlo simulations. Joint analysis of experimental and simulated spectra allowed quantifying the significance of radiative trapping. Results show the significant role of radiative trapping when predicting radiative heat transfer in the TCC-III engine.

42 ENGINEERING↗
DOE OSTI2022

A comprehensive model to capture electrical discharge and spark channel evolution during spark-ignition processes

Emissions reduction through engine efficiency improvements is a priority for automakers who have turned to unconventional engine operation such as highly dilute, boosted, and stratified charge. Given its importance to flame initiation and sustained turbulent flame propagation, reliable and accurate spark ignition models are necessary to design ignition systems that reduce cyclic indicated mean effective pressure (IMEP) variability and increase engine efficiency in these operation modes. In this paper, secondary electric circuit, short-circuit, blowout, and re-strike sub-models are added to the Lagrangian-Eulerian spark ignition (LESI) model to simulate electrical discharge and spark channel elongation in an inert cross-flow combustion vessel. First, the physics of the underlying sub-models are described and the governing equations discussed with the spark channel voltage expression playing a critical role. Then, the experimental and simulations setups are presented. The results section begins with the derivation of spark channel voltage from experimental results. Then, electrical discharge and spark channel elongation simulations in inert flows are carried out using LESI and compared against experimental results. Finally, the results validate the model's ability to accurately predict spark channel elongation, as well as the occurrence of short-circuits, blowouts, re-strikes, and end of discharge.

42 ENGINEERING↗
DOE OSTI2025

CASM Monte Carlo: Calculations of the thermodynamic and kinetic properties of complex multicomponent crystals

Monte Carlo techniques play a central role in statistical mechanics approaches that connect macroscopic thermodynamic and kinetic properties to the electronic structure of a material. This paper describes the implementation of Monte Carlo techniques for the study of multicomponent crystalline materials within the Clusters Approach to Statistical Mechanics (CASM) software suite, and demonstrates their use in model systems to calculate free energies and kinetic coefficients, study phase transitions, and construct phase diagrams from first principles. Many crystal structures are complex, with multiple sublattices occupied by differing sets of chemical species, along with the presence of vacancies or interstitial species. This imposes constraints on concentration variables, the form of thermodynamic potentials, and the values of kinetic transport coefficients. The framework used by CASM to formulate thermodynamic potentials and kinetic transport coefficients accounting for arbitrarily complex crystal structures is presented and demonstrated with examples of increasing complexity. Additionally, an overview of the capabilities of the CASM software specific to Monte Carlo methods is given, and a new CASM software package is introduced, casm-flow, which helps automate the setup, submission, management, and analysis of Monte Carlo simulations.

Cluster expansion↗