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,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Heat transfer in rocket combustion chambers

Complexities of liquid rocket engine heat transfer which involve the injector faceplate and film cooled walls are being investigated by computational analysis. A conjugate heat transfer analysis was used to describe localized heating phenomena associated with particular injector configurations and film coolant flows. These components were analyzed, and the analyses verified when appropriate test data were available. The component analyses are being synthesized into an overall flowfield/heat transfer model. A Navier-Stokes flow solver, the FDNS code, was used to make the analyses. Particular attention was given to the representation of the thermodynamic properties of the fluid streams. Unit flow models of specific coaxial injector elements have been developed and are being used to describe the flame structure near the injector faceplate.

Anderson, P.↗
NASA NTRS2005

Results of Aero/Acoustic Tests and Analytical Studies of a Two-Dimensional Eight-Lobe Mixer-Ejector Exhaust Nozzle at Takeoff Conditions

Hot flow aero-acoustic tests were conducted with Pratt & Whitney's High-Speed Civil Transport (HSCT) Mixer-Ejector Exhaust Nozzles by General Electric Aircraft Engines (GEAE) in the GEAE Anechoic Freejet Noise Facility (Cell 41) located in Evendale, Ohio. The tests evaluated the impact of various geometric and design parameters on the noise generated by a two-dimensional (2-D) shrouded, 8-lobed, mixer-ejector exhaust nozzle. The shrouded mixer-ejector provides noise suppression by mixing relatively low energy ambient air with the hot, high-speed primary exhaust jet. Additional attenuation was obtained by lining the shroud internal walls with acoustic panels, which absorb acoustic energy generated during the mixing process. Two mixer designs were investigated, the high mixing "vortical" and aligned flow "axial", along with variations in the shroud internal mixing area ratios and shroud length. The shrouds were tested as hardwall or lined with acoustic panels packed with a bulk absorber. A total of 21 model configurations at 1:11.47 scale were tested. The models were tested over a range of primary nozzle pressure ratios and primary exhaust temperatures representative of typical HSCT aero thermodynamic cycles. Static as well as flight simulated data were acquired during testing. A round convergent unshrouded nozzle was tested to provide an acoustic baseline for comparison to the test configurations. Comparisons were made to previous test results obtained with this hardware at NASA Glenn's 9- by 15-foot low-speed wind tunnel (LSWT). Laser velocimetry was used to investigate external as well as ejector internal velocity profiles for comparison to computational predictions. Ejector interior wall static pressure data were also obtained. A significant reduction in exhaust system noise was demonstrated with the 2-D shrouded nozzle designs.

Harrington, Douglas↗
NASA NTRS2019

Thermodynamic Analysis of Boundary Layer Species During Ablation

The focus of this work is on carbon clusters, hydrocarbons and other carbonaceous species which are constituents of pyrolysis gases injected into the boundary layer. Some of these molecules are observed to have absorption spectra in the VUV and UV region that match the emission spectra of atomic nitrogen and oxygen. Hence, they can potentially absorb the radiation impinging on the heat shield. Accurate thermodynamic data is not available for many of these species over the relevant temperature range for reentry. The objective of this work is to determine the thermochemical properties of potential radiation absorbing molecules using accurate ab initio quantum chemistry calculations. We have compiled a database of formation enthalpies, anharmonic vibrational frequencies and rotational constants for the ground and low-lying excited electronic states of these carbonaceous species. Using these thermochemical properties, we compute equilibrium mole fractions and observe that the mole fractions of some species of interest are significant, thus they could affect the radiative heat flux. Furthermore, comparisons of the formation enthalpies and mole fractions are made with data from literature.

heatshields↗
DOE OSTI2023

A multiscale phase field fracture approach based on the non-affine microsphere model for rubber-like materials

