Search NASA⌕ Search

SEARCH · Search NASA

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Investigation of Best-Practices and Computationally Inexpensive Radiative Exchange Models for Discrete Element Method Modeling of Aluminosilicate Particles in Concentrating Solar Power Environments

Chemically inert, aluminosilicate based particles have been investigated as both a thermal transport and sensible energy storage medium for concentrating solar power facilities. These particles will experience a wide range of operating temperatures (300-1000 K) and handling conditions (dense to dilute falling particle curtains, dense granular flows, or dense structures), requiring specially-designed and optimized infrastructures. The relative influence of collisional and frictional interactions between particles varies based on temperature-dependent particulate properties and greatly impacts the bulk, granular flow behavior. These underlying physics are captured using discrete element method modeling tools. However, this modeling method is computationally expensive as each particle position and interaction is tracked during the simulation. These modeling methods are further complicated by introducing temperature-dependent particle properties, high-temperature radiative exchange, and directional irradiation sources experienced by granular flows in concentrating solar power environments. In this study, coupled experimental and numerical slump testing of aluminosilicate particles was performed and computationally efficient radiative exchange models were evaluated to establish best-practices for discrete element method models for concentrating solar power environments. The three particle types investigated included Carbobead HSP 30 /60, Carbobead CP 30/60, and Granusil 4030. Existing modeling limitations and computationally-efficient multi-modal heat transfer models were evaluated using Aspherix®, a commercial discrete element method software. High-temperature (< 1073 K) slump testing of aluminosilicate particles was performed to investigate the deviation between experimentally-observed and numerically-predicted angles of repose introduced by computation-time reduction practices including the relaxation of the particle elastic modulus and coarse-graining. Coarse-graining is used to use a single modeled particle that is representative of a collection of smaller particles, decreasing the computational cost at the expense of geometric accuracy. Additionally, relaxation of the elastic modulus is used to reduce computational time at the expense of an increased, modeled particle overlap. Prior studies have determined that aluminosilicate particles retain a high elastic modulus at high temperatures (< 1073 K), requiring small simulation timesteps to ensure resolved contact forces resemble appropriate solid mechanics. A parametric study was performed to evaluate the influence of computation time improvements on the deviation between experimental and modeled angle of repose across high temperatures < 1073 K. Additionally, numerical case studies were performed on candidate particle systems at varying porosities and temperatures. These studies were performed to investigate the influence of computationally-efficient radiative-exchange modeling methods coupled to Aspherix® on modeled accuracy and computation time. The recently-developed distance-based approximation was evaluated in estimating radiative exchange between particles and participating surfaces located in close proximity. The distance based approximation was developed to use tabulated estimates of the radiative distribution factor between individual particles and surfaces in close proximity (< 40 particle radii). These methods were expanded to the aluminosilicate particles of interest, including the influence of particle size distributions. To capture radiative exchange between particles and surfaces not in close proximity (> 40 particle radii) and to capture the absorption of directional irradiation from concentrating solar resources, a volumetrically-averaged radiative distribution factor was calculated between the modeled granular flow and surfaces using Monte Carlo ray-tracing for participating media. Volume-averaged absorption and scattering coefficients were predicted using a volumetric discretization of the modeled domain with monodisperse approximations based on geometric optics and experimentally-determined scattering phase functions for aluminosilicate particles.

14 SOLAR ENERGY↗
DOE OSTI2024

Development and Experimental Optimization of High-Temperature Modeling Tools and Methods for Concentrated Solar Power Particle - Systems

