Engineering PapersSearch

SEARCH · Engineering Papers

Results for “radiative exchange 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 73 records · Page 4

Mars Propellant Liquefaction and Storage Performance Modeling using Thermal Desktop with an Integrated Cryocooler Model

NASAs current Mars architectures are assuming the production and storage of 23 tons of liquid oxygen on the surface of Mars over a duration of 500+ days. In order to do this in a mass efficient manner, an energy efficient refrigeration system will be required. Based on previous analysis NASA has decided to do all liquefaction in the propulsion vehicle storage tanks. In order to allow for transient Martian environmental effects, a propellant liquefaction and storage system for a Mars Ascent Vehicle (MAV) was modeled using Thermal Desktop. The model consisted of a propellant tank containing a broad area cooling loop heat exchanger integrated with a reverse turbo Brayton cryocooler. Cryocooler sizing and performance modeling was conducted using MAV diurnal heat loads and radiator rejection temperatures predicted from a previous thermal model of the MAV. A system was also sized and modeled using an alternative heat rejection system that relies on a forced convection heat exchanger. Cryocooler mass, input power, and heat rejection for both systems were estimated and compared against sizing based on non-transient sizing estimates.

fluid modeling

Mars Propellant Liquefaction Modeling in Thermal Desktop

NASAs current Mars architectures are assuming the production and storage of 23 tons of liquid oxygen on the surface of Mars over a duration of 500+ days. In order to do this in a mass efficient manner, an energy efficient refrigeration system will be required. Based on previous analysis NASA has decided to do all liquefaction in the propulsion vehicle storage tanks. In order to allow for transient Martian environmental effects, a propellant liquefaction and storage system for a Mars Ascent Vehicle (MAV) was modeled using Thermal Desktop. The model consisted of a propellant tank containing a broad area cooling loop heat exchanger integrated with a reverse turbo Brayton cryocooler. Cryocooler sizing and performance modeling was conducted using MAV diurnal heat loads and radiator rejection temperatures predicted from a previous thermal model of the MAV. A system was also sized and modeled using an alternative heat rejection system that relies on a forced convection heat exchanger. Cryocooler mass, input power, and heat rejection for both systems were estimated and compared against sizing based on non-transient sizing estimates.

Cryogenics

The problem of cooling the cold Io torus

Two models are developed for the transport of ions inward from Io's orbit, the first using the assumption that radial diffusion is the dominant transport mechanism while the second uses a combination of diffusive and convective transport. The models include thermal and number density transport, radiation, ionization and pickup of local neutrals, recombination, charge exchange, and Coulomb interactions. It is found that pure diffusive transport can account for the dramatic inward depletion of the torus only by invoking recombination or by postulating a massive increase in the production rate of torus ions sometime prior to the Voyager mission. Recombination is determined to be far too slow to be the cause of the observed density decrease inside of Io. The model combining convection and diffusion is shown to reasonably match the data, but only if a diffusion coefficient 100 times less than that derived from Pioneer observations is used. In addition, it is shown that the Pioneer derived diffusion rate combined with Voyager temperature and density measurements indicate a large nonradiative sink of energy in the inner torus.

Richardson, J. D.

Determination of the cross-sectional temperature distribution and boiling limitation of a heat pipe

A computer model is developed and verified which is capable of determining the cross-sectional temperature distribution within a heat pipe with an attached radiator fin; such heat pipes would be plugged into contact heat exchangers designed to carry heat from a space station habitation module to the radiator elements through a centralized fluid loop. The model can furnish information for determining the susceptibility of the monogroove heat pipe to boiling, as well as the location and magnitude of that boiling.

Peterson, G. P.

The liquid droplet radiator - An ultralightweight heat rejection system for efficient energy conversion in space

