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 433 records · Page 24

Shock-tube pyrolysis of acetylene - Sensitivity analysis of the reaction mechanism for soot formation

The impact of thermodynamic parameters on the sensitivity of model predictions of soot formation by shock-tube pyrolysis of acetylene were assessed analytically. The pyrolysis process was treated as having three components: initiation, the initial pyrolysis stages; cyclization, formation of larger molecules and radicals and small aromatic molecules; and polymerization, further growth of aromatic rings. Rate equations are reviewed for each component. Thermodynamic effects were assessed by varying the C2H-H and C2H3-H bond energies and the Ct-(Ct) group additivity value. Any change in the C2H-H bond energy had a significant impact on the temperature and the maximum amount of the soot yield. The findings underscore the necessity of using accurate thermodynamic data for modeling high-temperature chemical kinetics.

Frenklach, M.↗

Non-LTE line-blanketed model atmospheres of hot stars. 1: Hybrid complete linearization/accelerated lambda iteration method

A new munerical method for computing non-Local Thermodynamic Equilibrium (non-LTE) model stellar atmospheres is presented. The method, called the hybird complete linearization/accelerated lambda iretation (CL/ALI) method, combines advantages of both its constituents. Its rate of convergence is virtually as high as for the standard CL method, while the computer time per iteration is almost as low as for the standard ALI method. The method is formulated as the standard complete lineariation, the only difference being that the radiation intensity at selected frequency points is not explicity linearized; instead, it is treated by means of the ALI approach. The scheme offers a wide spectrum of options, ranging from the full CL to the full ALI method. We deonstrate that the method works optimally if the majority of frequency points are treated in the ALI mode, while the radiation intensity at a few (typically two to 30) frequency points is explicity linearized. We show how this method can be applied to calculate metal line-blanketed non-LTE model atmospheres, by using the idea of 'superlevels' and 'superlines' introduced originally by Anderson (1989). We calculate several illustrative models taking into accont several tens of thosands of lines of Fe III to Fe IV and show that the hybrid CL/ALI method provides a robust method for calculating non-LTE line-blanketed model atmospheres for a wide range of stellar parameters. The results for individual stellar types will be presented in subsequent papers in this series.

Hubeny, I.↗

A Variational Assimilation Method for Satellite and Conventional Data: Model 2 (version 1)

The Model II variational data assimilation model is the second of the four variational models designed to blend diverse meteorological data into a dynamically constrained data set. Model II differs from Model I in that it includes the thermodynamic equation as the fifth dynamical constraint. Thus, Model II includes all five of the primative equations that govern atmospheric flow for a dry atmosphere.

Achtemeier, Gary L.↗

Dynamic Adjustment of Model Parameters Using Assimilated Boundary Layer Height

The NASA Global Modeling and Assimilation Office is exploring the assimilation of boundary layer height (PBLH) derived from a variety of observations, including radiosondes, GNSS-RO, space- and ground-based lidar, and radar wind profilers. Previous work using the Goddard Earth Observing System (GEOS) has shown that within the lower troposphere, the representation of thermodynamic structure is strongly affected by model parameterizations. In assimilation experiments, parameterizations can compensate for analysis thermodynamic increments and reduce the observation information carried forward in time. To mitigate this effect, the present study explores a dynamic parameter adjustment based on the PBLH analysis field. The PBLH analysis increment is used to adjust model parameters that influence the depth of parameterized boundary layer mixing in both stable and convective regimes. Replay (nudging) experiments and cycled data assimilation experiments are used to evaluate the approach. The parameter adjustment is shown to reduce model compensation of the analysis thermodynamic tendencies and bring model parameterized mixing depths into closer agreement with observed PBLH.

Nathan Arnold↗

Seasonal meridional energy balance and thermal structure of the atmosphere of Uranus - A radiative-convective-dynamical model

A model is presented for the thermodynamics of the seasonal meridional energy balance and thermal structure of the Uranian atmosphere. The model considers radiation and small-scale convection, and dynamical heat fluxes due to large-scale baroclinic eddies. Phase oscillations with a period of 0.5 Uranian year are discerned in the total internal power and global enthalpy storage. The variations in the identity of the main transport agent with the magnitude of the internal heat source are discussed. It is shown that meridional heat transport in the atmosphere is sufficient to lower seasonal horizontal temperature contrasts below those predicted with radiative-convection models.

