Search NASA⌕ Search

SEARCH · Search NASA

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 127 records · Page 7

Time-Dependent Photoionization Modeling of Warm Absorbers: High-Resolution Spectra and Response to Flaring Light Curves

Time-dependent photoionization modeling of warm absorber (WA) outflows in active galactic nuclei can play an important role in understanding the interaction between WAs and the central black hole. The WA may be out of the equilibrium state because of the variable nature of the central continuum. In this paper, with the help of time dependent photoionization modeling, we study how the WA gas properties change with time and how it reacts to changing radiation fields. Incorporating a flaring incident light curve, we investigate the behavior of WAs using a photoionization code that simultaneously and consistently solves the time-dependent equations of level population, heating and cooling, and radiative transfer. We simulate the physical processes in the gas clouds, such as ionization, recombination, heating, cooling, and the transfer of ionizing radiation through the cloud and present high-resolution time-resolved absorption spectra. We demonstrate that time-dependent radiative transfer is important and the calculations that omit this effect quantitatively and systematically underestimate the absorption. Time-dependent photoionization models provide crucial insights into the characteristics of WAs and can be used to constrain their density and spatial distribution.

Dev R. Sadaula↗

3D Radiative Transfer in Eta Carinae: Application of the SimpleX Algorithm to 3D SPH Simulations of Binary Colliding Winds

Eta Carinae is an ideal astrophysical laboratory for studying massive binary interactions and evolution, and stellar wind-wind collisions. Recent three-dimensional (3D) simulations set the stage for understanding the highly complex 3D flows in Eta Car. Observations of different broad high- and low-ionization forbidden emission lines provide an excellent tool to constrain the orientation of the system, the primary's mass-loss rate, and the ionizing flux of the hot secondary. In this work we present the first steps towards generating synthetic observations to compare with available and future HST/STIS data. We present initial results from full 3D radiative transfer simulations of the interacting winds in Eta Car. We use the SimpleX algorithm to post-process the output from 3D SPH simulations and obtain the ionization fractions of hydrogen and helium assuming three different mass-loss rates for the primary star. The resultant ionization maps of both species constrain the regions where the observed forbidden emission lines can form. Including collisional ionization is necessary to achieve a better description of the ionization states, especially in the areas shielded from the secondary's radiation. We find that reducing the primary's mass-loss rate increases the volume of ionized gas, creating larger areas where the forbidden emission lines can form. We conclude that post processing 3D SPH data with SimpleX is a viable tool to create ionization maps for Eta Car.

3D radiative transfer simulations↗

Radiative transfer model for remote sensing of suspended sediments in water

A Monte Carlo simulation model of radiative transfer in turbid water is discussed. The model can be used to calculate characteristics of the backscattered signal from an illuminated body of water as a function of the turbidity level and spectral properties of the suspended particulates. The dependence of remote sensing applications on the concentration and spectral properties of sediments in the environmental waters is considered in terms of the model. Attention is directed to the effects of various inputs for the volume-scattering function on backscattered radiance from natural waters. The wavelength dependence of single scattering albedo is investigated.

Ghovanlou, A. H.↗

Analysis and Testing of High Temperature Fibrous Insulation for Reusable Launch Vehicles

Analytical models were developed to model the heat transfer through high-temperature fibrous insulation used in metallic thermal protection systems on reusable launch vehicles. The optically thick approximation was used to simulate radiation heat transfer through the insulation. Different models for gaseous conduction and solid conduction in the fibers, and for combining the various modes of heat transfer into a local, volume-averaged, thermal conductivity were considered. The governing heat transfer equations were solved numerically, and effective thermal conductivities were calculated from the steady-state results. An experimental apparatus was developed to measure the apparent thermal conductivity of insulation subjected to pressures, temperatures and temperature gradients representative of re-entry conditions for launch vehicles. The apparent thermal conductivity of an alumina fiber insulation was measured at nominal densities of 24, 48 and 96 kg/cu m. Data were obtained at environmental pressures from 10(exp 4) to 760 torr, with the insulation cold side maintained at room temperature and its hot side temperature varying up to 1000 C. The experimental results were used to evaluate the analytical models. The best analytical model resulted in effective thermal conductivity predictions that were within 8% of experimental results.