Rubber-like materials have a broad scope of applications due to their unique properties like high stretchability and increased toughness. Hence, computational models for simulating their fracture behavior are paramount for designing them against failures. In this study, the phase field fracture approach is integrated with a multiscale polymer model for predicting the fracture behavior in elastomers. At the microscale, damaged polymer chains are modeled to be made up of a number of elastic chain segments pinned together. Using the phase field approach, the damage in the chains is represented using a continuous variable. Both the bond stretch internal energy and the entropic free energy of the chain are assumed to drive the damage, and the advantages of this assumption are expounded. A framework for utilizing the non-affine microsphere model for damaged systems is proposed here by considering the minimization of a hypothetical undamaged free energy, ultimately connecting the chain stretch to the macroscale deformation gradient. At the macroscale, a thermodynamically consistent formulation is derived in which the total dissipation is assumed to be mainly due to the rupture of molecular bonds. Using a monolithic scheme, the proposed model is numerically implemented and the resulting three-dimensional simulation predictions are compared with existing experimental data. The capability of the model to qualitatively predict the propagation of complex crack paths and quantitatively estimate the overall fracture behavior is verified. Additionally, the effect of the length scale parameter on the predicted fracture behavior is studied for an inhomogeneous system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A generalized and adaptable tensor-contraction-based cluster expansion formalism for multicomponent solids

Density functional theory (DFT)-based simulations of materials have first-principles accuracy, but are very computationally expensive. For simulating various properties of multi-component alloys, the cluster expansion (CE) technique has served as the standard workaround to improve computational efficiency. However, the standard CE technique is difficult to extend to exotic and/or low-symmetry lattices, often implemented via iteration over particular cluster types, which must be enumerated per lattice structure. In this work, we introduce the tensor cluster expansion (TCE), implemented in the open-source code tce-lib, which maps correlation functions to mixed tensor contractions, eliminating the need to iterate over cluster types and additionally making the calculation of correlation functions well-suited for massively parallel architectures like GPUs. We show that local interaction energies are an immediate consequence of the TCE formalism, yielding nearly $\mathcal{O}$(1) energy difference calculations. We then use this formalism to fit CE models for the TaW and CoNiCrFeMn systems, and use these models to respectively compute the enthalpy of mixing curve and Cowley short-range order parameters, showing excellent agreement with ground truth data.

Cluster expansion↗
NASA NTRS2026

Optimization of a 10–Step Reduced CH₄/O₂ Combustion Mechanism for RDRE Applications

A 10-step, 9-species reduced mechanism for methane/oxygen combustion is evaluated and optimized for use in rotating detonation rocket engine (RDRE) simulations at an operating pressure of 10 bar. The mechanism is assessed against GRI-Mech 3.0 for ignition delay times, adiabatic equilibrium temperature and species, Chapman-Jouguet (CJ) detonation parameters, and laminar flame speed. A sequential Nelder–Mead optimization procedure is applied to the rate constants of reactions R1, R4, and R5, targeting ignition delay agreement over the equivalence ratio range φ = 1.0–1.4 and temperature range T = 1400–2000 K. The optimized mechanism achieves ignition delay agreement within 28% of GRI-Mech 3.0 for the temperature range relevant to RDRE operation (τ ≥ 5 μs). Adiabatic equilibrium temperatures and CJ detonation parameters are reproduced within 0.15% of GRI-Mech 3.0 values, confirming thermodynamic consistency of the mechanism. Laminar flame speed predictions deviate significantly from GRI-Mech 3.0, which is attributed to the absence of HO₂ chemistry in the reduced mechanism; this limitation does not affect RDRE simulation fidelity since laminar flame propagation is of secondary importance in RDREs.

Computational Fluid Dynamics↗
DOE OSTI2021

Chiral hydrodynamics in strong external magnetic fields

We construct the general hydrodynamic description of (3+1)-dimensional chiral charged (quantum) fluids subject to a strong external magnetic field with effective field theory methods. We determine the constitutive equations for the energy-momentum tensor and the axial charge current, in part from a generating functional. Furthermore, we derive the Kubo formulas which relate two-point functions of the energy-momentum tensor and charge current to 27 transport coefficients: 8 independent thermodynamic, 4 independent non-dissipative hydrodynamic, and 10 independent dissipative hydrodynamic transport coefficients. Five Onsager relations render 5 more transport coefficients dependent. We uncover four novel transport effects, which are encoded in what we call the shear-induced conductivity, the two expansion-induced longitudinal conductivities and the shear-induced Hall conductivity. Remarkably, the shear-induced Hall conductivity constitutes a novel non-dissipative transport effect. As a demonstration, we compute all transport coefficients explicitly in a strongly coupled quantum fluid via holography.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A thermodynamically consistent finite strain phase field approach to ductile fracture considering multi-axial stress states

