Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “example”

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,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Microbubble dynamics and heat transfer in boiling droplets

Dissipating large heat fluxes from a surface is critically important in numerous industrial and natural applications. Boiling based spray cooling and surface texturing are two of the most promising methods being investigated to address this problem. Although our understanding on these topics has significantly improved over past decades, critical gaps remain in the knowledgebase stymieing the realization of their full potential. As an example, while bubble growth in pool boiling have been investigated in detail, comparatively little is known about how the bubbles evolve inside boiling drops. In the present work, we have investigated for the first time, the microbubble dynamics inside water droplets boiling on superhydrophilic textured substrates using high-speed X-ray phase contrast imaging (XRPCI). Our observations show that the transient bubble density variation follows similar characteristics irrespective of the texture spacing at a given surface temperature. For an example microstructure, we found that the number of discrete bubbles on the surface decreases as temperature is increased although their growth rate increases. Here, we observe that bubble growth is highly non-uniform during the lifetime of a drop on the surface. Initially, bubbles grow under diffusion-limited regime, but at later times they grow as ~t 1.45 due to combined effects of coalescence and evaporation. In some conditions, we found that bubbles shrink dramatically af- ter the initial growth spurt presumably due to severe quenching of the surface, and migration of bubbles on the surface. Using the bubble sizes, for the first time we analyzed the heat flux removed by a single bubble and also by all the bubbles at a given time. We find that the highest dissipation through latent- heat component (~600 W/cm 2 ) occurs just in the beginning and thereafter it decreases. We expect that our findings and the analysis would guide further work on the topic and will aid in the overarching goal of engineering surfaces that are more efficient in boiling heat transfer.

42 ENGINEERING↗
DOE OSTI2023

Training material models using gradient descent algorithms

High temperature design requires accurate constitutive models to describe material inelastic deformation and failure behavior. Oftentimes, calibrating accurate models devolves into the problem of fitting the model parameters against experimental test data. Here, we present the pyopmat package, an open source framework for calibrating constitutive models against experiment data subjected to various loading conditions using machine learning techniques. The package calculates the exact gradient of the model response with respect to the parameters using a combination of automatic differentiation and the adjoint method. Given this exact gradient, we compare the performance of several gradient-based optimization techniques in fitting realistic constitutive models against data. Here, we demonstrate the efficiency and accuracy of our package through example problems using both synthetic data, generated using known parameter sets, under monotonic and cyclic loading conditions and also with an example applying the techniques developed here to actual high temperature creep-fatigue test data.

36 MATERIALS SCIENCE↗
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 OSTI2020

Persistent compositions of non-stoichiometric compounds with low bulk diffusivity: A theory and application to Nb3Sn superconductors

Non-stoichiometric compounds may develop a composition gradient when they are formed by reactive diffusional processes. This paper reports an interesting phenomenon that in compounds with low bulk diffusivities, which rely mainly on grain boundary diffusion for their growth, the final bulk compositions may be far from equilibrium, with a very low bulk diffusivity leading to fixed (persistent) compositions – one such example is Nb 3 Sn, a superconductor. We investigated the microchemistry at the reactive interface using atom probe tomography to clarify the diffusion reaction mechanism for this low-bulk-diffusivity case and thus propose a theory for what determines the compound composition, using Nb 3 Sn as an example for concreteness. Using certain approximations, we derive an explicit analytical equation that illustrates what factors determine its composition profile. We compare our model with the known facts of Nb 3 Sn and see good agreement. In particular, this model predicts that internal oxidation may lead to higher Sn contents than conventional, non-oxidized Nb 3 Sn. Our measurements show that this is indeed true, and that the internally-oxidized Nb 3 Sn also has higher upper critical fields, achieving up to 28.2 T at 4.2 K. We discuss the general applicability of this model to non-stoichiometric compounds with low bulk diffusivity, and propose it as a tool to help in the design and processing of such materials for compositional control.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A generalized approximate control variate framework for multifidelity uncertainty quantification