Daryabeigi, Kamran↗

Reinterpretation of Mariner 9 IRIS data on the basis of a simulation of radiative-conductive convective transfer in the dust laden Martian atmosphere

Time dependent vertical distributions of atmospheric temperature and static stability were determined by a radiative-convective-conductive heat transfer model attuned to Mariner 9 IRIS radiance data. Of particular interest were conditions of both the dust-laden and dust-free atmosphere in the middle latitudes on Mars during the late S.H. summer season. The numerical model simulates at high spatial and temporal resolution (52 atmospheric and 30 subsurface levels; with a time-step of 7.5 min.) the heat transports in the ground-atmosphere system. The algorithm is based on the solution of the appropriate heating rate equation which includes radiative, molecular-conductive and convective heat transfer terms. Ground and atmosphere are coupled by an internal thermal boundary condition.

Pallman, A. J.↗

Simulations of microphysical, radiative, and dynamical processes in a continental-scale forest fire smoke plume

The impact of a large forest fire smoke plume on atmospheric processes is studied through a numerical model of meteorology, aerosols, and radiative transfer. The simulated smoke optical depths at 0.63-micron wavelength are in agreement with analyses of satellite data and show values as high as 1.8. The smoke has an albedo of 35 percent, or more than double the clear-sky value, and cools the surface by as much as 5 K. An imaginary refractive index, n sub im, of 0.01 yields results which closely match the observed cooling, single scattering albedo, and the Angstrom wavelength exponent. An n exp im of 0.1, typical of smoke from urban fires, produces 9 K cooling. Coagulation causes the geometric mean radius by number to increase from the initial value of 0.08 micron to a final value of 0.15 micron, while the specific extinction and absorption increase by 40 and 25 percent, respectively.

Westphal, Douglas L.↗

Incorporating Endmember Variability into Spectral Mixture Analysis Through Endmember Bundles

Variation in canopy structure and biochemistry induces a concomitant variation in the top-of-canopy spectral reflectance of a vegetation type. Hence, the use of a single endmember spectrum to track the fractional abundance of a given vegetation cover in a hyperspectral image may result in fractions with considerable error. One solution to the problem of endmember variability is to increase the number of endmembers used in a spectral mixture analysis of the image. For example, there could be several tree endmembers in the analysis because of differences in leaf area index (LAI) and multiple scatterings between leaves and stems. However, it is often difficult in terms of computer or human interaction time to select more than six or seven endmembers and any non-removable noise, as well as the number of uncorrelated bands in the image, limits the number of endmembers that can be discriminated. Moreover, as endmembers proliferate, their interpretation becomes increasingly difficult and often applications simply need the aerial fractions of a few land cover components which comprise most of the scene. In order to incorporate endmember variability into spectral mixture analysis, we propose representing a landscape component type not with one endmember spectrum but with a set or bundle of spectra, each of which is feasible as the spectrum of an instance of the component (e.g., in the case of a tree component, each spectrum could reasonably be the spectral reflectance of a tree canopy). These endmember bundles can be used with nonlinear optimization algorithms to find upper and lower bounds on endmember fractions. This approach to endmember variability naturally evolved from previous work in deriving endmembers from the data itself by fitting a triangle, tetrahedron or, more generally, a simplex to the data cloud reduced in dimension by a principal component analysis. Conceptually, endmember variability could make it difficult to find a simplex that both surrounds the data cloud and has vertices that are realistic endmember spectra with reflectances between 0 and 1. In this paper, we create endmember bundles and bounding fraction images for an AVIRIS subscene simulated with a plant canopy radiative transfer model. The simulated subscene is spatially patterned after a subscene from the AVIRIS image acquired August, 1993 over La Copita, Texas. In addition, for comparison, we performed a traditional unmixing with image endmembers.

Bateson, C. Ann↗

