Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermodynamics of computing”

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 325 records · Page 18

A rapid method for the computation of equilibrium chemical composition of air to 15000 K

A rapid computational method has been developed to determine the chemical composition of equilibrium air to 15000 K. Eleven chemically reacting species, i.e., O2, N2, O, NO, N, NO+, e-, N+, O+, Ar, and Ar+ are included. The method involves combining algebraically seven nonlinear equilibrium equations and four linear elemental mass balance and charge neutrality equations. Computational speeds for determining the equilibrium chemical composition are significantly faster than the often used free energy minimization procedure. Data are also included from which the thermodynamic properties of air can be computed. A listing of the computer program together with a set of sample results are included.

Prabhu, Ramadas K.↗

User's manual for the one-dimensional hypersonic experimental aero-thermodynamic (1DHEAT) data reduction code

A FORTRAN computer code for the reduction and analysis of experimental heat transfer data has been developed. This code can be utilized to determine heat transfer rates from surface temperature measurements made using either thin-film resistance gages or coaxial surface thermocouples. Both an analytical and a numerical finite-volume heat transfer model are implemented in this code. The analytical solution is based on a one-dimensional, semi-infinite wall thickness model with the approximation of constant substrate thermal properties, which is empirically corrected for the effects of variable thermal properties. The finite-volume solution is based on a one-dimensional, implicit discretization. The finite-volume model directly incorporates the effects of variable substrate thermal properties and does not require the semi-finite wall thickness approximation used in the analytical model. This model also includes the option of a multiple-layer substrate. Fast, accurate results can be obtained using either method. This code has been used to reduce several sets of aerodynamic heating data, of which samples are included in this report.

Hollis, Brian R.↗

Effect of iodides on thermal behavior and phase partitioning in LiCl-KCl

Liquid-fueled molten salt reactors (MSRs) are designed to operate with fissile materials and, ultimately, fission products dissolved in the primary molten salt coolant. Understanding the speciation and transport of iodine—a high-yield fission product—is essential because this element’s accidental release poses significant environmental concerns due to its capacity to be readily absorbed by the human thyroid gland. Here, we report the impact of iodide species (LiI and KI) on phase transitions, phase distribution, and phase stability in LiCl–KCl-eutectic salt mixtures. The study employed a combination of computational and experimental techniques, including thermodynamic FactSage calculations, differential scanning calorimetry, and high-temperature X-ray diffraction. The results indicate that the presence of iodide (10–25 wt%) significantly alters the melting behavior of the LiCl–KCleutectic system. Adding 10 wt% LiI has a more-pronounced effect than 10 wt% KI, as LiI converts to KI, leading to formation of LiCl, thereby, altering the LiCl-KCl ratio which significantly affects the melting temperature of the mixture. Furthermore, the evolution of crystalline structure, solid-fraction composition, and the dynamics of mixed-halide solid–liquid partitioning as a function of temperature indicate the potential for selective iodide separation from chloride-salt mixtures via solid–liquid separation techniques. Overall, the presented findings provide valuable insights that are beneficial for the design and operation of MSRs, as well as for the safe handling and effective processing of used nuclear fuel using advanced pyrochemical techniques.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Integrated modeling methodology for ash agglomeration in poly-disperse fluidized beds using particle population framework

This article discusses a unique agglomeration modeling methodology developed based on binary collisions to combine the effects of heterogeneity in ash chemistry and granular physics. A simple population balance is defined to find changes in the particle size distribution (PSD) of a fluidized bed. Thermodynamic equilibrium calculations and a computational fluid dynamics (CFD) code, are used to obtain hydrodynamic parameters. A method to calculate and use a distribution of collision frequencies in fluidized beds, corresponding to a poly-disperse particle size distribution, was developed in order to incorporate the particle level heterogeneities. A distribution of collision frequencies obtained for poly-dispersed particles showed a three orders of magnitude higher collision frequency amongst the smaller particles than the coarser ones, at the initiation of agglomeration. Ash agglomeration occurred when the slag amount (binder) was less than 10 wt% at temperatures <850 °C for Pittsburgh No. 8 coal.

