Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase equilibria”

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

Evaluation of thermal processing and properties of 422 martensitic stainless steel for replacement of 4140 steel in diesel engine pistons

The thermal and mechanical properties of martensitic stainless steel 422 were evaluated for suitability as a drop-in replacement for 4140 steel in next generation heavy-duty diesel engine (HDDE) pistons. The time and temperature of the austenitization and tempering steps were studied to achieve optimum materials performance in piston applications, including the balance of thermal and mechanical properties and resistance to long-term thermal aging. Reducing the tempering temperature from 700 to 600 °C caused a substantial increase in elevated temperature strength from 25 to 600 °C, but had no significant influence on thermal conductivity, suggesting that thermal conductivity in 422 is dominated largely by composition and distribution of alloying elements and mostly independent of the sub-grain structure size and precipitate size. Compared to the current HDDE piston alloy 4140, 422 exhibits substantially higher elevated temperature strength and lower thermal conductivity, the latter which will cause 422 to operate at higher temperatures in pistons, possibly requiring a piston redesign to take advantage of the improved high temperature strength of 422. Piston material selection and alloy design strategies with potential to mitigate some of the shortcomings of martensitic stainless steels, such as 422, as drop-in replacements are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Thermophysical Properties of NaCl–UCl 3 –PuCl 3 Molten Salts: A Combined Computational and Experimental Study

Actinide-bearing molten salts for use as fuels are an essential part of next generation molten salt reactors. Yet, numerous multicomponent salt mixtures are underdeveloped or have not been investigated. Here, this study, based on a combination of experimental and modeling techniques, is dedicated to determining and understanding a variety of properties of the ternary system of NaCl–UCl 3 –PuCl 3 , which represents a scenario for burnup of NaCl–UCl 3 fuel, at two compositions (∼10 and 5 mol % PuCl 3 in eutectic NaCl–UCl 3 pseudobinary) and a range of temperatures. Evaluation of the heat flow and mass loss data showed the 0.61NaCl–0.30UCl 3 –0.09PuCl 3 salt had a melting temperature of 551 ± 5 °C. Two additional thermal effects were observed occurring at approximately 410 and 494 °C. The transition occurring at 410 °C may be due to the presence of oxide in the salt. Extrapolation of thermodynamic data indicates the transition occurring at 494 °C is due to the formation of a liquid phase. Experimental testing determined the density of this system is a linear function of temperature and can be represented by the equation ρ = 4.014–0.0010T(°C), R 2 = 0.992. Additionally, by using atomistic modeling, we found that increasing the PuCl 3 content from 5 to 10 mol % led to the formation of larger Pu 3+ clusters and slower transport of ions.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Modeling of non-equilibrium partitioning in direct metal deposited copper–iron concentrated alloys

The direct metal deposition (DMD) additive manufacturing process produces high cooling rates within a small melt pool and can lead to high amounts of solute trapping. These high cooling rates limit diffusion and lead to the formation of non-equilibrium phases. In this work, we utilize a numerical model to calculate the degree of solute trapping, defined as non-equilibrium partitioning. A theoretical case with overall composition fixed to 50Cu-50Fe at.%. was performed to observe the influence of increasing solidification rates. We then simulate DMD of equimolar Cu-Fe powder printed on mild steel substrate and the calculated non-equilibrium phase compositions were consistent with experimental observations reported earlier on this alloy composition. For a single deposited track, cooling rates are high enough to yield significant solute trapping. Here, the degree of solute trapping is highest near the free surface and has a gradient that correlates with the cooling rate gradient.

36 MATERIALS SCIENCE↗
DOE OSTI2025

CALPHAD-based Bayesian optimization to accelerate alloy discovery for high-temperature applications

Two crucial properties influencing the performance of high-temperature alloys are coefficient of thermal expansion (CTE) and phase constitution. It is desirable to have alloys with low CTE, which reduces CTE mismatch with the surface oxide and the likelihood of oxide spallation. Reducing the amount of brittle intermetallic phases such as Sigma (σ) enhances alloy ductility and processability. Here, we propose a multi-objective Bayesian Optimization (BO) model to simultaneously minimize CTE (at an operational temperature of 1150 °C) and T σ (temperature when the Sigma phase completely dissolves in the metal matrix), properties which are obtained from high-throughput CALculation of PHAse Diagrams (CALPHAD). The model successfully identifies several alloys with CTE ≤ 2 × 10 –5 /K and T σ ≤ 500 °C by exploring just 7% of the nickel–chromium–cobalt–aluminum–iron (Ni–Cr–Co–Al–Fe) composition space. Such multi-objective alloy design frameworks can be used to inform additive manufacturing experiments and accelerate alloy discovery for high-temperature energy applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The crystal structure of feitknechtite (β-MnOOH) and a new MnOOH polymorph