A novel, open-source radiative modeling toolset was developed to extend the functionality of particle-based modeling software (e.g. discrete element method (DEM)) to environmental conditions relevant to concentrated solar power applications. This toolset was optimized for deployment on desktop workstations instead of high-performance computing systems, to render such tools more accessible to the research community. Both particle-based modeling and radiative exchange modeling are computationally expensive and often require specialized programming expertise, making these methods cumbersome to use. Recent developments in DEM software by DCS Computing have greatly reduced these challenges, providing a graphical-user-interface based platform and modeling optimization for desktop workstations, HPCs, and cloud computing. The University of Dayton leveraged the experience of DCS Computing in developing a user-friendly, open-source radiative heat transfer expansion for DEM modeling. The University of Dayton DEM+ radiative modeling toolset was developed using a combination of fundamental experimental measurements, modeling, and simplified flow experiments over a range of temperatures and flow conditions. The toolset provides researchers with access to multiple radiative models including an accelerated Monte-Carlo Ray Tracing (application agnostic, highly computationally expensive), an expanded database of distance-based approximations (application limited, computationally light), and a weighted blending of the two methods capable of achieving over 90% reduction in computation time with equivalent accuracy compared to Monte-Carlo Ray Tracing. Through a graphical user interface, users can customize the radiative models to match their desired accuracy and available computational resources, improving access to particle based modeling for the research community. Ceramic sintered bauxite proppants were used in modeling and experimentally as a baseline. Both the radiative heat transfer and flow properties for particulate systems were investigated at elevated temperatures up to 800 °C. The major accomplishments for this work include a verified, open-source radiative modeling toolset to be distributed amongst the research community and the fabrication of three small-scale test facilities to investigate particle behavior and tune DEM flow properties for operation up to 800 °C. The findings have been shared with the research community via conference modeling workshops, deployment of the tools in DCS Computing Aspherix®, and open-source access to the developed radiative modeling tool. The development of next-generation CSP facilities and thermal energy storage systems based on ceramic particles requires providing access to computationally efficient and accurate modeling tools. Particles will experience a wide range of environments (20-800 °C) and handling conditions (dilute curtains or dense packing), requiring specially designed and optimized equipment. Optimizing solid particle physics models and establishing best-practices for particle modeling in CSP environments will assist researchers with designing optimized equipment, accelerating the deployment of more economically-competitive CSP facilities.

14 SOLAR ENERGY↗
NASA NTRS1991

Stagnation point nonequilibrium radiative heating and the influence of energy exchange models

A nonequilibrium radiative heating prediction method has been used to evaluate several energy exchange models used in nonequilibrium computational fluid dynamics methods. The radiative heating measurements from the FIRE II flight experiment supply an experimental benchmark against which different formulations for these exchange models can be judged. The models which predict the lowest radiative heating are found to give the best agreement with the flight data. Examination of the spectral distribution of radiation indicates that despite close agreement of the total radiation, many of the models examined predict excessive molecular radiation. It is suggested that a study of the nonequilibrium chemical kinetics may lead to a correction for this problem.

Hartung, Lin C.↗
NASA NTRS1991

Stagnation Point Nonequilibrium Radiative Heating and the Influence of Energy Exchange Models

A nonequilibrium radiative heating prediction method has been used to evaluate several energy exchange models used in nonequilibrium computational fluid dynamics methods. The radiative heating measurements from the FIRE II flight experiment supply an experimental benchmark against which different formulations for these exchange models can be judged. The models which predict the lowest radiative heating are found to give the best agreement with the flight data. Examination of the spectral distribution of radiation indicates that despite close agreement of the total radiation, many of the models examined predict excessive molecular radiation. It is suggested that a study of the nonequilibrium chemical kinetics may lead to a correction for this problem.

Hartung, Lin C.↗
NASA NTRS1989

Comparison of the measured and predicted response of the Earth Radiation Budget Experiment active cavity radiometer during solar observations

The Earth Radiation Budget Experiment consists of an array of radiometric instruments placed in earth orbit by the National Aeronautics and Space Administration to monitor the longwave and visible components of the earth's radiation budget. Presented is a dynamic electrothermal model of the active cavity radiometer used to measure the earth's total radiative exitance. Radiative exchange is modeled using the Monte Carlo method and transient conduction is treated using the finite element method. Also included is the feedback circuit which controls electrical substitution heating of the cavity. The model is shown to accurately predict the dynamic response of the instrument during solar calibration.

Mahan, J. R.↗
NASA NTRS2010

Orion Active Thermal Control System Dynamic Modeling Using Simulink/MATLAB