42 ENGINEERING↗

Solvent Isotherms and Structural Transitions in Nanoparticle Superlattice Assembly

Here, we introduce a Molecular Theory for Compressible Fluids (MOLT-CF) that enables us to compute free energies and other thermodynamic functions for nanoparticle superlattices with any solvent content, including the dry limit. Quantitative agreement is observed between MOLT-CF and united-atom molecular dynamics simulations performed to assess the reliability and precision of the theory. Among other predictions, MOLT-CF shows that the amount of solvent within the superlattice decreases approximately linearly with its vapor pressure and that in the late stages of drying, solvent-filled voids form at lattice interstitials. Applied to single-component superlattices, MOLT-CF predicts fcc-to-bcc Bain transitions for decreasing vapor pressure and for increasing ligand length, both in agreement with experimental results. We explore the stability of other single-component phases and show that the C14 Frank-Kasper phase, which has been reported in experiments, is not a global free-energy minimum. Implications for precise assembly and prediction of multicomponent nanoparticle systems are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Alkane Dehydrogenation and H/D Exchange by a Cationic Pincer-Ir(III) Hydride: Cooperative C–H Addition and β -H Elimination Modes Induce Anomalous Selectivity

We report that the cationic iridium complex ( iPr PCP)IrH + catalyzes the transfer-dehydrogenation of alkanes to give alkenes and hydrogen isotope exchange (HIE) of alkanes and arenes. Contrary to established selectivity trends found for C–H activation by transition metal complexes, strained cycloalkanes, including cyclopentane, cycloheptane, and cyclooctane, undergo C–H addition much more readily than n-alkanes, which in turn are much more reactive than cyclohexane. Aromatic C–H bonds also undergo H/D exchange much less rapidly than those of the strained cycloalkanes, but much more favorably than cyclohexane. The order of reactivity toward dehydrogenation correlates qualitatively with the reaction thermodynamics, but the magnitude is much greater than can be explained by thermodynamics. Accordingly, the cycloalkenes corresponding to the strained cycloalkanes undergo hydrogenation much more readily than cyclohexene, despite the less favorable thermodynamics of such hydrogenations. Here, computational (DFT) studies allow rationalization of the origin of reactivity and the unusual selectivity. Specifically, the initial C–H addition is strongly assisted by β-agostic interactions, which are particularly favorable for the strained cycloalkanes. Subsequent to α-C–H addition, the H atom of the β-agostic C–H bond is transferred directly to the hydride ligand of ( iPr PCP)IrH + to give a dihydrogen ligand.

02 PETROLEUM↗

Observation of laser ablation of silicon as a function of pulse length at constant fluence via time-resolved x-ray spectroscopy

We investigate the ablation of silicon as a function of laser pulse length at a constant fluence using time-resolved x-ray spectroscopy data obtained from OMEGA EP experiments at the University of Rochester's Laboratory for Laser Energetics. Our targets consisted of three-layer planar structures composed of Si (50 μm), Cu (25 μm), and SiO 2 (500 μm) layers. The Si layer was irradiated by a 351-nm laser with varying pulse widths of 250 ps, 500 ps, 1 ns, and 10 ns while maintaining a constant fluence of ~27.9 kJ/cm 2 . Electron temperatures and densities of the ablated plasma were determined by analyzing the time-resolved x-ray spectroscopy data through a comparison of experimental measurements with synthetic results obtained from Si atomic calculations in a steady state and non-local thermodynamic equilibrium. These calculations were computed using PrismSPECT. Additionally, radiation-hydrodynamics simulations with FLASH are used to generate simulated plasma-density and plasma-temperature profiles, which are then compared with the experimental measurements. Our analyses reveal that increasing the laser pulse length at a constant fluence results in a decrease in electron temperatures and densities. Furthermore, the longer pulses with lower intensities lead to deeper ablation regions before reaching the peak ablation but lower ionization balances in the silicon layer. Furthermore, these findings emphasize the critical role of laser pulse length in plasma ablation and shock generation for laser-impulse studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Computational investigation of water glasses using machine-learning potentials

