Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermodynamic simulation”

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 199 records · Page 11

Clouds-radiation interactions in a general circulation model - Impact upon the planetary radiation balance

Simultaneously conducted observations of the earth radiation budget and the cloud amount estimates, taken during the June 1979 - May 1980 Nimbus 7 mission were used to show interactions between the cloud amount and raidation and to verify a long-term climate simulation obtained with the latest version of the NCAR Community Climate Model (CCM). The parameterization of the radiative, dynamic, and thermodynamic processes produced the mean radiation and cloud quantities that were in reasonable agreement with satellite observations, but at the expense of simulating their short-term fluctuations. The results support the assumption that the inclusion of the cloud liquid water (ice) variable would be the best mean to reduce the blinking of clouds in NCAR CCM.

Smith, Laura D.↗

Dynamical Perturbation of the Stratosphere by a Pyrocumulonimbus Injection of Carbonaceous Aerosols

The Pacific Northwest pyrocumulonimbus Event (PNE) took place in British Columbia during the evening and night-time hours between the 12th and 13th of August 2017. Several pyroconvective clouds erupted in this occasion, and released in the upper troposphere/lower stratosphere unprecedented amounts of carbonaceous aerosols (300 ktn). Only a few years later, an even larger pyroCb injection took place over Australia. This event, named the Australian New Year (ANY) event, injected up to 1100 ktn of aerosol between December 29th 2019 and January 4th 2020. Such large injections of carbonaceous aerosol modify the stratospheric radiative budgets, locally perturbing stratospheric temperatures and winds. In this study, we use the Goddard Earth Observing System Chemistry Climate Model (GEOS CCM) to study the pertubations on the stratospheric meteorology induced by an aerosol injection of the magnitude of the PNE. Our simulations include the radiative interactions of aerosols, so that their impact on temperatures and winds are explicitly simulated. We show how the presence of the carbonaceous aerosols from the pyroCb causes the formation and maintenance of a synoptic scale stratospheric anticyclone. We follow this disturbance considering the potential vorticity anomaly and the brown carbon aerosol loading and we describe its dynamical and thermodynamical structure and its evolution in time. The analysis presented here shows that the simulated anticyclone undergoes daily expansion-compression cycles governed by the radiative heating, which are directly related to the vertical motion of the plume, and that the aerosol radiative heating is essential in maintaining the anticyclone itself.

Pyrocumulonimbus↗

Sulfur Oxidation and Contrail Precursor Chemistry

Sulfuric acid (H2SO4), formed in commercial aircraft operations via fuel-S (goes to) SO2 (goes to) SO3 (goes to) H2SO4 plays an important role in the formation of contrails. It is believed that the first step occurs inside the combustor, the second step in the engine exit nozzle, and the third step in the exhaust plume. Thus, measurements of the sulfur oxidation rates are critical to the understanding of contrail formation. Field measurements of contrails formed behind commercial aircraft indicate that significantly greater conversion of fuel-bound sulfur to sulfate aerosol occurs than can be explained by our current knowledge of contrail physics and chemistry. The conversion of sulfur from S(IV) to S(VI) oxidation state, required for sulfate aerosol formation, is thermodynamically favored for the conditions that exist within jet engines but is kinetically disfavored. The principal reaction pathway is O+SO2+M (goes to) SO3+M. The rates of this reaction have never been measured in the temperature and pressure regimes available to aircraft operation. In the first year (FY02) of this project, we performed a series of experiments to elucidate the rate information for the O+SO2+M (goes to) SO3+M reaction. The work performed is described following the proposed work plan. Because we used the H2/O2 system for an O-atom source and rate coefficients were obtained via computer simulation, construction of a reaction mechanism and either recalculation or estimation of thermodynamic properties of H(x)SO(y) species are described first.

DeWitt, Kenneth J.↗

A Novel Approach for Modeling Chemical Reaction in Generalized Fluid System Simulation Program