Friedson, James↗

A satellite-borne radar wind sensor (RAWS)

Modeling global atmospheric circulations and forecasting the weather would improve if worldwide information on winds aloft were available. Accurate prediction of weather is important to agriculture, shipping, air traffic, and many other fields. Global system models of climate are of great importance. Current global atmospheric models use pressure measurements and thermodynamic properties to calculate the effects of wind for use in Numerical Weather Prediction (NWP) models. Inputs to the NWP models are temperature, pressure and wind velocities at different heights. Clearly direct wind measurements could significantly improve the NWP model performance. The RAdar Wind Sounder (RAWS) program at the University of Kansas is a study of the feasibility and the trade-offs in the design of a space-based radar system to measure wind vectors. This can be done by measuring the Doppler shift of cloud and rain returns from three or more points and calculating the components of the wind vector. The RAWS study to date uses the candidate system selected after preliminary study of frequencies and sensitivities. Two frequencies chosen, 10 and 35 GHz, allow higher sensitivity for clouds and more penetration for rain. The past year was devoted to modeling the signal-to-noise ratio (SNR) achievable for the two frequencies. The determination of SNR versus cloud penetration depth used a cloud backscattering and attenuation model in the appropriate radar equation. Calculations assumed reasonable losses in reception and transmission, in addition to the atmospheric attenuation. We discovered that ice clouds provide a higher SNR than previously calculated, but some water clouds give lower SNRs than we calculated before. One of the primary issues in the SNR calculation was the choice of the drop size distribution. Although Xin used several distributions (e.g., log normal, Khrigian and Mazin), this year we used the Deirmendjian cloud model. SNR versus cloud penetration plots were generated to validate the candidate system. Rain, which appears in the cloud models at the lower altitudes, provides ample SNR, as do the higher clouds composed of ice particles. However, in some cloud situations we found the sensitivity for the clouds was marginal or inadequate. At 35 GHz, two of the cloud models characterized by 1 to 2 g/cu m of water content at altitudes extending from 150 to 1500 meters, produced a sufficient SNR. Other models, however, with water contents ranging from 0.5 to 4 g/cu m and altitudes up to 4000 meters, exhibit SNR of -3 to -23 dB, largely because of attenuation in the upper cloud layers. These results coupled with the lower SNR at 10 GHz, led to an investigation of alternate frequencies. The rain present beneath these clouds provides adequate SNR at 10 GHz, and in most cases, at GHz.

Moore, Richard K.↗

Applying Advanced Ground-Based Remote Sensing in the Southeast Asian Maritime Continent to Characterize Regional Proficiencies in Smoke Transport Modeling

This work describes some of the most extensive ground-based observations of the aerosol profile collected in Southeast Asia to date, highlighting the challenges in simulating these observations with a mesoscale perspective. An 84-h WRF Model coupled with chemistry (WRF-Chem) mesoscale simulation of smoke particle transport at Kuching, Malaysia, in the southern Maritime Continent of Southeast Asia is evaluated relative to a unique collection of continuous ground-based lidar, sun photometer, and 4-h radiosonde profiling. The period was marked by relatively dry conditions, allowing smoke layers transported to the site unperturbed by wet deposition to be common regionally. The model depiction is reasonable overall. Core thermodynamics, including landsea-breeze structure, are well resolved. Total model smoke extinction and, by proxy, mass concentration are low relative to observation. Smoke emissions source products are likely low because of undersampling of fires in infrared sun-synchronous satellite products, which is exacerbated regionally by endemic low-level cloud cover. Differences are identified between the model mass profile and the lidar profile, particularly during periods of afternoon convective mixing. A static smoke mass injection height parameterized for this study potentially influences this result. The model does not resolve the convective mixing of aerosol particles into the lower free troposphere or the enhancement of near-surface extinction from nighttime cooling and hygroscopic effects.

Aerosols; Air quality; Lidars/Lidar observations; ↗

Automotive Gas Turbine Power System-Performance Analysis Code

