Search NASA⌕ Search

SEARCH · Search NASA

Results for “Thermodynamic Modeling”

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 523 records · Page 29

Polynomial approximations of thermodynamic properties of arbitrary gas mixtures over wide pressure and density ranges

Computer programs for flow fields around planetary entry vehicles require real-gas equilibrium thermodynamic properties in a simple form which can be evaluated quickly. To fill this need, polynomial approximations were found for thermodynamic properties of air and model planetary atmospheres. A coefficient-averaging technique was used for curve fitting in lieu of the usual least-squares method. The polynomials consist of terms up to the ninth degree in each of two variables (essentially pressure and density) including all cross terms. Four of these polynomials can be joined to cover, for example, a range of about 1000 to 11000 K and 0.00001 to 1 atmosphere (1 atm = 1.0133 x 100,000 N/m sq) for a given thermodynamic property. Relative errors of less than 1 percent are found over most of the applicable range.

Allison, D. O.↗

A Generalized Fluid System Simulation Program to Model Flow Distribution in Fluid Networks

This paper describes a general purpose computer program for analyzing steady state and transient flow in a complex network. The program is capable of modeling phase changes, compressibility, mixture thermodynamics and external body forces such as gravity and centrifugal. The program's preprocessor allows the user to interactively develop a fluid network simulation consisting of nodes and branches. Mass, energy and specie conservation equations are solved at the nodes; the momentum conservation equations are solved in the branches. The program contains subroutines for computing "real fluid" thermodynamic and thermophysical properties for 33 fluids. The fluids are: helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, parahydrogen, water, kerosene (RP-1), isobutane, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, R-11, R-12, R-22, R-32, R-123, R-124, R-125, R-134A, R-152A, nitrogen trifluoride and ammonia. The program also provides the options of using any incompressible fluid with constant density and viscosity or ideal gas. Seventeen different resistance/source options are provided for modeling momentum sources or sinks in the branches. These options include: pipe flow, flow through a restriction, non-circular duct, pipe flow with entrance and/or exit losses, thin sharp orifice, thick orifice, square edge reduction, square edge expansion, rotating annular duct, rotating radial duct, labyrinth seal, parallel plates, common fittings and valves, pump characteristics, pump power, valve with a given loss coefficient, and a Joule-Thompson device. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. This paper also illustrates the application and verification of the code by comparison with Hardy Cross method for steady state flow and analytical solution for unsteady flow.

Majumdar, Alok↗

Laser Induced Aluminum Surface Breakdown Model

Laser powered propulsion systems involve complex fluid dynamics, thermodynamics and radiative transfer processes. Based on an unstructured grid, pressure-based computational aerothermodynamics; platform, several sub-models describing such underlying physics as laser ray tracing and focusing, thermal non-equilibrium, plasma radiation and air spark ignition have been developed. This proposed work shall extend the numerical platform and existing sub-models to include the aluminum wall surface Inverse Bremsstrahlung (IB) effect from which surface ablation and free-electron generation can be initiated without relying on the air spark ignition sub-model. The following tasks will be performed to accomplish the research objectives.

Chen, Yen-Sen↗

Approach to Modeling Boundary Layer Ingestion Using a Fully Coupled Propulsion-RANS Model

Airframe-propulsion integration concepts that use boundary layer ingestion have the potential to reduce aircraft fuel burn. One concept that has been recently explored is NASA's Starc-ABL aircraft configuration, which offers the potential for 12% mission fuel burn reduction by using a turbo-electric propulsion system with an aft-mounted electrically driven boundary layer ingestion propulsor. This large potential for improved performance motivates a more detailed study of the boundary layer ingestion propulsor design, but to date, analyses of boundary layer ingestion have used uncoupled methods. These methods account for only aerodynamic effects on the propulsion system or propulsion system effects on the aerodynamics, but not both simultaneously. This work presents a new approach for building fully coupled propulsive-aerodynamic models of boundary layer ingestion propulsion systems. A 1D thermodynamic cycle analysis is coupled to a RANS simulation to model the Starc-ABL aft propulsor at a cruise condition and the effects variation in propulsor design on performance are examined. The results indicates that both propulsion and aerodynamic effects contribute equally toward the overall performance and that the fully coupled model yields substantially different results compared to uncoupled. The most significant finding is that boundary layer ingestion, while offering substantial fuel burn savings, introduces throttle dependent aerodynamics effects that need to be accounted for. This work represents a first step toward the multidisciplinary design optimization of boundary layer ingestion propulsion systems.