The Generalized Fluid System Simulation Program (GFSSP) is a computer code developed at NASA Marshall Space Flight Center for analyzing steady state and transient flow rates, pressures, temperatures, and concentrations in a complex flow network. The code, which performs system level simulation, can handle compressible and incompressible flows as well as phase change and mixture thermodynamics. Thermodynamic and thermophysical property programs, GASP, WASP and GASPAK provide the necessary data for fluids such as helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, water, a hydrogen, isobutane, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, several refrigerants, nitrogen trifluoride and ammonia. The program which was developed out of need for an easy to use system level simulation tool for complex flow networks, has been used for the following purposes to name a few: Space Shuttle Main Engine (SSME) High Pressure Oxidizer Turbopump Secondary Flow Circuits, Axial Thrust Balance of the Fastrac Engine Turbopump, Pressurized Propellant Feed System for the Propulsion Test Article at Stennis Space Center, X-34 Main Propulsion System, X-33 Reaction Control System and Thermal Protection System, and International Space Station Environmental Control and Life Support System design. There has been an increasing demand for implementing a combustion simulation capability into GFSSP in order to increase its system level simulation capability of a liquid rocket propulsion system starting from the propellant tanks up to the thruster nozzle for spacecraft as well as launch vehicles. The present work was undertaken for addressing this need. The chemical equilibrium equations derived from the second law of thermodynamics and the energy conservation equation derived from the first law of thermodynamics are solved simultaneously by a Newton-Raphson method. The numerical scheme was implemented as a User Subroutine in GFSSP.

Sozen, Mehmet↗

Molecular Dynamics Simulations of Liquid and Polymer Electrolytes for Energy Storage Devices

Advancing beyond current lithium-ion technology is necessary in order to enable energy storage devices for electric airplanes. Electrolyte stability is a key limiting factor, yet the design of improved electrolytes remains a formidable challenge. Molecular dynamics (MD) simulations are a powerful tool for studying electrolytes, since they can be used to evaluate structural, thermodynamic, and transport properties, and can provide molecular-level detail often inaccessible to experimental techniques. Our computational materials groups at the NASA Ames Research Center has developed models and methods to accurately simulate both liquid and polymer electrolytes.We report the results from atomistic MD simulations of several electrolyte materials, with lithium salts dissolved in ionic liquids, dimethoxyethane (DME), and polyethylene oxide (PEO). For improved accuracy, we employ polarizable models, where each atom is given an environment-dependent atomic dipole. The simulations accurately predict bulk transport properties, including viscosity, diffusion, and ionic conductivity, in quantitative agreement with available experimental data. Moreover, the simulations provide important insights into the solvation structure of the lithium ions.We also report the results from coarse-grained MD simulations of polyanion electrolytes. In order to more efficiently capture the longer length- and time-scales of these systems, we employ a generic bead-spring model. These simulations provide important insight into how the polymer chain architecture and ionic interaction strengths affect the ionic aggregation behavior and cation dynamics. Despite the simplicity of the model, the simulations yield qualitative agreement with experimental data for similar systems.

Abbott, Lauren J.↗

Controls-structures-interaction dynamics during RCS control of the Orbiter/SRMS/SSF configuration

During the assembly flights of the Space Station Freedom (SSF), the Orbiter will either dock with the SSF and retract to the final berthed position, or will grapple the SSF using the Shuttle Remote Manipulator System (SRMS) and maneuver the SRMS coupled vehicles to their final berthed position. The SRMS method is expected to take approximately one to one and a half hours to complete and require periodic attitude corrections by either the Orbiter or the SSF reaction control system (RCS) or continuous control by a control moment gyro (CMG) system with RCS desaturation as required. Free drift of the attached vehicles is not currently thought to be acceptable because the desired system attitude will quickly deteriorate due to unbalanced gravity gradient and aerodynamic torques resulting in power generation problems, thermodynamic control problems, and communications problems. This paper deals with the simulation and control of the SRMS during trunnion/latch interaction dynamics and during RCS maneuvers. The SRMS servo drive joints have highly non-linear elastic characteristics which tend to degrade sensitive control strategies. In addition the system natural frequencies are extremely low and depend on the drive joint deflections and SRMS geometric position. The lowest mean period of oscillation for the Orbiter/SRMS/SSF(MB6) system in brakes hold mode positioned near the final berthed position is approximately 120 seconds. A detailed finite element model of the SRMS has been developed and used in a newly developed SRMS systems dynamics simulation to investigate the non-linear transient response dynamics of the Orbiter/SRMS/SSF systems. The present SRMS control strategy of brakes only recommended by the Charles Draper Labs is contrasted with a robust controller developed by the authors. The robust controller uses an optimal inear quadratic regulator (LQR) to optimally place the closed-loop poles of a multivariable continuous-time system within the common region of an open sector with the sector angle plus or minus 45 degrees from the negative real axis, and the left-hand side of a parallel to the imaginary axis in the complex s-plane. This guarantees that the critical damping ratio for the desired control modes is equal to or in excess of 0.707. The matrix sign function is used for solving the Riccati equations which appear in the controller design procedure. Fast and stable algorithms have recently been developed for the computation of the matrix sign function. Simulation results are given which demonstrate the potential CSI involvement for the current SRMS control system and the proposed control system.