An open cycle gas turbine numerical modelling code suitable for thermodynamic performance analysis (i.e. thermal efficiency, specific fuel consumption, cycle state points, working fluid flowrates etc.) of automotive and aircraft powerplant applications has been generated at the NASA Lewis Research Center's Power Technology Division. The use this code can be made available to automotive gas turbine preliminary design efforts, either in its present version, or, assuming that resources can be obtained to incorporate empirical models for component weight and packaging volume, in later version that includes the weight-volume estimator feature. The paper contains a brief discussion of the capabilities of the presently operational version of the code, including a listing of input and output parameters and actual sample output listings.

Juhasz, Albert J.↗

Multivariate Error Covariance Estimates by Monte-Carlo Simulation for Assimilation Studies in the Pacific Ocean

One of the most difficult aspects of ocean state estimation is the prescription of the model forecast error covariances. The paucity of ocean observations limits our ability to estimate the covariance structures from model-observation differences. In most practical applications, simple covariances are usually prescribed. Rarely are cross-covariances between different model variables used. Here a comparison is made between a univariate Optimal Interpolation (UOI) scheme and a multivariate OI algorithm (MvOI) in the assimilation of ocean temperature. In the UOI case only temperature is updated using a Gaussian covariance function and in the MvOI salinity, zonal and meridional velocities as well as temperature, are updated using an empirically estimated multivariate covariance matrix. Earlier studies have shown that a univariate OI has a detrimental effect on the salinity and velocity fields of the model. Apparently, in a sequential framework it is important to analyze temperature and salinity together. For the MvOI an estimation of the model error statistics is made by Monte-Carlo techniques from an ensemble of model integrations. An important advantage of using an ensemble of ocean states is that it provides a natural way to estimate cross-covariances between the fields of different physical variables constituting the model state vector, at the same time incorporating the model's dynamical and thermodynamical constraints as well as the effects of physical boundaries. Only temperature observations from the Tropical Atmosphere-Ocean array have been assimilated in this study. In order to investigate the efficacy of the multivariate scheme two data assimilation experiments are validated with a large independent set of recently published subsurface observations of salinity, zonal velocity and temperature. For reference, a third control run with no data assimilation is used to check how the data assimilation affects systematic model errors. While the performance of the UOI and MvOI is similar with respect to the temperature field, the salinity and velocity fields are greatly improved when multivariate correction is used, as evident from the analyses of the rms differences of these fields and independent observations. The MvOI assimilation is found to improve upon the control run in generating the water masses with properties close to the observed, while the UOI failed to maintain the temperature and salinity structure.

Borovikov, Anna↗

A Simplified Model for Detonation Based Pressure-Gain Combustors

A time-dependent model is presented which simulates the essential physics of a detonative or otherwise constant volume, pressure-gain combustor for gas turbine applications. The model utilizes simple, global thermodynamic relations to determine an assumed instantaneous and uniform post-combustion state in one of many envisioned tubes comprising the device. A simple, second order, non-upwinding computational fluid dynamic algorithm is then used to compute the (continuous) flowfield properties during the blowdown and refill stages of the periodic cycle which each tube undergoes. The exhausted flow is averaged to provide mixed total pressure and enthalpy which may be used as a cycle performance metric for benefits analysis. The simplicity of the model allows for nearly instantaneous results when implemented on a personal computer. The results compare favorably with higher resolution numerical codes which are more difficult to configure, and more time consuming to operate.

Paxson, Daniel E.↗

Single and Multi-Node Modeling of Direct, Submerged, and Self-Pressurization of A Cryogenic Propellant Tank Using Nodal Tools