This paper presents dynamic modeling of the crew exploration vehicle (Orion) active thermal control system (ATCS) using Simulink (Simulink, developed by The MathWorks). The model includes major components in ATCS, such as heat exchangers and radiator panels. The mathematical models of the heat exchanger and radiator are described first. Four different orbits were used to validate the radiator model. The current model results were compared with an independent Thermal Desktop (TD) (Thermal Desktop, PC/CAD-based thermal model builder, developed in Cullimore & Ring (C&R) Technologies) model results and showed good agreement for all orbits. In addition, the Orion ATCS performance was presented for three orbits and the current model results were compared with three sets of solutions- FloCAD (FloCAD, PC/CAD-based thermal/fluid model builder, developed in C&R Technologies) model results, SINDA/FLUINT (SINDA/FLUINT, a generalized thermal/fluid network-style solver ) model results, and independent Simulink model results. For each case, the fluid temperatures at every component on both the crew module and service module sides were plotted and compared. The overall agreement is reasonable for all orbits, with similar behavior and trends for the system. Some discrepancies exist because the control algorithm might vary from model to model. Finally, the ATCS performance for a 45-hr nominal mission timeline was simulated to demonstrate the capability of the model. The results show that the ATCS performs as expected and approximately 2.3 lb water was consumed in the sublimator within the 45 hr timeline before Orion docked at the International Space Station.

Wang, Xiao-Yen J.↗
NASA NTRS1981

A radiometric model of an earth radiation budget radiometer optical system with diffuse-specular surfaces

The Earth Radiation Budget Experiment (ERBE) is to fly on NASA's Earth Radiation Budget Satellite (ERBS) and on NOAA F and NOAA G. Large spatial scale earth energy budget data will be derived primarily from measurements made by the ERBE nonscanning instrument (ERBE-NS). A description is given of a mathematical model capable of simulating the radiometric response of any of the ERBE-NS earth viewing channels. The model uses a Monte Carlo method to accurately account for directional distributions of emission and reflection from optical surfaces which are neither strictly diffuse nor strictly specular. The model computes radiation exchange factors among optical system components, and determines the distribution in the optical system of energy from an outside source. Attention is also given to an approach for implementing the model and results obtained from the implementation.

Luther, M. R.↗
DOE OSTI2021

Progress on Computational Modeling of Water-Based NSTF (FY2021)

This report documents the FY21 progress and achievements made in the computational analyses of the water-based NSTF. Both system-level and high-fidelity Computational Fluid Dynamics (CFD) analyses were performed to gain a complete understanding of the complex flow and heat transfer phenomena in natural convection systems. The progress on the waterbased NSTF experimental testing is summarized in a companion report (ANL-ART-230). As a continuation of progress from previous years, in FY21 the RELAP5 model of the NSTF was first updated to include a representative heated cavity that incorporated convection and conduction means of heat transfer. Two different cavity models were developed, one that models the radiation heat exchange among all surfaces in one enclosure (Cavity Model 1), and another that employs multiple radiation enclosures with each considering the surfaces at one axial level (Cavity Model 2). This updated RELAP5 deck was benchmarked with single-phase test data to tune form loss coefficients of the elbows and tees, along with heat transfer coefficients within the cavity and off the external insulation panels. This tuned model was then found to accurately predict fluid temperatures within 2.5% and both system and riser flow rates within 6% during single-phase, steady-state simulations. For transient simulations, there is a slight over prediction of the flow rate, fluid temperatures, and heater temperatures during the heat-up period, which indicate that the transient heat loss may be underestimated. Following, the capabilities of RELAP5-3D were examined, which feature a conduction enclosure model that is not present in the current version of RELAP5-MOD3.3 used for all previous NSTF simulations. As an initial comparison of the two codes, single-phase simulations were run that used similar input decks. The results indicate that the single-phase flow predicted by the two codes was very similar, however RELAP5-3D appears to predict a higher system flow during the two-phase transient when compared to the flow rates predicted by RELAP5-MOD3.3. Furthermore, RELAP5-3D run did not capture the density wave oscillations or geysering phenomena, which are large oscillation flow instabilities observed in the experimental facility. To further the prediction capabilities for capturing two-phase instabilities, Options 55, 57, 58, 61, which are implemented by default in RELAP5-MOD3.3 but not in RELAP5-3D, were investigated. When enabling Option 61, RELAP5-3D agrees well with MOD3.3 in predicting the total system flow, while other options have negligible effect on the predicted total system flow by RELAP5-3D. However, none of these options enables RELAP5-3D to produce the same flow instabilities. For future steps, efforts are currently underway to investigate RELAP5-3D’s different behavior in predicting low-pressure two-phase instabilities. STAR-CCM+ was used to perform CFD analysis of the water-based test facility and provide simulation results for comparisons to experimental data. Efforts first focused on describing the characteristics of flow within the heated cavity. A study was performed on whether possible leakage flow between the two sides of the cavity could substantially impact and explain the results. An additional second assumption was examined which defined adiabatic side walls, so physical insulation structures were added to the models to allow conduction across this outer insulation. The results suggest that potential porosity of the insulation could be a major factor in the temperature distribution within the cavity and should be considered in any studies where it is desirable to know the cavity wall temperatures to higher accuracy.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2000