A heat rejection system for space is described which uses a recirculating free stream of liquid droplets in place of a solid surface to radiate waste heat. By using sufficiently small droplets (less than about 100 micron diameter) of low vapor pressure liquids (tin, tin-lead-bismuth eutectics, vacuum oils) the radiating droplet sheet can be made many times lighter than the lightest solid surface radiators (heat pipes). The liquid droplet radiator (LDR) is less vulnerable to damage by micrometeoroids than solid surface radiators, and may be transported into space far more efficiently. Analyses are presented of LDR applications in thermal and photovoltaic energy conversion which indicate that fluid handling components (droplet generator, droplet collector, heat exchanger, and pump) may comprise most of the radiator system mass. Even the unoptimized models employed yield LDR system masses less than heat pipe radiator system masses, and significant improvement is expected using design approaches that incorporate fluid handling components more efficiently. Technical problems (e.g., spacecraft contamination and electrostatic deflection of droplets) unique to this method of heat rejection are discussed and solutions are suggested.

Mattick, A. T.

The Liquid Droplet Radiator - an Ultralightweight Heat Rejection System for Efficient Energy Conversion in Space

A heat rejection system for space is described which uses a recirculating free stream of liquid droplets in place of a solid surface to radiate waste heat. By using sufficiently small droplets ( 100 micron diameter) of low vapor pressure liquids the radiating droplet sheet can be made many times lighter than the lightest solid surface radiators (heat pipes). The liquid droplet radiator (LDR) is less vulnerable to damage by micrometeoroids than solid surface radiators, and may be transported into space far more efficiently. Analyses are presented of LDR applications in thermal and photovoltaic energy conversion which indicate that fluid handling components (droplet generator, droplet collector, heat exchanger, and pump) may comprise most of the radiator system mass. Even the unoptimized models employed yield LDR system masses less than heat pipe radiator system masses, and significant improvement is expected using design approaches that incorporate fluid handling components more efficiently. Technical problems (e.g., spacecraft contamination and electrostatic deflection of droplets) unique to this method of heat rejectioon are discussed and solutions are suggested.

A T Mattick

Global transport of atmospheric smoke following a major nuclear exchange

A fully interactive zonally symmetric general circulation model is used to simulate the transport and evolution of a massive smoke cloud that is likely to form in the atmosphere following a major nuclear war. The presence of such a cloud significantly alters the simulated circulation and the subsequent transport of the smoke particles themselves. While the model indicates a tendency for interhemispheric exchange, the most immediate effect is that the radiatively active particles are carried into the stratosphere by the generation of strong vertical motions and intense convection regardless of their initial injection altitudes.

Haberle, R. M.

A LaGrangian Perspective of the Surface Energy Budget Contributions to Sea Ice Evolution in the ‘New Arctic’

Recent Arctic warming outpaces global warming and is characterized by dramatic changes in sea ice (record sea ice loss, thinner and more seasonal ice cover), warming SSTs, and a lengthening melt season. These changes impact the surface energy budget and the sensitivity of sea ice to radiative forcing, namely an increased sensitivity of thinner ice to local radiative forcings and a weaker buffering of surface-atmosphere exchanges of heat and moisture. However, these processes are a large source of uncertainty in global climate models with prior studies finding that poor representations of downwelling longwave radiation and sea ice albedo lead to model spread in sea ice conditions. A limitation of prior work is examining sea ice-atmosphere interactions through an Eulerian framework. This study utilizes a sea ice parcel dataset in a LaGrangian framework, allowing for tracking of surface energy budget-sea ice interactions over the lifespan of individual sea ice parcels; this framework also allows for the quantification of the sea ice response to episodic weather events along with seasonal drivers, including the drivers of rapid growth and melt events. First, lag-lead relationships between relevant surface energy budget terms and sea ice thickness changes are explored to understand the atmospheric conditions preceding rapid growth and melt. The impact of episodic events- moisture intrusion, warm air outbreaks, cyclones- on sea ice growth are also investigated by stratifying first-year parcels by sea ice thickness at the end of the growth season and analyzing the effect of atmospheric variability (number and intensity of episodic events) on parcel thickness change. Lastly, parcel growth and melt and the associated radiative anomalies are compared to snowfall events/snow thickness to determine how snowfall modifies the interaction between the local surface energy budget and sea ice mass balance.

