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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

How Frequent Will the Rarest Daily Rainfall Records of Hurricane Ida’s Remnants Be in the Future?

Abstract Gaining continued insights into the impact of global warming on the occurrence of hurricane-associated intense record downpours is essential for building climate resilient communities. This study investigates projected future changes in extreme rainfall over the Northeast United States, as represented by extreme daily amounts during Hurricane Ida in 2021. We used historical control simulations of Weather Research and Forecasting (WRF) Model generated from 40 years of weather events (1980–2014, 12 km) forced by the fifth generation European Centre for Medium-Range Weather Forecasts atmospheric reanalysis. These simulations are thermodynamically modified (2060–2100) via an imposed warming for the high-emission scenario of shared socioeconomic pathway (SSP585) from a range of general circulation models. Ground observations from the Global Historical Climatology Network (1950–2014) and WRF simulations (historical, 1980–2014, and future, 2060–2100) are integrated into a nonstationary generalized extreme value (GEV) framework to assess the frequency of Ida’s heaviest daily rain rates under the SSP585 scenario. Results show that Ida’s daily maximum rainfall recorded at different observation locations was higher than the single highest September daily maximum observed (1950–2014) for 5 out of 17 stations (∼30% of the stations). Ida-like extreme daily rain rates are projected to be, on average, more than 2 times more likely to occur at the end of the century in the simulations (with some regions as high as 5 times). This work demonstrates that integrating a high-resolution atmospheric model’s present-day and thermodynamically modified future simulations along with ground observations, within a nonstationary statistical framework, is crucial for understanding changing characteristics of extreme weather events. Significance Statement Daily scale extreme precipitation is expected to become more frequent and severe, as evidenced by observations and model simulations. While it is important to investigate how these intensifying heavy rainfall events affect current engineering standards, fewer studies have contextualized how warming impacts the most extreme rainfall from a single storm event relative to historical heavy downpours. In this study, we focused on the daily extreme rainfall associated with the extratropical transition of Hurricane Ida (2021), particularly over the northeastern United States—some of which exceeded the commonly used hydrologic design criteria for a 100-yr storm. Using a high-resolution atmospheric model simulation, we investigated how continued warming may influence the frequency of such daily rain rates. Under a high-emission scenario, these events are projected to become up to 5 times more likely at the end of the twenty-first century.

Dollan, Ishrat J↗
DOE OSTI2025

A chain stretch-based gradient-enhanced model for damage and fracture in elastomers

Similar to quasi-brittle materials, it has been recently shown that elastomers can exhibit a macroscopically diffuse damage zone that accompanies the fracture process. In this study, we introduce a stretch-based gradient-enhanced damage (GED) model that allows the fracture to localize and also captures the development of a physically diffuse damage zone. This capability contrasts with the paradigm of the phase field method for fracture, where a sharp crack is numerically approximated in a diffuse manner. Capturing fracture localization and diffuse damage in our approach is achieved by considering nonlocal effects that encompass network topology, heterogeneity, and imperfections. These considerations motivate the use of a statistical damage function dependent upon the nonlocal deformation state. From this model, fracture toughness is realized as an output. While GED models have been classically utilized for damage modeling of structural engineering materials (e.g., concrete), they face challenges when trying to capture the cascade from damage to fracture, often leading to damage zone broadening (de Borst and Verhoosel, 2016). This deficiency contributed to the popularity of the phase-field method over the GED model for elastomers and other quasi-brittle materials. Other groups have proceeded with damage-based GED formulations that prove identical to the phase-field method (Lorentz et al., 2012), but these inherit the aforementioned limitations. To address this issue in a thermodynamically consistent framework, we implement two modeling features (a nonlocal driving force bound and a simple relaxation function) specifically designed to capture the evolution of a physically meaningful damage field and the simultaneous localization of fracture, thereby overcoming a longstanding obstacle in the development of these nonlocal strain- or stretch-based approaches. Here, we discuss several numerical examples to understand the features of the approach at the limit of incompressibility, and compare them to the phase-field method as a benchmark for the macroscopic response and fracture energy predictions.

