SEARCH · Search NASA

Results for “Enthalpy method”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS2000

Effects of Space Environment on Flow and Concentration During Directional Solidification

A study of directional solidification of a weak binary alloy (specifically, Bi - 1 at% Sn) based on the fixed grid single domain approach is being undertaken. The enthalpy method is used to solve for the temperature field over the computational domain including both the solid and liquid phases; latent heat evolution is treated with the aid of an effective specific heat coefficient. A source term accounting for the release of solute into the liquid during solidification has been incorporated into the solute transport equation. The vorticity-stream function formulation is used to describe thermosolutal convection in the liquid region. In this paper we numerically investigate the effects of g-jitter on directional solidification. A background gravity of 1 micro-g has been assumed, and new results for the effects of periodic disturbances over a range of amplitudes and frequencies on solute field and segregation have been presented.

Benjapiyaporn, C.↗

NASA NTRS2000

Effect of g-jitter on Directional Solidification of a Binary Alloy

A study of directional solidification of a weak binary alloy (specifically, Bi - 1 at% Sn) based on the fixed grid single domain approach is being undertaken. The enthalpy method is used to solve for the temperature field over the computational domain including both the solid and liquid phases; latent heat evolution is treated with the aid of an effective specific heat coefficient. A source term accounting for the release of solute into the liquid during solidification has been incorporated into the solute transport equation. The vorticity-stream function formulation is used to describe thermosolutal convection in the liquid region. In this paper we present a numerical simulation of g-jitter: the small, rapid fluctuations in gravitational acceleration which may be experienced in an orbiting space vehicle. A background gravity of 1 micro-g has been assumed, and new results for the effects of orientation angle of the periodic disturbances over a range of amplitudes and frequencies on solute field and segregation have been presented.

Santiviriyapanich, P.↗

DOE OSTI2022

Experimental measurements and mathematical modeling of cold plate for aviation thermal management

Herein, this study, which has been motivated by the recent applications of the cold plate device in aviation thermal management, reports on physics-based mathematical models derived from the conservation laws of mass, momentum, and energy, and empiricism-based models. One of the objectives of the present work is to report on an elaborate and successful experimental work carried out on an additively manufactured device for the purpose of rigorously validating the numerical predictions. The excellent agreement between the numerical predictions and measured performance provides much needed confidence in the implementation, in the software package, of the offset-strip fin passage correlations, as well as in the software implementation of user-defined wavy fin correlations for aerospace heat exchangers and cold plates operating with ram air at Reynolds numbers below 8000. The contributions of this work can also be found in the development of a new and accurate thermal-hydraulic analysis procedure, referred to in this paper as plate-fin analogy. Results from this procedure are compared with those from thermal resistance network. The comparative study in this paper of the bulk and discrete enthalpy flux method is also new, as is the relative assessment of four off-set strip fin thermal-hydraulic models.

42 ENGINEERING↗

DOE OSTI2023

A general non-Fourier Stefan problem formulation that accounts for memory effects

The Stefan problem is the classical model of a melting phase change. In heterogeneous systems, such phase changes can exhibit non-Fourier (anomalous) behaviors, where the advance of the melt interface does not follow the expected time scaling. These situations can be modeled by replacing the derivatives, in the governing partial differential equations, with fractional order derivatives. In particular, replacing the time derivatives leads to non-Fourier models that account for memory effects in the system. In this work, by using appropriate time convolution integrals, a general thermodynamic balance statement for melting phase problems, explicitly accounting for memory effects, is developed. From this balance, a general model formulation applicable to problems involving melting over a temperature range (i.e., a mushy region) is derived. A key component in this model is the representation of memory effects through the use of fractional derivative based constitutive models of the enthalpy and heat flux. Further, on shrinking the mushy region to a single isotherm, a general sharp interface melting model is obtained. Here, in contrast to the classic Stefan problem, the fractional derivatives induce a natural regularization, such that the constitutive models for enthalpy and heat flux are continuous at the melt interface; a result confirmed through numerical simulation. To further support the theoretical findings, a physical example of a non-Fourier Stefan problem is presented. Overall the development and results in this paper underscore the importance of explicitly relating the development of fractional calculus models to the appropriate thermodynamic balance statements.