Modeling Specular Exchange Between Concentric Cylinders in a Radiative Shielded Furnace

The objective of this research is to develop and validate mathematical models to characterize the thermal performance of a radiative shielded furnace, the University of Alabama in Huntsville (UAH) Isothermal Diffusion Oven. The mathematical models are validated against experimental data obtained from testing the breadboard oven in a terrestrial laboratory environment. It is anticipated that the validation will produce math models capable of predicting the thermal performance of the furnace over a wide range of operating conditions, including those for which no experimental data is available. Of particular interest is the furnace core temperature versus heater power parametric and the transient thermal response of the furnace. Application to a microgravity environment is not considered, although it is conjectured that the removal of any gravity dependent terms from the math models developed for the terrestrial application should yield adequate results in a microgravity environment. The UAH Isothermal Diffusion Oven is designed to provide a thermal environment that is conducive to measuring the diffusion of high temperature liquid metals. In addition to achieving the temperatures required to melt a sample placed within the furnace, reducing or eliminating convective motions within the melt is an important design consideration [1]. Both of these influences are reflected in the design of the furnace. Reducing unwanted heat losses from the furnace is achieved through the use of low conductivity materials and reflective shielding. As evidenced by the highly conductive copper core used to house the sample within the furnace, convective motions can be greatly suppressed by providing an essentially uniform thermal environment. An oven of this design could ultimately be utilized in a microgravity environment, presumably as a experiment payload. Such an application precipitates other design requirements that limit the resources available to the furnace such as power, mass, volume, and possibly even time. Through the experimental and numerical results obtained, the power requirements and thermal response time of the breadboard furnace are quantified.

Schunk, Richard Gregory↗
NASA NTRS1980

The exchange of infrared radiative energy in the troposphere

On the basis of an atmospheric condition of a six-layer general circulation model of the National Center for Atmospheric Research, the spatial and the spectral partition of the long-wave radiation energy exchange in the model atmosphere is investigated and its implications are discussed. It is found that most of the energy is exchanged between adjacent layers and through the wings of strong lines. The strong line approximation is found to be a good approximation for water vapor transmittance. Cooling to space, however, is found to be a poor approximation for the total cooling rate. Both the strong line and cooling to space approximation are discussed in the content of the energy exchange findings.

Wu, M.-L. C.↗
NASA NTRSDate not supplied

Run Time Improvement Efforts for the Roman Space Telescope Thermal Analysis