Patrick Taylor

Some effects of global dust storms on the atmospheric circulation of Mars

A zonally symmetric primitive equation modelled on a sphere is used to numerically simulate the Martian atmosphere's response to various dust loads, as well as the ability of its meridional circulation to transport dust globally, where the circulation is driven by heating due to the absorption of solar and IR radiation by dust and CO2, in addition to sensible heat exchange with the ground. A preliminary experiment shows the model distribution of winds and temperature to compare favorably with zonally averaged values from a general circulation model. Experiments simulating the evolution of global dust storms show that dust is effectively transported by the zonal mean circulation, which rapidly intensifies as the dust spreads, and that the basic structure of the circulation is relatively insensitive to details, being mostly dependent on the heating of the tropical and subtropical atmosphere.

Haberle, R. M.

Tropical Convection and Climate Processes in a Cumulus Ensemble Model

Local convective-radiative equilibrium states of the tropical atmosphere are determined by the following external forcing: 1) Insolation, 2) Surface heat and moisture exchanges (primarily radiation and evaporation), 3) Heating and moistening induced by large-scale circulation. Understanding the equilibrium states of the tropical atmosphere in different external forcing conditions is of vital importance for studying cumulus parameterization, climate feedbacks, and climate changes. We extend our previous study using the Goddard Cumulus Ensemble (GCE) Model which resolves convective-radiative processes more explicitly than global climate models do. Several experiments are carried out under fixed insolation and sea surface temperature. The prescribed SST consists of a uniform warm pool (29C) surrounded by uniform cold SST (26C). The model produces "Walker"-type circulation with the ascending branch of the model atmosphere more humid than the descending part, but the vertically integrated temperature does not show a horizontal gradient. The results are compared with satellite measured moisture by SSM/I (Special Sensor Microwave/Imager) and temperature by MSU in the ascending and descending tropical atmosphere. The vertically integrated temperature and humidity in the two model regimes are comparable to the observed values in the tropics.

Sui, Chung-Hsiung

Backscatter of solar resonance radiation. I.

Calculation of the angular dependence of the intensity of solar Lyman alpha resonantly scattered from neutral interstellar hydrogen which has penetrated the solar system. A simple model which essentially neglects temperature effects, but which includes gravity, radiation pressure, photoionization and charge exchange is used for this calculation. The results are then compared with the observations. The resonant scattering of He I 584 A is also treated.

Johnson, H. E.

The radiation budget of a Cirrus layer deduced from simultaneous aircraft observations and model calculations

Several aircraft were employed during the FIRE Cirrus IFO in order to make nearly simultaneous observations of cloud properties and fluxes. A segment of the flight data collected on 28 October 1988 during which the NASA Ames ER-2 overflew the NCAR King Air was analyzed. The ER-2 flew at high altitude making observations of visible and infrared radiances and infrared flux and cloud height and thickness. During this segment, the King Air flew just above the cloud base making observations of ice crystal size and shape, local meteorological variables, and infrared fluxes. While the two aircraft did not collect data exactly coincident in space and time, they did make observations within a few minutes of each other. For this case study, the infrared radiation balance of the cirrus layer is of primary concern. Observations of the upwelling 10 micron radiance, made from the ER-2, can be used to deduce the 10 micron optical depth of the layer. The upwelling broadband infrared flux is also measured from the ER-2. At the same time, the upwelling and downwelling infrared flux at the cloud base is obtained from the King Air measurements. Information on cloud microphysics is also available from the King Air. Using this data in conjunction with atmospheric temperature and humidity profiles from local radiosondes, the necessary inputs for an infrared radiative transfer model can be developed. Infrared radiative transfer calculations are performed with a multispectral two-stream model. The model fluxes at the cloud base and at 19 km are then compared with the aircraft observations to determine whether the model is performing well. Cloud layer heating rates can then be computed from the radiation exchange.

Ackerman, Thomas P.