42 ENGINEERING↗

DOE OSTI2024

Non-Electricity Based Renewable Fuels: Theory and Computation for Solar Thermochemical Hydrogen

Dominated by photovoltaics and wind, current renewable energy sources generate mostly electricity, but 80% of the global final energy consumption occurs in form of fuels. Therefore, direct solar fuel generation would be a major breakthrough for the energy transition. Solar thermochemical hydrogen (STCH) is one of the very few potential routes towards scalable renewable fuels, but currently suffers from lack of an oxide working material that could optimally perform energy conversion within the thermodynamic boundary conditions. Theory and computation can contribute in two distinct ways, through materials search and discovery, but also by providing detailed mechanistic models for specific systems so to advance our understanding of possible design strategies. To enable high-throughput materials screening, we developed a defect graph neural network (dGNN) machine learning approach,[1] which accelerates the prediction of defect formation energies by replacing the tedious density functional theory (DFT) supercell calculations for all possible defect sites. This approach enables high-throughput database screening of oxides, which was integrated with thermodynamic modeling to extract the reduction entropies as additional selection criterion for STCH. Once potential candidate materials are identified, detailed models can guide materials design by predicting performance characteristics. One challenge is to quantitatively predict thermochemical equilibria at high concentrations when the redox active defects start to interact with each other, thereby impeding the formation of additional defects. Introducing a model for the free energy of defect interaction, parametrized on the basis of DFT data, we simulated the complete STCH redox cycle for (Sr,Ce)MnO3 alloys, achieving near-quantitative agreement with experimental data.[2] The analysis of these simulations reveals how defect interactions diminish the reduction entropy and H2 yield, suggesting to include these interactions in design considerations. Finally, we revisit the popular van't Hoff method for analyzing reduction enthalpies and entropies. This method is not ideal, as it involves a temperature-dependent convolution of gas-phase and solid-state entropies, causing uncertainties in the same order of magnitude as the physical quantities of interest. To avoid this problem, we suggest a simple alternative approach which can be applied to experimental and simulated data alike.

first-principles calculations↗

DOE OSTI2021

Systems, methods and tools for subterranean electrochemical characterization and enthalpy measurement in geothermal reservoirs

The present disclosure is directed to systems, methods and tools that measure ionic concentrations and downhole enthalpy of a flowing geothermal fluid in real-time at high-temperature and pressure. The systems, methods and tools include measuring the concentration of selected naturally occurring ions found in the liquid phase of the geothermal fluid throughout the wellbore using novel electrochemical sensor technologies. The change in liquid-phase ion concentration will be used to calculate the proportion of liquid to steam and allow for accurate enthalpy measurements. The techniques and technologies described here can be applied to any application of electrochemical sensing in extreme environments.

Cieslewski, Grzegorz↗

DOE OSTI2025

Comparison of multi-stage air treatment process divided by the same temperature and enthalpy difference

The multi-stage air treatment system has been proposed recently, and lower grade chilled/hot water could be used and energy efficiency could be improved. However, it has not been studied which division method of air treatment processes has higher energy efficiency. In this study, the model to calculate the energy consumption of multi-stage air treatment process is introduced, and the effects of two division methods, i.e. multi-stage air treatment process divided by the same temperature difference (ST method) or same enthalpy difference (SE method) between inlet and outlet at each stage, under 9 different air inlet parameters in the 2-stage and 3-stage air treatment processes are analysed and compared. The results show that (1) the system energy consumption of the SE method is generally lower than that of the ST method; (2) there is generally a larger energy consumption reduction rate of SE method when the air relative humidity is 70% compared to relative humidity of 50% and 90%; (3) the difference between ST method and SE method is not great, so both methods can be used for the design of multi-stage treatment system although SE method is normally recommended.

Wang, Wentao↗

NASA NTRS1988

Enthalpies of a binary alloy during solidification

