Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal equilibrium”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Dynamic reaction-induced phase separation in tunable, adaptive covalent networks

A series of catalyst-free, room temperature dynamic bonds derived from a reversible thia-Michael reaction are utilized to access mechanically robust dynamic covalent network films. The equilibrium of the thiol addition to benzalcyanoacetate-based Michael-acceptors can be directly tuned by controlling the electron-donating/withdrawing nature of the Michael-acceptor. By modulating the composition of different Michael-acceptors in a dynamic covalent network, a wide range of mechanical properties and thermal responses can be realized. Additionally, the reported systems phase-separate in a process, coined dynamic reaction-induced phase separation (DRIPS), that yields reconfigurable phase morphologies and reprogrammable shape-memory behaviour as highlighted by the heat-induced folding of a predetermined structure.

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

Analytical benchmark solution for 1-D neutron transport coupled with thermal conduction and material expansion

In this work we present fully analytical solutions for a class of finite, homogeneous, 1-D slab benchmark problems with nonlinear temperature feedback effects. The proposed class of benchmarks include multiplicative 1-D neutron transport (limited to quasistatic S{sub 2} with μ = ±1) coupled with thermal conduction, convection, Doppler broadening, and expansion effects along the length of the slab. This class of benchmark models, along with the corresponding analytical solutions, are valuable for validating multiphysics analysis frameworks that support coupled neutronics/thermal/structural calculations. Analytical solutions for the benchmarks are obtained by introducing an ansatz that the equilibrium flux and temperature distributions in the slab have the same shape. Specific values for the total microscopic cross section, σ{sub t,0}, and conductive heat transfer coefficient, h, that satisfy the assumed ansatz are then determined. A discussion of the procedure for generating benchmark models and analytical solutions is provided, along with numerical results for an example set of model parameters and thoughts on practical applications of the benchmarks for multiphysics code validation. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Space-time generalization of mutual information

The mutual information characterizes correlations between spatially separated regions of a system. Yet, in experiments we often measure dynamical correlations, which involve probing operators that are also separated in time. Here, we introduce a space-time generalization of mutual information which, by construction, satisfies several natural properties of the mutual information and at the same time characterizes correlations across subsystems that are separated in time. In particular, this quantity, that we call the space-time mutual information, bounds all dynamical correlations. We construct this quantity based on the idea of the quantum hypothesis testing. As a by-product, our definition provides a transparent interpretation in terms of an experimentally accessible setup. We draw connections with other notions in quantum information theory, such as quantum channel discrimination. Finally, we study the behavior of the space-time mutual information in several settings and contrast its long-time behavior in many-body localizing and thermalizing systems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

The effect of equilibration time on Al uptake in C-S-H

Calcium silicate hydrate (C-S-H) is the main hydration product of Portland cements. The partial replacement of Portland cements by supplementary cementitious materials can result in C-S-H with lower Ca/Si ratios and more aluminum and alkalis. The effect of equilibration time on Al uptake in C-S-H was investigated using equilibration times from 7 days up to 3 years. Lower Al concentrations were measured in the solution after longer equilibration times. In addition, a higher uptake of Al in C-S-H was observed based on the decrease in the content of secondary phases. Little further decrease in Al concentrations was observed after 2 years and longer at low Ca/Si. At high Ca/Si no significant change in solution concentrations was observed after more than 3 month, while the destabilization of secondary phases continued up to 1 year, indicating that a (meta)stable equilibrium was reached faster at higher Ca/Si ratios.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The Application of Differential Scanning Calorimetry to Investigate Precipitation Behavior in Nickel-Base Superalloys Under Continuous Cooling and Heating Conditions

In this study, a suite of experimental tools and fast-acting, numerical-simulation techniques was used to quantify the precipitation behavior of three nickel-base superalloys: IN-100, LSHR, and 718. Experimental methods comprised differential scanning calorimetry (DSC) to establish the specific heat as a function of temperature and selected direct-resistance heating trials (using a Gleeble ® machine) to obtain samples for microstructural analysis. For the DSC experiments, each alloy was cooled at a prescribed constant rate (between 5 and 20 K/min) after an initial soak/equilibration in the high-temperature, single-phase (supersolvus) temperature regime. On-heating DSC trials beginning at ambient temperature were also performed on alloy 718 in three different starting conditions: super- δ -solvus solution treated and water quenched (denoted as ST), solution treated and aged (STA), and solution treated and overaged (STOA). DSC results, revealing the thermal signatures associated with the kinetics of precipitation of γ ' (IN-100, LSHR) or γ ' and γ " (718), were interpreted using a previously-developed fast-acting routine that treats concurrent nucleation, growth, coarsening, and dissolution. For these simulations, special attention was paid to various thermo-kinetic input parameters including equilibrium solvus-approach curves, bulk free energies of transformation, matrix-precipitate interface energies, and effective diffusivities. For the γ - γ ' superalloys (IN-100 and LSHR), estimates of precipitate volume fraction as a function of temperature from the specific-heat data revealed semi-quantitative agreement with simulation predictions. For the γ - γ '- γ " superalloy (718), simulation predictions of precipitate volume fractions were converted to specific heat as a function of temperature and showed semi-quantitative agreement with the direct measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Air separation via a two-step solar thermochemical cycle based on (Ba,La) x Sr 1-x FeO 3-δ : Thermodynamic analysis