We describe and analyze a variance reduction approach for Monte Carlo (MC) sampling that accelerates the estimation of statistics of computationally expensive simulation models using an ensemble of models with lower cost. These lower cost models — which are typically lower fidelity with unknown statistics — are used to reduce the variance in statistical estimators relative to a MC estimator with equivalent cost. We derive the conditions under which our proposed approximate control variate framework recovers existing multifidelity variance reduction schemes as special cases. We demonstrate that existing recursive/nested strategies are suboptimal because they use the additional low-fidelity models only to efficiently estimate the unknown mean of the first low-fidelity model. As a result, they cannot achieve variance reduction beyond that of a control variate estimator that uses a single low-fidelity model with known mean. However, there often exists about an order-of-magnitude gap between the maximum achievable variance reduction using all low-fidelity models and that achieved by a single low-fidelity model with known mean. We show that our proposed approach can exploit this gap to achieve greater variance reduction by using non-recursive sampling schemes. The proposed strategy reduces the total cost of accurately estimating statistics, especially in cases where only low-fidelity simulation models are accessible for additional evaluations. Several analytic examples and an example with a hyperbolic PDE describing elastic wave propagation in heterogeneous media are used to illustrate the main features of the methodology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Mesoscale informed parameter estimation through machine learning: A case-study in fracture modeling

Scale bridging is a critical need in computational sciences, where the modeling community has developed accurate physics models from first principles, of processes at lower length and time scales that influence the behavior at the higher scales of interest. However, it is not computationally feasible to incorporate all of the lower length scale physics directly into upscaled models. This is an area where machine learning has shown promise in building emulators of the lower length scale models, which incur a mere fraction of the computational cost of the original higher fidelity models. We demonstrate the use of machine learning using an example in materials science estimating continuum scale parameters by emulating, with uncertainties, complicated mesoscale physics. Additionally, we describe a new framework to emulate the fine scale physics, especially in the presence of microstructures, using machine learning, and showcase its usefulness by providing an example from modeling fracture propagation. Our approach can be thought of as a data-driven dimension reduction technique that yields probabilistic emulators. Our results show well-calibrated predictions for the quantities of interests in a low-strain simulation of fracture propagation at the mesoscale level. Furthermore, on average, we achieve ~10% relative errors on time-varying quantities like total damage and maximum stresses. Successfully replicating mesoscale scale physics within the continuum models is a crucial step towards predictive capability in multi-scale problems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Sylvester-preconditioned adaptive-rank implicit time integrators for advection-diffusion equations with variable coefficients

Here, we consider the adaptive-rank integration of multi-dimensional time-dependent advection-diffusion partial differential equations (PDEs) with variable coefficients. We employ a standard finite-difference method for spatial discretization coupled with high-order diagonally implicit Runge-Kutta temporal schemes. The discrete equation is a generalized Sylvester equation (GSE), which we solve with a projection-based adaptive-rank algorithm structured around two key strategies: (i) constructing dimension-wise subspaces using a novel atypical extended Krylov strategy, and (ii) efficiently solving the basis coefficient matrix with a preconditioned GMRES solver. The low-rank decomposition is performed in 2D using SVD and with high-order SVD (HOSVD) in 3D to represent the tensor in a compressed Tucker format. For d-dimensional problems (here, d = 2 or 3), the computational complexity and memory storage of the approach are found numerically to scale as and $\mathscr{O}(Nr^2) + \mathscr{O} (r^{d+1})$ and $\mathscr{O}(Nr) + \mathscr{O} (r^{d})$, respectively, with the one-dimensional resolution and the maximal rank during the Krylov iteration (which we find to be largely independent of on our numerical examples). We present numerical examples that illustrate the advertised properties of the algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Microchannel reactive distillation for the conversion of aqueous ethanol to ethylene

Here we demonstrate the proof-of-concept for microchannel reactive distillation for alcohol-to-jet application: combining ethanol/water separation and ethanol dehydration in one unit operation. Ethanol is first distilled into the vapor phase, converted to ethylene and water, and then the water co-product is condensed to shift the reaction equilibrium. Process intensification is achieved through rapid mass transfer—ethanol stripping from thin wicks using novel microchannel architectures—leading to lower residence time and improved separation efficiency. Energy savings are realized with integration of unit operations. For example, heat of condensing water can offset vaporizing ethanol. Furthermore, the dehydration reaction equilibrium shifts towards completion by immediate removal of the water byproduct upon formation while maintaining aqueous feedstock in the condensed phase. For aqueous ethanol feedstock (40% w ), 71% ethanol conversion with 91% selectivity to ethylene was demonstrated at 220 °C, 600 psig, and 0.28 h -1 wt hour space velocity. 2.7 stages of separation were also demonstrated, under these conditions, using a device length of 8.3 cm. This provides a height equivalent of a theoretical plate (HETP), a measure of separation efficiency, of ~3.3 cm. By comparison, conventional distillation packing provides an HETP of ~30 cm. Thus, 9.1× reduction in HETP was demonstrated over conventional technology, providing a means for significant energy savings and an example of process intensification. Finally, preliminary process economic analysis indicates that by using microchannel reactive distillation technology, the operating and capital costs for the ethanol separation and dehydration portion of an envisioned alcohol-to-jet process could be reduced by at least 35% and 55%, respectively, relative to the incumbent technology, provided future improvements to microchannel reactive distillation design and operability are made.