System Mass Variation and Entropy Generation in 100k We Closed-Brayton-Cycle Space Power Systems

State-of-the-art closed-Brayton-cycle (CBC) space power systems were modeled to study performance trends in a trade space characteristic of interplanetary orbiters. For working-fluid molar masses of 48.6, 39.9, and 11.9 kg/kmol, peak system pressures of 1.38 and 3.0 MPa and compressor pressure ratios ranging from 1.6 to 2.4, total system masses were estimated. System mass increased as peak operating pressure increased for all compressor pressure ratios and molar mass values examined. Minimum mass point comparison between 72 percent He at 1.38 MPa peak and 94 percent He at 3.0 MPa peak showed an increase in system mass of 14 percent. Converter flow loop entropy generation rates were calculated for 1.38 and 3.0 MPa peak pressure cases. Physical system behavior was approximated using a pedigreed NASA Glenn modeling code, Closed Cycle Engine Program (CCEP), which included realistic performance prediction for heat exchangers, radiators and turbomachinery.

Barrett, Michael J.

System Mass Variation and Entropy Generation in 100-kWe Closed-Brayton-Cycle Space Power Systems

State-of-the-art closed-Brayton-cycle (CBC) space power systems were modeled to study performance trends in a trade space characteristic of interplanetary orbiters. For working-fluid molar masses of 48.6, 39.9, and 11.9 kg/kmol, peak system pressures of 1.38 and 3.0 MPa and compressor pressure ratios ranging from 1.6 to 2.4, total system masses were estimated. System mass increased as peak operating pressure increased for all compressor pressure ratios and molar mass values examined. Minimum mass point comparison between 72 percent He at 1.38 MPa peak and 94 percent He at 3.0 MPa peak showed an increase in system mass of 14 percent. Converter flow loop entropy generation rates were calculated for 1.38 and 3.0 MPa peak pressure cases. Physical system behavior was approximated using a pedigreed NASA Glenn modeling code, Closed Cycle Engine Program (CCEP), which included realistic performance prediction for heat exchangers, radiators and turbomachinery.

Barrett, Michael J.

Geothermal Heat Pump System Showcase: Short-Term Validation of Borehole Heat Exchanger Performance from Field Data to Numerical Modeling: Preprint

Since 2011, a geothermal heat pump (GHP) system has been operating to provide space heating and cooling for the Solar Radiation and Research Laboratory building at the National Laboratory of the Rockies (NLR) in Golden, Colorado. The system consists of 23 vertical boreholes, each extending to a depth of 300 ft (91 m), connected to 11 water-to-air heat pump units and four circulation pumps. Between fiscal years 2023 and 2025, additional power meters and temperature sensors were retrofitted to support detailed system performance assessment and model development. This study presents preliminary monitoring results and the development of an initial numerical model of the borehole heat exchanger field. The model incorporated site-specific geometry, ground thermal properties derived from thermal response tests, and ambient temperatures, and simulated system behavior over a representative operating day in September. Model predictions of outlet temperatures were compared against corresponding field measurements. Results showed that modeling initialized with a simplified linear subsurface temperature gradient presents systematic discrepancies in outlet temperature, whereas incorporating depth-resolved borehole temperature measurements for initialization yields substantially improved agreement with observations. The findings highlight the sensitivity of short-term predictive modeling to the representation of initial subsurface thermal conditions and underscore the value of high-resolution field measurements for model calibration and validation. These preliminary results inform ongoing efforts to extend the modeling framework to longer time horizons and to refine monitoring and modeling strategies that support the design guidance and operational optimization of GHP systems in research and commercial buildings.

15 GEOTHERMAL ENERGY

A model of the planetary boundary layer over a snow surface