A two-step solar thermochemical cycle was considered for air separation to produce N2 based on (Ba,La) x Sr 1-x FeO 3-δ perovskite reduction/oxidation (redox) reactions for A-site fractions of 0 ≤ x ≤ 0.2. The cycle steps encompassed (1) thermal reduction and O 2 release via concentrated solar input and (2) re-oxidation with air to uptake O 2 and produce high-purity N 2 . Thermogravimetry at temperatures between 400 and 1100 °C in atmospheres of 0.005 to 90% O 2 /Ar at 1 bar was performed to measure equilibrium nonstoichiometries. The compound energy formalism was applied to model redox thermodynamics for both Ba 2+ and La 3+ substitution. Non-linear regression was used to determine the empirical parameters based on the thermogravimetric measurements. The model was used to define partial molar reaction enthalpies and entropies and predicted equilibrium oxygen nonstoichiometry as functions of oxide stoichiometry, site fraction, temperature, and O 2 partial pressure. The thermodynamic analysis showed the materials are appealing for air separation at temperatures below 800 °C.

42 ENGINEERING↗
DOE OSTI2021

Impact of collisionality on turbulence in the edge of tokamak plasma using 3D global simulations

Collisionality is one of the key parameters in determining turbulent transport in the plasma edge, regulating phenomena such as 'shoulder formation', separation of scale lengths in the scrape-off layer (SOL), turbulence damping and zonal flow dynamics. Understanding its role is therefore of primary importance for future reactors like ITER. Getting reliable predictions and a better characterization of plasma flow properties when varying collisionality remains, however, a critical challenge for the simulations. This paper focuses on the impact of varying collisionality in a non-isothermal three-dimensional fluid model of the plasma edge. A high field side limited configuration encompassing open and closed magnetic field lines with parameters typical of a medium-sized tokamak is considered. The present model can consistently account for the variations of collisionality and its impact on both the parallel resistivity η ∥ and the ion and electron parallel thermal conductivities χ ∥e,i . Details on mean flow and turbulence properties are given. Changing collisionality leads to significant changes in the flow properties both on the mean and fluctuating quantities. In particular, lowering collisionality decreases the size of coherent structures, the fluctuation levels of turbulence, and steepens the density and temperature equilibrium profiles around the separatrix leading to a global reduction of the turbulent transport. The SOL width is observed to increase with collisionality, eventually resulting in the disappearance of the scale lengths separation between near and far SOL, consistently with previous experimental observations. At low collisionality, where the presence of narrow feature is well-established, a contribution of heat conduction increases up to compete with heat convection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Designing a validation experiment for radio frequency condensation

Abstract Theoretical studies have suggested that nonlinear effects can lead to ‘radio frequency (RF) condensation’, where an initially broad current profile can coalesce in islands when they reach sufficient width. In suitable conditions, RF condensation can ‘self-focus’ the driven current to the center of an island, improving stabilization efficiency and reducing control complexity. In unsuitable conditions, the effect can prematurely deplete the RF energy before it reaches the island center, impairing stabilization. It is predicted that the RF condensation effect can significantly impact reactor-scale tokamaks. This paper presents a set of simulations investigating the conditions under which RF condensation might be encountered in present-day tokamaks. For concreteness, the calculations use equilibrium reconstructions for two shots from DIII-D and AUG. The Current Condensation Amid Magnetic Islands (OCCAMI) simulation code has been used for this investigation. The code takes as its input a numerically specified axisymmetric EFIT equilibrium solution, and it perturbatively constructs a 3D field with an island embedded at the appropriate rational surface. In the OCCAMI code, the GENRAY code is used for ray tracing and for calculating the power deposition along a ray trajectory, and GENRAY is coupled self-consistently to a solution of the thermal diffusion equation in the island. The simulation results described in the paper illuminate the conditions required for experimental validation of the theory of RF condensation. The simulations also provide an explanation of why the effect was not noticed in experiments prior to the publication of theoretical papers on the subject.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Phonon-assisted formation of an itinerant electronic density wave