Schliesing, J. A.↗

The role of sea ice dynamics in global climate change

The topics covered include the following: general characteristics of sea ice drift; sea ice rheology; ice thickness distribution; sea ice thermodynamic models; equilibrium thermodynamic models; effect of internal brine pockets and snow cover; model simulations of Arctic Sea ice; and sensitivity of sea ice models to climate change.

Hibler, William D., III↗

COLD-SAT Dynamic Model Computer Code

COLD-SAT Dynamic Model (CSDM) computer code implements six-degree-of-freedom, rigid-body mathematical model for simulation of spacecraft in orbit around Earth. Investigates flow dynamics and thermodynamics of subcritical cryogenic fluids in microgravity. Consists of three parts: translation model, rotation model, and slosh model. Written in FORTRAN 77.

Bollenbacher, G.↗

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↗

Improving Boundary Layer Data Assimilation Using Observation Data from Multiple Observing Systems in the NASA GEOS System

The Planetary Boundary Layer (PBL) is a complex interface that mediates energy and moisture exchanges between the Earth’s surface and atmosphere. Accurate simulation and observation of PBL characteristics, such as PBL height and thermodynamic structure, have proven to be a challenge. In our latest efforts, we have focused on improving PBL thermodynamic structure using data from multiple observing systems in the Goddard Earth Observing System (GEOS), developed by the NASA Global Modeling and Assimilation Office (GMAO). We present strategies and results from assimilating PBL height data derived from radiosondes, GNSS radio occultation, space-based lidar (CALIPSO, CATS, IceSat-2), ground-based lidar (MPLNET), and radar wind profilers, including a novel global PBL height analysis dataset. We also discuss the impacts of better representing capping inversions by using PBL height data together with other observations in GEOS through adjustments to the background error covariance. Long-term statistics of the impact of assimilating and utilizing PBL height data in GEOS are presented. In addition, we explore an innovative approach to assimilate GNSS-RO refractivity data in the lower troposphere.

Eun-Gyeong Yang↗

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↗

Surface Emissivity Effects on Thermodynamic Retrieval of IR Spectral Radiance

The surface emissivity effect on the thermodynamic parameters (e.g., the surface skin temperature, atmospheric temperature, and moisture) retrieved from satellite infrared (IR) spectral radiance is studied. Simulation analysis demonstrates that surface emissivity plays an important role in retrieval of surface skin temperature and terrestrial boundary layer (TBL) moisture. NAST-I ultraspectral data collected during the CLAMS field campaign are used to retrieve thermodynamic properties of the atmosphere and surface. The retrievals are then validated by coincident in-situ measurements, such as sea surface temperature, radiosonde temperature and moisture profiles. Retrieved surface emissivity is also validated by that computed from the observed radiance and calculated emissions based on the retrievals of surface temperature and atmospheric profiles. In addition, retrieved surface skin temperature and emissivity are validated together by radiance comparison between the observation and retrieval-based calculation in the window region where atmospheric contribution is minimized. Both simulation and validation results have lead to the conclusion that variable surface emissivity in the inversion process is needed to obtain accurate retrievals from satellite IR spectral radiance measurements. Retrieval examples are presented to reveal that surface emissivity plays a significant role in retrieving accurate surface skin temperature and TBL thermodynamic parameters.

Zhou, Daniel K.↗