Modeling Mg II During Solar Flares. I. Partial Frequency Redistribution, Opacity, and Coronal Irradiation

The Interface Region Imaging Spectrograph has routinely observed the flaring Mg II near-ultraviolet (NUV) spectrum, offering excellent diagnostic potential and a window into the location of energy deposition. A number of studies have forward-modeled both the general properties of these lines and specific flare observations. Generally these have forward-modeled radiation via post-processing of snapshots from hydrodynamic flare simulations through radiation transfer codes. There has, however, not been a study of how the physics included in these radiation transport codes affects the solution. A baseline setup for forward-modeling Mg II in flares is presented and contrasted with approaches that add or remove complexity. It is shown for Mg II that (1) partial frequency distribution (PRD) is still required during flare simulations despite the increased densities; (2) using full angle-dependent PRD affects the solution but takes significantly longer to process a snapshot; (3) including Mg I in non-LTE (NLTE) results in negligible differences to the Mg II lines but does affect the NUV quasi-continuum; (4) only hydrogen and Mg II need to be included in NLTE; (5) ideally the nonequilibrium hydrogen populations, with nonthermal collisional rates, should be used rather than the statistical equilibrium populations; (6) an atom consisting of only the ground state, h and k upper levels, and continuum level is insufficient to model the resonance lines; and (7) irradiation from a hot, dense flaring transition region can affect the formation of Mg II. We discuss modifications to the RH code allowing straightforward inclusion of the transition region and coronal irradiation in flares.

Kerr, Graham S.↗

Prospects of Detecting Nonthermal Protons in Solar Flares via Lyman Line Spectroscopy: Revisiting the Orrall-Zirker Effect

Solar flares are efficient particle accelerators, with a substantial fraction of the energy released manifesting as nonthermal particles. While the role that nonthermal electrons play in transporting flare energy is well studied, the properties and importance of nonthermal protons are rather less well understood. This is in large part due to the paucity of diagnostics, particularly at the lower-energy (deka-keV) range of nonthermal proton distributions in flares. One means to identify the presence of deka-keV protons is by an effect originally described by Orrall & Zirker. In the Orrall–Zirker effect, nonthermal protons interact with ambient neutral hydrogen, and via charge exchange produce a population of energetic neutral atoms (ENAs) in the chromosphere. These ENAs subsequently produce an extremely redshifted photon in the red wings of hydrogen spectral lines. We revisit predictions of the strength of this effect using modern interaction cross sections, and numerical models capable of self-consistently simulating the flaring nonequilibrium ionization stratification, and the nonthermal proton distribution (and, crucially, their feedback on each other). We synthesize both the thermal and nonthermal emission from Ly α and Ly β, the most promising lines that may exhibit a detectable signal. These new predictions are weaker and more transient than prior estimates, but the effects should be detectable in fortuitous circumstances. We degrade the Ly β emission to the resolution of the Spectral Imaging of the Coronal Environment (SPICE) instrument on board Solar Orbiter, demonstrating that though likely difficult, it should be possible to detect the presence of nonthermal protons in flares observed by SPICE.

Solar flares↗

VIEW: A modification of the RAVFAC view factor program for use with the NASTRAN thermal analyzer

The NASTRAN thermal analyzer will include the capability to perform complete thermal analyses on structures. One of the inputs to NASTRAN required to simulate radiative heat transfer between surfaces will be the view factors (also called shape factors, form factors, configuration factors) between those surfaces. The purpose of the VIEW program is to compute these view factors and produce appropriate RADMTX and RADLST output to be used as NASTRAN bulk data.

Puccinelli, E. F.↗

Atmospheric effect on spectral signature - Measurements

In order to study the spectral atmospheric effect on remote sensing, measurements of the spectral signature of surface cover were conducted during hazy conditions. Simultaneous measurements of the aerosol optical thickness and its vertical distribution were carried out. The upward radiance in high flight were also simulated by radiative transfer computations, and the results were compared with the measurements. It is suggested that the radiances over dark areas can be used to derive the aerosol optical thickness. Combined with climatological information, the derived optical thickness can be used to perform corrections of the atmospheric effect.