Phase field models for ductile fracture have gained significant attention in the last two decades due to their ability in implicitly tracking the nucleation and propagation of cracks. However, most crack phase field formulations for elastoplastic solids focus only on the effects of plastic deformation, and do not consider the different multi-axial stress states that may arise in practical designs. In this work, a thermodynamically consistent phase field approach coupled with finite strain plasticity, considering multi-axial stress states is presented. In order to account for the coupling between plasticity and stress states, the Stress-Weighted Ductile Fracture Model (SWDFM) is utilized. The SWDFM represents a criterion for predicting ductile crack initiation under both monotonic and cyclic loadings based on histories of an internal plastic variable, stress triaxiality, and the Lode angle parameter. The excellent performance of the SWDFM for predicting ductile crack initiation motivates for its incorporation into a phase field approach for predicting both crack initiation and propagation through degradation of the fracture toughness. Moreover, based on the second law of thermodynamics, exact requirements are imposed on the rate at which the fracture toughness can evolve. Here, a novel function for degrading the plastic yield surface during the evolution of damage is introduced. This function, in line with experimental observations, leads to an accumulation of plastic deformation in damaged regions of a solid, and avoids numerical instabilities arising from concentrations of large plastic deformations in severely damaged regions. For validating the proposed model, results of computational simulations are compared to data from selected tests considering different multi-axial stress states. Comparisons of the numerical results with data from laboratory experiments demonstrate the capabilities of the proposed framework.

42 ENGINEERING↗
DOE OSTI2021

Adsorbate Free Energies from DFT-Derived Translational Energy Landscapes

Adsorption free energies are fundamental to surface chemistry and catalysis. Standard models combine some assumed analytical form of the translational potential energy surface, often parameterized against density functional theory (DFT) calculations, with an analytical expression for the resultant translational densities of states (DOS), free energy, and entropy. Here we compare the performance of such models against numerical evaluations of the DOS and thermodynamic functions derived from solutions to the translational Schrodinger equation. We compare results for a translational potential energy surface (PES) derived from nudged eleastic band calculations with those obtained from adsorbate rastering across a series of monatomic (O, S, C, N, H) and polyatomic (NH x ) adsorbates on (100) Pt and Au facets. We find that analytical models as commonly parameterized have mixed performance for describing the translational PES and that the consequences for computed free energies are modest but potentially significant in microkinetic models. Here, numerical solutions are possible for modest to no additional computational cost over analytical models and thus should be considered when reliable free energy estimates are needed or translational potential energy surfaces are available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

The R -Process Alliance: Exploring the cosmic scatter among ten r -process sites with stellar abundances

Context. The astrophysical origin of the rapid neutron-capture process (r-process), responsible for producing roughly half of the elements heavier than iron, remains uncertain. Detailed chemical signatures from the oldest, most metal-poor stars, which act as fossil records of the earliest nucleosynthesis events, can be used to identify the dominant r-process sites. Aims. We present a homogeneous chemical abundance analysis of ten r-process element-enhanced stars. These old and metal-poor stars are strongly enriched in r-process elements with minimal contamination from other nucleosynthetic sources. By focusing on this chemically pure sample, we aim to investigate intrinsic variations in the r-process abundance patterns and explore their implications for the nature and potential diversity of r-process sites. Methods. We performed a detailed chemical abundance analysis of high-resolution, high-signal-to-noise spectra. For each star, we inspected over 1400 individual absorption lines using a combination of equivalent width measurements and spectral synthesis. The analysis was conducted under the assumption of 1D local thermodynamic equilibrium and employing the MOOG radiative transfer code. Results. We derived abundances for 54 chemical species, including 29 neutron-capture (n-capture) elements, covering the full mass range of the r-process abundance pattern. A kinematic analysis reveals that stars likely originated from ten kinematically distinct systems. Based on this assumption, we used the sample to probe the maximum variation expected from ten independent r-process nucleosynthesis events and computed the intrinsic dispersion of each element relative to Zr and Eu for the light and heavy r-process elements, respectively. This exercise resulted in a remarkably low cosmic scatter across the ten r-process sites enriching these stars; for the rare earth and third peak elements, for example, we find σ [La/Eu] = 0.08 and σ [Os/Eu] = 0.11 dex, while the scatter between light and heavy elements, σ [Zr/Eu] , is slightly higher at 0.18 dex. Conclusions. The elemental abundance patterns across the ten independent r-process sites show remarkably small cosmic dispersions. This minimal dispersion suggests a high degree of uniformity in r-process yields across diverse astrophysical environments.