Boundary Layer Ingestion↗

The 2006 Cape Canaveral Air Force Station Range Reference Atmosphere Model Validation Study and Sensitivity Analysis to the National Aeronautics and Space Administration's Space Shuttle

The 2006 Cape Canaveral Air Force Station (CCAFS) Range Reference Atmosphere (RRA) is a statistical model summarizing the wind and thermodynamic atmospheric variability from surface to 70 kin. Launches of the National Aeronautics and Space Administration's (NASA) Space Shuttle from Kennedy Space Center utilize CCAFS RRA data to evaluate environmental constraints on various aspects of the vehicle during ascent. An update to the CCAFS RRA was recently completed. As part of the update, a validation study on the 2006 version was conducted as well as a comparison analysis of the 2006 version to the existing CCAFS RRA database version 1983. Assessments to the Space Shuttle vehicle ascent profile characteristics were performed to determine impacts of the updated model to the vehicle performance. Details on the model updates and the vehicle sensitivity analyses with the update model are presented.

Burns, Lee↗

Enabling the Assimilation of CrIS Shortwave Infrared Observations in Global NWP at NOAA. Part I: Background and Methods

Radiance observations from Earth-observing satellites have a significant positive impact on numerical weather prediction (NWP) forecasts, but some spectral regions are not fully exploited. Observations from hyperspectral infrared (IR) sounders in the longwave region (650-1100 cm -1 ), for instance, are routinely assimilated in many NWP models, but observations in the shortwave region (2155-2550 cm -1 ) are not. Each of these regions provides information on the temperature structure of the atmosphere, but the shortwave IR (SWIR) region is considered challenging to assimilate due to noise equivalent delta temperature (NEDT) that is highly variable depending on scene brightness temperature and to phenomena that are difficult to model, like non-Local Thermodynamic Equilibrium (NLTE) and solar reflectance. With recent advances in small-satellite technology, SWIR temperature sounders may provide an agile and cost-effective complement to the current constellation of IR sounders. Therefore, a better understanding of the use and impact of SWIR observations in data assimilation for NWP is warranted. In part one of this study, as presented here, the amount of unique information (as determined by Empirical Orthogonal Decomposition (EOD)) made available to a data assimilation system by Cross-track Infrared Sounder (CrIS) SWIR observations is reviewed, recent advancements to the Community Radiative Transfer Model (CRTM) for the simulation of CrIS shortwave radiances are tested, and enhancements to NOAA’s Global Data Assimilation System (GDAS) for the assimilation of CrIS SWIR observations are implemented and evaluated. Part two of this study, which seeks to assess the value of assimilating shortwave IR observations in global NWP, is also introduced.

Erin Jones↗

Nozzle flow with vibrational nonequilibrium