Studies suggest that feitknechtite (β-MnOOH) is a prevalent, and perhaps necessary, intermediate phase during the synthesis of birnessite-like phases, the abiotic oxidation of Mn2+, and the transformation of biogenic hexagonal phyllomanganates to more complex Mn oxides in laboratory and natural systems. Researchers have generally described feitknechtite as consisting of pyrochroite-like (or cadmium iodide-like) Mn-O octahedral layers, but a detailed crystal structure has not been reported. For this work, we used TEM/SAED and powder XRD and Rietveld refinements to derive the unit cell and, for the first time, report a complete structure description for feitknechtite (β-MnOOH). Rietveld refinements were also completed for three natural feitknechtite/hausmannite samples, and time-resolved synchrotron XRD experiments were used to follow the thermal transformation of feitknechtite to hausmannite. Additionally, we identified and report the structure for a second, and perhaps novel, MnOOH polymorph (proposed designation ε-MnOOH), mixed with the synthetic feitknechtite, that is similar to β-MnOOH but with a different layer stacking.

58 GEOSCIENCES↗
DOE OSTI2022

Thermodynamic modeling of adsorption at the liquid-solid interface

Adsorptive separation techniques are significantly energy efficient in comparison to conventional thermal separation techniques such as distillation. Despite extensive research and development activities undertaken for mixed-gas adsorption, the use of adsorption techniques for the separation of multicomponent liquid mixtures is still limited. A major barrier is the lack of accurate adsorption thermodynamic models, which form the scientific foundation of process simulation of such systems, making the translation to industrial scale challenging. In this work, we have rigorously computed the surface excess of adsorption for six binary liquid mixtures on silica gel at 303 K using the frameworks of the generalized Langmuir isotherm (gL) and the adsorbed solution theory (AST). The six binary liquid mixtures were formed by the pair-wise combinations of four components: benzene, 1,2-dichloroethane, cyclohexane, and n-heptane. We have based our calculations by considering simultaneous equilibria of three phases: saturated vapor phase, bulk liquid phase, and adsorbed phase. Further, the composition of the corresponding saturated vapor phase was estimated by the Nonrandom Two-Liquid activity coefficient model and experimental vapor-liquid equilibria data. The activity coefficients of the adsorbed phase, the central issue of multicomponent adsorption thermodynamics, were calculated using the adsorption Nonrandom Two-Liquid activity coefficient model. Devoid of simplifying assumptions, gL and AST provide rigorous thermodynamic frameworks for adsorption equilibria of multicomponent liquid mixtures.

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

Speciation and Organic Phase Structure in Nitric Acid Extraction with Trioctylamine

Understanding chemical speciation and intermolecular interactions in multicomponent liquids is essential to understanding their phase and chemical equilibria, which underpin chemical separation processes, including solvent extraction. Here we report on the extraction of nitric acid from its aqueous solutions into organic solutions of trioctylamine (TOA) in toluene, investigated with spectroscopic, X-ray scattering, and computational tools to understand molecular speciation in the organic phase and its relationship with the nanoscale structure of the organic phase. Trends in acid and water extraction clearly show two and three regimes, respectively, indicating different stoichiometric relationships, but speciation of HNO 3 , water, and amine in these regimes is not apparent. 1 H NMR of the organic phase shows that there are at least two distinct acidic protons in the organic phase while ATR-FTIR results show that the organic phase with excess acid extraction is a mixture of trioctylammonium-nitrate ion pairs (TOAHNO 3 ), and undissociated HNO 3 molecules. Comparison with DFT-computed IR spectra show that the chain-like configurations of TOAHNO 3 HNO 3 H 2 O are favored over TOAHNO 3 H 2 OHNO 3 , i.e., direct interaction between the nitrate and HNO 3 molecules is more favored compared to a water-mediated interaction. SAXS of the organic phases were modeled as sums of Ornstein-Zernike (O-Z) scattering and a prepeak feature in the higher Q region that corresponds to extractant packing. The extraction of undissociated HNO 3 by the ion pairs leads to an increased X-ray scattering contrast in the organic phase without any significant change in the correlation length. These results show that the organic phase nanostructure is more sensitive to the concentration of TOAHNO 3 and is relatively unaffected by excess acid extraction. Finally, these findings will enable a molecular understanding of the mechanisms behind metal extraction from acidic media with basic extractants.

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

