Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Radiative transfer simulations”

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 901 records · Page 50

Initial full core SMR simulations with NekRS

This document describes the completion of a recent milestone of the ExaSMR program concerning "Full core simulations". The ExaSMR project is developing tools for coupling Monte Carlo (MC) radiation transport solvers to a computational fluid dynamics (CFD) solver. Work is based on the Shift MC, OpenMC MC, and Nek5000/NekRS CFD codes. As part of this milestone, a novel set of pin-resolved CFD full-core simulations have been performed for the first time. These simulations represent a significant increase in capability in what is now possible with CFD on pre-Exascale systems. The simulations have been performed with the spectral element solvers NekRS on the supercomputer Summit. Multiple simulations campaigns have been conducted: (1) LES simulations in bare bundle, (2) RANS simulations in bare bundles, (3) RANS simulations with momentum sources and (4) Conjugate heat transfer calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

CRTM Microwave Cloud Scattering Lookup Tables and Radar Simulator

Microwave observations play a very important role in improving the weather forecasts. Although these observations are routinely assimilated into NWP models in clear-sky conditions, assimilation of all-sky microwave observations is very limited. Two main factors contributing to this limitation are inaccuracy in the input cloud and hydrometeor profiles used as input to the radiative transfer model and also error in scattering calculations performed by the radiative transfer model itself. The Mie theory is used by many fast RT models to estimate the optical properties of single particles. The Mie theory assumes spherical shapes for ice or snow particles with mixture of air and ice. However, hydrometeors scattering radiation at microwave frequencies have different shapes, sizes, and orientations. Therefore, using Mie theory to determine their optical properties leads to large uncertainties in all-sky radiative transfer calculations. The discrete dipole approximation (DDA) which approximates the optical properties of large objects in terms of discrete dipoles has shown promise in calculating the scattering properties of particles with different shapes in the microwave frequencies. This presentation focuses on recent advancements in the CRTM scattering calculations for frozen hydrometeors in the microwave frequencies using the DDA technique. In addition to using stand-alone CRTM calculations using collocated ATMS and reanalysis profiles, the data assimilation experiments conducted using the NOAA FV3GFS forecast system are used to evaluate the scattering improvements. Additionally, the backscattering information from the DDA database was used to implement a radar simulator into CRTM. The radar operator takes advantage of CRTM different modules to calculate clouds absorption and scattering properties. In addition to the forward model both adjoint and tangent linear of the radar simulator are implemented and evaluated as well. The radar simulator is currently being tested within the JEDI/GEOS data assimilation framework to facilitate the assimilation of radar measurements such as CloudSat CPR and GPM DPR into the NASA GEOS model.

Isaac Moradi↗

The Uncertainty of Biomass Estimates from Modeled ICESat-2 Returns Across a Boreal Forest Gradient

The Forest Light (FLIGHT) radiative transfer model was used to examine the uncertainty of vegetation structure measurements from NASA's planned ICESat-2 photon counting light detection and ranging (LiDAR) instrument across a synthetic Larix forest gradient in the taiga-tundra ecotone. The simulations demonstrate how measurements from the planned spaceborne mission, which differ from those of previous LiDAR systems, may perform across a boreal forest to non-forest structure gradient in globally important ecological region of northern Siberia. We used a modified version of FLIGHT to simulate the acquisition parameters of ICESat-2. Modeled returns were analyzed from collections of sequential footprints along LiDAR tracks (link-scales) of lengths ranging from 20 m-90 m. These link-scales traversed synthetic forest stands that were initialized with parameters drawn from field surveys in Siberian Larix forests. LiDAR returns from vegetation were compiled for 100 simulated LiDAR collections for each 10 Mg · ha(exp -1) interval in the 0-100 Mg · ha(exp -1) above-ground biomass density (AGB) forest gradient. Canopy height metrics were computed and AGB was inferred from empirical models. The root mean square error (RMSE) and RMSE uncertainty associated with the distribution of inferred AGB within each AGB interval across the gradient was examined. Simulation results of the bright daylight and low vegetation reflectivity conditions for collecting photon counting LiDAR with no topographic relief show that 1-2 photons are returned for 79%-88% of LiDAR shots. Signal photons account for approximately 67% of all LiDAR returns, while approximately 50% of shots result in 1 signal photon returned. The proportion of these signal photon returns do not differ significantly (p greater than 0.05) for AGB intervals greater than 20 Mg · ha(exp -1). The 50m link-scale approximates the finest horizontal resolution (length) at which photon counting LiDAR collection provides strong model fits and minimizes forest structure uncertainty in the synthetic Larix stands. At this link-scale AGB greater than 20 Mg · ha(exp -1) has AGB error from 20-50% at the 95% confidence level. These results suggest that the theoretical sensitivity of ICESat-2 photon counting LiDAR measurements alone lack the ability to consistently discern differences in inferred AGB at 10 Mg · ha(exp -1) intervals in sparse forests characteristic of the taiga-tundra ecotone.