"Observatory thermal models for large, complex missions, such as the Wide Field InfraRed Survey Telescope (WFIRST) mission, produce an immense amount of data to be processed. Configuration management of the model throughout the project life cycle has mainly focused on which versions of the subsystem models form the current observatory level configuration. However, the results produced by the model are not nearly as well The Roman Space Telescope (RST), formerly known as the Wide Field InfraRed Survey Telescope, is the next great astrophysics observatory mission to follow the James Webb Space Telescope with a planned launch in 2026. As a large scale, flagship mission for NASA with challenging wave front error stability requirements, a single model approach for both thermal discipline analysis and thermo-optical distortion analysis has been used since the early days of the project. In alleviating the need to maintain two separate models for different analysis types, it imposes run time penalties on the thermal analysis with a large model with significant radiation heat exchange. Throughout the lifecycle of the project, the component models have steadily grown in size, resulting in a continuous growth of the overall observatory model with each update and consequently a considerable increase in the model run time. While ongoing efforts to reduce run time are continuously investigated, previous efforts had primarily focused on timestep size and total simulation time to reach quasi-equilibrium. More recently, studies were performed on the total number of radiation couplings (radks) included in the model, which has a nearly linear impact on run time, but increases exponentially with node count. As standard practice for spacecraft analysis, small radks were excluded from the temperature solution based on the assumption that their interchange/view factors have a negligible impact on heat flow. Four approaches were investigated to reduce the model run time while minimizing the impact on accuracy: (1) the Equivalent Radiation Network node, (2) Progressive Radk Inclusion as solution proceeds, (3) Targeted Radk Filtering for critical/non critical areas, and lastly (4) Representation of culled radks with Backloads. Furthermore, the investigation of model run time also revealed that cold cases took noticeably longer to run than hot cases; the root computational inefficiencies were explored along with the computation penalty of linearization of the external radks and recalculation of temperature dependent linear couplings at each timestep. This paper outlines the details of each of the above approaches and their impact on run time and model accuracy.

Thermal Analysis↗
NASA NTRSDate not supplied

Run Time Improvement Efforts for the Roman Space Telescope Thermal Analysis

Observatory thermal models for large, complex missions, such as the Wide Field InfraRed Survey Telescope (WFIRST) mission, produce an immense amount of data to be processed. Configuration management of the model throughout the project life cycle has mainly focused on which versions of the subsystem models form the current observatory level configuration. However, the results produced by the model are not nearly as well The Roman Space Telescope (RST), formerly known as the Wide Field InfraRed Survey Telescope, is the next great astrophysics observatory mission to follow the James Webb Space Telescope with a planned launch in 2026. As a large scale, flagship mission for NASA with challenging wave front error stability requirements, a single model approach for both thermal discipline analysis and thermo-optical distortion analysis has been used since the early days of the project. In alleviating the need to maintain two separate models for different analysis types, it imposes run time penalties on the thermal analysis with a large model with significant radiation heat exchange. Throughout the lifecycle of the project, the component models have steadily grown in size, resulting in a continuous growth of the overall observatory model with each update and consequently a considerable increase in the model run time. While ongoing efforts to reduce run time are continuously investigated, previous efforts had primarily focused on timestep size and total simulation time to reach quasi-equilibrium. More recently, studies were performed on the total number of radiation couplings (radks) included in the model, which has a nearly linear impact on run time, but increases exponentially with node count. As standard practice for spacecraft analysis, small radks were excluded from the temperature solution based on the assumption that their interchange/view factors have a negligible impact on heat flow. Four approaches were investigated to reduce the model run time while minimizing the impact on accuracy: (1) the Equivalent Radiation Network node, (2) Progressive Radk Inclusion as solution proceeds, (3) Targeted Radk Filtering for critical/non critical areas, and lastly (4) Representation of culled radks with Backloads. Furthermore, the investigation of model run time also revealed that cold cases took noticeably longer to run than hot cases; the root computational inefficiencies were explored along with the computation penalty of linearization of the external radks and recalculation of temperature dependent linear couplings at each timestep. This paper outlines the details of each of the above approaches and their impact on run time and model accuracy.

Thermal Analysis↗
DOE OSTI2023

Coupled Experimental and High-Temperature Discrete-Element Method Modeling Studies of Aluminosilicate Particle Handling in Concentrated Solar Power Environments