The purpose of the paper is to present a method of calculating the enthalpy of a dendritic alloy during solidification. The enthalpies of the dendritic solid and interdendritic liquid of alloys of the Pb-Sn system are evaluated, but the method could be applied to other binaries, as well. The enthalpies are consistent with a recent evaluation of the thermodynamics of Pb-Sn alloys and with the redistribution of solute in the same during dendritic solidification. Because of the heat of mixing in Pb-Sn alloys, the interdendritic liquid of hypoeutectic alloys (Pb-rich) of less than 50 wt pct Sn has enthalpies that increase as temperature decreases during solidification.

Poirier, D. R.↗

NASA NTRS1999

Mid-Field Sonic Boom Extrapolation Methodology

In the design cycle of low-boom airplanes, sonic boom prediction must be accurate and efficient. The classical linear method, Whitham's F-function theory, has been widely applied to predict sonic boom signatures. However, linear theory fails to capture the nonlinear effects created by large civil transport. Computational fluid dynamics (CFD) has been used successfully to predict sonic boom signals at the near and mid fields. Nevertheless, it is computationally expansive in airplane design runs. In the present study, the method of characteristics is used to predict sonic boom signals in an efficient fashion. The governing equations are the axisymmetric Euler's equations with constant enthalpy. Since the method solves Euler's equations, it captures more nonlinear effects than the classical Whitham's F-function technique. Furthermore, the method of characteristics is an efficient marching scheme for initial value problems. In this study, we will first review the current CFD extrapolation technique and the work previously done in sonic boom extrapolation. Then, we will introduce the governing equations and the method of characteristics. Finally, we will show that the present method yields the same accurate results as previous CFD techniques, but with higher efficiency.

Cheung, Samson↗

NASA NTRSDate not supplied

Using Thermally Imperfect and Ideal Gas Equations for Determining Conditions in Slotted Throat Wind Tunnels

The derivation and application of two back pressure controlling methods for setting the dynamic pressure, Mach number, and Reynolds number of subsonic wind tunnels with a slotted-throat test section are discussed. The first method uses a pressure difference to determine the tunnel dynamic pressure via the compressible Bernoulli equations and is used in the ideal gas path in the CFD code. The second method uses an enthalpy difference in the flow to determine the velocity in the test section and can be applied to multispecies, calorically imperfect and calorically perfect gas flow simulations. For a tunnel operating point close to standard atmosphere conditions, Mach number, pressure, and temperature along the tunnel centerline were very closely matched between the various thermodynamic equation assumption controller methods.

Ground Testing↗

DOE OSTI2024

Assessing Gas-Phase Ion Reactivity of 50 Elements with NO and the Direct Application for 239 Pu in Complex Matrices Using ICP-MS/MS

Understanding the reactivity of metal cations with various reaction gases in ICP-MS/MS is important to determine the best gas to use for a given analyte/interference pair. In this study, nitric oxide (NO) was investigated as the reaction gas following previous experimental designs. The reactions with 50 elements were investigated to examine periodic trends in reactivity, validate theoretical modeling of reaction enthalpies as a method to screen reactant gases, and provide a baseline for potential in-line gas separation methods. ICP-MS/MS studies involving actinides are typically limited to Th, U and Pu, with analyses of Np and Am rarely reported in the literature. To date, only two previous methods have investigated the use of NO in ICP-MS/MS analyses. To showcase the utility of NO, a method was developed to measure 239 Pu in the presence of environmental matrix constituent and other actinides, like what could be expected from post-detonation debris, with no chemical separation prior to analysis. 239 Pu + was reacted to form 239 Pu 16 O + , eliminating interferences derived from the sample matrix by measuring 239 Pu + intensity at m/z = 255 ( 239 Pu 16 O + ). To validate NO for 238 U 1 H + interference removal in environmental matrices, standard reference materials were diluted to 1 mg per gram of solution and spiked to 0.05 pg/g of 239 Pu and 1 µg/g 238 U (Pu/U = 5 x 10 -8 ). Measured 239 Pu concentrations were within 6% of the spiked value. In conclusion, these results demonstrate that reliable 239 Pu measurements can be made at levels relevant to nuclear forensics without the need for extensive chemical matrix separation prior to analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Influence of Carboxymethyl Cellulose as a Thickening Agent for Glauber’s Salt-Based Low Temperature PCM