Elastomers↗
DOE OSTI2020

Tuning the quantumness of simple Bose systems: A universal phase diagram

We present a comprehensive theoretical study of the phase diagram of a system of many Bose particles interacting with a two-body central potential of the so-called Lennard-Jones form. First-principles path-integral computations are carried out, providing essentially exact numerical results on the thermodynamic properties. The theoretical model used here provides a realistic and remarkably general framework for describing simple Bose systems ranging from crystals to normal fluids to superfluids and gases. The interplay between particle interactions on the one hand and quantum indistinguishability and delocalization on the other hand is characterized by a single quantumness parameter, which can be tuned to engineer and explore different regimes. Taking advantage of the rare combination of the versatility of the many-body Hamiltonian and the possibility for exact computations, we systematically investigate the phases of the systems as a function of pressure (P) and temperature (T), as well as the quantumness parameter. Here, we show how the topology of the phase diagram evolves from the known case of 4 He, as the system is made more (and less) quantum, and compare our predictions with available results from mean-field theory. Possible realization and observation of the phases and physical regimes predicted here are discussed in various experimental systems, including hypothetical muonic matter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Recent Advances of PyROS: A Pyomo Solver for Nonconvex Two-Stage Robust Optimization in Process Systems Engineering

This poster highlights uncertainty and technical risk reduction capabilities in CCSI2, with a focus on robust optimization. It presents recent advances of the two-stage robust optimization (RO) solver PyROS and applications to advanced energy systems optimization. To demonstrate the computational performance and reliability of PyROS, a benchmarking study on a library of over 8,500 small-scale RO problems is presented. Further, PyROS is used to obtain robust system designs of a MEA-based CO2 absorber under uncertainty in the thermodynamic property models for a variety of CO2 capture rate threshold requirements. Overall, the results demonstrate that the PyROS solver, including recent extensions to multi-stage RO settings, provides a reliable avenue to optimize the design and operation of advanced energy systems subject to various sources of parametric uncertainty.

Sherman, Jason↗
DOE OSTI2024

Recent Advances of PyROS: A Pyomo Solver for Nonconvex Two-Stage Robust Optimization in Process Systems Engineering

This poster highlights uncertainty and technical risk reduction capabilities in CCSI2, with a focus on robust optimization. It presents recent advances of the two-stage robust optimization (RO) solver PyROS and applications to advanced energy systems optimization. To demonstrate the computational performance and reliability of PyROS, a benchmarking study on a library of over 8,500 small-scale RO problems is presented. Further, PyROS is used to obtain robust system designs of a MEA-based CO2 absorber under uncertainty in the thermodynamic property models for a variety of CO2 capture rate threshold requirements. Overall, the results demonstrate that the PyROS solver, including recent extensions to multi-stage RO settings, provides a reliable avenue to optimize the design and operation of advanced energy systems subject to various sources of parametric uncertainty.

Sherman, Jason↗
DOE OSTI2025

Effect of MgSO 4 addition on alkali sulphates induced hot corrosion of a β-NiAl coating

The effect of MgSO 4 addition on alkali sulphates induced hot corrosion of a β-NiAl coating was studied by performing a series of tests with deposits of Na 2 SO 4 , Na 2 SO 4 -20 %K 2 SO 4 and Na 2 SO 4 -12 %K 2 SO 4 -35 %MgSO 4 (mol.% for all) at 700 °C in air with addition of SO 2 . For elucidating the observed corrosion phenomena, various salt-oxide/NiAl powder mixtures were exposed to the same hot corrosion test conditions. Two pseudo-binary phase diagrams of the Na, K, Mg-sulphate systems were calculated using an in-house developed thermodynamic database. Addition of MgSO 4 to Na 2 SO 4 -20 %K 2 SO 4 deposit salt inhibited the hot corrosion attack during 24 h exposure, which is related to consumption of aggressive K 2 SO 4 by the formation of solid K 2 Mg 2 (SO 4 ) 3 . However, extending the exposure with Na 2 SO 4 -12 %K 2 SO 4 -35 %MgSO 4 deposit to 100 h resulted in severe corrosion attack. The latter observation is explained by a reaction between K 2 Mg 2 (SO 4 ) 3 and NiAl causing formation of a Mg, Ni and Al containing spinel accompanied by release of the aggressive K 2 SO 4 into the liquid phase.

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

