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 1,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Modeling the internal combustion engine

A flexible and computationally economical model of the internal combustion engine was developed for use on large digital computer systems. It is based on a system of ordinary differential equations for cylinder-averaged properties. The computer program is capable of multicycle calculations, with some parameters varying from cycle to cycle, and has restart capabilities. It can accommodate a broad spectrum of reactants, permits changes in physical properties, and offers a wide selection of alternative modeling functions without any reprogramming. It readily adapts to the amount of information available in a particular case because the model is in fact a hierarchy of five models. The models range from a simple model requiring only thermodynamic properties to a complex model demanding full combustion kinetics, transport properties, and poppet valve flow characteristics. Among its many features the model includes heat transfer, valve timing, supercharging, motoring, finite burning rates, cycle-to-cycle variations in air-fuel ratio, humid air, residual and recirculated exhaust gas, and full combustion kinetics.

Zeleznik, F. J.↗
NASA NTRS1990

Investigation of mechanical and thermal properties of microwave-sintered lunar simulant materials using 2.45 GHz radiation

The mechanical and thermal properties of lunar simulant material were investigated. An alternative method of examining thermal shock in microwave-sintered lunar samples was researched. A computer code was developed that models how the fracture toughness of a thermally shocked lunar simulant sample is related to the sample hardness as measured by a micro-hardness indentor apparatus. This technique enables much data to be gathered from a few samples. Several samples were sintered at different temperatures and for different times at the temperatures. The melting and recrystallization characteristics of a well-studied binary system were also investigated to see if the thermodynamic barrier for the nucleation of a crystalline phase may be affected by the presence of a microwave field. The system chosen was the albite (sodium alumino silicate) anorthite system (calcium alumino silicate). The results of these investigations are presented.

Meek, T. T.↗
NASA NTRS2015

Numerical Modeling of Self-Pressurization and Pressure Control by Thermodynamic Vent System in a Cryogenic Tank

This paper presents a numerical model of a system-level test bed - the multipurpose hydrogen test bed (MHTB) using Generalized Fluid System Simulation Program (GFSSP). MHTB is representative in size and shape of a fully integrated space transportation vehicle liquid hydrogen (LH2) propellant tank and was tested at Marshall Space Flight Center (MSFC) to generate data for cryogenic storage. GFSSP is a finite volume based network flow analysis software developed at MSFC and used for thermo-fluid analysis of propulsion systems. GFSSP has been used to model the self-pressurization and ullage pressure control by Thermodynamic Vent System (TVS). A TVS typically includes a Joule-Thompson (J-T) expansion device, a two-phase heat exchanger, and a mixing pump and spray to extract thermal energy from the tank without significant loss of liquid propellant. Two GFSSP models (Self-Pressurization & TVS) were separately developed and tested and then integrated to simulate the entire system. Self-Pressurization model consists of multiple ullage nodes, propellant node and solid nodes; it computes the heat transfer through Multi-Layer Insulation blankets and calculates heat and mass transfer between ullage and liquid propellant and ullage and tank wall. TVS model calculates the flow through J-T valve, heat exchanger and spray and vent systems. Two models are integrated by exchanging data through User Subroutines of both models. The integrated models results have been compared with MHTB test data of 50% fill level. Satisfactory comparison was observed between test and numerical predictions.

Majumdar, Alok↗
NASA NTRS2018

Tetrahedral-Mesh Simulations of Shock-Turbulence Interaction

Despite decades of development of unstructured mesh methods, direct numerical simulations (DNS) of turbulent flows are still predominantly performed on structured or unstructured hexahedral meshes with high-order finite-difference methods, weighted essentially nonoscillatory (WENO) schemes, or hybrid schemes formed by their combinations. Tetrahedral meshes offer easy mesh generation and adaptation around complex geometries and the potential of an orientation-free grid that would benefit the isotropic nature of small-scale dissipation, as well as the solution accuracy of intermediate scales. To advance the state of the art of unstructured-mesh simulation capabilities for shock/turbulence interaction, DNS using pure tetrahedral meshes are carried out with the space-time conservation element, solution element (CESE) method in this research. By its design, the CESE method is constructed based on a non-dissipative scheme and is a genuinely multidimensional numerical framework that is free from the use of an approximate Riemann-solver. The numerical framework also provides the ability to add numerical dissipation (the nondissipative scheme acts as the reference state like that of the reversible state in thermodynamics) when needed (with justification from mathematics/physics). The above-mentioned features along with the CESE method's consistent shock-capturing approach and strong enforcement of flux conservation in spacetime offers a novel method to accurately simulate turbulent flows and their interaction with shocks using tetrahedral meshes. Two canonical problems, namely, isotropic turbulence interaction with a normal shock and a Mach 2.9 turbulent boundary layer flow over a 24deg compression corner are investigated in this study. Computational results show reasonably good agreement with experimental data and results from structured-mesh, high-order simulations available in the literature. Successful validation of these canonical problems demonstrated here paves the way for future high-fidelity supersonic flow simulations involving complex-geometries.