A model of the planetary boundary layer over a snow surface has been developed. It contains the vertical heat exchange processes due to radiation, conduction, and atmospheric turbulence. Parametrization of the boundary layer is based on similarity functions developed by Hoffert and Sud (1976), which involve a dimensionless variable, dependent on boundary-layer height and a localized Monin-Obukhov length. The model also contains the atmospheric surface layer and the snowpack itself, where snowmelt and snow evaporation are calculated. The results indicate a strong dependence of surface temperatures, especially at night, on the bursts of turbulence which result from the frictional damping of surface-layer winds during periods of high stability, as described by Businger (1973). The model also shows the cooling and drying effect of the snow on the atmosphere, which may be the mechanism for air mass transformation in sub-Arctic regions.

Halberstam, I.

Enthalpy-Based Thermal Evolution of Loops: II. Improvements to the Model

This paper further develops the zero-dimensional (0D) hydrodynamic coronal loop model "Enthalpy-based Thermal Evolution of Loops" (EBTEL) originally proposed by Klimchuk et al (2008), which studies the plasma response to evolving coronal heating. It has typically been applied to impulsive heating events. The basis of EBTEL is the modelling of mass exchange between the corona and transition region and chromosphere in response to heating variations, with the key parameter being the ratio of transition region to coronal radiation. We develop new models for this parameter that now include gravitational stratification and a physically motivated approach to radiative cooling. A number of examples are presented, including nanoflares in short and long loops, and a small flare. It is found that while the evolution of the loop temperature is rather insensitive to the details of the model, accurate tracking of the density requires the inclusion of our new features. In particular, we are able to now obtain highly over-dense loops in the late cooling phase and decreases to the coronal density arising due to stratification. The 0D results are compared to a 1D hydro code (Hydrad). The agreement is acceptable, with the exception of the flare case where some versions of Hydrad can give significantly lower densities. This is attributed to the method used to model the chromosphere in a flare. EBTEL is suitable for general use as a tool for (a) quick-look results of loop evolution in response to a given heating function and (b) situations where the modelling of hundreds or thousands of elemental loops is needed. A single run takes a few seconds on a contemporary laptop.

Cargill, P. J.

Sensitivity of the Assimilated Ozone in the UTLS to Model and Data Selection Changes

This presentation will discuss the sensitivity of assimilated ozone fields in the upper troposphere and lower stratosphere (UTLS) to a number of factors, focusing mainly on aspects of data selection and the prediction model. This is important, because assimilation represents an attempt to construct our best estimates of the true ozone field; however, inaccuracies in the UTLS ozone distribution translate into an uncertainty in factors such as the calculated radiative forcing of climate or the inferred stratosphere-troposphere exchange (STE) of ozone. The 3D ozone data assimilation system, from NASA's Global Modeling and Assimilation Office (GMAO), combines observations of total ozone column and stratospheric profiles with predictions from an off-line, parameterized chemistry and transport model (pCTM) to produce six-hourly, global analyses. The first experiments discussed assimilate ozone retrievals from the Earth-Probe Total Ozone Mapping Spectrometer (EPTOMS) and stratospheric profiles from the Solar Backscatter UltraViolet/2 (SBUV/2) instrument. The SBUV/2 ozone data have a coarse vertical resolution, with increased uncertainty below the ozone maximum, and TOMS provides only total ozone columns. Thus, the assimilated ozone profiles in the UTLS region are only weakly constrained by the incoming SBUV and TOMS data. Consequently, the assimilated ozone distribution should be sensitive to changes in inputs to the statistical analysis scheme. Sensitivity studies have been conducted to examine the responses to TOMS and SBUV/2 data selection, modifications of the forecast and observation error covariance models, and the model formulation (turning off chemistry or using different wind analyses in the pCTM). The second set of experiments includes an additional data type: ozone retrieved from infrared limb-emission by MIPAS on Envisat. These data offer not only improved vertical resolution in the stratosphere, but also give measurements in the polar night. Comparisons of the assimilated ozone fields from both sets of experiments with independent observations, primarily ozone sondes, are used to determine the impact of each of these changes. It is shown that many of the changes have a significant impact on the UTLS ozone estimates. Implications for interpretation of STE and radiative forcing of climate are discussed.

Pawson, Steven