Astronomy and AstroPhysics↗
DOE OSTI2024

Final Scientific/Technical Report Rapid Design and Manufacturing of High-Performance Materials for Turbine Blades

This research demonstrated the concept of carbide precipitation-strengthened refractory high entropy alloys (RHEA). The advantage of a precipitation strengthened alloy is all phases are in thermodynamic equilibrium promoting microstructural stability, and consequently retention of properties at elevated temperatures. Additionally, as with any precipitation strengthened – or age hardened alloy, components can be heat-treated after manufacturing to manipulate the microstructure and optimize properties for performance. This is an advantage of the precipitation strengthened alloys over composites and ceramics materials, where the microstructure and properties are essentially fixed upon the initial materials manufacturing stage. High throughput (HT), multi-scale computer modeling was used to identifying novel RHEA compositions with desired characteristics needed for precipitation strengthening. The results showed that carbides precipitated and consequently, the strength of the alloys (measured in compression) increased after heat-treatment, which is the desired effect. The project also demonstrated the feasibility of producing articles from RHEA by additive manufacturing (AM).

36 MATERIALS SCIENCE↗
DOE OSTI2024

Rapid Design and Manufacturing of High-Performance Materials for Turbine Blades

This research demonstrated the concept of carbide precipitation-strengthened refractory high entropy alloys (RHEA). The advantage of a precipitation strengthened alloy is all phases are in thermodynamic equilibrium promoting microstructural stability, and consequently retention of properties at elevated temperatures. Additionally, as with any precipitation strengthened – or age hardened - alloy, components can be heat-treated after manufacturing to manipulate the microstructure and optimize properties for performance. This is an advantage of the precipitation strengthened alloys over composites and ceramics materials, where the microstructure and properties are essentially fixed upon the initial materials manufacturing stage. High throughput (HT), multi-scale computer modeling was used to identifying novel RHEA compositions with desired characteristics needed for precipitation strengthening. Designs were validated by producing small ingots via are melting. The results showed that carbides precipitated and consequently, the strength of the alloys (measured in compression) increased after heat-treatment, which is the desired effect. The project also demonstrated the feasibility of producing articles from RHEA by additive manufacturing (AM). Electron beam melting (EBM) and laser direct energy deposition (L-DED) additive manufacturing methods, were explored with various processing parameters were interrogated for both methods. Sound (crack-free and dense) precipitation strengthened RHEA samples were produced via EDM AM, demonstrating the feasibility of the concept.

36 MATERIALS SCIENCE↗
NASA NTRS1980

The multispecies modeling of the premixed, laminar steady-state ozone flame

Species dependent kinetic, transport and thermodynamic coefficients were employed in a one dimensional model of the premixed, laminar, steady state ozone flame. Convenient expressions for these coefficients are reported. They are based on independent measurements, no arbitrary parameters are used. The governing equations are solved using a relaxation technique and the partial differential equation package, PDECOL. Species and temperature profiles and the burning velocities are found over the range of initial ozone mole fraction of 0.25 to 1.00. The computed burning velocities are no more than 30% greater than the measurements of Streng and Grosses. Comparison with the computed results of Warnatz shows agreement within + or - 12%, even though quite different expressions for some of the kinetic coefficients were used. These differences are most obvious in the atomic oxygen and temperature profiles at an initial ozone mole fraction of unity.

Heimerl, J. M.↗
NASA NTRS1987

A computational study of the flowfield surrounding the Aeroassist Flight Experiment vehicle