LiDAR↗

Mission Simulation of Space Lidar Measurements for Seasonal and Regional CO2 Variations

Results of mission simulation studies are presented for a laser-based atmospheric [82 sounder. The simulations are based on real-time carbon cycle process modeling and data analysis. The mission concept corresponds to the Active Sensing of [82 over Nights, Days, and Seasons (ASCENDS) recommended by the US National Academy of Sciences Decadal Survey of Earth Science and Applications from Space. One prerequisite for meaningful quantitative sensor evaluation is realistic CO2 process modeling across a wide range of scales, i.e., does the model have representative spatial and temporal gradients? Examples of model comparison with data will be shown. Another requirement is a relatively complete description of the atmospheric and surface state, which we have obtained from meteorological data assimilation and satellite measurements from MODIS and [ALIPS0. We use radiative transfer model calculations, an instrument model with representative errors ' and a simple retrieval approach to complete the cycle from "nature" run to "pseudo-data" CO2, Several mission and instrument configuration options are examined/ and the sensitivity to key design variables is shown. We use the simulation framework to demonstrate that within reasonable technological assumptions for the system performance, relatively high measurement precision can be obtained, but errors depend strongly on environmental conditions as well as instrument specifications. Examples are also shown of how the resulting pseudo - measurements might be used to address key carbon cycle science questions.

Kawa, Stephan↗

Extending radiative transfer models by use of Bayes rule

This paper presents a procedure that extends some existing radiative transfer modeling techniques to problems in atmospheric science where curvature and layering of the medium and dynamic range and angular resolution of the signal are important. Example problems include twilight and limb scan simulations. Techniques that are extended include successive orders of scattering, matrix operator, doubling, Gauss-Seidel iteration, discrete ordinates and spherical harmonics. The procedure for extending them is based on Bayes' rule from probability theory.

Whitney, C.↗

Evaluating Multispectral Snowpack Reflectivity With Changing Snow Correlation Lengths

This study investigates the sensitivity of multispectral reflectivity to changing snow correlation lengths. Matzler's ice-lamellae radiative transfer model was implemented and tested to evaluate the reflectivity of snow correlation lengths at multiple frequencies from the ultraviolet (UV) to the microwave bands. The model reveals that, in the UV to infrared (IR) frequency range, the reflectivity and correlation length are inversely related, whereas reflectivity increases with snow correlation length in the microwave frequency range. The model further shows that the reflectivity behavior can be mainly attributed to scattering rather than absorption for shallow snowpacks. The largest scattering coefficients and reflectivity occur at very small correlation lengths (approximately 10(exp -5 m) for frequencies higher than the IR band. In the microwave range, the largest scattering coefficients are found at millimeter wavelengths. For validation purposes, the ice-lamella model is coupled with a multilayer snow physics model to characterize the reflectivity response of realistic snow hydrological processes. The evolution of the coupled model simulated reflectivities in both the visible and the microwave bands is consistent with satellite-based reflectivity observations in the same frequencies. The model results are also compared with colocated in situ snow correlation length measurements (Cold Land Processes Field Experiment 2002-2003). The analysis and evaluation of model results indicate that the coupled multifrequency radiative transfer and snow hydrology modeling system can be used as a forward operator in a data-assimilation framework to predict the status of snow physical properties, including snow correlation length.

Kang, Do Hyuk↗

SSME Propellant Path Leak Detection

The primary objective of this phase of the investigation is the experimental validation of techniques for detecting and analyzing propellant path external leaks which have a high probability of occurring on the SSME. The selection of candidate detection methods requires a good analytic model for leak plumes which would develop from external leaks and an understanding of radiation transfer through the leak plume. One advanced propellant path leak detection technique is obtained by using state-of-art technology of infrared (IR) thermal imaging systems combined with computer, digital image processing and expert systems for the engine protection. The feasibility of the IR leak plume detection will be evaluated on subscale simulated laboratory plumes to determine sensitivity, signal to noise, and general suitability for the application. The theoretical analysis was undertaken with the objective of developing and testing simple, easy-to-use models to predict the amount of radiation coming from a radiation source, background plate (BP), which can be absorbed, emitted and scattered by the gas leaks.

Roger Crawford↗

Prediction of Ablation Rates from Solid Surfaces Exposed to High Temperature Gas Flow

A mathematical model and a solution algorithm is developed to study the physics of high temperature heat transfer and material ablation and identify the problems associated with the flow of hydrogen gas at very high temperatures and velocities through pipes and various components of Nuclear Thermal Rocket (NTR) motors. Ablation and melting can be experienced when the inner solid surface of the cooling channels and the diverging-converging nozzle of a Nuclear Thermal Rocket (NTR) motor is exposed to hydrogen gas flow at temperatures around 2500 degrees Kelvin and pressures around 3.4 MPa. In the experiments conducted on typical NTR motors developed in 1960s, degradation of the cooling channel material (cracking in the nuclear fuel element cladding) and in some instances melting of the core was observed. This paper presents the results of a preliminary study based on two types of physics based mathematical models that were developed to simulate the thermal-hydrodynamic conditions that lead to ablation of the solid surface of a stainless steel pipe exposed to high temperature hydrogen gas near sonic velocities. One of the proposed models is one-dimensional and assumes the gas flow to be unsteady, compressible and viscous. An in-house computer code was developed to solve the conservations equations of this model using a second-order accurate finite-difference technique. The second model assumes the flow to be three-dimensional, unsteady, compressible and viscous. A commercial CFD code (Fluent) was used to solve the later model equations. Both models assume the thermodynamic and transport properties of the hydrogen gas to be temperature dependent. In the solution algorithm developed for this study, the unsteady temperature of the pipe is determined from the heat equation for the solid. The solid-gas interface temperature is determined from an energy balance at the interface which includes heat transfer from or to the interface by conduction, convection, radiation, and ablation. Two different ablation models are proposed to determine the heat loss from the solid surface due to the ablation of the solid material. Both of them are physics based. Various numerical simulations were carried out using both models to predict the temperature distribution in the solid and in the gas flow, and then predict the ablation rates at a typical NTR motor hydrogen gas temperature and pressure. Solid mass loss rate per foot of a pipe was also calculated from these predictions. The results are presented for fully developed turbulent flow conditions in a sample SS pipe with a 6 inch diameter.

Akyuzlu, Kazim M.↗

Cryogenic System Preliminary Design for a 0.5m-Long, Conduction-Cooled Nb3Sn Undulator Magnet Prototype

There are several NbTi superconducting undulator (SCU) magnets currently in operation at the Advanced Photon Source (APS) at Argonne National Laboratory (ANL). The development on Nb3Sn-based superconducting undulator magnets at APS is underway due to the potential to further enhance the performance of the SCUs. Superconducting undulator magnets need to keep temperature gradients minimized in order to retain thermal and operating current margin. We have designed the thermal links for efficient heat conduction using 3D finite element analysis (FEA) simulation in COMSOL Multiphysics software, which was later used for the calculation of the temperature distribution across a 0.5 m long, conduction-cooled, Nb3Sn undulator magnet prototype that includes both conductive heat transfer and radiative heating components. We have modelled the evolution of the thermal properties of the magnet winding as well as other cold parts during cool-down from ambient temperature with the operation of a SHI RDK-415D cryocooler to examine the estimated time that is needed for the cooldown as well as the baseline temperature we could achieve. A key result was that a maximum coil delta T = 0.11 K along the designed thermal links and a temperature range of 3.42 K-3.53 K for the winding were predicted at steady state. A pair of G10 support rods for carrying the magnet as well as a pair of current leads made of metals were also designed and optimized for the testing system, targeting for an operation temperature of 4 K and a coil current of 950 A DC.

Zhang, Danlu↗

Variance Reduction within Implicit Monte Carlo Thermal Radiation Transport using the Local Importance Function Transform [Slides]

System of equations derived by Fleck and Cummings to model time-dependent thermal radiative transfer (TRT) problems: Solution of radiation specific intensity (I), non-linearly coupled to material internal energy (U); Absorption and re-emission approximated by effective scattering events. Useful for high-energy density physics simulations: i.e., astrophysics, inertial confinement fusion (ICF). Whereas deterministic methods are fully discretized and free of statistical noise, Monte Carlo methods allow for dynamic sampling of the phase space at the price of statistical noise; Monte Carlo methods are also characterized by a slow $\frac{1}{√Ν}$ convergence rate.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Numerical studies of collapsing interstellar clouds

Numerical simulation of the structure and evolution of interstellar clouds was initiated. Steps were taken toward an integrated treatment of the dynamical, thermal, and chemical processes entering model calculations. A detailed study was made of radiative transfer in molecular lines to allow model predictions to be tested against empirical data. The calculations have successfully reproduced and explained several observed cloud properties, including abundances of complex molecular species and the apparent depletion of CO in dense cores.

Bodenheimer, P. H.↗

Accretion disk coronae in high-luminosity systems

We present the results of self-consistent models of Compton-heated accretion disk coronae. The models are calculated using a new method for computing monochromatic radiative transfer n two dimensions. The method splits the radiation into direct and scattered components. The direct radiation is computed by calculating the optical depth along rays, while transfer of the scattered radiation is approximated by flux-limited diffusion. The resulting code agrees with more accurate treatments to within 50%, and is highly efficient, making it practical for use in large hydrodynamic simulations. The coronal models are used to confirm the results of earlier work, and to extend it to higher luminosities. In contrast to earlier work, which found the outer disks to be shadowed by the inner corona at high luminosities, we find our results to form an almost continuous extension of the models at lower luminosities. This is due to the presence of multiply scattered radiation, which acts to partially offset the loss of direct radiation from the central source. Although the analytic methods derived at lower luminosities cannot be used to derive the coronal structure for L/L(sub Edd) approx. greater than 0.1, the results of the models are amenable to semiempirical fits. We also discuss possible observational consequences of the results for coronal veiling and line fluorescence from the disk.

Murray, Stephen D.↗

Passive Thermal Coating Observatory Operating in Low-Earth Orbit (PATCOOL) – Cubesat Design to Test Passive Thermal Coatings in Space

The PATCOOL is a NASA sponsored, University of Florida developed 3U Cu-beSat meant to investigate the feasibility of using a cryogenic selective surface coating as a new, more efficient way of passively cooling components in space. Initial tests on the ground demonstrate that this coating should provide a much higher reflectance of the Sun’s irradiant power than any existing coating, while still providing far-infrared power emission. The ultimate validation of this tech-nology requires on-orbit testing. PATCOOL hosts a 4-sample housing, with the samples shaped as thin cylinders (coin-like). Two samples are coated with state-of-the-art material, while the other pair uses the new coating to be evaluated. The temperatures of all samples during the mission (minimum 72 hours of data col-lection) are measured via thermistors. The samples are connected via thin Kevlar strings to the housing, to minimize heat transfer. The housing is designed to shield the samples from Earth’s thermal radiation, and the CubeSat is attitude stabilized and controlled via a gravity gradient boom, magnetorquers and a reaction wheel set. Thermal Desktop simulations show PATCOOL’s ability to thermally isolate the samples from heat exchanges other than with Sun and deep space, thanks to its thermal design and the chosen attitude profile.

Ojeda, Carlos↗

PAssive Thermal Coating Observatory Operating in Low Earth Orbit (PATCOOL)

The PATCOOL is a NASA sponsored, University of Florida developed 3U CubeSat meant to investigate the feasibility of using a cryogenic selective surface coating as a new, more efficient way of passively cooling components in space. Initial tests on the ground demonstrate that this coating should provide a much higher reflectance of the Sun's irradiant power than any existing coating, while still providing far-infrared power emission. The ultimate validation of this technology requires on-orbit testing. PATCOOL hosts a 4-sample housing, with the samples shaped as thin cylinders (coin-like). Two samples are coated with state-of-the-art material, while the other pair uses the new coating to be evaluated. The temperatures of all samples during the mission (minimum 72 hours of data collection) are measured via thermistors. The samples are connected via thin Kevlar strings to the housing, to minimize heat transfer. The housing is designed to shield the samples from Earth's thermal radiation, and the CubeSat is attitude stabilized and controlled via a gravity gradient boom, magnetorquers and a reaction wheel set. Thermal Desktop simulations show PATCOOL's ability to thermally isolate the samples from heat exchanges other than with Sun and deep space, thanks to its thermal design and the chosen attitude profile

Ojeda, Carlos↗

FORWARD AND INVERSE L-BAND RADIATIVE TRANSFER MODELING OVER THE DRY CHACO, USING SMOS OBSERVATIONS, LAND SURFACE MODELING AND IN SITU DATA

Passive microwave L-band remote sensing is well known for its sensitivity to surface soil moisture over land.The signal is also affected by other dynamic variables such as vegetation and soil salinity. In this research, L-band microwave observations of the Soil Moisture Ocean Salinity (SMOS) mission are used to explore soil surface salinity, moisture and vegetation, in the Argentinean Dry Chaco, an area with possible emerging dryland salinity. A Radiative Transfer Model (RTM) with inclusion of a correction to the soil’s dielectric constant for salinity was used in forward and inverse mode, using either in situ data or Catchment Land Surface Model (CLSM) simulations as RTM input. The forward analysis pointed out shortcomings in the modeled soil moisture and soil surface temperature estimates. The impact of salinity on forward simulations over the Dry Chaco was limited.The RTM inversion using 10 years of SMOS brightness temperature observations resulted in realistic estimates of vegetation and roughness, both with and withouto ptimizing a correction term for the dielectric constant in terms of salinity equivalents. However,the retrieval of the correction term to the dielectric constant was very uncertain and not representative of soil surface salinity, but rather of open water, texture uncertainty or soil moisture bias.

L-band↗

Electromagnetic Time and Wavelength Domain Predictions for Accreting Binary Black Holes

Accreting massive binary black holes in gaseous environments are keymulti-messenger sources for the LISA mission. We will present some ofthe latest predictions of the observational signatures of accreting binary black holes across the electromagnetic spectrum from the fieldand provide some of our own. In particular, we will provide a briefsummary of the progress made in our group to understand these systemstheoretically using high-performance simulations of generalrelativistic magnetohydrodynamics (MHD) and dynamic spacetimes.Simulating these systems realistically is challenging asradiation-coupled MHD must be considered over large dynamic ranges inspace and time. Constraining our view to the circumbinary disk regionfirst, we will show how the evolution duration, accretion disk size,mass ratio, and thermodynamics model affect the structure andvariability of the accretion flow. We particularly emphasize how theseparameters influence the over-density feature, or "lump", which orbitsthe binary near the edge of the cavity, since it is responsible formost of the electromagnetic emission's variability---a key signatureof a system being a binary. Extending to smaller length scales, wewill report on simulations following accretion all the down to theevent horizons so that we may begin to investigate how black hole spinaffects mini-disk dynamics, accretion rate, and jet power. We willalso report on our radiative transfer predictions for how black holespin alters the spectral energy distribution of these systems.Wrapping up, we will itemize the remaining key theoreticaltechnologies and directions needed to be accomplished to betterunderstand these sources, and assess strategies for searching for themin the coming decade.

Scott C Noble↗

Simulation of the present-day atmospheric ozone, odd nitrogen, chlorine and other species using a coupled 2-D model in isentropic coordinates

A zonally averaged two-dimensional model for chemical species in the atmosphere (Olaguer et al., 1990), with coupled dynamics, radiation transfer, and chemistry, is used to study the present-day atmospheric ozone, odd nitrogen, chlorine, other species. The model utilizes all zonally averaged physical equations of momentum, energy, and mass and determines self-consistently the advective and diffusive tranasport parameters from the temperature specific to the period of observation. A comparison of the model results with observations showed good agreement.

Yang, H.↗

Brightness temperature simulations for the physical interpretation of AMSU moisture channels

The use of the Advanced Microwave Sounding Unit (AMSU) moisture imagery for the dynamic, kinematic, and synoptic interpretations requires the availability of idealized 'soundings' of temperature, water vapor, and cloudiness as input to a radiative transfer model for determining how atmospheric conditions control the top-of-atmosphere brightness temperatures. This paper uses idealized vertical configurations of atmospheric water vapor and liquid cloud water to illustrate their effects on the AMSU moisture channels. Calculations with various cloud-drop size distribution functions (SDFs) showed that, at least for the specific SDFs used in this work, the Rayleigh approximation is invalid for mode radii greater than 30-40 microns, for frequencies above 89.0 GHz. It is also invalid above about 20 microns for frequencies above 157 GHz. It is concluded that AMSU moisture channels might not be suitable for detecting convective instability over land using image-based approaches.

Muller, Bradley M.↗