10 SYNTHETIC FUELS↗
DOE OSTI2024

Reactive flash sintering and characterization of bulk high entropy nitrides

Over the past decade, numerous high-entropy ceramics have been synthesized, often displaying attractive properties. However, the study on facile preparation of bulk high entropy nitrides (HEN) are limited, despite its broad potential applications. This research demonstrates for the first time rapid fabrication (within ∼6 min) of bulk high-entropy nitrides, especially (Al 0.17 Nb 0.17 Ta 0.17 Ti 0.32 Zr 0.17 )N, from binary nitride powder mixtures using a highly efficient reactive flash sintering (RFS) technique. X-ray diffraction (XRD) shows the HENs from RFS are near single-phase solid solutions with a rock salt crystal structure, while in situ synchrotron study carried out during RFS captured in real time the formation of HEN, which was preserved upon cooling, suggesting thermodynamic stability of the HEN phase, even up to extreme pressure (∼35.6 GPa). Microscopic analyses using SEM, STEM, and EDS reveal decent uniformity for HEN with no obvious segregation of elements, even to submicron scale. Some properties of the obtained bulk HENs are consistent with expectations. For example, their hardness and bulk modulus are close to estimates based on rule-of-mixture (ROM) values from the constituent binary nitrides. Meanwhile, some other measured properties seem to show surprises. For example, the fracture toughness for the HENs (e.g., 7.81 ± 1.40 MPa•m 1/2 or higher) turns out to be more than double of the expected ROM estimates. The significantly improved fracture toughness is attributed to the observed nano-layered structure of the HENs, despite the HEN’s cubic crystal structure and high hardness. In addition, the oxidation resistance shows improvement up till ∼800°C, possibly due to Ta doping that suppress oxygen vacancy formation in the oxide shell, while the 5-metal HEN of (Al 0.17 Nb 0.17 Ta 0.17 Ti 0.32 Zr 0.17 )N displays superconductivity (T c of ∼5–7 K from magnetism and resistivity measurements, slightly lower than ROM estimate), despite insulating property of starting AlN. Furthermore, future study combining experimental investigation using larger samples to confirm the observed increase in fracture toughness and oxidation resistance, theoretical modeling at different length scale, and more detailed structural/chemical characterization, especially at the atomic scale, are all needed to fully understand the inter-relationships between composition, processing, structure, and novel properties for these HENs and the development of related new materials for different applications.

Flash sintering↗
DOE OSTI2020

Comparison of four learning-based methods for predicting groundwater redox status