Venkatachari, Balaji Shankar↗
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 OSTI2021

Shape-Selective Synthesis of Alkylcyclopentenyl Cations in Zeolites and Spectroscopic Distinction of Constitutional Isomers

Alkylcyclopentenyl cations belong to the long-lived intermediates that make up the “hydrocarbon pool” during the catalytic conversion of methanol on zeolites, and recent works show that such cations contribute to olefin and aromatics formation. From liquid phase chemistry, two types of alkylcyclopentenyl cations are known and distinguished by the substituent at the central carbon (C-2) of the allylic system: a more stable type with a methyl group at the C-2 and a less stable type with hydrogen at the C-2. The following three linked objectives are pursued: (i) IR spectroscopic distinction of the different substitution patterns at the allylic system of alkylcyclopentenyl cations, (ii) the role of the zeolite framework in determining the substitution pattern, and (iii) identification of alkylcyclopentenyl cations from precursors relevant to methanol-to-olefins conversion. UV–vis and IR spectroscopy are applied in situ to characterize alkylcyclopentenyl cations produced by adsorption of a pentaalkylcyclopentadiene or by adsorption and thermally induced cyclization of 2,6-dimethyl-2,4,6-octatriene. Prior knowledge of electronic spectra is combined with DFT-computed and experimental IR spectra to establish the frequency of the C–H vibration of the C-2 hydrogen-substituted type, and a characteristic red shift of the asymmetric allylic stretching vibration (Δν ≈ -20 to -30 cm –1 ) after replacing hydrogen by methyl at the central carbon of the allylic system. Although DFT demonstrates that both types of ions fit into medium- and large-pore zeolites and that the C-2 methyl-substituted type is thermodynamically favored even in the pores of the considered zeolites, formation of alkylcyclopentenyl cations is found by UV–vis and IR spectra to be shape-selective. The bulkier C-2 methyl-substituted type is detected in large-pore zeolites (MOR, BEA) and in the intersections of medium-pore zeolites (MFI), whereas in channels of medium size (TON), the less bulky C-2 hydrogen-substituted type is exclusively formed. Finally, the type of ion formed and its overall size are dictated by the zeolite framework and are independent of the precursor; the same type of alkylcyclopentenyl cation as found through cyclization of dimethyloctatriene could be generated from ethene.

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

The Martini 3 Lipidome: Expanded and Refined Parameters Improve Lipid Phase Behavior

Lipid membranes are central to cellular life. Complementing experiments, computational modeling has been essential in unraveling complex lipid-biomolecule interactions, crucial in both academia and industry. The Martini model, a coarse-grained force field for efficient molecular dynamics simulations, is widely used to study membrane phenomena but has faced limitations, particularly in capturing realistic lipid phase behavior. Here, we present refined Martini 3 lipid models with a mapping scheme that distinguishes lipid tails that differ by just two carbon atoms, enhancing the structural resolution and thermodynamic accuracy of model membrane systems including ternary mixtures. The expanded Martini lipid library includes thousands of models, enabling simulations of complex and biologically relevant systems. These advancements establish Martini as a robust platform for lipid-based simulations across diverse fields.

Lipids↗
DOE OSTI2021

Epitaxial nanotwinned metals and alloys: synthesis-twin structure–property relations

Extensive studies have been performed to understand the structure–property relationship in nanotwinned (NT) metals, where the average twin thickness or twin spacing is less than 100 nm. NT metals may have high strength, ductility, good electrical conductivity and corrosion resistance and have promising applications in microelectronics, structural materials and coatings for gas and oil industry. Most prior studies focus on the computational and experimental investigations of NT metals with low stacking fault energy (SFE) and often in polycrystalline forms. In this highlight, we review recent findings on epitaxial growth of NT metals and alloys with low and high SFEs deposited on single crystal substrates, and focus primarily on advancing our comprehension towards tailoring the microstructures and properties of NT metals via extrinsic factors, including the synthesis technique, parametric control and growth template, and intrinsic factors, such as SFE and thermodynamic, kinetic alloying effects. The review intends to illuminate how to synthesize and design epitaxial NT structures to achieve desirable properties.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Thermalization from quantum entanglement: Jet simulations in the massive Schwinger model

We investigate the emergence of thermalization in a quantum-field-theoretic model mimicking the production of jets in QCD: the massive lattice Schwinger model coupled to external sources. Specifically, we compute the expectation values of local operators as functions of time and compare them to their thermal counterparts, quantify the overlap between the evolving density matrix and the thermal one, and compare the dynamics of the energy-momentum tensor to predictions from relativistic hydrodynamics. Through these studies, we find that the system approaches thermalization at late times and elucidate the mechanisms by which quantum entanglement drives thermalization in closed field-theoretic systems. Our results show how thermodynamic behavior emerges in real time from unitary quantum dynamics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1999