The pressurization of cryogenic propellant tanks will always be an important process so long as cryogenic liquids are being considered as fuel sources or used for other in-space applications. Pressure control of the tank ullage is necessary for managing propellant flowrates to an engine or a receiver tank, and modeling of the process is used to predict the pressurant requirements and the amount of propellant boiloff. Direct ullage pressurization is the more traditional approach to tank pressurization, as the physics are straight-forward, and ample test (flight) data have been collected and analyzed over the past several decades. Submerged injection pressurization is an alternate method for tank pressurization and has been shown to reduce pressurant requirements, subcool the propellant, and reduce the risk of ullage collapse. Additionally, the pressurant gas entering the ullage is usually much colder when using the submerged pressurization approach, resulting in reduced propellant boiloff. These benefits are at the expense of vaporizing a small percentage of the propellent. Both tank pressurization methods are viable options for current and future space missions, and it is important to have the capability of analyzing the tank ullage conditions for both approaches. Our previous work has demonstrated the development of a Generalized Fluid System Simulation Program (GFSSP) model, which contains a thermodynamic equilibrium heat and mass transfer subroutine capable of effectively analyzing both direct and submerged pressurization systems [1-2]. This subroutine has most recently been enhanced to include the non-equilibrium effect of pressurant dissolution into the propellant. To date the ullage has always been represented as a single node, and although the simulated single-node temperatures have good comparison with the volume-averaged ullage temperatures computed from test data, the physics of the thermal stratification in the ullage were never captured, and adjustment factors in the model were required. The purpose of this paper is to introduce the development of a multi-node ullage model using GFSSP and to discuss the improvements of the simulated ullage temperature distribution and its resulting effects on ullage heat transfer processes. Test data from the Cryogenic Propellant Storage and Transfer Engineering Developmental Unit (CPST EDU) was used for model validation. For additional comparison, a Thermal Desktop (TD) model was also developed to analyze the CPST EDU direct ullage pressurization tests using both a single node and multi-node approach. The model includes the direct pressurant line, vent line, fill/drain line, and a TD FloCAD Compartment. The TD FloCAD Compartment is employed to represent the liquid and ullage as single volumes inside the tank, to include a liquid/vapor interface, and to generate network level objects such as lumps (analogous to nodes in GFSSP), paths, and ties between the fluid and thermal elements. An established heat load on the model tank was leveraged from a pre-existing higher-fidelity model correlated to CPST EDU test data.

pressurization↗

Equilibrium gas flow computations. I - Accurate and efficient calculation of equilibrium gas properties

This paper treats the accurate and efficient calculation of thermodynamic properties of arbitrary gas mixtures for equilibrium flow computations. New improvements in the Stupochenko-Jaffe model for the calculation of thermodynamic properties of diatomic molecules are presented. A unified formulation of equilibrium calculations for gas mixtures in terms of irreversible entropy is given. Using a highly accurate thermo-chemical data base, a new, efficient and vectorizable search algorithm is used to construct piecewise interpolation procedures with generate accurate thermodynamic variable and their derivatives required by modern computational algorithms. Results are presented for equilibrium air, and compared with those given by the Srinivasan program.

Liu, Yen↗

Generalized Fluid System Simulation Program

A general purpose program implemented on a computer analyzes steady state and transient flow in a complex fluid network, modeling phase changes, compressibility, mixture thermodynamics and external body forces such as gravity and centrifugal force. A preprocessor provides for the inter- active development of a fluid network simulation having nodes and branches. Mass, energy, and specie conservation equations are solved at the nodes, and momentum conservation equations are solved in the branches. Contained herein are subroutines for computing "real fluid" thermodynamic and thermophysical properties for 12 fluids, and a number of different source options are provided for model- ing momentum sources or sinks in the branches. 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. Application and verification of this invention are provided through an example problem, which demonstrates that the predictions of the present invention compare most reasonably with test data.

Alok Kumar Majumdar↗

Effects of Increasing the Category Resolution of the Sea Ice Thickness Distribution in a Coupled Climate Model on Arctic and Antarctic Sea Ice Mean State

Many modern sea ice models used in global climate models represent the subgrid-scale heterogeneity in sea ice thickness with an ice thickness distribution (ITD), which improves model realism by representing the significant impact of the high spatial heterogeneity of sea ice thickness on thermodynamic and dynamic processes. Most models default to five thickness categories. However, little has been done to explore the effects of the resolution of this distribution (number of categories) on sea-ice feedbacks in a coupled model framework and resulting representation of the sea ice mean state. Here, we explore this using sensitivity experiments in CESM2 with the standard 5 ice thickness categories and 15 ice thickness categories. Increasing the resolution of the ITD in a run with preindustrial climate forcing results in substantially thicker Arctic sea ice year-round. Analyses show that this is a result of the ITD influence on ice strength. With 15 ITD categories, weaker ice occurs for the same average thickness, resulting in a higher fraction of ridged sea ice. In contrast, the higher resolution of thin ice categories results in enhanced heat conduction and bottom growth and leads to only somewhat increased winter Antarctic sea ice volume. The spatial resolution of the ICESat-2 satellite mission provides a new opportunity to compare model outputs with observations of seasonal evolution of the ITD in the Arctic (ICESat-2; 2018–2021). Comparisons highlight significant differences from the ITD modeled with both runs over this period, likely pointing to underlying issues contributing to the representation of average thickness.

