Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Forming mechanism of equilibrium and non-equilibrium metallurgical phases in dissimilar aluminum/steel (Al–Fe) joints

Abstract Forming metallurgical phases has a critical impact on the performance of dissimilar materials joints. Here, we shed light on the forming mechanism of equilibrium and non-equilibrium intermetallic compounds (IMCs) in dissimilar aluminum/steel joints with respect to processing history (e.g., the pressure and temperature profiles) and chemical composition, where the knowledge of free energy and atomic diffusion in the Al–Fe system was taken from first-principles phonon calculations and data available in the literature. We found that the metastable and ductile (judged by the presently predicted elastic constants) Al 6 Fe is a pressure ( P ) favored IMC observed in processes involving high pressures. The MoSi 2 -type Al 2 Fe is brittle and a strong P -favored IMC observed at high pressures. The stable, brittle η-Al 5 Fe 2 is the most observed IMC (followed by θ-Al 13 Fe 4 ) in almost all processes, such as fusion/solid-state welding and additive manufacturing (AM), since η-Al 5 Fe 2 is temperature-favored, possessing high thermodynamic driving force of formation and the fastest atomic diffusivity among all Al–Fe IMCs. Notably, the ductile AlFe 3 , the less ductile AlFe, and most of the other IMCs can be formed during AM, making AM a superior process to achieve desired IMCs in dissimilar materials. In addition, the unknown configurations of Al 2 Fe and Al 5 Fe 2 were also examined by machine learning based datamining together with first-principles verifications and structure predictions. All the IMCs that are not P- favored can be identified using the conventional equilibrium phase diagram and the Scheil-Gulliver non-equilibrium simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Rapid depletion analysis of flowing-pebble reactor systems at equilibrium using SCALE

Several high-temperature gas-cooled reactor concepts (and more recently, salt-cooled designs such as the fluoride salt-cooled high-temperature reactor) feature core designs employing continuously circulating fuel pebbles. These reactor designs permit both continuous online refueling of fuel elements as well as higher overall achievable discharge burnups. However, rapid calculation of time-dependent fuel isotopic inventories proves challenging for this class of dynamic systems with current analysis tools. While iterative approaches employing coupled neutron transport have been developed to solve this issue, rapid depletion analysis techniques are needed to calculate time-dependent inventories for individual pebbles and batches (and thus the construction of full- core inventory at equilibrium). We propose a depletion analysis strategy for this type of system for cores at equilibrium. Drawing upon previous neutronic analysis of the PBMR-400 equilibrium core, we demonstrate the viability of developing collapsed one-group cross section libraries suitable for performing rapid depletion analyses with SCALE. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Free-boundary MRxMHD equilibrium calculations using the stepped-pressure equilibrium code

Here the stepped-pressure equilibrium code (SPEC) (Hudson et al 2012 Phys. Plasmas 19 , 112 502) is extended to enable free-boundary multi-region relaxed magnetohydrodynamic (MRxMHD) equilibrium calculations. The vacuum field surrounding the plasma inside an arbitrary 'computational boundary', , is computed, and the virtual-casing principle is used iteratively to compute the normal field on so that the equilibrium is consistent with an externally produced magnetic field. Recent modifications to SPEC are described, such as the use of Chebyshev polynomials to describe the radial dependence of the magnetic vector potential, and a variety of free-boundary verification calculations are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Air Exchange Rate Impact on Activity Equilibrium Factors and Inhalation Fractional Equilibrium Factors for Rn, Xe, Kr, Ar, Ne, and Their Progeny in Vapor Intrusion, Risk, and Dose Models