Modeling Nanosecond-Pulsed Spark Discharge and Flame Kernel Evolution

Dilute combustion, either using exhaust gas recirculation or with excess air, is considered a promising strategy to improve the thermal efficiency of internal combustion engines. However, the dilute air-fuel mixture, especially under intensified turbulence and high-pressure conditions, poses significant challenges for ignitability and combustion stability, which may limit the attainable efficiency benefits. In-depth knowledge of the flame kernel evolution to stabilize ignition and combustion in a challenging environment is crucial for effective engine development and optimization. To date, a comprehensive understanding of ignition processes that result in the development of fully predictive ignition models usable by the automotive industry does not yet exist. Spark-ignition consists of a wide range of physics that includes electrical discharge, plasma evolution, joule-heating of gas, and flame kernel initiation and growth into a self-sustainable flame. In this study, an advanced approach is proposed to model spark-ignition energy deposition and flame kernel growth. To decouple the flame kernel growth from the electrical discharge, a nanosecond-pulsed high-voltage discharge is used to trigger spark-ignition in an optically accessible small ignition test vessel with a quiescent mixture of air and methane. Initial conditions for the flame kernel, including its thermodynamic state and species composition, are derived from a plasma-chemical equilibrium calculation. The geometric shape and dimension of the kernel are characterized using a multi-dimensional thermal plasma solver. Here, the proposed modeling approach is evaluated using a high-fidelity computational fluid dynamics procedure to compare the simulated flame kernel evolution against flame boundaries from companion Schlieren images.

42 ENGINEERING↗
DOE OSTI2023

On the Thermodynamic Condition for Adsorption Azeotropes

Adsorption azeotropy is a common phenomenon in mixed-gas adsorption equilibria. Since the first literature report of binary adsorption azeotropes in 1933, numerous studies investigated the thermodynamic conditions for adsorption azeotropes but failed to reach definitive conclusions. Based on the generalized Langmuir isotherm model for multicomponent adsorption equilibria which takes into account the vacant site as part of the adsorbed phase, this study presents the thermodynamic condition for adsorption azeotropes as derived from the generalized Langmuir isotherm to be the equality of the ratios of adsorbed phase activity coefficient γi and adsorption equilibrium constant $K$$^{o}_{i}$ for the two adsorbates in the binary 1–2 adsorption system, i.e., γ 1 /$K$$^{o}_{1}$ = γ 2 /$K$$^{o}_{2}$. This adsorption azeotropic condition is analogous to the vapor–liquid equilibrium azeotropic condition, i.e., the equality of the products of liquid phase activity coefficient and saturation vapor pressure Pisat for the two components, i.e., γ 1 $P$$^{sat}_{1}$ = γ 2 $P$$^{sat}_{2}$. We validated the thermodynamic condition for adsorption azeotropes with 14 azeotrope-forming adsorption systems in this study. As a result, we investigated the effects of the pressure, temperature, and adsorbed phase nonideality on azeotrope formation.

42 ENGINEERING↗
DOE OSTI2021

The critical role of synthesis conditions on small polaron carrier concentrations in hematite—A first-principles study