Madison M. Smith↗

Thermodynamic Data to 20,000 K For Monatomic Gases

This report contains standard-state thermodynamic functions for 50 gaseous atomic elements plus deuterium and electron gas, 51 singly ionized positive ions, and 36 singly ionized negative ions. The data were generated by the NASA Lewis computer program PAC97, a modified version of PAC91 reported in McBride and Gordon. This report is being published primarily to document part of the data currently being used in several NASA Lewis computer programs. The data are presented in tabular and graphical format and are also represented in the form of least-squares coefficients. The tables give the following data as functions of temperature : heat capacity, enthalpy, entropy Gibbs energy, enthalpy of formation, and equilibrium constant. A brief discussion and a comparison of calculated results are given for several models for calculating ideal thermodynamic data for monatomic gases.

Gordon, Sanford↗

Session 6: Dynamic Modeling and Systems Analysis

These presentations cover some of the ongoing work in dynamic modeling and dynamic systems analysis. The first presentation discusses dynamic systems analysis and how to integrate dynamic performance information into the systems analysis. The ability to evaluate the dynamic performance of an engine design may allow tradeoffs between the dynamic performance and operability of a design resulting in a more efficient engine design. The second presentation discusses the Toolbox for Modeling and Analysis of Thermodynamic Systems (T-MATS). T-MATS is a Simulation system with a library containing the basic building blocks that can be used to create dynamic Thermodynamic Systems. Some of the key features include Turbo machinery components, such as turbines, compressors, etc., and basic control system blocks. T-MAT is written in the Matlab-Simulink environment and is open source software. The third presentation focuses on getting additional performance from the engine by allowing the limit regulators only to be active when a limit is danger of being violated. Typical aircraft engine control architecture is based on MINMAX scheme, which is designed to keep engine operating within prescribed mechanical/operational safety limits. Using a conditionally active min-max limit regulator scheme, additional performance can be gained by disabling non-relevant limit regulators

control↗

Retrieval of wind temperature and pressure from single Doppler radar and a numerical model

A 4D data assimilation algorithm to obtain 3D wind and thermodynamic fields from radial velocities and a numerical model is derived. It is a blend of the continuous updating technique of Charney et al. (1969) and a diagnostic pressure and temperature retrieval technique suggested by Gal-Chen (1978). This model is tested only against model-generated data which are viewed as 'data observed from a real atmosphere'. A thermal bubble developing in a dry, neutral environment is simulated. The thermal recovery technique is found to be a necessary procedure to obtain successful data assimilation results. The present assimilation method is found to be capable of reducing the 'observational error' and making the assimilation run converge toward the control run. Without data insertion, errors in the nonobserved wind component are large.

Liou, Yu-Chieng↗

Tensegrity: the architectural basis of cellular mechanotransduction

Physical forces of gravity, hemodynamic stresses, and movement play a critical role in tissue development. Yet, little is known about how cells convert these mechanical signals into a chemical response. This review attempts to place the potential molecular mediators of mechanotransduction (e.g. stretch-sensitive ion channels, signaling molecules, cytoskeleton, integrins) within the context of the structural complexity of living cells. The model presented relies on recent experimental findings, which suggests that cells use tensegrity architecture for their organization. Tensegrity predicts that cells are hard-wired to respond immediately to mechanical stresses transmitted over cell surface receptors that physically couple the cytoskeleton to extracellular matrix (e.g. integrins) or to other cells (cadherins, selectins, CAMs). Many signal transducing molecules that are activated by cell binding to growth factors and extracellular matrix associate with cytoskeletal scaffolds within focal adhesion complexes. Mechanical signals, therefore, may be integrated with other environmental signals and transduced into a biochemical response through force-dependent changes in scaffold geometry or molecular mechanics. Tensegrity also provides a mechanism to focus mechanical energy on molecular transducers and to orchestrate and tune the cellular response.

Review↗