This research concerns the modeling and numerical solutions of the coupled system of compressible Navier-Stokes equations in cylindrical coordinates under conditions of equilibrium and nonequilibrium thermodynamics. The problem considered was the modeling of a high temperature diatomic gas N2 flowing through a converging-diverging high expansion nozzle. The problem was modeled in two ways. The first model uses a single temperature with variable specific heats as functions of this temperature. For the second model we assume that the various degrees of freedom all have a Boltzmann distribution and that there is a continuous redistribution of energy among the various degrees of freedom as the gas passes through the nozzle. Each degree of freedom is assumed to have its own temperature and, consequently, each system state can be characterized by these temperatures. This suggests that formulation of a second model with a vibrational degree of freedom along with a rotational-translation degree of freedom, each degree of freedom having its own temperature. Initially the vibrational degree of freedom is excited by heating the gas to a high temperature. As the high temperature gas passes through the nozzle throat there is a sudden drop in temperature along with a relaxation time for the vibrational degree of freedom to achieve equilibrium with the rotational-translation degree of freedom. That is, we assume that the temperature change upon passing through the throat is so great that the changes in the vibrational degree of freedom occur at a much slower pace and consequently lags behind the rotational-translational energy changes. This lag results in a finite relaxation time. In this context the term nonequilibrium is used to denote the fact that the energy content of the various degrees of freedom are characterized by two temperatures. We neglect any chemical reactions which could also add nonequilibrium effects. We develop the energy equations for the nonequilibrium model from first principles. The resulting equations, which model the nozzle flow, can be expressed in various forms. In most forms the resulting equations are coupled systems of nonlinear partial differential equations subject to certain boundary conditions. To solve the resulting coupled system of nonlinear partial differential equations, several numerical techniques were investigated: (1) the explicit MacCormack method, (2) the explicit-implicit MacCormack method, (3) the method of operator splitting, (4) factorization schemes, and (5) the Steger-Warming scheme.

Heinbockel, J. H.↗

Thermodynamic effects on developed cavitation

The results of an investigation of thermodynamic effects are presented. Distributions of temperature and pressure in a developed cavity were measured for zero- and quarter-caliber ogives. A semiempirical entrainment theory was developed to correlate the measured temperature depression in the cavity. This theory correlates the maximum temperature depression expressed in dimensionless form as the Jakob number in terms of the dimensionless numbers of Nusselt, Reynolds, Froude, and Peclet, and dimensionless cavity length, L/D. The results show that in general, the temperature depression increases with L/D and temperature and the cavitation number based on measured cavity pressure is a function of L/D for a given model contour, independent of the thermodynamic effect.

Holl, J. W.↗

The 2006 Kennedy Space Center Range Reference Atmosphere Model Validation Study and Sensitivity Analysis to the Performance of the National Aeronautics and Space Administration's Space Shuttle Vehicle

The Kennedy Space Center (KSC) Range Reference Atmosphere (RRA) is a statistical model that summarizes wind and thermodynamic atmospheric variability from surface to 70 km. The National Aeronautics and Space Administration's (NASA) Space Shuttle program, which launches from KSC, utilizes the KSC RRA data to evaluate environmental constraints on various aspects of the vehicle during ascent. An update to the KSC RRA was recently completed. As part of the update, the Natural Environments Branch at NASA's Marshall Space Flight Center (MSFC) conducted a validation study and a comparison analysis to the existing KSC RRA database version 1983. Assessments to the Space Shuttle vehicle ascent profile characteristics were performed by JSC/Ascent Flight Design Division to determine impacts of the updated model to the vehicle performance. Details on the model updates and the vehicle sensitivity analyses with the update model are presented.

Burns, Lee↗

Loci/STREAM Sharp Interface Model Validation of Cryogenic Propellant Tank Self-Pressurization

