Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “engineering thermodynamics”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Seeding the explosion of a high-current-density conductor in a controlled manner through the addition of micron-scale surface defects

Inhomogeneities in a current-carrying conductor promote non-uniform heating and expansion through the complex feedback between current density, electrical resistivity, Ohmic heating, temperature, and hydrodynamics. Three-dimensional-magnetohydrodynamic (3D-MHD) simulations suggest that μm-scale resistive inclusions or voids seed local overheating and through hydrodynamic explosion generate continuously growing craters which become several times larger than the initial perturbation. The ejected mass is the genesis of an electrothermally driven plasma filament which develops at lower current than plasmas on uniform surfaces adjacent to the defect. This result suggests that 1D or even 2D treatments are largely inadequate for detailed prediction of plasma formation. To test computational predictions, z-pinch experiments driven to 1 MA studied ultra-high-purity aluminum rods which were then machined to include pairs of quasi-hemispherical voids or “engineered defects (ED)” on the current-carrying surface. ED are the dominant current-density perturbation and reproducibly drive local overheating which can be compared with 3D-MHD simulation. Data from high-resolution-gated imagers of visible surface emissions confirm many simulation predictions, including the surface topography of local overheating, and the propensity for neighboring ED to prematurely source plasmas which then connect to form a plasma filament. Results also provide conditional support of theory which suggests heating similarity; that is, heating is independent of ED size for geometrically scaled ED.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Diffusion coefficient and nucleation density studies on electrochemical deposition of aluminum from chloroaluminate ionic liquid electrolytes

In this study, aluminum was electrodeposited from ionic liquids comprising a melt of alkyl imidazolium chlorides and aluminum chloride (AlCl 3 ). The ionic liquids utilized in this study were 1-ethyl-3-methyl-imidazolium chloride (EMIC)-AlCl 3 , 1-butyl-3-methylimidazolium chloride (BMIC)-AlCl 3 , and 1-hexyl-3-methyl-imidazolium chloride (HMIC)-AlCl 3 at the AlCl 3 mole fraction of 0.667 (molar ratio of 1:2). The electrochemical behavior of chloroaluminate species in the three ionic liquids was investigated by cyclic voltammetry (CV) and chronoamperometry (CA) techniques at different temperatures. The cyclic voltammograms indicated that the reduction of Al 2 Cl 7 – species to metallic aluminum followed a diffusion-controlled phenomenon. Moreover, even at a less negative applied potential, the higher current density was obtained for EMIC-AlCl 3 ionic liquid. That indicates that EMIC-AlCl 3 favors less energy consumption during the electrodeposition. The chronoamperometric analysis revealed that the onset of aluminum deposition from such ionic liquids proceeds via a three-dimensional instantaneous nucleation process. The concentrations of Al 2 Cl 7 – ions in the ionic liquids were calculated based on the thermodynamic data, which were 2520, 2271, and 2035 mol m –3 for EMIC-AlCl 3 , BMIC-AlCl 3 , and HMIC-AlCl 3 , respectively. The diffusion coefficient (D) values of Al 2 Cl 7 – species in such ionic liquids were also calculated at various temperatures. The D values determined from the CA technique at 363 K are 2.16 × 10 –11 , 1.03 × 10 –11 , and 0.87 × 10 –11 m 2 s –1 for EMIC-AlCl 3 , BMICAlCl 3 , and HMIC-AlCl 3 , respectively. In addition, the D value increased as temperature increased and decreased as the hydrocarbon group in the ionic liquid increased. As a result, the calculated grain sizes of nucleation range from 2 to 4 μm for these ionic liquids and are in good agreement with experimental data obtained from SEM micrographs.

42 ENGINEERING↗
DOE OSTI2020

Primal interface debonding formulation for finite strain isotropic plasticity