Abstract Electronic instabilities drive ordering transitions in condensed matter. Despite many advances in the microscopic understanding of the ordered states, a more nuanced and profound question often remains unanswered: how do the collective excitations influence the electronic order formation? Here, we experimentally show that a phonon affects the spin density wave (SDW) formation after an SDW-quench by femtosecond laser pulses. In a thin film, the temperature-dependent SDW period is quantized, allowing us to track the out-of-equilibrium formation path of the SDW precisely. By exploiting its persistent coupling to the lattice, we probe the SDW through the transient lattice distortion, measured by femtosecond X-ray diffraction. We find that within 500 femtoseconds after a complete quench, the SDW forms with the low-temperature period, directly bypassing a thermal state with the high-temperature period. We argue that a wavevector-matched phonon launched by the quench changes the formation path of the SDW through the dynamic pinning of the order parameter.

42 ENGINEERING↗
DOE OSTI2020

A mathematical model of solid-state dewetting of barium thin films on W (112)

Solid state dewetting occurs when a thin solid film is heated. The temperature of dewetting depends on the thickness of the film; dewetting can be observed in the range of 1/3 to 1/2 of the bulk melting temperature. While remaining solid, the film behaves in a manner similar to liquids dewetting and agglomerating to forming islands or droplets. One of the possible mechanisms is the conversion of a metastable thin film geometry into a more stable form. Heating the metastable film gives the film atoms higher mobility, and the films retract, dewetting the surface. This atomic motion can be restricted due to surface anisotropy. Herein we present in situ emission microscopy observations of barium thin films deposited onto W (112) by thermal evaporation. From the modeling viewpoint, the evolution of the film in this system could be divided in four stages: (i) the nucleation and growth of the thin film as a simply connected region; (ii) formation of droplets/islands/hillocks; (iii) nucleation of holes; (iv) evolution of the components of the disconnected film to their equilibrium state. We present a continuum model that is qualitatively consistent with the evolution of the film observed at the initial stage of the experiment and discuss the later stages of the evolution of surface structures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Eigenmode analysis of the sheared-flow Z-pinch

Experiments have demonstrated that a Z-pinch can persist for thousands of times longer than the growth time of global magnetohydrodynamic (MHD) instabilities such as the m=0 sausage and m=1 kink modes. These modes have growth times on the order of ta=a/vi, where vi is the ion thermal speed and a is the pinch radius. Axial flows with duz/dr ≲ vi/a have been measured during the stable period, and the commonly accepted theory is that this amount of shear is sufficient to stabilize these modes as predicted by numerical studies using the ideal MHD equations. However, these studies only consider specific equilibrium profiles that typically have a modest magnitude for the logarithmic pressure gradient, qP≡d ln P/d ln r, and may not represent experimental conditions. Linear stability of the sheared-flow Z-pinch is studied here via a direct eigen-decomposition of the matrix operator obtained from the linear ideal MHD equations. Several equilibrium profiles with a large variation of qP are examined. Considering a practical range of k, 1/3 ≲ ka ≲ 10, it is shown that the shear required to stabilize m=0 modes can be expressed as duz/dr≥Cγ0/(ka)α. Here, γ0=γ0(ka) is the profile-specific growth rate in the absence of shear, which scales approximately with |qP|. Both C and α are profile-specific constants, but C is order unity and α≈1. It is further demonstrated that even a large value of shear, duz/dr=3vi/a, is not sufficient to provide linear stabilization of the m=1 kink mode for all profiles considered. This result is in contrast to the currently accepted theory predicting stabilization at much lower shear, duz/dr=0.1vi/a, and suggests that the experimentally observed stability cannot be explained within the linear ideal-MHD model.

Angus, J. R. (ORCID:0000000314740002)↗
DOE OSTI2022

HEAT TRANSFER ASPECTS OF TWO PISTON SPLAT QUENCHING