This paper documents a model validation effort for the prediction of cryogenic propellant tank selfpressurization physics using the Loci/STREAM sharp interface CFD model developed at NASA MSFC. Cryogenic fluid management (CFM) applications for long duration spaceflight missions require accurate modeling of heat transfer and thermodynamics. Therefore, detailed validation of those physics is important to ensure simulation accuracy. In the past, the MSFC Propulsion Fluid Dynamics Branch has validated and used the Loci/STREAM Volume of Fluid (VOF) module extensively for propellant slosh applications, which are dominated by fluid motion. For cases where the interface gas/liquid interface is essentially static, the Loci/STREAM sharp interface was developed for more rapid and efficient CFM analysis. In the present study, Loci/STREAM sharp interface simulations were compared with the NASA K-site self-pressurization experiment of a flight weight, partially full, LH2 tank in normal gravity. This study investigated several parametric effects and made a number of important observations on the performance of the Loci/STREAM sharp interface method: - First, a homogeneous model was derived to verify the CFD model. The homogeneous model assumes both gas and liquid phases are uniform and are at saturation temperature. The predicted pressure rise rates at 3.5 W/m2 and 2.0 W/m2 with a 49% fill level are all lower than the experimental measurement. This is expected due to the observed temperature stratification near the gas-liquid interface during testing. - The k-w and k-e turbulence models predict comparable pressurization rates, which are much lower than the experiment measurement. On the other hand, the pressure values and pressurization rates from both models are almost the same as those of the homogeneous model. The current turbulence model is too diffusive and is not capable of capturing the thermal stratification in the liquid and ullage. - When using the laminar flow model, the pressurization rate is found to be very close to the experimental value when a quasi-steady flow field is allowed to develop by continuously venting the tank prior to pressurization as was done during testing. On the other hand, without initial venting, the rate is not constant. Additionally, the development time to reach a quasi-steady flow field differed from the experiment due to the lack of turbulent mixing away from the gas-liquid interface that is expected in reality. - Using the laminar flow model and appropriate initial conditioning, the pressurization rate is within 15% of the experimental value for 49% fill level at 3.5 W/m 2 and is within 18% for 2.0 W/m 2 heat loading. - Using the laminar flow model and appropriate initial conditioning, the pressurization rate is within 18% of the experimental value for a 29% fill level and within 15% for an 83% fill level.

H. Q. Yang↗

An analysis of the high-latitude thermospheric wind pattern calculated by a thermospheric general circulation model. I - Momentum forcing

A diagnostic processor (DP) was developed for analysis of hydrodynamic and thermodynamic processes predicted by the NCAR thermospheric general circulation model (TGCM). The TGCM contains a history file on the projected wind, temperature and composition fields at each grid point for each hour of universal time. The DP assimilates the history file plus ion drag tensors and drift velocities, specific heats, coefficients of viscosity, and thermal conductivity and calculates the individual forcing terms for the momentum and energy equations for a given altitude. Sample momentum forcings were calculated for high latitudes in the presence of forcing by solar radiation and magnetospheric convection with a 60 kV cross-tail potential, i.e., conditions on Oct. 21, 1981. It was found that ion drag and pressure forces balance out at F region heights where ion drift velocities are small. The magnetic polar cap/auroral zone boundary featured the largest residual force or net acceleration. Diurnal oscillations were detected in the thermospheric convection, and geostrophic balance was dominant in the E layer.

Killeen, T. L.↗

Interannual Variability in Surface LW Fluxes Over the Tropical Oceans As Seen in ISCCP-FD and GEWEX SRB Data Sets