In this work, a framework is developed for modeling ductile damage of nonlinear materials whose plastic deformation is characterized using rate independent classical plasticity. This method relies on the assumption that the free energy can be decomposed into elastic, plastic and damage parts. A thermodynamically consistent method is derived which satisfies the second law of thermodynamics in the Clausius–Duhem inequality form. The dissipation associated with plasticity takes place in the domain only, while damage dissipation is localized to the interface. The method is developed using Variational Multiscale ideas to obtain definitions of the interface fluxes within a primal formulation analogous to the Discontinuous Galerkin method, which ensures weakly vanishing interface gap prior to reaching a damage initiation criterion. The local nonlinear problem to calculate both plastic deformation gradient and damage variable follows an incremental approach similar to classical plasticity return mapping algorithm. This elastoplastic damage formulation is developed for material undergoing finite strain, and it naturally accommodates a trapezoidal traction separation law (TSL) whose shape can be varied to model either ductile interface behavior or brittle interface behavior. The formulation's performance is assessed through modeling a patch test and a compact tension specimen.

42 ENGINEERING↗
DOE OSTI2020

Effect of Molecular Structure of Quinones and Carbon Electrode Surfaces on the Interfacial Electron Transfer Process

Quinones can undergo thermodynamically reversible proton-coupled electron transfer reactions and are being applied as electroactive compounds in aqueous organic batteries. However, the electrochemical reversibility of these compounds is affected not only by their molecular structure but also by the properties of a carbon-based electrode surface. This study combines experimental and theoretical approaches to understand this dependence. Specifically, we study the electron transfer kinetics of two synthesized quinone derivatives and two commercially available ones with a glassy carbon, a highly ordered pyrolytic graphite, and a high-edge-density graphite electrode (HEDGE). The electrochemical reversibility is notably improved on the HEDGE, which shows a higher density of defects and presents oxygenated functional groups at its surface. The electron transfer kinetics are controlled by adsorbed species onto the HEDGE. Molecular dynamics simulation and quantum mechanics calculations suggest defects with oxygen-containing functional groups, such as C-O and C$=$O, on HEDGE surfaces drive the interaction with the functional groups of the molecules, during physisorption from van der Waals forces. The presence of sulfonic acid side groups and a greater number of aromatic rings in the molecular structure may contribute to a higher stabilization of quinone derivatives on HEDGEs. We propose that high-performance carbon-based electrodes can be obtained without catalysts for organic batteries, by the engineering of carbon-based surfaces with edge-like defects and oxygenated functional groups.

25 ENERGY STORAGE↗
DOE OSTI2021

Best practices in catalysis: A perspective

Catalysis, from its roots in petrochemical refining and conversion, has emerged as a transdisciplinary field that now encompasses synthesis of materials and molecules that enable applications in energy conversion and storage, environmental remediation, medicine, plastics, and fertilizer production, among numerous others. A handful of disciplines can claim relevance and success over such an extended period of time and continue to claim a preeminent role in defining the state-of-the-art in science and technology. Syncretic and rapid advancements in formulation and spectroscopic characterization of materials and molecules useful as catalysts, high-level density functional and molecular orbital theory calculations, and a strong foundation in concepts of physical chemistry, thermodynamics, and chemical kinetics offer new and abundant opportunities at the present day for addressing the grand challenge of controlling chemical transformations using catalysis. Here, we examine what we have learned of concepts that underpin heterogeneous catalysis but more importantly, how we learned to archive our knowledge in context of a set of best practices and standards that have emerged in the course of our learnings—ones we seek to highlight herein. Our perspective emphasizes concepts in synthesis, characterization, kinetics, and theory, because these four elements combined enable description of molecular acts that happen on surfaces and how fast they occur. In revisiting best practices in heterogeneous catalysis in these sub-fields and in authorship and peer review we aspire to augment clarity, reproducibility, and rigor in the science and practice of catalysis.

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

Data-driven particle dynamics: Structure-preserving coarse-graining for emergent behavior in non-equilibrium systems