A symmetric total variation diminishing (STVD) algorithm has been applied to the solution of the three-dimensional hypersonic flowfield surrounding the Aeroassist Flight Experiment (AFE) vehicle. Both perfect-gas and chemical nonequilibrium models have been used. The perfect-gas flows were computed at two different Reynolds numbers, including a flight trajectory point at maximum dynamic pressure, and on two different grids. Procedures for coupling the solution of the species continuity equations with the Navier-Stokes equations in the presence of chemical nonequilibrium are reviewed and tested on the forebody of the AFE and on the complete flowfield assuming noncatalytic wall and no species diffusion. Problems with the STVD algorithm unique to flows with variable thermodynamic properties (real gas) are identified and algorithm modifications are suggested. A potential heating problem caused by strong flow impingement on the nozzle lip in the near wake at 0-deg angle of attack has been identified.

Gnoffo, Peter A.↗
NASA NTRS1992

Heat transfer in rocket engine combustion chambers and nozzles

The complexities of liquid rocket engine heat transfer which involve the injector faceplate and regeneratively and film cooled walls are being investigated by computational analysis. A conjugate heat transfer analysis will be used to describe localized heating phenomena associated with particular injector configurations and coolant channels and film coolant dumps. These components are being analyzed, and the analysis verified with appropriate test data. Finally, the component analysis will be synthesized into an overall flowfield/heat transfer model. The FDNS code is being used to make the component analyses. Particular attention is being given to the representation of the thermodynamic properties of the fluid streams and to the method of combining the detailed models to represent overall heating. Unit flow models of specific coaxial injector elements have been developed and will be described. Film cooling simulations of film coolant flows typical of the subscale Space Transportation Main Engine (STME) being experimentally studied by Pratt and Whitney have been made, and these results will be presented. Other film coolant experiments have also been simulated to verify the CFD heat transfer model being developed. The status of the study and its relevance as a new design tool are covered. Information is given in viewgraph form.

Anderson, P. G.↗
NASA NTRS1993

An upwind, kinetic flux-vector splitting method for flows in chemical and thermal non-equilibrium

We have developed new upwind kinetic difference schemes for flows with non-equilibrium thermodynamics and chemistry. These schemes are derived from the Boltzmann equation with the resulting Euler schemes developed as moments of the discretized Boltzmann scheme with a locally Maxwellian velocity distribution. Splitting the velocity distribution at the Boltzmann level is seen to result in a flux-split Euler scheme and is called Kinetic Flux Vector Splitting (KFVS). Extensions to flows with finite-rate chemistry and vibrational relaxation is accomplished utilizing nonequilibrium kinetic theory. Computational examples are presented comparing KFVS with the schemes of Van Leer and Roe for a quasi-one-dimensional flow through a supersonic diffuser, inviscid flow through two-dimensional inlet, and viscous flow over a cone at zero angle-of-attack. Calculations are also shown for the transonic flow over a bump in a channel and the transonic flow over an NACA 0012 airfoil. The results show that even though the KFVS scheme is a Riemann solver at the kinetic level, its behavior at the Euler level is more similar to the existing flux-vector splitting algorithms than to the flux-difference splitting scheme of Roe.

Eppard, W. M.↗
NASA NTRS1993

Heat transfer in rocket engine combustion chambers and nozzles

Complexities of liquid rocket engine heat transfer which involve the injector faceplate and regeneratively and film cooled walls are being investigated by computational analysis. A conjugate heat transfer analysis will be used to describe localized heating phenomena associated with particular injector configurations and coolant channels and film coolant dumps. These components are being analyzed, and the analyses verified with appropriate test data. Finally, the component analyses will be synthesized into an overall flowfield/heat transfer model. The FDNS code is being used to make the component analyses. Particular attention is being given to the representation of the thermodynamic properties of the fluid streams and to the method of combining the detailed models to represent overall heating. Unit flow models of specific coaxial injector elements have been developed and will be described. Since test data from the NLS development program are not available, new validation heat transfer data have been sought. Suitable data were obtained from a Rocketdyne test program on a model hydrocarbon/oxygen engine. Simulations of these test data will be presented. Recent interest in the hybrid motor have established the need for analyses of ablating solid fuels in the combustion chamber. Analysis of a simplified hybrid motor will also be presented.

Anderson, P. G.↗