One notable aspect of Earth s climate is that although the planet appears to be very close to radiative balance at top-of-atmosphere (TOA), the atmosphere itself and underlying surface are not. Profound exchanges of energy between the atmosphere and oceans, land and cryosphere occur over a range of time scales. Recent evidence from broadband satellite measurements suggests that even these TOA fluxes contain some detectable variations. Our ability to measure and reconstruct radiative fluxes at the surface and at the top of atmosphere is improving rapidly. In this work we will evaluate two recently released estimates of radiative fluxes, focusing primarily on surface estimates. The International Satellite Cloud Climatology Project FD radiative flux profiles are available from mid-1 983 to near present and have been constructed by driving the radiative transfer physics from the Goddard Institute for Space Studies (GISS) global model with ISCCP clouds and TOVS (TIROS Operational Vertical Sounder)thermodynamic profiles. Fu!l and clear sky SW and LW fluxes are produced. A similar product from the NASA/GEWEX Surface Radiation Budget Project using different radiative flux codes and thermodynamics from the NAS/Goddard Earth Observing System (GEOS-1) assimilation model makes a similar calculation of surface fluxes. However this data set currently extends only through 1995. Significant differences in both interannual variability as well as trends are found between among these data sets. For radiative fluxes these differences are traced to TOVS thermodynamic soundings used to drive the ISCCP-FD calculations. Errors in near surface temperature and precipitable water cascade into ISCCP upward and downward IR flux components, demonstrably affecting interannual variability. Revised estimates of clear-sky fluxes over ocean are made using statistical algorithms and water vapor from the (SSM/I) Special Sensor Microwave Imager. These calculations show strong near-surface water vapor feedback over the tropical oceans in association with SST changes. However, it is also shown that ISCCP longwave cloud forcing, common to both the ISCCP-FD and GEWEX SRB retrievals, is the main driver of a long-term decrease in net LW flux to the surface during the near-20 year period covered by these revised estimates.

Robertson, F. R.↗

An internal variable assisted formulation of moisture sorption models for polymers and composites

The purpose of the paper is to demonstrate that constitutive equations for moisture sorption and swelling can be obtained with an internal variable method which is rooted in the theory of irreversible thermodynamics. Realism is injected into the model through the concept of free volume. The power of the method is illustrated by the solution of the formidable nonlinear boundary value problem for diffusion with swelling-stress dependent diffusivity.

Hiel, C. C.↗

Conceptual Feasibility Study of the Hyperloop Vehicle for Next-Generation Transport

The Hyperloop concept is proposed as a faster, cheaper alternative to high-speed rail and traditional short-haul aircraft. It consists of a passenger pod traveling through a tube under light vacuum while being propelled and levitated by a combination of permanent and electro-magnets. The concept addresses NASA's research thrusts for growth in demand, sustainability, and technology convergence for high-speed transport. Hyperloop is a radical departure from other advanced aviation concepts, however it remains an aeronautics concept that tackles the same strategic goals of low-carbon propulsion and ultra-effcient vehicles. System feasibility was investigated by building a multidisciplinary vehicle sizing model that takes into account aerodynamic, thermodynamic, structures, electromagnetic, weight, and mission analyses. The sizing process emphasized the strong coupling between the two largest systems: the tube and the passenger pod. The model was then exercised to examine Hyperloop from a technical and cost perspective. The structural sizing analysis of the travel tube demonstrates potential for signi cant capital cost reductions by considering an underwater route. Examination of varying passenger capacity indicates that the system can be operated with a wide range of passenger loads without significant change in operating expenses. Lastly, a high-level sizing study simulated variations in tube area, pressure, pod speed, and passenger capacity showing that there is a tube pressure that minimizes operating energy usage. The value of this optimal tube pressure is highly sensitive to numerous design details. These combined estimates of energy consumption, passenger throughput, and mission analyses all support Hyperloop as a faster and cheaper alternative to short-haul flights. The tools and expertise used to quantify these results also demonstrate how traditional aerospace design methods can be leveraged to handle the complex and coupled design process. Much of the technology development required for the Hyperloop is shared with next-generation aircraft. Furthermore, the substantial public interest and active commercial development make it an ideal candidate as an aircraft technology driver and test bed.

Hyperloop Vehicle↗

Numerical Predictions of Damage and Failure in Carbon Fiber Reinforced Laminates Using a Thermodynamically-Based Work Potential Theory

A thermodynamically-based work potential theory for modeling progressive damage and failure in fiber-reinforced laminates is presented. The current, multiple-internal state variable (ISV) formulation, referred to as enhanced Schapery theory (EST), utilizes separate ISVs for modeling the effects of damage and failure. Consistent characteristic lengths are introduced into the formulation to govern the evolution of the failure ISVs. Using the stationarity of the total work potential with respect to each ISV, a set of thermodynamically consistent evolution equations for the ISVs are derived. The theory is implemented into a commercial finite element code. The model is verified against experimental results from two laminated, T800/3900-2 panels containing a central notch and different fiber-orientation stacking sequences. Global load versus displacement, global load versus local strain gage data, and macroscopic failure paths obtained from the models are compared against the experimental results.