The molecular origins of water’s anomalous properties have long been a subject of scientific inquiry. The liquid–liquid phase transition hypothesis, which posits the existence of distinct low-density and high-density liquid states separated by a first-order phase transition terminating at a critical point, has gained increasing experimental and computational support and offers a thermodynamically consistent framework for many of water’s anomalies. However, experimental challenges in avoiding crystallization near the postulated liquid–liquid critical point have focused attention to water’s canonical glassy states: low-density and high-density amorphous ice. Here, we use two Deep Potential machine-learning models, trained on the Strongly Constrained and Appropriately Normed density functional and the highly accurate Many-Body Polarizable potential, to conduct an investigation of water’s glassy phenomenology based on quantum mechanical calculations. Despite not being explicitly trained on amorphous ices, both models accurately capture the structure and transformation of the water glasses, including their interconversion along different thermodynamic paths. Isobaric quenching of liquid water at various pressures generates a continuum of intermediate amorphous ices and density fluctuations increase near the liquid–liquid critical pressure. The glass transition temperatures of the amorphous ices produced at different pressures exhibit two distinct branches, corresponding to low-density and high-density amorphous ice behaviors, consistent with experiment and the liquid–liquid transition hypothesis. Extrapolating transformation pressures from isothermal compressions to experimental compression rates brings our simulations into excellent agreement with data. Our findings demonstrate that machine-learning potentials trained on equilibrium phases can effectively model nonequilibrium glassy behavior and pave the way for studying long-timescale, out-of-equilibrium processes with quantum mechanical accuracy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Uranium Bed Design Parameters for Tritium Plants Supporting Fusion Reactors

Uranium hydride is commonly used to store hydrogen or its isotopes in a solid state. The Self-Assaying Tritium Accountancy and Containment Unit for ITER (STACI) is a 5.2-kg bed of depleted uranium (dU) capable of holding up to 33 mol of hydrogen or its isotopes. Here we summarize the data analysis of past experimental campaigns with STACI, with the aim of describing the kinetics and thermodynamics of the hydriding process. Computed tomography imaging was performed on STACI both before and after its experimental campaign, and a high degree of swelling was observed in the dU. Literature on studies in regard to the swelling of large (multikilogram) quantities of uranium hydride for storage applications was not identified during this study. Data from the experimental campaign, as well as data on the formation reaction, are presented. The authors hope to create an analytical model of STACI based on these data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Origin of mean-field behavior in an elastic Ising model

Simple elastic models of spin-crossover compounds are known empirically to exhibit classical critical behavior. Here, we demonstrate how the long-range interactions responsible for this behavior arise naturally upon integrating out mechanical fluctuations of such a model. A mean-field theory applied to the resulting effective Hamiltonian quantitatively accounts for both thermodynamics and kinetics observed in computer simulations, including a barrier to magnetization reversal that grows extensively with system size. For nanocrystals, which break translational symmetry, a straightforward extension of mean-field theory yields similarly accurate results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multimodal Ligand Non-Innocence for Reduction Catalysis (Final Technical Report)

The goal of this project was to develop multimodal ligand non-innocence in late transition metal complexes of aminophenolate ligands and to apply this ligand-centered reactivity to small-molecule reduction reactions relevant to energy-conversion chemistry. Specific aims were, 1) Prepare new coordination complexes of aminophenolate ligands and characterize them using a battery of structural, spectroscopic, and computational methods. 2) Benchmark the thermodynamics and kinetics of e - , H + , H • , and H - transfer reactivity to establish the baseline reactivity of new metal complexes containing aminophenolate ligands. 3) Apply new aminophenolate complexes to catalytic reduction reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Venus - Atmospheric evolution.

Thermodynamic equilibrium composition of Venus atmosphere computed from atomic composition deduced spectroscopically

SPECTROSCOPY↗

Skylab program CSM verification analysis report

The application of the SINDA computer program for the transient thermodynamic simulation of the Apollo fuel cell/radiator system for the limit condition of the proposed Skylab mission is described. Results are included for the thermal constraints imposed upon the Pratt and Whitney fuel cell power capability by the Block 2 EPS radiator system operating under the Skylab fixed attitude orbits.