Process Modeling of CO 2 Absorption with Monoethanolamine Aqueous Solutions Using Rotating Packed Beds

A first-principle process simulation model is presented for the chemical absorption of carbon dioxide (CO 2 ) with monoethanolamine (MEA) aqueous solutions using rotating packed beds (RPB). Built on a proven rate-based packed bed absorber model, the RPB model rigorously simulates the phase and chemical equilibria at the vapor-liquid interface, the heat and mass transfer across the gas and liquid films, the fast reactions between MEA and CO 2 in the liquid film, and the RPB hydraulics. Here, estimation of the mass transfer rate across the liquid film is central to accurate simulation of the CO 2 absorption process with MEA aqueous solutions. We show that the literature lab-scale RPB data for CO 2 removal efficiency can be satisfactorily correlated by introducing a correction factor for the effective packing surface area predicted by the Onda correlation. Given the validated RPB model, we further show that, among the gas-phase mass transfer coefficient, the liquid-phase mass transfer coefficient, and the reaction rate constant for the reaction between amine and CO 2 , the reaction rate constant is the controlling step with the highest potential to enhance the CO 2 absorption performance in RPB.

42 ENGINEERING↗
DOE OSTI2024

Shape equilibria of vesicles with rigid planar inclusions

Motivated by recent studies of two-phase lipid vesicles possessing 2D solid domains integrated within a fluid bilayer phase, we study the shape equilibria of closed vesicles possessing a single planar, circular inclusion.

Chemistry↗
DOE OSTI2022

Effect of injected flux and current temporal phasing on self-organization in the HIT-SI3 experiment

The HIT-SI3 device at the University of Washington uses three oscillating inductive helicity injectors to form and sustain spheromak plasma equilibria. By adjusting the temporal phase of the injector waveforms with respect to each other, the toroidal spectrum of the imposed perturbations can be controlled. Using a recently implemented GPU-based control system, the available mode spectra were explored experimentally by scanning the space of relative injector phasing. In this space, significant variation in the toroidal mode spectrum ($n$ = 1, 2, 3) of the perturbations was observed. Additionally, variation in characteristics of driven equilibria was also observed, including a ≈ $30$% range in toroidal current gain ($I$ $Φ$ / $I$ $Inj$ ). Experimental results are compared with both a composite-equilibria and nonlinear dynamic model, including extended MHD simulations using the NIMROD code and composite Taylor state equilibria computed using the PSI-Tet code. In conclusion, qualitative agreement is seen with the nonlinear models, but not with composite-equilibria models, suggesting the use of nonlinear models to better capture observed plasma dynamics and provide predictive use for future experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Machine learning the metastable phase diagram of covalently bonded carbon

Abstract Conventional phase diagram generation involves experimentation to provide an initial estimate of the set of thermodynamically accessible phases and their boundaries, followed by use of phenomenological models to interpolate between the available experimental data points and extrapolate to experimentally inaccessible regions. Such an approach, combined with high throughput first-principles calculations and data-mining techniques, has led to exhaustive thermodynamic databases (e.g. compatible with the CALPHAD method), albeit focused on the reduced set of phases observed at distinct thermodynamic equilibria. In contrast, materials during their synthesis, operation, or processing, may not reach their thermodynamic equilibrium state but, instead, remain trapped in a local (metastable) free energy minimum, which may exhibit desirable properties. Here, we introduce an automated workflow that integrates first-principles physics and atomistic simulations with machine learning (ML), and high-performance computing to allow rapid exploration of the metastable phases to construct “metastable” phase diagrams for materials far-from-equilibrium. Using carbon as a prototypical system, we demonstrate automated metastable phase diagram construction to map hundreds of metastable states ranging from near equilibrium to far-from-equilibrium (400 meV/atom). We incorporate the free energy calculations into a neural-network-based learning of the equations of state that allows for efficient construction of metastable phase diagrams. We use the metastable phase diagram and identify domains of relative stability and synthesizability of metastable materials. High temperature high pressure experiments using a diamond anvil cell on graphite sample coupled with high-resolution transmission electron microscopy (HRTEM) confirm our metastable phase predictions. In particular, we identify the previously ambiguous structure of n -diamond as a cubic-analog of diaphite-like lonsdaelite phase.

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