Exposure to the radioactive noble gasses, especially radon, is of high concern and poses a significant risk to humans in an indoor air environment as the second leading cause of lung cancer in the United States. To evaluate and minimize the risks posed by these gasses, it is important to understand their radiological and physical properties. The EPA’s Radon Vapor Intrusion Screening Level (RVISL) calculator calculates indoor air RVISLs based on target working levels (WLs), target excess lifetime cancer risk (ELCR), and annual dose limits for the actinon (Rn-219), thoron (Rn-220), and radon (Rn-222) decay series. The RVISLs are based on inhalation and submersion in gas cloud exposure routes for residential and commercial settings. The RVISLs are analogous to preliminary remediation goals (PRGs) and dose compliance concentrations (DCCs), where the isotope-specific values are in units of activity concentration (activity per unit volume). If the concentration of a parent isotope of radon or its progeny is found to exceed the RVISL, then further action to ensure cleanup of the contaminant may be necessary. In residential and commercial settings, the RVISLs will vary based on the air exchange rate present. The EPA’s Radionuclide PRG and DCC Calculators also assess the risk/dose from noble gases in the air due to household use of water like showering. In this study, a computational method in MATLAB was developed to determine the impact of the air exchange rate on the activity equilibrium factor (A eq ) and the inhalation fractional equilibrium factor (F eq ). Both factors are values that reflect the equilibrium concentrations of progeny to their parent in the air. These factors have a direct impact on the RVISL, PRG, and DCC calculations of WL, ELCR, and annual dose, respectively. This study builds on a previous report that only focused on actinon, thoron, and radon by revisiting the original A eq and F eq calculation methods, as well as including the values for the Rn-207, Rn-209, Rn-210, Rn-211, Rn-215, Rn-216, Rn-217, Rn-218, Rn-223, Ne-24, Ar-42, Ar-43, Ar-44, Kr-74, Kr-75, Kr-76, Kr-77, Kr-88, Kr-89, Xe-120, Xe-121, Xe-122, Xe-123, Xe-135m, and Xe-138 decay chains, which are not currently available in literature. The EPA’s RVISL calculator will be updated to include the new A eq and F eq values for the actinon, thoron, and radon decay chains, while the rest of the calculators will incorporate all the new A eq and F eq values as appropriate.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2020

Comparison of equilibrium climate sensitivity estimates from slab ocean, 150-year, and longer simulations

We compare equilibrium climate sensitivity (ECS) estimates from pairs of long (= 800-year) control and abruptly quadrupled CO2 simulations with shorter (150, 300 year) coupled atmosphere-ocean simulations and Slab Ocean Models (SOM). Consistent with previous work, ECS estimates from shorter coupled simulations based on annual averages for years 1-150 underestimate those from SOM (-8% ± 13%) and long (-14% ± 8%) simulations. Analysis of only years 21-150 improved agreement with SOM (-2% ± 14%) and long (-8% ± 10%) estimates. Use of pentadal averages for years 51-150 results in improved agreement with long simulations (-4% ± 11%). While ECS estimates from current generation US models based on SOM and coupled annual averages of years 1-150 range from 2.6°C to 5.3°C, estimates based longer simulations of the same models range from 3.2°C to 7.0°C. Such variations between methods argues for caution in comparison and interpretation of ECS estimates across models.

Dunne, John P.↗
DOE OSTI2026

Equilibrium and non-equilibrium effects in high pressure phase transformations of carbon

The behavior of carbon in the range 1–100 GPa and 1–10 kK is central to problems in planetary interiors, inertial confinement fusion targets, and high-pressure synthesis of carbon-based materials, but experiments in this regime are difficult and often provide only indirect constraints on phase behavior. As a result, phase boundary loci, structure, and limits of metastability at high pressure remain uncertain. In this work, machine-learning enhanced atomistic simulations are used to address this knowledge gap. We determine the melt line up to 100 GPa, the graphite-diamond phase boundary up to the melt line, and analyze structure of the coexisting phases. We show that the coexisting liquid evolves smoothly with pressure without evidence for a first-order liquid–liquid transition. Orientation-resolved graphite melting simulations indicate that basal-plane interfaces develop a dewetting layer and undergo layer-by-layer melting, producing kinetic hysteresis and an apparent orientation dependence of the melt line. Non-equilibrium quenches from the melt are used to construct a kinetically limiting graphite–diamond phase boundary for rapid quenches from above the melt line, and show that graphite is metastable at pressures of up to ≈ 25 GPa. These results provide bounds on equilibrium and metastable behavior in carbon relevant for interpreting high-pressure experiments and for designing synthesis pathways to specific carbon microstructures.