Schaefer, J. L.↗

Program Gives Data On Physical Properties Of Hydrogen

TAB II computer program provides values of thermodynamic and transport properties of hydrogen in useful format. Also, provides values for equilibrium hydrogen and para-hydrogen. Program fast, moderately accurate, and operates over wide ranges of input variables. Written in FORTRAN 77.

Roder, H. M.↗

Three CLIPS-based expert systems for solving engineering problems

We have written three expert systems, using the CLIPS PC-based expert system shell. These three expert systems are rule based and are relatively small, with the largest containing slightly less than 200 rules. The first expert system is an expert assistant that was written to help users of the ASPEN computer code choose the proper thermodynamic package to use with their particular vapor-liquid equilibrium problem. The second expert system was designed to help petroleum engineers choose the proper enhanced oil recovery method to be used with a given reservoir. The effectiveness of each technique is highly dependent upon the reservoir conditions. The third expert system is a combination consultant and control system. This system was designed specifically for silicon carbide whisker growth. Silicon carbide whiskers are an extremely strong product used to make ceramic and metal composites. The manufacture of whiskers is a very complicated process. which to date. has defied a good mathematical model. The process was run by experts who had gained their expertise by trial and error. A system of rules was devised by these experts both for procedure setup and for the process control. In this paper we discuss the three problem areas of the design, development and evaluation of the CLIPS-based programs.

Parkinson, W. J.↗

Interpreting Assemblages with Titanite (Sphene): It Does not have to be Greek to You.

Assemblages with titanite, pyroxene(s), olivine, ilmenite, magnetite, and quartz can be used to constrain the intensive and compositional variables that operate during crystallization. Such assemblages are relatively rare in metamorphic rocks, but they are more common in igneous rocks and more frequently reported in plutonic than volcanic rocks. We used the program QUILF, enhanced with thermodynamic data for titanite, to compute stable reactions among titanite (CaTiOSiO4), Fe-Mg-Ti ilmenite and magnetite (hereafter ilmenite and magnetite), Ca-Mg-Fe pyroxenes and olivine, and quartz, and to evaluate some of the factors that control titanite stability. Calculations at 1, 3, and 6 Kbar and 650, 850, 1100 0 C, in the system CaO - MgO - FeO Fe2O3 - TiO2 - SiO2, suggest that the reactions: Augitc + Ilmenite = Titanite + Magnetite + Quartz and Augite + Ilmenite + Quartz = Titanite + Orthopyroxene, impose well defined fugacity of O2, alpha(sub SiO2), , and compositional restrictions to the assemblages: (1) Titanite + Magnetite + Quartz, (2) Titanite + Orthopyroxene, (3) Augite + Ilmenite, and consequently titanite stability. From our calculations in this system we can draw the following general conclusions: (1) The assemblage Titanite + Magnetite + Quartz is always a good indicator of relatively high fugacity of O2, and it is likely more common in relatively iron-rich bulk compositions and for decreasing temperature and pressure conditions. (2) At high temperatures (>= 650 C) titanite is not stable in quartz-saturated rocks that contain the assemblage Orthopyroxene + Augite + Ilmenite + Magnetite. (3) In quartz-saturated rocks the coexistence of titanite and magnetite with either orthopyroxene or olivine requires a confluence of conditions relating bulk composition, fugacity of O2, and slow cooling. Thus, such assemblages must be rare. (4) Regardless of T and fugacity of O2 conditions, and bulk-composition, titanite is not stable in quartz-absent rocks that contain Olivine + Orthopyroxene + Augite + Ilmenite + Magnetite. Decreasing temperature and pressure conditions appear to favor titanite crystallization, thus, it is not unsurprising that titanite is frequently observed in slowly cooled rocks, albeit, in association with amphibole. We argue that the titanite + amphibole association is likely favored by high water activity, regardless of oxygen fugacity. Because water activity increases during crystallization of a pluton, the association titanite + amphibole, and consequently titanite, is likely to be more common in plutonic rocks than in volcanic rocks.

Xirouchakis, Dimitrios M.↗