Achieving highly efficient energy conversion with transition metal oxides necessitates overcoming conductivity limitations due to the formation of small polarons. Detailed understanding of the interplay among intrinsic defects, dopants, and electron polarons can help devise strategies for achieving higher carrier concentrations, therefore improving carrier conductivity. This work employs first-principles calculations to reliably predict electron polaron concentrations in a prominent polaronic oxide, hematite (Fe 2 O 3 ), by resolving interactions between charged defects and electron polarons and keeping charge neutrality condition among all charged species. This work addresses that both V O and Fe i can be primary donors in undoped hematite depending on the synthesis conditions, such as synthesis temperature and oxygen partial pressure, despite the fact that V O owns an extremely high ionization energy compared to k B T. Furthermore, from calculations of a plethora of n-type dopants (group IV and V elements), we find that Ti, Ge, Sb, and Nb are able to raise electron polaron concentrations in hematite significantly without considering dopant clustering. However, the magnitude of electron polaron concentration increase would be smaller if the dopant has a high tendency of clustering, such as Ti. Here we reveal the critical role of synthesis conditions on tuning electron polaron concentrations of both undoped and doped hematite. Our theoretical analysis provides important insights and general design principles for engineering more conductive polaronic oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Predicting the viability of beta-lactamase: How folding and binding free energies correlate with beta-lactamase fitness

One of the long-standing holy grails of molecular evolution has been the ability to predict an organism’s fitness directly from its genotype. With such predictive abilities in hand, researchers would be able to more accurately forecast how organisms will evolve and how proteins with novel functions could be engineered, leading to revolutionary advances in medicine and biotechnology. In this work, we assemble the largest reported set of experimental TEM-1 β-lactamase folding free energies and use this data in conjunction with previously acquired fitness data and computational free energy predictions to determine how much of the fitness of β-lactamase can be directly predicted by thermodynamic folding and binding free energies. We focus upon β-lactamase because of its long history as a model enzyme and its central role in antibiotic resistance. Based upon a set of 21 β-lactamase single and double mutants expressly designed to influence protein folding, we first demonstrate that modeling software designed to compute folding free energies such as FoldX and PyRosetta can meaningfully, although not perfectly, predict the experimental folding free energies of single mutants. Interestingly, while these techniques also yield sensible double mutant free energies, we show that they do so for the wrong physical reasons. We then go on to assess how well both experimental and computational folding free energies explain single mutant fitness. We find that folding free energies account for, at most, 24% of the variance in β-lactamase fitness values according to linear models and, somewhat surprisingly, complementing folding free energies with computationally-predicted binding free energies of residues near the active site only increases the folding-only figure by a few percent. This strongly suggests that the majority of β-lactamase’s fitness is controlled by factors other than free energies. Overall, our results shed a bright light on to what extent the community is justified in using thermodynamic measures to infer protein fitness as well as how applicable modern computational techniques for predicting free energies will be to the large data sets of multiply-mutated proteins forthcoming.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Harnessing metastability for grain size control in multiprincipal element alloys during additive manufacturing

Abstract Controlling microstructure in fusion-based metal additive manufacturing (AM) remains a significant challenge due to the many parameters that directly impact solidification condition. Multiprincipal element alloys (MPEAs), also known as high entropy alloys, offer a vast compositional space to design for microstructural engineering due to their chemical complexity and exceptional properties. Here, we use the FeMnCoCr system as a model platform for exploring alloy design in MPEAs for AM. By exploiting the decreasing stability of the face-centered cubic phase with increasing Mn content, we achieve notable grain refinement and breakdown of epitaxial columnar grain growth. We employ a multifaceted approach encompassing thermodynamic modeling, operando synchrotron X-ray diffraction, multiscale microstructural characterization, and mechanical testing to gain insight into the solidification physics and its ramifications on the resulting microstructure of FeMnCoCr MPEAs. This work aims toward tailoring desirable grain sizes and morphology through targeted manipulation of phase stability, thereby advancing microstructure control in AM applications.

Wakai, Akane↗
DOE OSTI2022

Acetyl-CoA Synthesis Through a Bicyclic Carbon-Fixing Pathway in Gas-Fermenting Bacteria