Lyu, Yanjun [Department of Materials Science and E↗
DOE OSTI2021

Phase Equilibrium of Water with Hexagonal and Cubic Ice Using the SCAN Functional

Machine learning models are rapidly becoming widely used to simulate complex physicochemical phenomena with ab initio accuracy. Here, we use one such model as well as direct density functional theory (DFT) calculations to investigate the phase equilibrium of water, hexagonal ice (Ih), and cubic ice (Ic), with an eye toward studying ice nucleation. The machine learning model is based on deep neural networks and has been trained on DFT data obtained using the SCAN exchange and correlation functional. We use this model to drive enhanced sampling simulations aimed at calculating a number of complex properties that are out of reach of DFT-driven simulations and then employ an appropriate reweighting procedure to compute the corresponding properties for the SCAN functional. This approach allows us to calculate the melting temperature of both ice polymorphs, the driving force for nucleation, the heat of fusion, the densities at the melting temperature, the relative stability of ices Ih and Ic, and other properties. We find a correct qualitative prediction of all properties of interest. In some cases, quantitative agreement with experiment is better than for state-of-the-art semiempirical potentials for water. Furthermore, our results also show that SCAN correctly predicts that ice Ih is more stable than ice Ic.

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

On the Equilibrium Electrostatic Potential and Light-Induced Charge Redistribution in Halide Perovskite Structures

Lead halide perovskites are semiconductor materials which are employed as nonintentionally doped absorbers inserted between two selective carrier transport layers (SCTL), realizing a p-i-n or n-i-p heterojunction. In our study, we have developed and investigated a lateral device, based on methylammonium lead iodide (MAPbI3) in which the p-i-n heterojunction develops in the horizontal direction. Our research suggests that the effective doping level in the MAPbI3 film should be very low, below 1012 cm-3. Along the vertical direction, this doping level is not enough to screen the electric field of the buried heterojunction with the SCTL. The perovskite work function is therefore affected by the work function of the SCTL underneath. From drift-diffusion simulations, we show that intrinsic perovskite-SCTL structures develop mV range surface photovoltages (SPVs) under continuous illumination. However, perovskite-SCTL structures can develop SPVs of hundreds of mV, as confirmed by our measurements. We therefore analyzed the compatibility between low doping and low defect densities in the perovskite layer and such high SPV values using numerical modeling. It is shown that these high SPV values could originate from electronic processes due to large band offsets in the buried perovskite-SCTL heterojunctions, or at the SCTL-transparent conductive oxide (TCO) buried heterojunction. However, such electronic processes can hardly explain the long SPV persistence after switching off the illumination.

defects↗
DOE OSTI2022

The roles of kinematic constraint and diffusion in non-equilibrium solid state phase transformations of Ti-6Al-4V

A solid state phase transformation of Ti-6Al-4V was studied using high speed in situ x-ray diffraction measurements made during rapid cooling of a cold metal transfer arc weld bead deposited onto a water cooled substrate. Analysis of body centered cubic (BCC) and hexagonal close packed (HCP) lattices revealed an abrupt, nonlinear shift in the lattice parameters of both phases just after the HCP phase had nucleated. Postmortem transmission electron microscopy confirmed that V diffusion was mostly suppressed during cooling. Together, these results indicate that at this cooling rate of approximately 10 4 K/s, which is representative of cooling rates of many additive manufacturing and welding processes, kinematic coherency of the BCC–HCP interfaces gives rise to the anomalous lattice expansion and contraction behaviors of both phases during the initial nucleation and growth stages of (mostly) martensitic transformation from BCC to HCP; the role of diffusion in such lattice anomalies is shown to be minimal.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Quantum criticality and topology in non-equilibrium systems

A key goal of condensed matter physics research is to identify new phases of matter, and to understand the universal features of the phase transitions between them. In the past decade, physicists have uncovered a wealth of new phases with interesting surface properties, exemplified by the theoretical prediction and subsequent experimental discovery of topological insulators and superconductors. Traditional condensed matter systems are usually in a thermal equilibrium state and typically at very low temperature. Very recently, experimental advances have sparked interest in the non-equilibrium setting. Non-equilibrium systems can host new phases and phenomena with no equilibrium counterpart, and could also enable robust ways to build quantum memory devices to store and manipulate quantum information in a coherent manner. These phases and phenomena are inherently “dynamical”: they are described not by changes in the arrangement or structure of the constituent particles, but instead marked by sharp distinctions in how the particles move and exchange energy or quantum information. The discovery of robust non-equilibrium phases raises many fundamental questions: Can we develop a systematic theory of states of matter and of dynamical transitions between such states? How can such states be realized and probed experimentally? The main goal of this project was to explore the emergence of topological phases and quantum criticality (two cornerstones of modern condensed matter physics in equilibrium) in such non-equilibrium quantum systems. Specific goals included (1) using tensor networks to efficiently represent non-equilibrium states of matter and their phase transitions; (2) studying and designing new probes for periodically driven systems; and (3) developing analytic and numerical tools to analyze non-equilibrium topological phase transitions. Advances in these directions were achieved using novel techniques appropriate to study the non-equilibrium dynamics in many-body quantum systems combining strong interactions and randomness. Taken together, these results provide a new conceptual framework for understanding the emergence of quantum critical and topological properties in quantum systems far from thermal equilibrium.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Experimental and theoretical determinations of hydrogen isotopic equilibrium in the system CH 4 —H 2 —H 2 O from 3 to 200°C

The stable isotopic composition of methane (CH 4 ) is commonly used to fingerprint natural gas origins. Over the past 50 years, there have been numerous proposals that both microbial and thermogenic CH 4 can form in or later attain hydrogen isotopic equilibrium with water (H 2 O) and carbon isotopic equilibrium with carbon dioxide (CO 2 ). Evaluation of such proposals requires knowledge of the equilibrium fractionation factors between CH 4 and H 2 O or CO 2 at the temperatures where microbial and thermogenic CH 4 form in or are found in the environment, which is generally less than 200°C. Experimental determinations of these fractionation factors are only available above 200°C, requiring extrapolation of these results beyond the calibrated range or the use of theoretical calculations at lower temperatures. Here, we provide a calibration of the equilibrium hydrogen isotopic fractionation factor for CH 4 and hydrogen gas (H 2 ) ( D α CH4(g)–H2(g) ) based on experiments using γ-Al 2 O 3 and Ni catalysts from 3 to 200°C. Results were regressed as a 2 nd order polynomial of 1000 × ln D α CH4(g)–H2(g) vs. 1/T (K -1 ) yielding: 1000 × l n D α C H 4 ( g ) - H 2 ( g ) = 3.5317 × 10 7 T 2 + 2.7749 × 10 5 T - 179.48 We combine this calibration with previous experimental determinations of hydrogen isotope equilibrium between H 2 , H 2 O(g), and H 2 O(l) and we provide an interpolatable experimental calibration of 1000 × ln D α CH4(g)–H2O(l) from 3 to 200°C. Our resulting 4th order polynomial is the following equation: 1000 × l n D α C H 4 ( g ) - H 2 O l = - 7.9443 × 10 12 T 4 + 8.7772 × 10 10 T 3 - 3.4973 × 10 8 T 2 + 5.4398 × 10 5 T - 382.05 At 3°C, the value from our calibration differs by 93‰ relative to what would be calculated based on the extrapolation of the only experimental calibration currently available to temperatures below its calibrated range (lowest temperature of 200°C; Horibe and Craig, 1995). We additionally provide new theoretical estimates of hydrogen isotopic equilibrium between CH 4 (g), H 2 (g), and H 2 O(g) and carbon isotopic equilibrium between CH 4 (g) and CO 2 (g) using Path Integral Monte Carlo (PIMC) calculations. Our PIMC calculations for hydrogen isotopic equilibrium between CH 4 and H 2 agree 1:1 with our experiments. Finally, we compile carbon and hydrogen isotopic measurements of CH 4 , CO 2 , and H 2 O from various environmental systems and compare observed differences between carbon and hydrogen isotopes to those expected based on isotopic equilibrium. We find that isotopic compositions of some microbial gases from marine sedimentary, coalbed, and shale environments are consistent with those expected for CH 4 H 2 O(l) hydrogen and CH 4 CO 2 carbon isotopic equilibrium. In contrast, microbial terrestrial and pure culture gases are not consistent with both CH 4 H 2 O(l) hydrogen and CH 4 CO 2 carbon isotopic equilibrium. Overall, these results are explained qualitatively using previously developed conceptual models that link free energy gradients available to microorganisms to the degree that their enzymes can promote isotope-exchange reactions between CH 4 , CO 2 , and H 2 O.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