Flow of GE90 Turbofan Engine Simulated

The objective of this task was to create and validate a three-dimensional model of the GE90 turbofan engine (General Electric) using the APNASA (average passage) flow code. This was a joint effort between GE Aircraft Engines and the NASA Lewis Research Center. The goal was to perform an aerodynamic analysis of the engine primary flow path, in under 24 hours of CPU time, on a parallel distributed workstation system. Enhancements were made to the APNASA Navier-Stokes code to make it faster and more robust and to allow for the analysis of more arbitrary geometry. The resulting simulation exploited the use of parallel computations by using two levels of parallelism, with extremely high efficiency.The primary flow path of the GE90 turbofan consists of a nacelle and inlet, 49 blade rows of turbomachinery, and an exhaust nozzle. Secondary flows entering and exiting the primary flow path-such as bleed, purge, and cooling flows-were modeled macroscopically as source terms to accurately simulate the engine. The information on these source terms came from detailed descriptions of the cooling flow and from thermodynamic cycle system simulations. These provided boundary condition data to the three-dimensional analysis. A simplified combustor was used to feed boundary conditions to the turbomachinery. Flow simulations of the fan, high-pressure compressor, and high- and low-pressure turbines were completed with the APNASA code.

Veres, Joseph P.↗
DOE OSTI2022

Investigation of Rare Earth Elements Mobility Using Molecular Dynamics

Rare Earth Elements (REE) are critical many applications such as electronics, permanent magnets, agricultural etc. Effective separation of REEs is a challenge due to their similarities in chemical and physical properties. Current separation processes are costly and environmentally detrimental. Therefore, research on REEs mobility in water-based solvents is necessary to design effective separation processes. In this research, the key kinetic and thermodynamic parameters, including diffusion coefficient and hydration shell are determined under the influence of external electric field. The diffusion coefficient of metal ions increases due to electric field in pure water.

99 GENERAL AND MISCELLANEOUS↗
NASA NTRS1974

Free energies of formation of WC and WzC and the thermodynamic properties of carbon in solid tungsten

The activity of carbon in the two-phase regions - W + WC and W + W2C was obtained from the carbon content of iron rods equilibrated with mixtures of metal plus carbide powders. From this activity data the standard free energies of formation of WC and W2C were calculated. The temperature of the invariant reaction W2C = W + WC was fixed at 1570 + or - 5K. Using available solubility data for C in solid W, the partial molar free energy of C in the dilute solid solution was also calculated. The heat of solution of C in W, and the excess entropy for the interstitial solid solution, were computed, assuming that the carbon atoms reside in the octahedral interstices of bcc W.

Gupta, D. K.↗
DOE OSTI2019

Model systems for screening and investigation of lithium metal electrode chemistry and dendrite formation

The use of lithium metal as an electrode for electrochemical energy storage will provide a significant impact on practical energy storage technology. Unfortunately, the use of lithium metal is plagued with challenging chemical problems. Specifically, the formation of a solid electrolyte interphase layer and the nucleation and growth of lithium dendrites: both must be addressed and controlled in order to achieve a practically useable pure lithium metal electrode. Currently sophisticated experimental techniques and computationally expensive simulations are being used to probe these problems but these methods are arduous and time consuming which delays timely evaluation and insight into the rapidly changing field of advanced energy storage. Here, we report the use of DFT simulations of lithium nanoclusters to investigate and explore lithium metal chemistry with inexpensive computational methods to gain greater insight into electrochemical reductions and the nucleation and growth of dendrites. DME, LiTFSI, and LiFSI reduction energetics and structures with electrode effects from lithium metal are reported providing better physical description of the absolute reduction potential characterization. The electronic structure of the lithium nanoclusters were used to investigate the nucleation and growth of lithium dendrites from an ab initio perspective. The results demonstrate that kinetic processes have more control over non uniform growth than thermodynamic processes. Based on this information, a non ab initio model was created in Matlab that shows the initial stages of dendrite nucleation considering approximately 2000 atoms.

25 ENERGY STORAGE↗
DOE OSTI2022

Quantum membrane phases in synthetic lattices of cold molecules or Rydberg atoms

Here, we calculate properties of dipolar interacting ultracold molecules or Rydberg atoms in a semisynthetic three-dimensional configuration—one synthetic dimension plus a two-dimensional real-space optical lattice or periodic microtrap array—using the stochastic Green's function quantum Monte Carlo method. Through a calculation of thermodynamic quantities and appropriate correlation functions, along with their finite-size scalings, we show that there is a second-order transition to a low-temperature phase in which two-dimensional “sheets” form in the synthetic dimension of internal rotational or electronic states of the molecules or Rydberg atoms, respectively. Simulations for different values of the interaction V, which acts between atoms or molecules that are adjacent both in real and synthetic space, allow us to compute a phase diagram. We find a finite-temperature transition at sufficiently large V as well as a quantum phase transition—a critical value V c below which the transition temperature vanishes.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1994