Knowing the location where groundwater denitrification occurs, or by proxy the groundwater redox status (oxic, mixed, and anoxic), is valuable information for assessing and managing potential agricultural land-use impacts on freshwater quality. We compare the efficacy of supervised (Linear Discriminant Analysis LDA; Boosted Regression Trees, BRT; and Random Forest, RF) and unsupervised (Modified Self-Organizing Map, MSOM) learning-based methods to predict groundwater redox status in the agriculturally dominated Tasman, Waikato, and Wellington regions of New Zealand. Thresholds applied to regional groundwater-quality samples provide redox status variables and learn heuristics constrained by these variables and applied to spatial factors (climate, elevation, geologic, hydrology soils, and well depth) identify optimal sets of regional predictor variables. A split-sample approach is used to train and test the learning methods ability to predict redox status using the optimal predictor variables. Overall, the supervised methods demonstrate a prediction bias toward oxic conditions and inability to perform statistically well when using independent regional data; for example, consider kappa statistics for BRT (Tasman: 0.42, Waikato: 0.38, Wellington: 0.17), RF (Tasman: 0.42, Waikato: 0.47, Wellington: 0.17 and LDA (Tasman: 0.46, Waikato: 0.32, Wellington: 0.17). By contrast, the unsupervised method performs statistically well when predicting oxic, mixed, and anoxic conditions and corresponding depths when using independent regional data; for example, consider MSOM kappa statistics for Tasman: 0.78, Waikato: 0.80, Wellington: 0.76. The unsupervised learning method provides the added benefits of being (1) able to combine predictions into 3D regional anoxic probability plots for interpreting the spatial influence of paleosols and groundwater flowpaths on redox status, and (2) readily extended to map 3D redox status across New Zealand and other countries despite data bias and sparsity.

machine learning algorithms, Random Forest, booste↗
DOE OSTI2023

Even spheres as joint spectra of matrix models

The Clifford spectrum is a form of joint spectrum for noncommuting matrices. This theory has been applied in photonics, condensed matter and string theory. In applications, the Clifford spectrum can be efficiently approximated using numerical methods, but this only is possible in low dimensional example. In this paper we examine the higher-dimensional spheres that can arise from theoretical examples. We also describe a constructive method to generate five real symmetric almost commuting matrices that have a K-theoretical obstruction to being close to commuting matrices. For this, we look to matrix models of topological electric circuits.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Reassessing early-age strength development of high-volume fly ash concretes for precast buildings

Increasing beneficial use of fresh or landfilled fly ash as a replacement for Portland cement can be more challenging for the construction of precast buildings or similar applications requiring rapid strength development. Therefore, the framework presented in this paper aims to reassess high-volume fly ash concretes but in the context of facilitating more sustainable precast buildings. More specifically, the framework was used to characterize strength development of concrete mixes with a target minimum 24-hour compressive strength of 24.1 MPa (3500 psi), selected as an example strength development metric to demonstrate the framework, and comprised of 40% Class C, Class F, and landfilled (harvested) fly ash – as a high-volume replacement of Type III or Type IL cement. High-early strength was driven by optimized dosages of commercial grade gypsum and accelerating admixtures, in addition to optimal aggregate packing and mix proportioning strategies. Early-age mechanical properties including compressive strength, modulus of rupture, and modulus of elasticity were reevaluated within 24 hours of batching with respect to common precast production demands. Simple data analyses were then used to highlight cases where currently accepted design provisions for the aforementioned properties are either overly-conservative or unconservative with respect to test data. Furthermore, the framework and demonstration of example mixes presented herein aim to promote confidence for using larger fractions of fresh or landfilled fly ashes for precast buildings to further enhance environmental benefits without sacrificing pertinent early-age structural performance.

42 ENGINEERING↗
DOE OSTI2020

UQpy: A general purpose Python package and development environment for uncertainty quantification

In this paper, we present the UQpy software toolbox, an open-source Python package for general uncertainty quantification (UQ) in mathematical and physical systems. The software serves as both a user-ready toolbox that includes many of the latest methods for UQ in computational modeling and a convenient development environment for Python programmers advancing the field of UQ. The paper presents an introduction to the software's architecture and existing capabilities, divided in the code in a set of modules centered around different UQ tasks such as sampling methods, generation of random processes and random fields, probabilistic inverse modeling, reliability analysis, surrogate modeling, and active learning. The paper also highlights the importance of the RunModel module, which is used to drive simulations in the uncertainty analyses performed in UQpy. This module conveniently allows the user to define computational models directly in Python, or to run simulations from a third-party software in serial or in parallel. To illustrate the various capabilities, two examples are tracked throughout the paper and analyzed repeatedly for various UQ tasks. The first is a Python model solving a nonlinear structural dynamics problem, used to illustrate UQpy's capabilities in sampling and forward propagation of high dimensional random vectors (stochastic processes), and probabilistic inference. The second model is a third-party Abaqus finite element model solving the thermomechanical response of a beam structure. This example is used to illustrate UQpy's capabilities in variance reduction sampling techniques, reliability analysis, surrogate modeling and active learning techniques.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Compositionally complex perovskite oxides: Discovering a new class of solid electrolytes with interface-enabled conductivity improvements

Compositionally complex ceramics (CCCs), including high-entropy ceramics, offer a vast, unexplored compositional space for materials discovery. Herein, we propose and demonstrate strategies for tailoring CCCs via a combination of non-equimolar compositional designs and control of grain boundaries (GBs) and microstructures. Using oxide solid electrolytes for all-solid-state batteries as an example, we have discovered a class of compositionally complex perovskite oxides (CCPOs) with improved lithium ionic conductivities beyond the limit of conventional doping. For example, we demonstrate that the ionic conductivity can be improved by >60% in (Li 0.375 Sr 0.4375 )(Ta 0.375 Nb 0.375 Zr 0.125 Hf 0.125 )O 3-δ compared with the (Li 0.375 Sr 0.4375 )(Ta 0.75 Zr 0.25 )O 3-δ (LSTZ) baseline. Furthermore, the ionic conductivity can be improved by another >70% via quenching, achieving >270% of the LSTZ. Notably, we demonstrate GB-enabled conductivity improvements via both promoting grain growth and altering GB structures through compositional designs and processing. In a broader perspective, this work suggests new routes for discovering and tailoring CCCs for energy storage and many other applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A General Strategy for Batch Development of High-Performance and Cost-Effective Sodium Layered Cathodes

High-performance and low-cost transition metal (TM) layered oxides using earth abundant elements are promising cathodes for Na-ion batteries. However, it is challenging to obtain desired materials because the large Na size, different Na occupations and various layer stacking sequences multiply the complication in determining the structure of a given composition and exacerbate uncertainty to the structure-property correlation. In this work, we use the attainment of desired NaxMnyNiZTM1-y-zO2-based cathode materials as an example to demonstrate a general roadmap for batch development of sodium layered cathodes towards practical applications. A synthesis phase diagram of NaxMnyNi1-yO2 was created for pre-screening and rational selection of the platform material of P2/O3-structured Na0.85Mn0.6Ni0.4O2. Cationic potential was leveraged in elemental substitution to further promote the material structural stability and electrochemical performance. Several cost-effective O3 and P2/O3 hybrid cathode materials have been obtained, all of which demonstrate excellent performance. In particular, the Na0.85Mn0.5Ni0.4Ti0.1O2 delivers a high specific capacity of ~130 mAh/g between 2-4 V and 91% retention after 500 cycles. The work discovers multiple materials as high-performance and cost-effective Na-ion battery cathodes and offers critical guidance to the rational design of future layered cathode materials.

Xiao, Biwei↗
DOE OSTI2024

Lie–Poisson Neural Networks (LPNets): Data-based computing of Hamiltonian systems with symmetries

An accurate data-based prediction of the long-term evolution of Hamiltonian systems requires a network that preserves the appropriate structure under each time step. Every Hamiltonian system contains two essential ingredients: the Poisson bracket and the Hamiltonian. Hamiltonian systems with symmetries, whose paradigm examples are the Lie–Poisson systems, have been shown to describe a broad category of physical phenomena, from satellite motion to underwater vehicles, fluids, geophysical applications, complex fluids, and plasma physics. The Poisson bracket in these systems comes from the symmetries, while the Hamiltonian comes from the underlying physics. We view the symmetry of the system as primary, hence the Lie–Poisson bracket is known exactly, whereas the Hamiltonian is regarded as coming from physics and is considered not known, or known approximately. Using this approach, we develop a network based on transformations that exactly preserve the Poisson bracket and the special functions of the Lie–Poisson systems (Casimirs) to machine precision. We present two flavors of such systems: one, where the parameters of transformations are computed from data using a dense neural network (LPNets), and another, where the composition of transformations is used as building blocks (G-LPNets). We also show how to adapt these methods to a larger class of Poisson brackets. We apply the resulting methods to several examples, such as rigid body (satellite) motion, underwater vehicles, a particle in a magnetic field, and others. The methods developed in this paper are important for the construction of accurate data-based methods for simulating the long-term dynamics of physical systems.

Data-based modeling↗
DOE OSTI2023

G4CMP: Condensed matter physics simulation using the G$\scriptsize{\mathrm{EANT}}$4 toolkit

G4CMP simulates phonon and charge transport in cryogenic semiconductor crystals using the Geant4 toolkit. The transport code is capable of simulating the propagation of acoustic phonons as well as electron and hole charge carriers. Herein processes for anisotropic phonon propagation, oblique charge-carrier propagation, and phonon emission by accelerated charge carriers are included. The simulation reproduces theoretical predictions and experimental observations such as phonon caustics, heat-pulse propagation times, and mean charge-carrier drift velocities. In addition to presenting the physics and features supported by G4CMP, this report outlines example applications from the dark matter and quantum information science communities. These communities are applying G4CMP to model and design devices for which the energy transported by phonons and charge carriers is germane to the performance of superconducting instruments and circuits placed on silicon and germanium substrates. The G4CMP package is available to download from GitHub: github.com/kelseymh/G4CMP.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