composite materials↗

Electron Energy Partition Across Interplanetary Shocks. III. Analysis

An analysis of model fit results of 15,210 electron velocity distribution functions (VDFs), observed within ±2 hr of 52 interplanetary (IP) shocks by the Wind spacecraft near 1 au, is presented as the third and final part on electron VDFs near IP shocks. The core electrons and protons dominate in the magnitude and change in the partialto-total thermal pressure ratio, with the core electrons often gaining as much or more than the protons. Only a moderate positive correlation is observed between the electron temperature and the kinetic energy change across the shock, while weaker, if any, correlations were found with any other macroscopic shock parameter. No VDF parameter correlated with the shock normal angle. The electron VDF evolves from a narrowly peaked core with flaring suprathermal tails in the upstream to either a slightly hotter core with steeper tails or much hotter flattop core with even steeper tails downstream of the weaker and strongest shocks, respectively. Both quasi-static and fluctuating fields are examined as possible mechanisms modifying the VDF, but neither is sufficient alone. For instance, flattop VDFs can be generated by nonlinear ion acoustic wave stochastic acceleration (i.e., inelastic collisions), while other work suggested they result from the combination of quasi-static and fluctuating fields. This three-part study shows that not only are these systems not thermodynamic in nature; even kinetic models may require modification to include things like inelastic collision operators to properly model electron VDF evolution across shocks or in the solar wind.

Solar wind↗

The extratropical 40-day oscillation in the UCLA general circulation model. Part 1: Atmospheric angular momentum

Variations in atmospheric angular momentum (AAM) are examined in a three-year simulation of the large-scale atmosphere with perpetual January forcing. The simulation is performed with a version of the University of California at Los Angeles (UCLA) general circulation model that contains no tropical Madden-Julian Oscillation (MJO). In addition, the results of three shorter experiments with no topography are analyzed. The three-year standard topography run contains no significant intraseasonal AAM periodicity in the tropics, consistent with the lack of the MJO, but produces a robust, 42-day AAM oscillation in the Northern Hemisphere (NH) extratropics. The model tropics undergoes a barotropic, zonally symmetric oscillation, driven by an exchange of mass with the NH extratropics. No intraseasonal periodicity is found in the average tropical latent heating field, indicating that the model oscillation is dynamically rather than thermodynamically driven. The no-mountain runs fail to produce an intraseasonal AAM oscillation, consistent with a topographic origin for the NH extratropical oscillation in the standard model. The spatial patterns of the oscillation in the 500-mb height field, and the relationship of the extratropical oscillation to intraseasonal variations in the tropics, will be discussed in Part 2 of this study.

Marcus, S. L.↗

Phase relations in the Fe-Ni-Cr-S system and the sulfidation of an austenitic stainless steel

The stability fields of various sulfide phases that form on Fe-Cr, Fe-Ni, Ni-Cr and Fe-Cr-Ni alloys were developed as a function of temperature and the partial pressure of sulfur. The calculated stability fields in the ternary system were displayed on plots of log P sub S sub 2 versus the conjugate extensive variable which provides a better framework for following the sulfidation of Fe-Cr-Ni alloys at high temperatures. Experimental and estimated thermodynamic data were used in developing the sulfur potential diagrams. Current models and correlations were employed to estimate the unknown thermodynamic behavior of solid solutions of sulfides and to supplement the incomplete phase diagram data of geophysical literature. These constructed stability field diagrams were in excellent agreement with the sulfide phases and compositions determined during a sulfidation experiment.

Jacob, K. T.↗