Kaufman, Yoram J.↗

Remote sensing of precipitation structures using combined microwave radar and radiometric techniques

After reviewing published multiparameter radar observations and their interpretation for microphysical retrieval, the paper presents radiative transfer model simulations performed using radar measurements and a cloud model. The qualitative and quantitative aspects of the remote sensing of precipitation are discussed with particular attention given to the potential of estimating scattering-based ice water path.

Vivekanandan, J.↗

Coupling Satellite and Ground-Based Instruments to Map Climate Forcing by Anthropogenic Aerosols

Climate forcing by anthropogenic aerosols is a significant but highly uncertain factor in global climate change. Only satellites can offer the global coverage essential to reducing this uncertainty; however, satellite measurements must be coupled with correlative, in situ measurements both to constrain the aerosol optical properties required in satellite retrieval algorithms and to provide chemical identification of aerosol sources. This grant funded the first two years of a three-year project which seeks to develop methodologies for combining spaceborne lidar with in-situ aerosol data sets to improve estimates of direct aerosol climate forcing. Progress under this two-year grant consisted in the development and deployment of a new in-situ capability for measuring aerosol 180' backscatter and the extinction-to-backscatter ratio. This new measurement capacity allows definitive lidar/in-situ comparisons and improves our ability to interpret lidar data in terms of climatically relevant quantities such as the extinction coefficient and optical depth. Measurements were made along the coast of Washington State, in Central Illinois, over the Indian Ocean, and in the Central Pacific. Thus, this research, combined with previous measurements by others, is rapidly building toward a global data set of extinction-to-backscatter ratio for key aerosol types. Such information will be critical to interpreting lidar data from the upcoming PICASSO-CENA, or P-C, satellite mission. Another aspect of this project is to investigate innovative ways to couple the lidar-satellite signal with targeted in-situ measurements toward a direct determination of aerosol forcing. This aspect is progressing in collaboration with NASA Langley's P-C lidar simulator and radiative transfer modeling by the University of Lille, France.

Charlson, Robert J.↗

Toward an Improved Single-Particle Model for Large Irregular Grains

To interpret remote spectral observations, scattering and absorption in a particulate surface are simulated via radiative transfer models. The standard model for this purpose among the planetary science community is the Hapke model. This model (like many others) uses two parameters to characterize the optical behavior of individual grains in a particulate surface, the single-scattering albedo omega and phase function p(g). These terms describe, respectively, the quantity and the angular distribution of light scattered by an individual grain. Unfortunately, these parameters are strictly optical. They can be rather difficult to interpret in terms of more interesting particle properties such as grain sizes, shapes, and compositions, that a remote sensing experiment might seek to discover. An equivalent slab approximation is typically used to relate omega to the grain size and optical constants of the material. This approach can mimic the wavelength-dependent absorption behavior of irregular grains, as long as the imaginary index kappa is much less than 1, the shape is equant, and the grain size D is much larger than the wavelength lambda. Unfortunately, the equivalent slab approach provides no information about p(g), which also has a strong dependence on optical constants and particle form.

Grundy, W. M.↗

A Model for Estimation of Rain Rate on Tropical Land from TRMM Microwave Imager Radiometer Observations