A thermodynamic theory of dynamic fragmentation

We present a theory of dynamic fragmentation of brittle materials based on thermodynamic arguments. We recover the expressions for average fragment size and number as originally derived by Grady. We extend the previous work by obtaining descriptions of fragment size distribution and compressibility change due to the fragmentation process. The size distribution is assumed to be proportional to the spectral power of the strain history and a sample distribution is presented for a fragmentation process corresponding to a constant rate strain history. The description of compressibility change should be useful in computational studies of fragmentation. These results should provide insight into the process of fragmentation of brittle materials from hypervelocity impact.

Ching H Yew↗
DOE OSTI2022

Tuning Enzyme Thermostability via Computationally Guided Covalent Stapling and Structural Basis of Enhanced Stabilization

Enhancing the thermostability of enzymes without impacting their catalytic function represents an important yet challenging goal in protein engineering and biocatalysis. We recently introduced a novel method for enzyme thermostabilization that relies on the computationally guided installation of genetically encoded thioether “staples” into a protein via cysteine alkylation with the noncanonical amino acid O-2-bromoethyl tyrosine (O2beY). Here, we demonstrate the functionality of an expanded set of electrophilic amino acids featuring chloroacetamido, acrylamido, and vinylsulfonamido side-chain groups for protein stapling using this strategy. Using a myoglobin-based cyclopropanase as a model enzyme, our studies show that covalent stapling with pchloroacetamido-phenylalanine (pCaaF) provides higher stapling efficiency and enhanced stability (thermodynamic and kinetic) compared to the other stapled variants and the parent protein. Interestingly, molecular simulations of conformational flexibility of the cross-links show that the pCaaF staple allows fewer energetically feasible conformers than the other staples, and this property may be a broader indicator of stability enhancement. Using this strategy, pCaaF-stapled variants with significantly enhanced stability against thermal denaturation (ΔT m ' = +27 °C) and temperature-induced heme loss (ΔT 50 = +30 °C) were obtained while maintaining high levels of catalytic activity and stereoselectivity. Crystallographic analyses of singly and doubly stapled variants provide key insights into the structural basis for stabilization, which includes both direct interactions of the staples with protein residues and indirect interactions through adjacent residues involved in heme binding. This work expands the toolbox of protein stapling strategies available for protein stabilization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1988

Hydrogen turbines for space power systems: A simplified axial flow gas turbine model

Hydrogen cooled, turbine powered space weapon systems require a relatively simple, but reasonably accurate hydrogen gas expansion turbine model. Such a simplified turbine model would require little computational time and allow incorporation into system level computer programs while providing reasonably accurate volume/mass estimates. This model would then allow optimization studies to be performed on multiparameter space power systems and provide improved turbine mass and size estimates for the various operating conditions (when compared to empirical and power law approaches). An axial flow gas expansion turbine model was developed for these reasons and is in use as a comparative bench mark in space power system studies at Sandia. The turbine model is based on fluid dynamic, thermodynamic, and material strength considerations, but is considered simplified because it does not account for design details such as boundary layer effects, shock waves, turbulence, stress concentrations, and seal leakage. Although the basic principles presented here apply to any gas or vapor axial flow turbine, hydrogen turbines are discussed because of their immense importance on space burst power platforms.

Hudson, Steven L.↗
NASA NTRSDate not supplied

Investigating Fifth Oxide Effect on CMXAS Glass Properties

Coatings for hot section aero turbine engine structural materials are subject to thermochemical degradation after interacting with calcium-magnesium-aluminosilicates (CMAS). Molten CMAS viscosity is the primary glass property of importance, correlating with coating infiltration and reaction depth. CMAS viscosity was measured using a spindle-dipped viscometer at temperatures between 1300-1550°C in 25°C increments. Experimental data were compared to Thermo-Calc and FactSage computational models, which showed good agreement with experimental trends. CMXAS compositions, where X denotes a fifth oxide, of Ca30.75-Mg9-X5-Al13-Si42.25 (single cation oxide mol%) for X = Fe2+, Fe3+, Gd3+, Hf4+, Na1+, Ti4+, Y3+, Yb3+, Zr4+ additions were used. Oxide precipitation was observed in X = Hf4+ , Zr4+ CMXAS specimens, indicating coating dissolution is an increased threat for T ≥ 1390°C. Molten CMAS and CMXAS viscosity measurements were shown to correlate with net cation field strength, offering an alternative in property prediction where thermodynamic data are unavailable.

Clark Avery Luckhardt↗