Multiscale systems are ubiquitous in science and technology, but are notoriously challenging to simulate as short spatiotemporal scales must be appropriately linked to emergent bulk physics. When expensive high-dimensional dynamical systems are coarse-grained into low-dimensional models, the entropic loss of information leads to emergent physics which are dissipative, history-dependent, and stochastic. To machine learn coarse-grained dynamics from time-series observations of particle trajectories, we propose a framework using the metriplectic bracket formalism that preserves these properties by construction; most notably, the framework guarantees discrete notions of the first and second laws of thermodynamics, conservation of momentum, and a discrete fluctuation-dissipation balance crucial for capturing non-equilibrium statistics. We introduce the mathematical framework abstractly before specializing to a particle discretization. As labels are generally unavailable for entropic state variables, we introduce a novel self-supervised learning strategy to identify emergent structural variables. We validate the method on benchmark systems and demonstrate its utility on two challenging examples: (1) coarse-graining star polymers at challenging levels of coarse-graining while preserving non-equilibrium statistics, and (2) learning models from high-speed video of colloidal suspensions that capture coupling between local rearrangement events and emergent stochastic dynamics. We provide open-source implementations in both PyTorch and LAMMPS, enabling large-scale inference and extensibility to diverse particle-based systems.

Computational Engineering, Finance, and Science (c↗
DOE OSTI2021

Octane Index Applicability over the Pressure-Temperature Domain

Modern boosted spark-ignition (SI) engines and emerging advanced compression ignition (ACI) engines operate under conditions that deviate substantially from the conditions of conventional autoignition metrics, namely the research and motor octane numbers (RON and MON). The octane index (OI) is an emerging autoignition metric based on RON and MON which was developed to better describe fuel knock resistance over a broader range of engine conditions. Prior research at Oak Ridge National Laboratory (ORNL) identified that OI performs reasonably well under stoichiometric boosted conditions, but inconsistencies exist in the ability of OI to predict autoignition behavior under ACI strategies. Instead, the autoignition behavior under ACI operation was found to correlate more closely to fuel composition, suggesting fuel chemistry differences that are insensitive to the conditions of the RON and MON tests may become the dominant factor under these high efficiency operating conditions. This investigation builds on earlier work to study autoignition behavior over six pressure-temperature (PT) trajectories that correspond to a wide range of operating conditions, including boosted SI operation, partial fuel stratification (PFS), and spark-assisted compression ignition (SACI). A total of 12 different fuels were investigated, including the Co-Optima core fuels and five fuels that represent refinery-relevant blending streams. It was found that, for the ACI operating modes investigated here, the low temperature reactions dominate reactivity, similar to boosted SI operating conditions because their PT trajectories lay close to the RON trajectory. Additionally, the OI metric was found to adequately predict autoignition resistance over the PT domain, for the ACI conditions investigated here, and for fuels from different chemical families. This finding is in contrast with the prior study using a different type of ACI operation with different thermodynamic conditions, specifically a significantly higher temperature at the start of compression, illustrating that fuel response depends highly on the ACI strategy being used.

02 PETROLEUM↗
DOE OSTI2022

Performance evaluation of a wrapped around condenser for heat pump water heater applications

In this report, accurate performance evaluation of a wrapped-around condenser for heat pump water heaters (WHPWH) is critical since the COP of the system depends heavily on thermal stratification and condenser design. An analytical, quasi-steady state heat balance method has been developed to determine the optimal spacing between adjacent channels, tube diameter, tube shape, and total refrigerant charge amount for the condenser. The heat transfer rate is compared among three approaches: the refrigerant thermodynamic model, the condenser-wrap fin model, and the analytical natural convection model on the inside of the tank for the three regions based on the refrigerant phase (e.g., super-heated, saturated, or subcooled). The heat transfer rate was predicted by a combination of slug-plug and annular flow for a D-shaped helical tube. Experimental data were used as boundary conditions for validating the model. The rate of heat transfer based on tube shape and tank wall temperatures was compared by CFD analysis. Parametric analysis indicates a tradeoff between refrigerant mass, pressure drop, tube diameter, and tube length to maximize the heat transfer rate. The model suggests the saturated region length can be extended by 400%, and the condenser pressure drop can be reduced by 23% with an optimal spacing pattern.

42 ENGINEERING↗
DOE OSTI2020

H2OI95: A Stand-Alone Fortran Code for Evaluating the IAPWS-95 Equation-of-State Model for Water (Rev. 1)

H2OI95 is a stand-alone Fortran code for evaluating the IAPWS-95 equation-of-state model for the thermodynamic properties of water. It also evaluates the corresponding thermochemical properties of water consistent with the CODATA recommendations. The IAPWS-95 model is based on a model equation for the dimensionless Helmholtz energy for which the primary variables are the inverse reduced temperature (τ = T cr /T ) and reduced density (δ = ρ/ρ cr ). Here T is the absolute temperature (K), ρ is density (kg/m³), and the subscript “ cr ” refers to the critical point of water (647.096 K and 22.064 MPa pressure in this model, for which ρ cr is 322 kg/m³). The code solves four basic types of problems, distinguished by the specified inputs: 1. Temperature (K) and density (ρ kg/m³) or reduced density (δ) 2. Temperature (K) and pressure (MPa). 3. Temperature (K) on the saturation (liquid-vapor equilibrium) curve 4. Pressure (MPa) on the saturation curve. Each type of problem is run using a corresponding input (text) file. All but the first type of problem require iteration. For example, to solve for desired temperature and pressure, the reduced density must be adjusted to give the desired pressure. Iteration is accomplished using the Newton-Raphson method, though the secant method is also used in solving the fourth type of problem. For the last three types of problems, H2OI95 has been used to conduct numerical studies of convergence and the problem of multiple numerical solutions, some of which are not physically valid. Obtaining valid results depends mainly on appropriate choice of starting values for the density. The default values used in H2OI95 appear to consistently lead to generally desired results. With modification (not addressed here), H2OI95 can be used to support SUPCRT92 and similar codes that compute chemical thermodynamic properties of species and reactions over a wide range of temperature and pressure (273.16-1273K and 0-1000 MPa). Version 1.1 of this software calculates an extended set of water properties, including some derivable solely from IAPWS-95, such as the isobaric thermal expansivity and the isothermal compressibility, and still other properties based on supplementary models sanctioned by IAPWS, including the dynamic viscosity, thermal conductivity, surface tension, and static dielectric constant, and various other properties (notably Debye-Hückel parameters) derivable from the dielectric constant and its partial derivatives with respect to temperature and pressure. A coding error in Version 1.0 was identified while testing the implementation of the viscosity model. This error is described in Appendix C. It mainly affects calculation of “higher order” properties such as compressibility and heat capacity near the critical point.

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

Variable-density effects in incompressible non-buoyant shear-driven turbulent mixing layers

The asymmetries that arise when a mixing layer involves two miscible fluids of differing densities are investigated using incompressible (low-speed) direct numerical simulations. The simulations are performed in the temporal configuration with very large domain sizes, to allow the mixing layers to reach prolonged states of fully turbulent self-similar growth. Imposing a mean density variation breaks the mean symmetry relative to the classical single-fluid temporal mixing layer problem. Unlike prior variable-density mixing layer simulations in which the streams are composed of the same fluids with dissimilar thermodynamic properties, the density variations are presently due to compositional differences between the fluid streams, leading to different mixing dynamics. Variable-density (non-Boussinesq) effects introduce strong asymmetries in the flow statistics that can be explained by the strongest turbulence increasingly migrating to the lighter fluid side as free-stream density difference increases. Interface thickness growth rates also reduce, with some thickness definitions particularly sensitive to the corresponding changes in alignment between density and streamwise velocity profiles. Additional asymmetries in the sense of statistical distributions of densities at a given position within the mixing layer reveal that fine scales of turbulence are preferentially sustained in lighter fluid, which also is where fastest mixing occurs. These effects influence statistics involving density fluctuations, which have important implications for mixing and more complicated phenomena that are sensitive to the mixing dynamics, such as combustion.

42 ENGINEERING↗
DOE OSTI2021

Interfacial and Confinement-Mediated Organization of Gas Hydrates, Water, Organic Fluids, and Nanoparticles for the Utilization of Subsurface Energy and Geological Resources

Harnessing the subsurface geologic environments in an efficient and environmentally sustainable manner is challenged by uncertainties associated with predicting the fate of fluids and sustaining porosity and permeability in subsurface geologic environments. Some of these uncertainties arise from confined and interfacially induced structures of fluids in subsurface geologic environments. The formation of gas hydrates, phase transitions of confined fluids, assembly and deposition of heavy hydrocarbons, and the agglomeration and fate of nanoparticles in confined environments are summarized in this review. Nanoscale confinement contributes to anisotropic structures and dynamics of fluids, which is the basis for anomalous phase transition thermodynamics, reactivity, transport, and geomechanical behavior. In this review, we discuss the structures of confined fluids and deviation in observed properties from bulk fluids. The factors influencing the structures of confined fluids can be generally divided into two groups: (a) pore characteristics including pore size, pore surface chemistry, and pore geometry and (b) confined fluid/solid characteristics such as molecular structure, concentrations, charges, pore filling, and presence of additives. Scientific advancements and knowledge gaps in our understanding of the structures of confined fluids and the associated differences in observed properties compared to bulk fluids are discussed. Here, the phenomena discussed in this review are of particular relevance to our efforts in harnessing the subsurface environments for a low carbon future by increasing the utilization of geothermal energy, using CO 2 as a working fluid, and storing CO 2 in subsurface geologic environments.

42 ENGINEERING↗
DOE OSTI2023

Achieving High Strength and Creep Resistance in Inconel 740H Superalloy through Wire-Arc Additive Manufacturing and Thermodynamic-Guided Heat Treatment

Inconel 740H superalloy is commonly used in advanced ultra-supercritical power plants since it possesses excellent strength and creep resistance. This study investigates the microstructure and mechanical properties of Inconel 740H superalloy fabricated using wire-arc additive manufacturing. The as-printed microstructure consisted of columnar γ grains with the Laves phase and (Nb, Ti)C carbides as secondary phases. The anisotropy in grain structure increased from the bottom to the top regions, while the hardness was highest in the middle portion of the build. To guide the post-heat treatment design, thermodynamic and kinetic simulations were employed to predict the temperature and time. Complete recrystallization with the Laves phase dissolution occurred throughout the build after homogenization at 1200 °C for 2 h. The peak hardness was achieved after aging at 760 °C for 12 h with the M 23 C 6 carbides decorating the grain boundaries and γ’ precipitates in the grain interior. The yield strength (655 MPa) and ductility (29.5%) in the post-heat treated condition exceeded the design targets (620 MPa, 20%). Stress rupture tests at 750 °C showed that the high-temperature performance was at par with the wrought counterparts. The fracture mode after rupture was identified to be intergranular with the presence of grain boundary cavities along with grain boundary sliding.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Non-intrusive data-driven model reduction for differential–algebraic equations derived from lifting transformations

In this paper we present a non-intrusive data-driven approach for model reduction of nonlinear systems. The approach considers the particular case of nonlinear partial differential equations (PDEs) that form systems of partial differential–algebraic equations (PDAEs) when lifted to polynomial form. Such systems arise, for example, when the governing equations include Arrhenius reaction terms (e.g., in reacting flow models) and thermodynamic terms (e.g., the Helmholtz free energy terms in a phase-field solidification model). Using the known structured form of the lifted algebraic equations, the approach computes the reduced operators for the algebraic equations explicitly, using straightforward linear algebra operations on the basis matrices. The reduced operators for the differential equations are inferred from lifted snapshot data using operator inference, which solves a linear least squares regression problem. The approach is illustrated for the nonlinear model of solidification of a pure material. The lifting transformations reformulate the solidification PDEs as a system of PDAEs that have cubic structure. The operators of the lifted system for this solidification example have affine dependence on key process parameters, permitting us to learn a parametric reduced model with operator inference. Numerical experiments show the effectiveness of the resulting reduced models in capturing key aspects of the solidification dynamics.

42 ENGINEERING↗
DOE OSTI2024

Kinetically controlled metal-elastomer nanophases for environmentally resilient stretchable electronics

Nanophase mixtures, leveraging the complementary strengths of each component, are vital for composites to overcome limitations posed by single elemental materials. Among these, metal-elastomer nanophases are particularly important, holding various practical applications for stretchable electronics. However, the methodology and understanding of nanophase mixing metals and elastomers are limited due to difficulties in blending caused by thermodynamic incompatibility. Here, we present a controlled method using kinetics to mix metal atoms with elastomeric chains on the nanoscale. We find that the chain migration flux and metal deposition rate are key factors, allowing the formation of reticular nanophases when kinetically in-phase. Moreover, we observe spontaneous structural evolution, resulting in gyrified structures akin to the human brain. The hybridized gyrified reticular nanophases exhibit strain-invariant metallic electrical conductivity up to 156% areal strain, unparalleled durability in organic solvents and aqueous environments with pH 2–13, and high mechanical robustness, a prerequisite for environmentally resilient devices.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A Neural Network Approach to Predict Gibbs Free Energy of Ternary Solid Solutions

Here, we present a data-centric deep learning (DL) approach using neural networks (NNs) to predict the thermodynamics of ternary solid solutions. We explore how NNs can be trained with a dataset of Gibbs free energies computed from a CALPHAD database to predict ternary systems as a function of composition and temperature. We have chosen the energetics of the FCC solid solution phase in 226 binaries consisting of 23 elements at 11 different temperatures to demonstrate the feasibility. The number of binary data points included in the present study is 102,000. We select six ternaries to augment the binary dataset to investigate their influence on the NN prediction accuracy. We examine the sensitivity of data sampling on the prediction accuracy of NNs over selected ternary systems. It is anticipated that the current DL workflow can be further elevated by integrating advanced descriptors beyond the elemental composition and more curated training datasets to improve prediction accuracy and applicability.

42 ENGINEERING↗
DOE OSTI2022

What is in a name: Defining “high entropy” oxides

High entropy oxides are emerging as an exciting new avenue to design highly tailored functional behaviors that have no traditional counterparts. Study and application of these materials are bringing together scientists and engineers from physics, chemistry, and materials science. The diversity of each of these disciplines comes with perspectives and jargon that may be confusing to those outside of the individual fields, which can result in miscommunication of important aspects of research. In this Perspective, we provide examples of research and characterization taken from these different fields to provide a framework for classifying the differences between compositionally complex oxides, high entropy oxides, and entropy stabilized oxides, which is intended to bring a common language to this emerging area. We highlight the critical importance of understanding a material’s crystallinity, composition, and mixing length scales in determining its true definition.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Practical and thermodynamic constraints on electromicrobially accelerated CO 2 mineralization

By the end of the century, tens of gigatonnes of CO 2 will need to be removed from the atmosphere every year to maintain global temperatures. Natural weathering of ultramafic rocks and subsequent mineralization reactions can convert CO 2 into ultra-stable carbonates. Although this will draw down all excess CO 2 , it will take thousands of years. CO 2 mineralization could be accelerated by weathering ultramafic rocks with biodegradable lixiviants. We show that if these lixiviants come from cellulosic biomass, this demand could monopolize the world’s biomass supply. We demonstrate that electromicrobial production technologies (EMP) that combine renewable electricity and microbial metabolism could produce lixiviants for as little as $\$200$ to $\$400$ per tonne at solar electricity prices achievable within the decade. We demonstrate that EMP could make enough lixiviants to sequester a tonne of CO 2 for less than $\$100$. This work highlights the potential of this approach and the need for extensive R&D.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Counterdiabatic Driving with Performance Guarantees

Counterdiabatic (CD) driving has the potential to speed up adiabatic quantum state preparation by suppressing unwanted excitations. However, existing approaches either require intractable classical computations or are based on approximations that do not have performance guarantees. We propose and analyze a nonvariational, system-agnostic CD expansion method and analytically show that it converges exponentially quickly in the expansion order. In finite systems, the required resources scale inversely with the spectral gap, which we argue is asymptotically optimal. To extend our method to the thermodynamic limit and suppress errors stemming from high-frequency transitions, we leverage finite-time adiabatic protocols. In particular, we show that a time determined by the quantum speed limit is sufficient to prepare the desired ground state, without the need to optimize the adiabatic trajectory. Numerical tests of our method on the quantum Ising chain show that our method can outperform state-of-the-art variational CD approaches.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