Over the tropical land regions observations of the 85 GHz brightness temperature (T(sub 85v)) made by the TRMM Microwave Imager (TMI) radiometer when analyzed with the help of rain rate (R(sub pR)) deduced from the TRMM Precipitation Radar (PR) indicate that there are two maxima in rain rate. One strong maximum occurs when T(sub 85) has a value of about 220 K and the other weaker one when T(sub 85v) is much colder approx. 150 K. Together with the help of earlier studies based on airborne Doppler Radar observations and radiative transfer theoretical simulations, we infer the maximum near 220 K is a result of relatively weak scattering due to super cooled rain drops and water coated ice hydrometeors associated with a developing thunderstorm (Cb) that has a strong updraft. The other maximum is associated with strong scattering due to ice particles that are formed when the updraft collapses and the rain from the Cb is transit2oning from convective type to stratiform type. Incorporating these ideas and with a view to improve the estimation of rain rate from existing operational method applicable to the tropical land areas, we have developed a rain retrieval model. This model utilizes two parameters, that have a horizontal scale of approx. 20km, deduced from the TMI measurements at 19, 21 and 37 GHz (T(sub 19v), T(sub 21v), T(sub 37v). The third parameter in the model, namely the horizontal gradient of brightness temperature within the 20 km scale, is deduced from TMI measurements at 85 GHz. Utilizing these parameters our retrieval model is formulated to yield instantaneous rain rate on a scale of 20 km and seasonal average on a mesoscale that agree well with that of the PR.

Prabhakara, C.↗

Higher-Order Finite Elements for Computing Thermal Radiation

Two variants of the finite-element method have been developed for use in computational simulations of radiative transfers of heat among diffuse gray surfaces. Both variants involve the use of higher-order finite elements, across which temperatures and radiative quantities are assumed to vary according to certain approximations. In this and other applications, higher-order finite elements are used to increase (relative to classical finite elements, which are assumed to be isothermal) the accuracies of final numerical results without having to refine computational meshes excessively and thereby incur excessive computation times. One of the variants is termed the radiation sub-element (RSE) method, which, itself, is subject to a number of variations. This is the simplest and most straightforward approach to representation of spatially variable surface radiation. Any computer code that, heretofore, could model surface-to-surface radiation can incorporate the RSE method without major modifications. In the basic form of the RSE method, each finite element selected for use in computing radiative heat transfer is considered to be a parent element and is divided into sub-elements for the purpose of solving the surface-to-surface radiation-exchange problem. The sub-elements are then treated as classical finite elements; that is, they are assumed to be isothermal, and their view factors and absorbed heat fluxes are calculated accordingly. The heat fluxes absorbed by the sub-elements are then transferred back to the parent element to obtain a radiative heat flux that varies spatially across the parent element. Variants of the RSE method involve the use of polynomials to interpolate and/or extrapolate to approximate spatial variations of physical quantities. The other variant of the finite-element method is termed the integration method (IM). Unlike in the RSE methods, the parent finite elements are not subdivided into smaller elements, and neither isothermality nor other unrealistic physical conditions are assumed. Instead, the equations of radiative heat transfer are integrated numerically over the parent finite elements by use of a computationally efficient Gaussian integration scheme.

Gould, Dana C.↗

A Novel Non-Intrusive Method to Resolve the Thermal-Dome-Effect of Pyranometers: Radiometric Calibration and Implications

Traditionally the calibration equation for pyranometers assumes that the measured solar irradiance is solely proportional to the thermopile's output voltage; therefore only a single calibration factor is derived. This causes additional measurement uncertainties because it does not capture sufficient information to correctly account for a pyranometer's thermal effect. In our updated calibration equation, temperatures from the pyranometer's dome and case are incorporated to describe the instrument's thermal behavior, and a new set of calibration constants are determined, thereby reducing measurement uncertainties. In this paper, we demonstrate why a pyranometer's uncertainty using the traditional calibration equation is always larger than a-few-percent, but with the new approach can become much less than 1% after the thermal issue is resolved. The highlighted calibration results are based on NIST-traceable light sources under controlled laboratory conditions. The significance of the new approach lends itself to not only avoiding the uncertainty caused by a pyranometer's thermal effect but also the opportunity to better isolate and characterize other instrumental artifacts, such as angular response and non-linearity of the thermopile, to further reduce additional uncertainties. We also discuss some of the implications, including an example of how the thermal issue can potentially impact climate studies by evaluating aerosol's direct-radiative effect using field measurements with and without considering the pyranometer's thermal effect. The results of radiative transfer model simulation show that a pyranometer's thermal effect on solar irradiance measurements at the surface can be translated into a significant alteration of the calculated distribution of solar energy inside the column atmosphere.

Ji, Qiang↗