This thesis investigates the heat transfer aspects of the two-piston splat quench system by modifying surface roughness and by using copper and stainless steel platens. In additive manufacturing, specifically selective laser melting (SLM), metal is subjected to a very rapid heating and cooling that produces cooling rates in the range from 105-106 °C/s. This rapid solidification produces microstructures that deviate far from equilibrium. Splat quenching is an experimental method that produces rapidly solidified structures under cooling rates that are comparable to those found in additive manufacturing. The splat quenching approach allows for more efficient testing of the rapid solidification behavior of a wide variety of alloys. In this thesis, two platen characteristics, surface roughness and material thermal conductivity, were examined to determine their role in controlling the heat transfer and the subsequent cooling rate during splat quenching of austenitic stainless steel 316L. To that objective, 5 surface finishes were produced using metallographic techniques, from 7046-203 Å Rms. Copper and 316L stainless steel platen material were used to reach different levels of thermal conductance. The combination of these 2 variables created 10 scenarios, which were examined using, backscatter imaging, secondary imaging, and electron back scatter diffraction. Cooling rates were estimated using cell sizing obtained through electrolytic etching and secondary imaging. This thesis determined that both surface modification and changes in material type had a statistically significant impact on the cell sizes observed. Where material type far outweighed surface roughness in its effect on cell size.

Roze, Jonathan↗
DOE OSTI2023

A statistical mechanics framework for polymer chain scission, based on the concepts of distorted bond potential and asymptotic matching

To design increasingly tough, resilient, and fatigue-resistant elastomers and hydrogels, the relationship between controllable network parameters at the molecular level (bond type, non-uniform chain length, entanglement density, etc.) to macroscopic quantities that govern damage and failure must be established. Many of the most successful constitutive models for elastomers have been rooted in statistical mechanical treatments of polymer chains. Typically, such constitutive models have used variants of the freely jointed chain model with rigid links. However, since the free energy state of a polymer chain is dominated by enthalpic bond distortion effects as the chain approaches its rupture point, bond extensibility ought to be accounted for if the model is intended to capture chain rupture. To that end, a new bond potential is supplemented to the freely jointed chain model (as derived in the u FJC framework of Buche and Silberstein (2021) and Buche et al. (2022)), which we have extended to yield a tractable, closed-form model of single chain behavior that should be amenable to continuum-level constitutive model development. Inspired by the asymptotically matched u FJC model response in both the low/intermediate chain force and high chain force regimes, a simple, quasi-polynomial bond potential energy function is derived. This bond potential exhibits harmonic behavior near the equilibrium state and anharmonic behavior for large bond stretches tending to a characteristic energy plateau (akin to the Lennard-Jones and Morse bond potentials). Using this bond potential, approximate yet highly-accurate analytical functions for bond stretch and chain force dependent upon chain stretch are established. Then, using this polymer chain model, a stochastic thermal fluctuation-driven chain rupture framework is developed. This framework is based upon a force-modified tilted bond potential that accounts for distortional bond potential energy, allowing for the derivation and subsequent calculation of the dissipated chain scission energy. Here, the cases of rate-dependent and rate-independent scission are accounted for throughout the rupture framework. The impact of Kuhn segment number on chain rupture behavior is also investigated. The model is fit to single chain mechanical response data collected from atomic force microscopy tensile tests for validation and to glean deeper insight into the molecular physics taking place. Due to their analytical nature, this polymer chain model and the associated rupture framework can in the future be implemented in finite element models accounting for fracture and fatigue in polydisperse elastomer networks.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Radiation-Enhanced Anion Transport in Hematite

The influence of radiation-induced (1 MeV energy H+ to ~0.1 dpa at 450°C), non-equilibrium point defect populations on mass transport is studied with an integrated campaign of experimental and theoretical methods. Using epitaxial thin films of hematite with embedded 18O tracer layers and nanoscale atom probe tomography measurements, it is shown that anion self-diffusion is enhanced by at least 2 orders of magnitude under irradiation compared to thermal diffusion alone. Complementary scanning transmission electron microscopy of vacuum annealed specimens, reveals associated microstructural changes in the oxide films, including local phase transformation to Fe3O4 and the development of nanoscale voids from vacancy coalescence. Point defect formation and migration energies were computed from density functional theory and applied within the context of chemical rate theory to analyze contributions from both interstitial and vacancy mechanisms to self-diffusion in thermal and irradiation conditions. Comparisons are made between calculated, literature and newly measured self-diffusion values, revealing good agreement on the magnitude of radiation-enhanced anion diffusion. Further, the model suggests a transition from vacancy to interstitialcy mechanisms at low temperatures and high oxygen activity, consistent with the varied activation energies reported from prior studies.

iron oxides, isotopic tracers, atom probe tomograp↗
DOE OSTI2021

The coupling and competition of crystallization and phase separation, correlating thermodynamics and kinetics in OPV morphology and performances