Minimizing the Electromechanical Stresses in Poloidal Field Coils by Optimizing their Numbers and Locations using FREDA Framework

Poloidal field (PF) and central solenoid (CS) coils play a crucial role in sustaining the equilibrium and preserving the shape of highly confined tokamak plasmas. Ensuring that PF coil current and mechanical stress stay within superconducting and structural limitations is an important check in the design assessment. Minimizing the PF coil currents and mechanical stresses influences reliability, cost, and performance. A free-boundary MHD equilibrium code—FreeGS is employed within the fusion reactor design and assessment (FREDA) whole facility modeling (WFM) framework to construct the plasma equilibrium based on the configuration and currents in the PF coils. Here, we present the capability of the FreeGS code to minimize the currents, forces, and electromagnetic stresses on the PF coils by optimizing their number, sizes, structures, and locations while maintaining an MHD stable plasma configuration with a large confinement factor. The workflow is initialized with a configuration of plasma parameters and coils’ locations from the 0-D tokamak build systems code in the FREDA framework. Then, FreeGS is called to calculate the initial equilibrium at the minimum total current in PF coils. Thereafter, FreeGS’s internal optimizer minimizes the currents and hoop and central forces on the PF coils while maintaining the reference equilibrium. Finally, the input configuration is updated with the optimized parameters for equilibria over the ramp-up phase of a burning-plasma operation. FREDA’s whole facility optimization capability, which includes all magnetic field coil systems, blanket, vacuum vessel (VV), first wall, divertor, etc., is under development and out of the scope for this study.

Hassan, Ehab [ORNL] (ORCID:0000000181060301)↗
DOE OSTI2021

MB-Fit: Software infrastructure for data-driven many-body potential energy functions

Many-body potential energy functions (MB-PEFs), which integrate data-driven representations of many-body short-range quantum mechanical interactions with physics-based representations of many-body polarization and long-range interactions, have recently been shown to provide high accuracy in the description of molecular interactions from the gas to the condensed phase. Here, we present MB-Fit, a software infrastructure for the automated development of MB-PEFs for generic molecules within the TTM-nrg (Thole-type model energy) and MB-nrg (many-body energy) theoretical frameworks. Besides providing all the necessary computational tools for generating TTM-nrg and MB-nrg PEFs, MB-Fit provides a seamless interface with the MBX software, a many-body energy and force calculator for computer simulations. Given the demonstrated accuracy of the MB-PEFs, particularly within the MB-nrg framework, we believe that MB-Fit will enable routine predictive computer simulations of generic (small) molecules in the gas, liquid, and solid phases, including, but not limited to, the modeling of quantum isomeric equilibria in molecular clusters, solvation processes, molecular crystals, and phase diagrams.

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

Assessment of diffuse-interface methods for compressible multiphase fluid flows and elastic-plastic deformation in solids

This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook, Subramaniam et al., and Adler and Lele, in which artificial diffusion terms are added to the individual phase mass fraction transport equations and are coupled with the other conservation equations. The second method is the gradient-form approach that is based on the quasi-conservative method of Shukla et al., in which the diffusion and sharpening terms (together called regularization terms) are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. The third approach is the divergence-form approach that is based on the fully conservative method of Jain et al., in which the regularization terms are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. In the present study, all three diffuse-interface methods are used in conjunction with a four-equation, multicomponent mixture model, in which pressure and temperature equilibria are assumed among the various phases. The primary objective of this work is to compare these three methods in terms of their ability to: maintain constant interface thickness throughout the simulation; conserve mass, momentum, and energy; and maintain accurate interface shape for long-time integration. The second objective of this work is to consistently extend these methods to model interfaces between solid materials with strength. To assess and compare the methods, they are used to simulate a wide variety of problems, including (1) advection of an air bubble in water, (2) shock interaction with a helium bubble in air, (3) shock interaction and the collapse of an air bubble in water, and (4) Richtmyer–Meshkov instability of a copper–aluminum interface. The current work focuses on comparing these methods in the limit of relatively coarse grid resolution, which illustrates the true performance of these methods. In conclusion, this is because it is rarely practical to use hundreds of grid points to resolve a single bubble or drop in large-scale simulations of engineering interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Generalization of thermodynamic Langmuir isotherm for mixed‐gas adsorption equilibria