Chemically inert, aluminosilicate based particles have been investigated as both a thermal transport and sensible energy storage medium for concentrated solar power facilities. These particles will experience a wide range of operating temperatures (300-1000 K) and handling conditions (dense to dilute falling particle curtains, dense granular flows, or dense structures), requiring specially-designed and optimized infrastructures. The relative influence of collisional and frictional interactions between particles varies based on temperature-dependent particulate properties and greatly impacts the bulk, granular flow behavior. These underlying physics are captured using discrete element method modeling tools. However, this modeling method is computationally expensive as each particle position and interaction is tracked during the simulation. These modeling methods are further complicated by introducing temperature-dependent particle properties, high-temperature radiative exchange, and directional irradiation sources experienced by granular flows in concentrated solar power environments. Coupled experimental and numerical studies of aluminosilicate particles in rotary kilns and dense particle curtains were performed for bulk temperatures up to 1073 K. The three particle types investigated included Carbobead HSP 30 /60, Carbobead CP 30/60, and Granusil 4030. Temperature, spatial, and velocity profile data were extracted from experimental runs using embedded K-type thermocouple probes and particle image velocimetry techniques. Experimental and numerical studies were compared using spatial temperature profiles, velocity fields, and shape profiles of the bulk, granular flows. Numerical models were developed using commercially available discrete element method modeling software, Aspherix®. Existing Aspherix® functionality was expanded by introducing coupled radiative exchange modeling tools. The laboratory-scale rotary kiln was developed to investigate the steady-state heat and mass transfer performance of aluminosilicate particles based on particle type, bulk handling temperature, and wall roughness. The rotational speed of the rotary kiln was varied to control the relative impact of collisional and frictional effects upon the granular flow behavior. Heat and mass transfer performance was categorized based on the Froude number and the observed flow regimes of slipping, rolling, cascading, and centrifuging. Coupled discrete element method modeling studies were used to evaluate the effects of temperature-dependent, particulate mechanical properties upon bulk flow behavior and upon the relative effects of radiative, advective, and/or conductive heat transfer. A high-temperature (< 1073 K) falling particle curtain was similarly fabricated to investigate the heat and mass transfer performance of aluminosilicate particles in particle handling situations dominated by inter-particle collisions. The impact of particle type, flow preheat temperatures (< 1073K), and bulk mass flow rates were investigated upon the particle curtain shape, temperature, and velocity profiles. Coupled discrete element method modeling studies were performed to evaluate the varying impact of temperature-dependent, particulate mechanical properties on the bulk flow behavior and the temperature profile of the particle curtain.

14 SOLAR ENERGY↗
NASA NTRS2001

Probabilities of radiation-induced inter- and intrachromosomal exchanges and their dependence on the DNA content of the chromosome

A biophysical model has been developed that is based on the assumptions that an interphase chromosome occupies a spherical territory and that chromosome exchanges are formed by the misrejoining of two DNA double-strand breaks induced within a defined interaction distance. The model is used to explain the relative frequencies of inter- and intrachromosomal exchanges and the relationship between radiation-induced aberrations in individual chromosomes and the DNA content of the chromosome. Although this simple model predicts a higher ratio of inter- to intrachromosomal exchanges for low-LET radiation than for high-LET radiation, as has been suggested by others, we argue that the comparison of the prediction of the model with experimental results is not straightforward. With the model, we also show that the probability of the formation of interchromosomal exchanges is proportional to the "surface area" of the chromosome domain plus a correction term. The correction term is small if the interaction distance is less than 1 microm for both low- and high-LET radiations.

NASA Center JSC↗
NASA NTRS1985

Interface shape in Czochralski grown crystals - Effect of conduction and radiation

A sequential modular computational scheme has been proposed for the calculation of temperature profiles in the melt and crystal including the melt-crystal interface shape and the crystal pulling rate for a fixed crystal radius. The heat transfer in the melt is assumed to occur by conduction only. The shape of melt-gas meniscus, described by Laplace-Young equation, is incorporated in the model and its effects on interface shape and pulling rate are examined. Further, the model incorporates the detailed radiation interaction among the various surfaces in the puller using the Gebhart enclosure theory which accounts for both direct and reflected radiation. The effects of various process parameters on the interface shape and the pulling rate have been investigated. The results of this study indicate that the shape of the melt-gas meniscus has a significant effect on the pulling rate and the interface shape in conformity with earlier studies. Further, it is shown that a simple radiation model (Stefan's model) is inadequate to model the radiation heat exchange in the Czochralski puller apparatus and the detailed Gebhart analysis is necessary for accurate calculation of both the temperature profile in the crystal as well as the interface shape.

Srivastava, R. K.↗