The active layer morphology transition of organic photovoltaics under non-equilibrium conditions are of vital importance in determining the device power conversion efficiency and stability; however, a general and unified picture on this issue has not been well addressed. Using combined in situ and ex situ morphology characterizations, morphological parameters relating to kinetics and thermodynamics of morphology evolution are extracted and studied in model systems under thermal annealing. The coupling and competition of crystallization and demixing are found to be critical in morphology evolution, phase purification and interfacial orientation. A unified model summarizing different phase diagrams and all possible kinetic routes is proposed. The current observations address the fundamental issues underlying the formation of the complex multi-length scale morphology in bulk heterojunction blends and provide useful morphology optimization guidelines for processing devices with higher efficiency and stability.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Measuring the density, viscosity, and surface tension of molten titanates using electrostatic levitation in microgravity

Rare earth and barium titanates are useful as ferroelectric, dielectric, and optical materials. Measurements of their thermophysical properties in the liquid state can help guide melt processing technologies for their manufacture and advance understanding of fragile liquids' behavior and glass formation. Here, we report the density, thermal expansion, viscosity, and surface tension of molten BaTi 2 O 5 , BaTi 4 O 9 , and 83TiO 2 -17RE 2 O 3 (RE = La or Nd). Measurements were made using electrostatic levitation and droplet oscillation techniques in microgravity, which provide access to quiescent liquid droplets and deep supercooling of 510–815 K below the equilibrium melting points. Densities were measured over 900–2400 K. Viscosities were similar for all four compositions, increasing from ~10 mPa s near 2100 K to ~30 mPa s near 1750 K. Surface tensions were 450–490 dyn cm -1 for the rare earth titanates and 383–395 dyn cm -1 for the barium titanates; surface tensions of all compositions had small or negligible temperature dependence over 1700–2200 K. In conclusion, for solids recovered after melt quenching, x-ray microtomography revealed the fracture mechanics in crystalline products and minimal internal porosity in glass products, likely arising from entrapped gas bubbles. Internal microstructures were generally similar for products processed either in microgravity or in a terrestrial aerodynamic levitator.

42 ENGINEERING↗
DOE OSTI2025

Plant Reload Optimization (prlo)

The PRLO framework is built on a modular and extensible architecture that tightly couples advanced evolutionary optimization algorithms with nuclear fuel depletion solvers (i.e., nuclear physics neutronics code). It supports exploring complex, high-dimensional design spaces constrained by user-specified operational, safety, and economic constraints. Objectives such as minimizing fresh fuel enrichment, flattening radial and axial power distributions, and maximizing discharge burnup are evaluated. PRLO’s equilibrium cycle optimization capability enables the identification of core configurations that maintain fuel cycle sustainability over extended planning horizons. Its integration with the RAVEN platform facilitates optimization of loading patterns or fuel shuffling schemes across multiple cycles. The interface with SIMULATE, a licensed industry-standard nodal code developed by Studsvik, ensures accurate neutronic and thermal-hydraulic feedback for reactor core design. PRLO’s automated workflow engine supports iterative design refinement, enabling utilities to streamline core design processes and meet evolving performance and regulatory targets.

Kim, Junyung [Idaho National Laboratory] (00090005↗
DOE OSTI2026

Metastable formation and disordering kinetics of body-centered orthorhombic CrNi 2 laths in a Cr-Ni binary alloy

Chromium (Cr) alloys combine low density with high-temperature strength but suffer from brittleness and rapid softening. Body-centered orthorhombic (BCO) CrNi 2 laths have been shown to improve Cr-alloy high-temperature strength retention, yet their thermal stability and transformation behavior remain unclear. CrNi 2 conventionally forms through long-range ordering from a face-centered cubic (FCC) phase. Using multiscale microscopy and neutron diffraction, we show that CrNi 2 instead nucleates from a body-centered cubic (BCC) matrix, in a binary Cr 85 Ni 15 alloy. Despite aging above the equilibrium ordering temperature, CrNi 2 persists for ∼220 h at 760 °C before fully disordering, indicating sluggish transformation kinetics. The formation and decomposition follow differing transformation pathways (BCC → CrNi 2 → FCC), revealing kinetic asymmetry not captured by equilibrium Cr-Ni phase diagrams. In conclusion, these findings redefine the transformation behavior of CrNi 2 , establish its metastable kinetic window, and suggest alloying strategies to stabilize CrNi 2 for precipitation strengthening of high-temperature Cr-alloys.

CrNi2 transformation kinetics↗