Gas-fermenting acetogens can upgrade one-carbon (C1) compounds (such as CO2 and CO) to the two-carbon (C2) metabolite acetyl coenzyme A (CoA) and convert sugar feedstocks to acetyl-CoA with minimal CO2 emissions. Fulfilling the biosynthetic potential of these microbes requires overcoming challenges in pathway engineering. Here we design a synthetic acetyl-CoA bi-cycle-in addition to the natural carbon-fixing pathways-for C2 metabolite synthesis. This pathway produces an acetyl-CoA by fixation of two CO2 equivalents via three functional modules acting in sequence: carbon fixation, gluconeogenesis and non-oxidative glycolysis. The pathway was examined by in silico thermodynamic and kinetic analyses. The prototypic pathway was implemented in a syngas-fermenting organism, Clostridium ljungdahlii DSM 13528, by expressing a heterologous phosphoketolase that can work with other native enzymes in the host acetogen. The carbon conversion pathway is possible under various growth conditions and is independent of the Wood-Ljungdahl pathway for the valorization of H2 and CO2. This study reports the improvement of carbon conversion using a reductive acetyl-CoA bi-cycle and the potential impact of redox homoeostasis in the acetogenic host for industrial applications of gas fermentation.

Acetyl-CoA↗
DOE OSTI2020

A Numerical Study on the Impact of Cathode Catalyst Layer Loading on the Open Circuit Voltage

The open-circuit voltage (OCV) is the measured potential when no current is drawn from a cell. Additionally, theoretically this potential is given by the equilibrium potential, which is determined from thermodynamics based on the operating conditions. In practice, the actual potential at zero current is measured to be approximately 300 mV lower than the theoretical value.

42 ENGINEERING↗
DOE OSTI2024

An MPMD approach coupling electromagnetic continuum mechanics approximations in ALEGRA

In this work, two complementary approximations for describing aspects of continuum electromagnetics in moving media are discussed: electroquasistatic and magnetoquasistatic. Each has been implemented in the finite element shock code ALEGRA for modeling dynamic electromechanical phenomena on typical engineering time scales, with fully integrated circuit coupling. The approximations can be obtained by consistent asymptotic balancing of Maxwell’s equations relative to timescales associated with magnetic diffusion, charge relaxation, and electromagnetic wave propagation. In ALEGRA, the electroquasistatic approximation is used for ferroelectric (FE) modeling, while the magnetoquasistatic approximation is used for magnetohydrodynamic (MHD) modeling. In this paper we introduce for the first time a detailed derivation of a useful quasi-steady “low-R m ” variant of the MHD approximation applicable for cases, such as with detonators, where the thermodynamic pressure arising from Joule heating dominates over magnetic forces. An additional purpose of this paper is to present a coupling mode using Multiple Program-Multiple Data (MPMD) message passing communication that allows the user to run 3D FE problems together with 2D and/or 3D MHD problems with the respective simulation domains coupled through a common circuit equation. The MPMD coupling capability is used here to model the dynamic coupling of a notional ferroelectric generator with an RP-87 exploding bridgewire detonator. The simulated bridgewire heats up and bursts under current generated by simulated depoling of the ferroelectric generator, as a demonstration of the MPMD capability.

42 ENGINEERING↗
DOE OSTI2025

Stress-dependent χ phase transformation in a Ni-based superalloy

The ongoing push to elevate operating temperatures in aerospace gas turbine engines – driven by goals of enhanced fuel efficiency and reduced CO 2 emissions – mandates advancements in the creep resistance of Ni- and Co-based superalloys, which are integral for critical engine components. This study elucidates the role of stress assisted localized phase transformations in the creep properties of these alloys. By leveraging chemo-mechanical coupling, self-healing γ′ precipitates are designed to immobilize planar defects, thereby increasing creep resistance. Employing advanced characterization techniques such as high-resolution Scanning Transmission Electron Microscopy (HR-STEM), in conjunction with atomistic simulations and thermodynamic calculations, novel deformation pathways facilitated by χ local phase transformation (LPT) strengthening have been uncovered; notably, the formation of χ nano-laths through microtwinning and superlattice intrinsic stacking fault (SISF) shearing. This study highlights critical insights into the compositional boundaries necessary for optimizing LPT strengthening while avoiding deleterious bulk formation of η/χ phases. These advancements will guide the design of new alloys maximizing high-temperature creep strength for advanced aerospace applications.