This work is focused on a novel, promising low temperature phase change material (PCM), based on the eutectic Glauber’s salt composition. To allow phase transition within the refrigeration range of temperatures of +5 °C to +12 °C, combined with a high repeatability of melting–freezing processes, and minimized subcooling, the application of three variants of sodium carboxymethyl cellulose (Na-CMC) with distinct molecular weights (700,000, 250,000, and 90,000) is considered. The primary objective is to optimize the stabilization of this eutectic PCM formulation, while maintaining the desired enthalpy level. Preparation methods are refined to ensure repeatability in mixing components, thereby optimizing performance and stability. Additionally, the influence of Na-CMC molecular weight on stabilization is examined through differential scanning calorimetry (DSC), T-history, and rheology tests. The PCM formulation of interest builds upon prior research in which borax, ammonium chloride, and potassium chloride were used as additives to sodium sulfate decahydrate (Glauber’s salt), prioritizing environmentally responsible materials. The results reveal that CMC with molecular weights of 250 kg/mol and 90 kg/mol effectively stabilize the PCM without phase separation issues, slowing crystallization kinetics. Conversely, CMC of 700 kg/mol proved ineffective due to the disruption of gel formation at its low gel point, hindering higher concentrations. Calculations of ionic concentration indicate higher Na ion content in PCM stabilized with 90 kg/mol CMC, suggesting increased ionic interactions and gel strength. A tradeoff is discovered between the faster crystallization in lower molecular weight CMC and the higher concentration required, which increases the amount of inert material that does not participate in the phase transition. After thermal cycling, the best formulation had a latent heat of 130 J/g with no supercooling, demonstrating excellent performance. This work advances PCM’s reliability as a thermal energy storage solution for diverse applications and highlights the complex relationship between Na-CMC molecular weight and PCM stabilization.

25 ENERGY STORAGE↗

NASA NTRS1992

Finite element simulations involving simultaneous multiple interface fronts in phase change problems

The present paper describes the simulation of phase change problems involving simultaneous multiple interface fronts employing the finite element method. Much of the past investigations employing finite elements have been restricted to primarily a single phase change situation. The existence of more than one phase, that is, the presence of multiple phase fronts poses certain challenges and further complications. However, the results provide a very interesting thermal behavior for this class of problems. In this paper, attention is focused on fixed grid methods and the trapezoidal family of one-step methods using the enthalpy formulations. Illustrative examples which handle simultaneous multiple fronts in phase change problems are presented.

Ouyang, Tianhong↗

NASA NTRS2014

Experimental and Theoretical Study of Thermodynamics of the Reaction of Titania and Water at High Temperatures

The transpiration method was used to determine the volatility of titanium dioxide (TiO2) in water vapor-containing environments at temperatures between 1473 and 1673 K. Water contents ranged from 0 to 76 mole % in oxygen or argon carrier gases for 20 to 250 hr exposure times. Results indicate that oxygen is not a key contributor to volatilization and the primary reaction for volatilization in this temperature range is: TiO2(s) + H2O(g) = TiO(OH)2(g). Data were analyzed with both the second and third law methods to extract an enthalpy and entropy of formation. The geometry and vibrational frequencies of TiO(OH)2(g) were computed using B3LYP density functional theory, and the enthalpy of formation was computed using the coupled-cluster singles and doubles method with a perturbative correction for connected triple substitutions [CCSD(T)]. Thermal functions are calculated using both a structure with bent and linear hydroxyl groups. Calculated second and third heats show closer agreement with the linear hydroxyl group, suggesting more experimental and computational spectroscopic and structural work is needed on this system.

equilibrium↗

NASA NTRS1961

A Technique for Determining the Nozzle-Flow Properties of Air in an Equilibrium, Nonequilibrium, or Frozen State