Comments on "The Sensitivity Study of Radiative-Convective Equilibrium in the Tropics with a Convective Resolving Model"

In general, there are two broad scientific objectives when using cloud resolving models (CRMs or cloud ensemble models-CEMs) to study tropical convection. The first one is to use them as a physics resolving models to understand the dynamic and microphysical processes associated with the tropical water and energy cycles and their role in the climate system. The second approach is to use the CRMs to improve the representation of moist processes and their interaction with radiation in large-scale models. In order to improve the credibility of the CRMs and achieve the above goals, CRMs using identical initial conditions and large-scale influences need to produce very similar results. Two CRMs produced different statistical equilibrium (SE) states even though both used the same initial thermodynamic and wind conditions. Sensitivity tests to identify the major physical processes that determine the SE states for the different CRM simulations were performed. Their results indicated that atmospheric horizontal wind is treated quite differently in these two CRMs. The model that had stronger surface winds and consequently larger latent and sensible heat fluxes from the ocean produced a warmer and more humid modeled thermodynamic SE state. In addition, the domain mean thermodynamic state is more unstable for those experiments that produced a warmer and more humid SE state. Their simulated wet (warm and humid) SE states are thermally more stable in the lower troposphere (from the surface to 4-5 km in altitude). The large-scale horizontal advective effects on temperature and water vapor mixing ratio are needed when using CRMs to perform long-term integrations to study convective feedback under specified large-scale environments. In addition, it is suggested that the dry and cold SE state simulated was caused by enhanced precipitation but not enough surface evaporation. We find some problems with the interpretation of these three phenomena.

Tao, W.-K.↗

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↗

Thermodynamic properties of a high pressure subcritical UF6/He gas volume (irradiated by an external source)

A computer simulation study concerning a compressed fissioning UF6 gas is presented. The compression is to be achieved by a ballistic piston compressor. Data on UF6 obtained with this compressor were incorporated in the simulation study. As a neutron source to create the fission events in the compressed gas, a fast burst reactor was considered. The conclusion is that it takes a neutron flux in excess of 10 to the 15th power n/sec sq cm to produce measurable increases in pressure and temperature, while a flux in excess of 10 to 19th power n/sq cm sec would probably damage the compressor.

Sterritt, D. E.↗

Hydrodynamic simulations of flare/surge events

A one-dimensional, hydrodynamic, time-dependent model which simulates the major observed dynamics of flare associated surges was developed. The thermodynamics, the surge mass, the time scales, the physical dimensions, and the velocities of typical surge events are reproduced. The surge is created by a sudden pressure increase at the top of the chromosphere, and this pressure pulse produces a disturbance which is followed, with a time-dependent numerical solution, as it propagates upward through the transition region and into the corona. The leading edge of the disturbance is a weak shock which has only a slight effect on the original transition region and coronal thermodynamics. The major effect occurs behind the shock where the temperature is decreased and the density is increased, and this cool, dense region moves upward then falls downward. After the material begins returning to the chromosphere, a second shock is formed which propagates upward, brings the infalling material to rest, and returns the atmosphere to hydrostatic equilibrium.

Steinolfson, R. S.↗

The interaction of a surface cold front with a prefrontal thermodynamically well-mixed boundary layer

The dynamics of a surface cold front as it proceeds from sea to shore and interacts with a thermodynamically well-mixed boundary layer are investigated by means of numerical simulations using a modified version of the two-dimensional model of Reeder and Smith (1987). The theoretical basis of the model is explained, and the results are presented graphically and characterized in detail. Good agreement is found with observational data on Australian summer cold fronts, and the transverse frontal circulation is shown to be significantly affected by diabatic heating.

Reeder, Michael J.↗

Durability of Carbon-Carbon Composites to a Simulated Martian Atmosphere

A new facility was designed and constructed to investigate the effects of the Martian atmosphere on high temperature radiator components. Preliminary experimental results from carbon-carbon composite radiator samples suggest that carbon is unstable in the simulated Martian atmosphere at a temperature of 600 C, consistent with previous thermodynamic modeling. However, more work is needed to identify the exact mechanism. The versatile design of the facility permits additional study of materials stability in the presence of other atmospheres.

Jaworske, D. A.↗