Abstract This work presents a comprehensive thermodynamic model for both pure component isotherms and mixed‐gas adsorption equilibria. A generalization of thermodynamic Langmuir isotherm, the proposed model assumes competitive adsorption of multiple adsorbates on adsorbent surface for mixed‐gas adsorption equilibria, and it applies an area‐based adsorption nonrandom two‐liquid activity coefficient model in the activity coefficient calculations for the adsorbate phase. The resulting generalized Langmuir isotherm properly captures both surface loading dependence and adsorbate phase composition dependence for mixed‐gas adsorption equilibria. The model is validated with accurate representations of gas adsorption equilibrium data for varieties of unary, binary, and ternary gas systems. The model results are further compared with those calculated from extended Langmuir isotherm and Ideal Adsorbed Solution Theory.

Hamid, Usman↗
DOE OSTI2024

A collision operator for describing dissipation in noncanonical phase space

The phase space of a noncanonical Hamiltonian system is partially inaccessible due to dynamical constraints (Casimir invariants) arising from the kernel of the Poisson tensor. When an ensemble of noncanonical Hamiltonian systems is allowed to interact, dissipative processes eventually break the phase space constraints, resulting in a thermodynamic equilibrium described by a Maxwell–Boltzmann distribution. However, the time scale required to reach Maxwell–Boltzmann statistics is often much longer than the time scale over which a given system achieves a state of thermal equilibrium. Examples include diffusion in rigid mechanical systems, as well as collisionless relaxation in magnetized plasmas and stellar systems, where the interval between binary Coulomb or gravitational collisions can be longer than the time scale over which stable structures are self-organized. Here, we focus on self-organizing phenomena over spacetime scales such that particle interactions respect the noncanonical Hamiltonian structure, but yet act to create a state of thermodynamic equilibrium. We derive a collision operator for general noncanonical Hamiltonian systems, applicable to fast, localized interactions. This collision operator depends on the interaction exchanged by colliding particles and on the Poisson tensor encoding the noncanonical phase space structure, is consistent with entropy growth and conservation of particle number and energy, preserves the interior Casimir invariants, reduces to the Landau collision operator in the limit of grazing binary Coulomb collisions in canonical phase space, and exhibits a metriplectic structure. We further show how thermodynamic equilibria depart from Maxwell–Boltzmann statistics due to the noncanonical phase space structure, and how self-organization and collisionless relaxation in magnetized plasmas and stellar systems can be described through the derived collision operator.

Boltzmann equation↗
DOE OSTI2023

Induced phase transformation in ionizable colloidal nanoparticles

Abstract Acid–base equilibria directly influence the functionality and behavior of particles in a system. Due to the ionizing effects of acid–base functional groups, particles will undergo charge exchange. The degree of ionization and their intermolecular and electrostatic interactions are controlled by varying the pH and salt concentration of the solution in a system. Although the pH can be tuned in experiments, it is hard to model this effect using simulations or theoretical approaches. This is due to the difficulty in treating charge regulation and capturing the cooperative effects in a colloidal suspension with Coulombic interaction. In this work, we analyze a suspension of ionizable colloidal particles via molecular dynamics (MD) simulations, along with Monte Carlo simulations for charge regulation (MC-CR) and derive a phase diagram of the system as a function of pH. It is observed that as pH increases, particles functionalized with acid groups change their arrangement from face-centered cubic (FCC) packing to a disordered state. We attribute these transitions to an increase in the degree of charge polydispersity arising from an increase in pH. Our work shows that charge regulation leads to amorphous solids in colloids when the mean nanoparticle charge is sufficiently high. Graphical abstract

36 MATERIALS SCIENCE↗