One of the problems associated with the development of high-enthalpy hypersonic test facilities hos been the determination of nozzle flow properties. This note presents a technique for determining the flow properties of air in an equilibrium, non-equilibrium, or frozen state using two test section measurements. A knowledge of two stagnation properties is required before the technique can be applied. The stagnation pressure can be measured, and reference 1 offers a method for determining stagnation enthalpy. Results of the calculations using the method of reference 1 are presented in reference 2 in a chart (chart 21) which can be readily used to determine the stagnation enthalpy from measurements of stagnation pressure, mass flow of air through the tunnel, and throat diameter.

Hayman, Lovick O.↗

NASA NTRS2024

Defining the Operational Envelope for Air Flows in the Miniature Arc-Jet Research Chamber (mARC II)

The second-generation 30 kW miniature Arc-jet Research Chamber (mARC II) at NASA Ames Research Center produces high enthalpy flows relevant for entry systems ground testing. The mARC II facility has recently undergone upgrades, including the installation of a new vacuum system to address the issues preventing it from maintaining underexpanded flow under test conditions. In this work, we present data obtained from an Integrated Systems Testing campaign and provide an initial assessment of arc-jet performance following the upgrades. Air is used as the working gas for the standard mARC II arc-heater configuration with two constrictor disks. Seven runs were investigated for five test conditions to assess the lowest achievable stagnation point heat fluxes for air flow rates of 0.15 or 0.25 g s −1 . The heat flux was measured using a water-cooled Gardon gauge (Ø4.76 mm, 3/16" hemispherical) at 70 mm from the nozzle exit plane. The new vacuum system produced test box pressures in the medium (fine) vacuum range (~0.03 torr, 4 Pa) prior to gas addition and successfully maintained underexpanded flow after gas addition. The upgrade yielded a ~4X reduction in heat flux relative to the previous system for the same set test conditions. We report the lowest heat fluxes measured in mARC II to date, ranging from 26 to 81 W cm −2 , for sonic flow enthalpies of 4–14 MJ kg −1 . Bulk enthalpies estimated using an energy balance method (EB 2 ) are reported for the first time using mARC II.Initial data suggests EB 2 generally estimates lower enthalpies than sonic flow methods for mARC II. Lastly, laminar axisymmetric Navier–Stokes simulations were performed using the NASA DPLR code. Numerical heat flux results show good agreement with experiments at low arc powers (3–15% difference at the minimum set arc current), but discrepancy increases with arc power (49% difference at the maximum set arc current).

arc-jet↗

NASA NTRS2024

Defining the Operational Envelope for Air Flows in the Miniature Arc-Jet Research Chamber (mARC II)

Pitot tube↗

DOE OSTI2022

An evaluation for geometries, formation enthalpies, and dissociation energies of diatomic and triatomic (C, H, N, O), NO 3 , and HNO 3 molecules from the PAW DFT method with PBE and optB88-vdW functionals

The structural geometries, formation enthalpies, and dissociation energies of all diatomic and triatomic molecules consisting of the four basic elements C, H, N, and/or O are calculated using the projector augmented wave density functional theory (DFT) method with the Perdew–Burke–Ernzerhof and optB88-vdW exchange-correlation functionals. The calculations are also extended to two larger molecules NO 3 and HNO 3 , which consist of four and five atoms, respectively. In total, 82 molecules or isomers are considered in the calculations. The geometric parameters including 42 bond lengths and 15 bond angles of these molecules calculated using the planewave DFT method are highly satisfactory, relative to the available experimental data. The error analysis is also performed for 49 formation enthalpies and 138 dissociation energies (including 51 atomization energies as well as the corresponding bond dissociation energies). The results are also compared with the previous data from various atomic-orbital-based methods for molecules and from similar or different planewave DFT methods for various solids and other molecules. This provides an informative and instructive evaluation especially for calculating the large-size material systems containing these small molecules as well as for developing the DFT methods further.

74 ATOMIC AND MOLECULAR PHYSICS↗