Egan, Ashton J. [Friedrich-Alexander University Er↗
DOE OSTI2020

Performance Characterization of a Natural Gas–Air Rotating Detonation Engine

An experimental study of a rotating detonation engine (RDE) operating with natural gas and air at elevated chamber pressures and air preheat temperatures was conducted to quantify its performance at conditions representative of land-based power generation gas turbine engines. Here, the thrust produced by the combustor was measured to characterize its work output potential. High-frequency pressure transducers and broadband chemiluminescence measurements of the flame provided information about the wave structure and dynamics. Analysis of common performance metrics demonstrated the necessity of normalizing any RDE performance parameter by the driving system potential, typically the reactant manifold pressure. Application of a thermodynamic performance model to a generic RDE identified the area ratio between the RDE exhaust and injection throats as the primary parameter affecting delivered pressure gain. The model was further applied to draw comparison with experimental measurements of net pressure gain for identical flow conditions. Only one of the two tested injector configurations followed the predicted trends, suggesting that performance of the second was governed by physical processes other than the reactant thermodynamics. Although an absolute pressure gain was not demonstrated, it is promising that the natural gas–air RDE delivered up to 90% of the theoretical performance.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2022

Thermodynamic Constraints on Electromicrobial Protein Production

Global consumption of protein is projected to double by the middle of the 21st century. However, protein production is one of the most energy intensive and environmentally damaging parts of the food supply system today. Electromicrobial production technologies that combine renewable electricity and CO 2 -fixing microbial metabolism could dramatically increase the energy efficiency of commodity chemical production. Here we present a molecular-scale model that sets an upper limit on the performance of any organism performing electromicrobial protein production. We show that engineered microbes that fix CO 2 and N 2 using reducing equivalents produced by H 2 -oxidation or extracellular electron uptake could produce amino acids with energy inputs as low as 64 MJ kg −1 , approximately one order of magnitude higher than any previous estimate of the efficiency of electromicrobial protein production. This work provides a roadmap for development of engineered microbes that could significantly expand access to proteins produced with a low environmental footprint.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

An experimental and numerical investigation of HD diesel engine DOC efficiency in oxidizing NO to NO 2

Reducing pollutant emissions from heavy-duty (HD) diesel engines is critical due to their significant environmental impact, particularly concerning NOx emissions. Understanding and optimizing modern diesel oxidation catalyst (DOC) and selective catalytic reduction (SCR) performance is essential for improving exhaust aftertreatment (EAT) system efficiency to meet stringent emissions regulations. The oxidation of nitric oxide (NO) to nitrogen dioxide (NO 2 ) in DOC plays a key role in improving SCR efficiency in reducing NO x . This study investigates the DOC performance in oxidizing NO to NO2 and its impact on the SCR efficiency of a 2021 MY Navistar E39 HD diesel engine. The influence of engine speed, load, exhaust gas temperature, and composition on DOC efficiency is experimentally investigated. The relationships between DOC inlet temperature, oxygen availability and NO 2 /NO x ratio at the DOC inlet are examined to better understand their effects on the overall DOC efficiency. The results indicate that DOC NO oxidation efficiency is highly dependent on exhaust temperature, with optimal oxidation occurring within a specific temperature range (275-350°C). Below this threshold, the chemical reactions are kinetically limited, while at higher temperatures, thermodynamic constraints reduce the efficiency of DOC in oxidizing NO to NO 2 . The experimental data further reveal that the NO 2 /NO x ratio peaks at medium loads before declining at higher loads due to reduced residence time and mass transfer effects. Additionally, the SCR NO x conversion efficiency is significantly influenced by the NO 2 /NO x ratio, achieving peak performance when the NO 2 /NO x ratio approaches 0.5. A DOC chemistry model was developed and validated against the experimental data to predict DOC oxidation behavior under various operating conditions. The findings of this study provide insights into the interdependencies between DOC and SCR performance, contributing to the optimization of SCR systems for optimized NOx reduction.

33 ADVANCED PROPULSION SYSTEMS↗