The coupling of radiative transfer to quasi 1-D flows with thermochemical nonequilibrium
NASA Technical Reports Server (NTRS)
Gokcen, Tahir; Park, Chul
1991-01-01
Quasi-one-dimensional nonequilibrium nozzle flows with coupled radiative transfer are considered. The strongly coupled formulation of radiation and flowfield leads to a governing set of integro-differential equations. A fully implicit numerical method using the full matrix inversion or block iteration methods is presented to solve these equations. The nonequilibrium gas model consists of two chemical species, molecular and atomic nitrogen. The thermodynamic state of the gas is described by two temperatures, translational-rotational and vibrational, and the thermal radiation is assumed to be governed by the vibrational temperature. In radiative transfer, gases are assumed to be absorbing and emitting, and a detailed spectral dependency of the absorption coefficient is prescribed for a non-gray gas. The numerical solutions of strongly radiating nonequilibrium flows are presented for both gray and non-gray gases.
NASA Astrophysics Data System (ADS)
Yang, Jun; Leconte, Jérémy; Wolf, Eric T.; Goldblatt, Colin; Feldl, Nicole; Merlis, Timothy; Wang, Yuwei; Koll, Daniel D. B.; Ding, Feng; Forget, François; Abbot, Dorian S.
2016-08-01
An accurate estimate of the inner edge of the habitable zone is critical for determining which exoplanets are potentially habitable and for designing future telescopes to observe them. Here, we explore differences in estimating the inner edge among seven one-dimensional radiative transfer models: two line-by-line codes (SMART and LBLRTM) as well as five band codes (CAM3, CAM4_Wolf, LMDG, SBDART, and AM2) that are currently being used in global climate models. We compare radiative fluxes and spectra in clear-sky conditions around G and M stars, with fixed moist adiabatic profiles for surface temperatures from 250 to 360 K. We find that divergences among the models arise mainly from large uncertainties in water vapor absorption in the window region (10 μm) and in the region between 0.2 and 1.5 μm. Differences in outgoing longwave radiation increase with surface temperature and reach 10–20 W m‑2 differences in shortwave reach up to 60 W m‑2, especially at the surface and in the troposphere, and are larger for an M-dwarf spectrum than a solar spectrum. Differences between the two line-by-line models are significant, although smaller than among the band models. Our results imply that the uncertainty in estimating the insolation threshold of the inner edge (the runaway greenhouse limit) due only to clear-sky radiative transfer is ≈10% of modern Earth’s solar constant (i.e., ≈34 W m‑2 in global mean) among band models and ≈3% between the two line-by-line models. These comparisons show that future work is needed that focuses on improving water vapor absorption coefficients in both shortwave and longwave, as well as on increasing the resolution of stellar spectra in broadband models.
NASA Astrophysics Data System (ADS)
Hawes, Frederick T.; Berk, Alexander; Richtsmeier, Steven C.
2016-05-01
A validated, polarimetric 3-dimensional simulation capability, P-MCScene, is being developed by generalizing Spectral Sciences' Monte Carlo-based synthetic scene simulation model, MCScene, to include calculation of all 4 Stokes components. P-MCScene polarimetric optical databases will be generated by a new version (MODTRAN7) of the government-standard MODTRAN radiative transfer algorithm. The conversion of MODTRAN6 to a polarimetric model is being accomplished by (1) introducing polarimetric data, by (2) vectorizing the MODTRAN radiation calculations and by (3) integrating the newly revised and validated vector discrete ordinate model VDISORT3. Early results, presented here, demonstrate a clear pathway to the long-term goal of fully validated polarimetric models.
Electronic-to-vibrational energy transfer efficiency in the O/1 D/-N2 and O/1 D/-CO systems
NASA Technical Reports Server (NTRS)
Slanger, T. G.; Black, G.
1974-01-01
With the aid of a molecular resonance fluorescence technique, which utilizes optical pumping from the v = 1 level of the ground state of CO by A 1 Pi-X 1 Sigma radiation, a study is made of the efficiency of E-V transfer from O(1 D) to CO. O(1 D) is generated at a known rate by O2 photodissociation at 1470 A in an intermittent mode, and the small modulation of the fluorescent signal associated with CO (v = 1) above the normal thermal background is interpreted in terms of E-V transfer efficiency. The CO (v = 1) lifetime in this system is determined mainly by resonance trapping of the IR fundamental band, and is found to be up to ten times longer than the natural radiative lifetime. For CO, (40 plus or minus 8)% of the O(1 D) energy is converted into vibrational energy. By observing the effect of N2 on the CO (v = 1) fluorescent intensity and lifetime, it is possible to obtain the E-V transfer efficiency for the system O(1 D)-N2 relative to that for O(1 D)-CO. The results indicate that the efficiency for N2 is (83 plus or minus 10)% of that for CO.
Thermal radiation heat transfer.
NASA Technical Reports Server (NTRS)
Siegel, R.; Howell, J. R.
1972-01-01
A comprehensive discussion of heat transfer by thermal radiation is presented, including the radiative behavior of materials, radiation between surfaces, and gas radiation. Among the topics considered are property prediction by electromagnetic theory, the observed properties of solid materials, radiation in the presence of other modes of energy transfer, the equations of transfer for an absorbing-emitting gas, and radiative transfer in scattering and absorbing media. Also considered are radiation exchange between black isothermal surfaces, radiation exchange in enclosures composed of diffuse gray surfaces and in enclosures having some specularly reflecting surfaces, and radiation exchange between nondiffuse nongray surfaces. The use of the Monte Carlo technique in solving radiant-exchange problems and problems of radiative transfer through absorbing-emitting media is explained.
Atmospheric Radiative Transfer
NASA Astrophysics Data System (ADS)
Perliski, Lori
Because radiative transfer cuts across many scientific disciplines with applications including remote sensing, climate, atmospheric chemistry, and photobiology, there is a need for comprehensive books on this subject that can appeal to a wide readership. While Atmospheric Radiative Transfer takes strides toward filling this niche by addressing a broad range of topics, it is dry reading and suffers from lack of detail. The book was based on a graduate-level course taught at the University of Sciences and Technologies in Lille, France, and indeed, the text reads much like an expanded outline perhaps derived from lecture notes.Part one deals with general radiative transfer, and part two covers Earth's radiation budget, the climate system, and remote sensing techniques. The radiative transfer equation and solutions for absorbing and scattering atmospheres are discussed as are the details of absorption, such as energy levels, line strengths, line intensities, equivalent widths, and weak- and strong-line limits.
NASA Astrophysics Data System (ADS)
Kalkofen, Wolfgang
2009-07-01
Preface; Introduction; Part I. Operator Perturbation: 1. Survey of operator perturbation methods W. Kalkofen; 2. Line formation in expanding atmospheres: multilevel calculations using approximate lambda operators W. R. Hamann; 3. Stellar atmospheres in non-LTE: model construction and line formation calculations using approximate lambda operators K. Werner; 4. Acceleration of convergence L. H. Auer; 5. Line formation in a time-dependent atmosphere W. Kalkofen; 6. Iterative solution of multilevel transfer problems Eugene H. Avrett and Rudolf Loeser; 7. An algorithm for the simultaneous solution of thousands of transfer equations under global constraints Lawrence S. Anderson; 8. Operator perturbation for differential equations W. Kalkofen; Part II. Polarised Radiation: 9. A gentle introduction to polarised radiative transfer David E. Rees; 10. Non-LTE polarised radiative transfer in special lines David E. Rees and Graham A. Murphy; 11. Transfer of polarised radiation using 4x4 matrices E. Landi Degli'Innocenti; 12. Radiative transfer in the presence of strong magnetic fields A. A. van Ballegooijen; 13. An integral operator technique of radiative transfer in spherical symmetry A. Peraiah; 14. Discrete ordinate matrix method M. Schmidt and R. Wehrse.
NASA Astrophysics Data System (ADS)
Kalkofen, Wolfgang
1988-01-01
Preface; Introduction; Part I. Operator Perturbation: 1. Survey of operator perturbation methods W. Kalkofen; 2. Line formation in expanding atmospheres: multilevel calculations using approximate lambda operators W. R. Hamann; 3. Stellar atmospheres in non-LTE: model construction and line formation calculations using approximate lambda operators K. Werner; 4. Acceleration of convergence L. H. Auer; 5. Line formation in a time-dependent atmosphere W. Kalkofen; 6. Iterative solution of multilevel transfer problems Eugene H. Avrett and Rudolf Loeser; 7. An algorithm for the simultaneous solution of thousands of transfer equations under global constraints Lawrence S. Anderson; 8. Operator perturbation for differential equations W. Kalkofen; Part II. Polarised Radiation: 9. A gentle introduction to polarised radiative transfer David E. Rees; 10. Non-LTE polarised radiative transfer in special lines David E. Rees and Graham A. Murphy; 11. Transfer of polarised radiation using 4x4 matrices E. Landi Degli'Innocenti; 12. Radiative transfer in the presence of strong magnetic fields A. A. van Ballegooijen; 13. An integral operator technique of radiative transfer in spherical symmetry A. Peraiah; 14. Discrete ordinate matrix method M. Schmidt and R. Wehrse.
1D problems of radiation pressure on elastic solids
NASA Astrophysics Data System (ADS)
Požar, Tomaž; Možina, Janez
2015-08-01
We treat the light-matter interaction due to radiation pressure in one dimension using the fundamental, nonrelativistic conservation principles of energy and momentum. Additionally, we assume that the center of mass-energy maintains the same uniform motion if the interaction takes place or not. Since we handle solids as elastic objects, the results are consistent with the principle of causality and agree with recent experimental observations. We analyze the problem of reflection of a light pulse from a fully-reflective mirror and show that its reflection gives rise to an elastic wave with a measurable amplitude and a correct Doppler shift of the reflected pulse. We also analyze the problem of light pulse transmission into an anti-reflection coated, non-dispersive and lossless dielectric, where an elastic wave may as well be accompanied by a mechanical wave escorting the light pulse. We show that the Balazs rigid box thought experiment can be also realized in elastic dielectrics where some of the energy of the incident light is transferred to the wave motion. It follows from our approach that the electromagnetic momentum of the light pulse in the dielectric acquires Abraham's form only when a single type of the mechanical waves accompanies the interaction.
Hoshino, Shoma; Nakano, Yukio; Araki, Mitsunori; Ishiwata, Takashi; Tsukiyama, Koichi
2016-06-01
We report the first observation of collision induced state-to-state energy transfer from the 2u ((1)D2) (v2u = 3-7) ion-pair state of I2 using a perturbation facilitated optical-optical double resonance technique through the c (1)Πg∼ B (3)Π(0) hyperfine mixed double-faced valence state as the intermediate state. The excitation of the 2u ((1)D2) state yielded the weak UV fluorescence from the wide range of vibrational levels in the nearby 2g ((1)D2) state. The vibrational distribution in the 2g ((1)D2) state derived by the Franck-Condon simulation of the UV fluorescence showed that the population in the 2u ((1)D2) state transfers mostly to the 2g ((1)D2) vibronic levels which are located energetically above the laser-prepared level. The radiative lifetimes and the self-quenching rate constants were determined to be 21.3 ± 0.1 and 44.6 ± 0.8 ns, and (1.30 ± 0.01) × 10(-9) and (2.26 ± 0.17) × 10(-9) cm(3) molecule(-1) s(-1) for the 2u ((1)D2) (v2u = 3) and 2g ((1)D2) (v2g = 5) states, respectively. The rate constant for the 2u ((1)D2) - 2g ((1)D2) collision induced state-to-state energy transfer was also evaluated to be (1.89 ± 0.01), (3.07 ± 0.07), and (3.77 ± 0.05) × 10(-10) cm(3) molecule(-1) s(-1) for the v2u = 3, 5, and 7 levels, respectively. The very large self-quenching cross sections for the ion-pair states of I2 could be explained by the harpoon mechanism. PMID:27165483
Utrecht Radiative Transfer Courses
NASA Astrophysics Data System (ADS)
Rutten, R. J.
2003-01-01
The Utrecht course ``The Generation and Transport of Radiation'' teaches basic radiative transfer to second-year students. It is a much-expanded version of the first chapter of Rybicki & Lightman's ``Radiative Processes in Astrophysics''. After this course, students understand why intensity is measured per steradian, have an Eddington-Barbier feel for optically thick line formation, and know that scattering upsets LTE. The text is a computer-aided translation by Ruth Peterson of my 1992 Dutch-language course. My aim is to rewrite this course in non-computer English and make it web-available at some time. In the meantime, copies of the Peterson translation are made yearly at Uppsala -- ask them, not me. Eventually it should become a textbook. The Utrecht course ``Radiative Transfer in Stellar Atmospheres'' is a 30-hour course for third-year students. It treats NLTE line formation in plane-parallel stellar atmospheres at a level intermediate between the books by Novotny and Boehm-Vitense, and Mihalas' ``Stellar Atmospheres''. After this course, students appreciate that epsilon is small, that radiation can heat or cool, and that computers have changed the field. This course is web-available since 1995 and is regularly improved -- but remains incomplete. Eventually it should become a textbook. The three Utrecht exercise sets ``Stellar Spectra A: Basic Line Formation'', ``Stellar Spectra B: LTE Line Formation'', and ``Stellar Spectra C: NLTE Line Formation'' are IDL-based computer exercises for first-year, second-year, and third-year students, respectively. They treat spectral classification, Saha-Boltzmann population statistics, the curve of growth, the FAL-C solar atmosphere model, the role of H-minus in the solar continuum, LTE formation of Fraunhofer lines, inversion tactics, the Feautrier method, classical lambda iteration, and ALI computation. The first two sets are web-available since 1998; the third will follow. Acknowledgement. Both courses owe much to previous
Spectrally Invariant Approximation within Atmospheric Radiative Transfer
NASA Technical Reports Server (NTRS)
Marshak, A.; Knyazikhin, Y.; Chiu, J. C.; Wiscombe, W. J.
2011-01-01
Certain algebraic combinations of single scattering albedo and solar radiation reflected from, or transmitted through, vegetation canopies do not vary with wavelength. These spectrally invariant relationships are the consequence of wavelength independence of the extinction coefficient and scattering phase function in vegetation. In general, this wavelength independence does not hold in the atmosphere, but in cloud-dominated atmospheres the total extinction and total scattering phase function vary only weakly with wavelength. This paper identifies the atmospheric conditions under which the spectrally invariant approximation can accurately describe the extinction and scattering properties of cloudy atmospheres. The validity of the assumptions and the accuracy of the approximation are tested with 1D radiative transfer calculations using publicly available radiative transfer models: Discrete Ordinate Radiative Transfer (DISORT) and Santa Barbara DISORT Atmospheric Radiative Transfer (SBDART). It is shown for cloudy atmospheres with cloud optical depth above 3, and for spectral intervals that exclude strong water vapor absorption, that the spectrally invariant relationships found in vegetation canopy radiative transfer are valid to better than 5%. The physics behind this phenomenon, its mathematical basis, and possible applications to remote sensing and climate are discussed.
Benchmarks and models for 1-D radiation transport in stochastic participating media
Miller, D S
2000-08-21
Benchmark calculations for radiation transport coupled to a material temperature equation in a 1-D slab and 1-D spherical geometry binary random media are presented. The mixing statistics are taken to be homogeneous Markov statistics in the 1-D slab but only approximately Markov statistics in the 1-D sphere. The material chunk sizes are described by Poisson distribution functions. The material opacities are first taken to be constant and then allowed to vary as a strong function of material temperature. Benchmark values and variances for time evolution of the ensemble average of material temperature energy density and radiation transmission are computed via a Monte Carlo type method. These benchmarks are used as a basis for comparison with three other approximate methods of solution. One of these approximate methods is simple atomic mix. The second approximate model is an adaptation of what is commonly called the Levermore-Pomraning model and which is referred to here as the standard model. It is shown that recasting the temperature coupling as a type of effective scattering can be useful in formulating the third approximate model, an adaptation of a model due to Su and Pomraning which attempts to account for the effects of scattering in a stochastic context. This last adaptation shows consistent improvement over both the atomic mix and standard models when used in the 1-D slab geometry but shows limited improvement in the 1-D spherical geometry. Benchmark values are also computed for radiation transmission from the 1-D sphere without material heating present. This is to evaluate the performance of the standard model on this geometry--something which has never been done before. All of the various tests demonstrate the importance of stochastic structure on the solution. Also demonstrated are the range of usefulness and limitations of a simple atomic mix formulation.
LRAT: Lightning Radiative Transfer
NASA Technical Reports Server (NTRS)
Phanord, Dieudonne D.
1993-01-01
In this report, we extend to cloud physics the work done for single and multiple scattering of electromagnetic waves. We consider the scattering of light, visible or infrared, by a spherical cloud represented by a statistically homogeneous ensemble of configurations of N identical spherical water droplets whose centers are uniformly distributed in its volume V. The ensemble is specified by the average number rho of scatterers in unit volume and by rho f(R) with f(R) as the distribution function for separations R of pairs. The incident light, vector-phi(sub 0) a plane electromagnetic wave with harmonic time dependence, is from outside the cloud. The propagation parameter kappa(sub 0) and the index of refraction eta(sub 0) determine physically the medium outside the distribution of scatterers. We solve the interior problem separately to obtain the bulk parameters for the scatterer equivalent to the ensemble of spherical droplets. With the interior solution or the equivalent medium approach, the multiple scattering problem is reduced to that of an equivalent single scatterer excited from outside illumination. A dispersion relation which determines the bulk propagation parameter K and the bulk index of refraction eta of the cloud is given in terms of the vector equivalent scattering amplitude vector-G and the dyadic scattering amplitude tilde-g of the single object in isolation. Based on this transfer model we will have the ability to consider clouds composed of inhomogeneous distribution of water and/or ice particles and we will be able to take into account particle size distributions within the cloud. We will also be able to study the effects of cloud composition (i.e., particle shape, size, composition, orientation, location) on the polarization of the single or the multiple scattered waves. Finally, this study will provide a new starting point for studying the problem of lightning radiative transfer.
An operator perturbation method for polarized line transfer. I. Non-magnetic regime in 1D media.
NASA Astrophysics Data System (ADS)
Faurobert-Scholl, M.; Frisch, H.; Nagendra, K. N.
1997-06-01
In this paper we generalize an Approximate Lambda Iteration (ALI) technique developed for scalar transfer problems to a vectorial transfer problem for polarized radiation. Scalar ALI techniques are based on a suitable decomposition of the Lambda operator governing the integral form of the transfer equation. Lambda operators for scalar transfer equations are diagonally dominant, offering thus the possibility to use iterative methods of the Jacobi type where the iteration process is based on the diagonal of the Lambda operator (Olson et al., 1986, JQSRT 35, 431). Here we consider resonance polarization, created by the scattering of an anisotropic radiation field, for spectral lines formed with complete frequency redistribution in a 1D axisymmetric medium. The problem can be formulated as an integral equation for a 2-component vector (Rees, 1978PASJ...30..455R) or, as shown by Ivanov (1995A&A...303..621I), as an integral equation for a (2x2) matrix source function which involves the same generalized Lambda operator as the vector integral equation. We find that this equation holds also in the presence of a weak turbulent magnetic field. The generalized Lambda operator is a (2x2) matrix operator. The element {11} describes the propagation of the intensity and is identical to the Lambda operator of non-polarized problems. The element {22} describes the propagation of the polarization. The off-diagonal terms weakly couple the intensity and the polarization. We propose a block Jacobi iterative method and show that its convergence properties are controlled by the propagator for the intensity. We also show that convergence can be accelerated by an Ng acceleration method applied to each element of the source matrix. We extend to polarized transfer a convergence criterion introduced by Auer et al. (1994A&A...292..599A) based on the grid truncation error of the converged solution.
DOPEX-1D2C: A one-dimensional, two-constraint radiation shield optimization code
NASA Technical Reports Server (NTRS)
Lahti, G. P.
1973-01-01
A one-dimensional, two-constraint radiation sheild weight optimization procedure and a computer program, DOPEX-1D2C, is described. The DOPEX-1D2C uses the steepest descent method to alter a set of initial (input) thicknesses of a spherical shield configuration to achieve a minimum weight while simultaneously satisfying two dose-rate constraints. The code assumes an exponential dose-shield thickness relation with parameters specified by the user. Code input instruction, a FORTRAN-4 listing, and a sample problem are given. Typical computer time required to optimize a seven-layer shield is less than 1/2 minute on an IBM 7094.
Radiative Transfer: Methods and Applications
NASA Astrophysics Data System (ADS)
Mayer, Bernhard; Emde, Claudia; Buras, Robert; Kylling, Arve
Solar and terrestrial radiation is the driver of atmospheric dynamics and chemistry and can be exploited by remote sensing algorithms to determine atmospheric composition. For this purpose, accurate radiative transfer models are needed. Here, a modern radiative transfer tool developed over many years at the Institute of Atmospheric Physics is explained. As an application, the remote sensing of cloud microphysics using the angular distribution of reflected solar radiance in the rainbow and backscatter glory is shown, with special emphasis on the polarization of radiation.
A 1-D radiative conductive model to study the SOIR/VEx thermal profiles
NASA Astrophysics Data System (ADS)
Mahieux, Arnaud; Erwin, Justin T.; Chamberlain, Sarah; Robert, Séverine; Carine Vandaele, Ann; Wilquet, Valérie; Thomas, Ian; Yelle, Roger V.; Bertaux, Jean-Loup
2015-04-01
SOIR is an infrared spectrometer on board Venus Express that probes the Venus terminator region since 2006. The measurements are taken on the morning and evening sides of the terminator, covering all latitudes from the North Pole to the South Pole. Its wavelength range - 2.2 to 4.3 μm - allows a detailed chemical inventory of the Venus atmosphere [1-5], such as CO2, CO, H2O, HCl, HF, SO2 and aerosols. CO2 is detected from 70 km up to 165 km, CO from 70 km to 140 km, and the minor species typically below 110 km down to 70 km. Number density profiles of these species are computed from the measured spectra. Temperature profiles are obtained while computing the spectral inversion of the CO2 spectra combined with the hydrostatic law [6]. These temperature measurements show a striking permanent temperature minimum (at 125 km) and a weaker temperature maximum (over 100-115 km). The time variability of the CO2 density profiles spans over two orders of magnitude, and a clear trend is seen with latitude. The temperature variations are also important, of the order of 35 K for a given pressure level, but the latitude variation are small. Miss-RT, a 1D radiative transfer model has been developed to reproduce the SOIR terminator profiles, derived from the Mars thermosphere code presented in [7]. This model has been expanded to better account for the CO2, CO, and O non-LTE radiative heating and cooling processes which have to be considered in the dense atmosphere of Venus. Radiative cooling by minor species detected by SOIR (e.g. HCl, SO2, and H2O) are found to be small in comparison to the 15 μm CO2 cooling. Aerosol cooling in the 60-90km altitude range may be important to the thermal balance. There is a good agreement between the 1D model temperature profile and the mean SOIR temperature profile. Further we can suggest parameters that can be adjusted to improve the agreement between the model and measurements. The remaining differences can be attributed to the atmosphere
Quantum propagation and confinement in 1D systems using the transfer-matrix method
NASA Astrophysics Data System (ADS)
Pujol, Olivier; Carles, Robert; Pérez, José-Philippe
2014-05-01
The aim of this article is to provide some Matlab scripts to the teaching community in quantum physics. The scripts are based on the transfer-matrix formalism and offer a very efficient and versatile tool to solve problems of a physical object (electron, proton, neutron, etc) with one-dimensional (1D) stationary potential energy. Resonant tunnelling through a multiple-barrier or confinement in wells of various shapes is particularly analysed. The results are quantitatively discussed with semiconductor heterostructures, harmonic and anharmonic molecular vibrations, or neutrons in a gravity field. Scripts and other examples (hydrogen-like ions and transmission by a smooth variation of potential energy) are available freely at http://www-loa.univ-lille1.fr/˜pujol in three languages: English, French and Spanish.
NASA Technical Reports Server (NTRS)
Barker, Howard W.; Kato, Serji; Wehr, T.
2012-01-01
The main point of this study was to use realistic representations of cloudy atmospheres to assess errors in solar flux estimates associated with 1D radiative transfer models. A scene construction algorithm, developed for the EarthCARE satellite mission, was applied to CloudSat, CALIPSO, and MODIS satellite data thus producing 3D cloudy atmospheres measuring 60 km wide by 13,000 km long at 1 km grid-spacing. Broadband solar fluxes and radiances for each (1 km)2 column where then produced by a Monte Carlo photon transfer model run in both full 3D and independent column approximation mode (i.e., a 1D model).
Gen Purpose 1-D Finite Element Network Fluid Flow Heat Transfer System Simulator
Energy Science and Technology Software Center (ESTSC)
1993-08-02
SAFSIM (System Analysis Flow Simulator) is a FORTRAN computer program to simulate the integrated performance of systems involving fluid mechanics, heat transfer, and reactor dynamics. SAFSIM provides sufficient versatility to allow the engineering simulation of almost any system, from a backyard sprinkler system to a clustered nuclear reactor propulsion system. In addition to versatility, speed and robustness are primary SAFSIM development goals. SAFSIM contains three basic physics modules: (1) a one-dimensional finite element fluid mechanicsmore » module with multiple flow network capability; (2) a one-dimensional finite element structure heat transfer module with multiple convection and radiation exchange capability; and (3) a point reactor dynamics module with reactivity feedback and decay heat capability. SAFSIM can be used for compressible and incompressible, single-phase, multicomponent flow systems.« less
Auroral resonance line radiative transfer
Gladstone, G.R. )
1992-02-01
A model is developed for simulating the two-dimensional radiative transfer of resonance line emissions in auroras. The method of solution utilizes Fourier decomposition of the horizontal dependence in the intensity field so that the two-dimensional problem becomes a set of one-dimensional problems having different horizontal wavenumbers. The individual one-dimensional problems are solved for using a Feautrier-type solution of the differential-integral form of the radiative transfer equation. In the limit as the horizontal wavenumber becomes much larger than the local line-center extinction coefficient, the scattering integral becomes considerably simplified, and the final source function is evaluated in closed form. The two-dimensional aspects of the model are tested against results for nonresonance radiative transfer studies, and the resonance line part of the model is tested against results of existing plane-parallel resonance line radiative transfer codes. Finally, the model is used to simulate the intensity field of O{sub I} 1,304{angstrom} for hard and soft auroras of various Gaussian horizontal widths. The results demonstrate the importance of considering the effects of two-dimensional radiative transfer when analyzing auroral resonance line data.
HELIOS-CR A 1-D radiation-magnetohydrodynamics code with inline atomic kinetics modeling
NASA Astrophysics Data System (ADS)
Macfarlane, J. J.; Golovkin, I. E.; Woodruff, P. R.
2006-05-01
HELIOS-CR is a user-oriented 1D radiation-magnetohydrodynamics code to simulate the dynamic evolution of laser-produced plasmas and z-pinch plasmas. It includes an in-line collisional-radiative (CR) model for computing non-LTE atomic level populations at each time step of the hydrodynamics simulation. HELIOS-CR has been designed for ease of use, and is well-suited for experimentalists, as well as graduate and undergraduate student researchers. The energy equations employed include models for laser energy deposition, radiation from external sources, and high-current discharges. Radiative transport can be calculated using either a multi-frequency flux-limited diffusion model, or a multi-frequency, multi-angle short characteristics model. HELIOS-CR supports the use of SESAME equation of state (EOS) tables, PROPACEOS EOS/multi-group opacity data tables, and non-LTE plasma properties computed using the inline CR modeling. Time-, space-, and frequency-dependent results from HELIOS-CR calculations are readily displayed with the HydroPLOT graphics tool. In addition, the results of HELIOS simulations can be post-processed using the SPECT3D Imaging and Spectral Analysis Suite to generate images and spectra that can be directly compared with experimental measurements. The HELIOS-CR package runs on Windows, Linux, and Mac OSX platforms, and includes online documentation. We will discuss the major features of HELIOS-CR, and present example results from simulations.
Radiative transfer in spherical atmospheres
NASA Technical Reports Server (NTRS)
Kalkofen, W.; Wehrse, R.
1984-01-01
A method for defining spherical model atmospheres in radiative/convective and hydrostatic equilibrium is presented. A finite difference form is found for the transfer equation and a matrix operator is developed as the discrete space analog (in curvilinear coordinates) of a formal integral in plane geometry. Pressure is treated as a function of temperature. Flux conservation is maintained within the energy equation, although the correct luminosity transport must be assigned for any given level of the atmosphere. A perturbed integral operator is used in a complete linearization of the transfer and constraint equations. Finally, techniques for generating stable solutions in economical computer time are discussed.
Transfer Matrix Approach to 1d Random Band Matrices: Density of States
NASA Astrophysics Data System (ADS)
Shcherbina, Mariya; Shcherbina, Tatyana
2016-08-01
We study the special case of n× n 1D Gaussian Hermitian random band matrices, when the covariance of the elements is determined by the matrix J=(-W^2triangle +1)^{-1} . Assuming that n≥ CW log W≫ 1 , we prove that the averaged density of states coincides with the Wigner semicircle law up to the correction of order W^{-1}.
Radiative Transfer Under Inhomogeneous Configurations
NASA Astrophysics Data System (ADS)
Bendicho, P. Fabiani
1998-06-01
We present, for the first time, three dimensional (3D) radiative transfer (RT) results with realistic atomic models (multilevel) and without using the local thermodinamical equilibrium approximation (non-LTE). We have developed a new code based on efficient iterative methods (Trujillo Bueno, and Fabiani Bendicho 1995; Fabiani Bendicho, Trujillo Bueno and Auer 1997) characterized by a very high convergence rate. With this 3D multilevel code and using a schematic atmospheric model we are able to demonstrate that one may need self-consistent multidimensional RT calculations in order to interpret high spatial resolution solar spectroscopic observations.
Radiative transfer in dusty nebulae
NASA Technical Reports Server (NTRS)
Dana, R. A.
1977-01-01
The effects of dust scattering on observable optical and infrared parameters, and the accuracy of approximate solutions were examined. The equation of radiative transfer in a static and homogeneous, but not necessarily uniform, distribution gas and dust around a central empty core with a point source of energy at its center was solved. The dust properties were characterized by a phenomenological extinction cross section, albedo and parameters describing the anisotropy of dust scattering. For ultraviolet photons, ionization equilibrium equations for the gas were solved, and for infrared photons a self-consistent dust temperature was calculated. Ray tracing was used to solve for the angular dependence of the intensity.
Radiative Transfer in Protoplanetary Disks
NASA Astrophysics Data System (ADS)
Graziani, L.; Aiello, S.; Belleni-Morante, A.; Cecchi-Pestellini, C.
2008-09-01
Abstract Protoplanetary disks are the precursors of planetary systems. All building materials needed to assembly the planetary systems are supplied by these reservoirs, including many organic molecules [1,2]. Thus, the physical and chemical properties in Protoplanetary disks set the boundary conditions for the formation and evolution of planets and other solar system bodies. In standard radiative scenario structure and chemistry of protoplanetary disks depend strongly on the nature of central star around which they formed. The dust temperature is manly set by the stellar luminosity, while the chemistry of the whole disk depends on the UV and X ray fluxes [3,4,6,8]. Therefore, a knowledge as accurate as possible of the radiative transfer (RT) inside disks is a prerequisite for their modelling. Actually, real disks are complex, stratified and inhomogeneous environments requiring a detailed dust mixture modelling and the ability to follow the radiation transfer across radial and vertical gradients. Different energetic processes as the mass accretion processes onto the star surface, the viscous dissipative heating dominating the midplane region, and the flared atmospheres radiation reprocessing, have a significant role in the disk structuring [4,5,8]. During the last 10 years many authors suggested various numerical and analytical techniques to resolve the disk temperature structure providing vertical temperature profiles and disk SED databases [4,6]. In this work we present the results of our semi analytical and numerical model solving the radiative transfer problem in two separate interesting disk regions: 1) Disk atmospheres at large radius, r > 10 AU. 2) Vertical disk structure over 1 < r < 10 AU and 10 < r < 100 AU. A simplified analytical approach based on P-N approximation [7] for a rectified disk surface (suitable for limited range of r) is compared and contrasted with a more accurate Monte Carlo integration [5]. Our code can handle arbitrary dust
Stochastic Radiative transfer and real cloudiness
Evans, F.
1995-09-01
Plane-parallel radiative transfer modeling of clouds in GCMs is thought to be an inadequate representation of the effects of real cloudiness. A promising new approach for studying the effects of cloud horizontal inhomogeneity is stochastic radiative transfer, which computes the radiative effects of ensembles of cloud structures described by probability distributions. This approach is appropriate because cloud information is inherently statistical, and it is the mean radiative effect of complex 3D cloud structure that is desired. 2 refs., 1 fig.
The role of solitons in charge and energy transfer in 1D molecular chains
NASA Astrophysics Data System (ADS)
Ivić , Zoran
1998-03-01
The idea that polarons and solitons could play the crucial role in the transport processes in biological structures, has been critically reexamined on the basis of the general theory of self-trapping phenomena. The criteria which enable one to determine conditions for the existence and stability of polarons and solitons and to determine their character, in dependence of the values of the basic physical parameters of the system, were formulated. Validity of the so-called Davydov's soliton model was discussed on the basis of these criteria. It was found that the original Davydov's proposal, based upon the idea of the soliton creation due to the single excitation (particle, vibron, etc.) self-trapping, cannot explain the intramolecular energy transfer in α-helix and acetanilide. However, Davydov theory is flexible enough to describe the single electron transfer in some systems (α-helix and acetanilide for example). In the many-particle systems, dressing effect, due to the quantum nature of phonons, may cause the creation of the bound states of the several excitons in the molecular chain. The possibility of creation of the soliton states of this type is discussed for the simple Fröhlich's one-dimensional model. The regions of the system parameter space where different mechanisms dominate the behaviour of such entities are characterized.
XCHEM-1D: A Heat Transfer/Chemical Kinetics Computer Program for multilayered reactive materials
Gross, R.J.; Baer, M.R.; Hobbs, M.L.
1993-10-01
An eXplosive CHEMical kinetics code, XCHEM, has been developed to solve the reactive diffusion equations associated with thermal ignition of energetic materials. This method-of-lines code uses stiff numerical methods and adaptive meshing to resolve relevant combustion physics. Solution accuracy is maintained between multilayered materials consisting of blends of reactive components and/or inert materials. Phase change and variable properties are included in one-dimensional slab, cylindrical and spherical geometries. Temperature-dependent thermal properties have been incorporated and the modification of thermal conductivities to include decomposition effects are estimated using solid/gas volume fractions determined by species fractions. Gas transport properties, including high pressure corrections, have also been included. Time varying temperature, heat flux, convective and thermal radiation boundary conditions, and layer to layer contact resistances have also been implemented.
BART: Bayesian Atmospheric Radiative Transfer fitting code
NASA Astrophysics Data System (ADS)
Cubillos, Patricio; Blecic, Jasmina; Harrington, Joseph; Rojo, Patricio; Lust, Nate; Bowman, Oliver; Stemm, Madison; Foster, Andrew; Loredo, Thomas J.; Fortney, Jonathan; Madhusudhan, Nikku
2016-08-01
BART implements a Bayesian, Monte Carlo-driven, radiative-transfer scheme for extracting parameters from spectra of planetary atmospheres. BART combines a thermochemical-equilibrium code, a one-dimensional line-by-line radiative-transfer code, and the Multi-core Markov-chain Monte Carlo statistical module to constrain the atmospheric temperature and chemical-abundance profiles of exoplanets.
NASA Astrophysics Data System (ADS)
González-Rodríguez, Pedro; Ilan, Boaz; Kim, Arnold D.
2016-06-01
We introduce the one-way radiative transfer equation (RTE) for modeling the transmission of a light beam incident normally on a slab composed of a uniform forward-peaked scattering medium. Unlike the RTE, which is formulated as a boundary value problem, the one-way RTE is formulated as an initial value problem. Consequently, the one-way RTE is much easier to solve. We discuss the relation of the one-way RTE to the Fokker-Planck, small-angle, and Fermi pencil beam approximations. Then, we validate the one-way RTE through systematic comparisons with RTE simulations for both the Henyey-Greenstein and screened Rutherford scattering phase functions over a broad range of albedo, anisotropy factor, optical thickness, and refractive index values. We find that the one-way RTE gives very good approximations for a broad range of optical property values for thin to moderately thick media that have moderately to sharply forward-peaked scattering. Specifically, we show that the error made by the one-way RTE decreases monotonically as the anisotropic factor increases and as the albedo increases. On the other hand, the error increases monotonically as the optical thickness increases and the refractive index mismatch at the boundary increases.
NASA Astrophysics Data System (ADS)
Gangopadhyay, Sunandan; Kulkarni, Shailesh
2008-01-01
We apply the method of Banerjee and Kulkarni [R. Banerjee and S. Kulkarni, Phys. Rev. D 77, 024018 (2008).PRVDAQ0556-282110.1103/PhysRevD.77.024018] to provide a derivation of Hawking radiation from the Garfinkle-Horowitz-Strominger (stringy) black hole which falls in the class of the most general spherically symmetric black holes (-g≠1) and also the nonextremal D1-D5 black hole using only covariant gravitational anomalies.
Gangopadhyay, Sunandan; Kulkarni, Shailesh
2008-01-15
We apply the method of Banerjee and Kulkarni [R. Banerjee and S. Kulkarni, Phys. Rev. D 77, 024018 (2008).] to provide a derivation of Hawking radiation from the Garfinkle-Horowitz-Strominger (stringy) black hole which falls in the class of the most general spherically symmetric black holes ({radical}(-g){ne}1) and also the nonextremal D1-D5 black hole using only covariant gravitational anomalies.
Radiative heat transfer in porous uranium dioxide
Hayes, S.L.
1992-12-01
Due to low thermal conductivity and high emissivity of UO{sub 2}, it has been suggested that radiative heat transfer may play a significant role in heat transfer through pores of UO{sub 2} fuel. This possibility was computationally investigated and contribution of radiative heat transfer within pores to overall heat transport in porous UO{sub 2} quantified. A repeating unit cell was developed to model approximately a porous UO{sub 2} fuel system, and the heat transfer through unit cells representing a wide variety of fuel conditions was calculated using a finite element computer program. Conduction through solid fuel matrix as wekk as pore gas, and radiative exchange at pore surface was incorporated. A variety of pore compositions were investigated: porosity, pore size, shape and orientation, temperature, and temperature gradient. Calculations were made in which pore surface radiation was both modeled and neglected. The difference between yielding the integral contribution of radiative heat transfer mechanism to overall heat transport. Results indicate that radiative component of heat transfer within pores is small for conditions representative of light water reactor fuel, typically less than 1% of total heat transport. It is much larger, however, for conditions present in liquid metal fast breeder reactor fuel; during restructuring of this fuel type early in life, the radiative heat transfer mode was shown to contribute as much as 10-20% of total heat transport in hottest regions of fuel.
1D Scaling with Ablation for K-Shell Radiation from Stainless Steel Wire Arrays
Giuliani, J. L.; Thornhill, J. W.; Dasgupta, A.; Davis, J.; Clark, R. W.; Jones, B.; Cuneo, M.; Coverdale, C. A.; Deeney, C.
2009-01-21
A 1D Lagrangian magnetohydrodynamic z-pinch simulation code is extended to include wire ablation. The plasma transport coefficients are calibrated to reproduce the K-shell yields measured on the Z generator for three stainless steel arrays of diameter 55 mm and masses ranging from 1.8 to 2.7 mg. The resulting 1D scaling model is applied to a larger SS array (65 mm and 2.5 mg) on the refurbished Z machine. Simulation results predict a maximum K-shell yield of 77 kJ for an 82 kV charging voltage. This maximum drops to 42 kJ at 75 kV charging. Neglecting the ablation precursor leads to a {approx}10% change in the calculated yield.
Critical ingredients of Type Ia supernova radiative-transfer modelling
NASA Astrophysics Data System (ADS)
Dessart, Luc; Hillier, D. John; Blondin, Stéphane; Khokhlov, Alexei
2014-07-01
We explore the physics of Type Ia supernova (SN Ia) light curves and spectra using the 1D non-local thermodynamic equilibrium (non-LTE) time-dependent radiative-transfer code CMFGEN. Rather than adjusting ejecta properties to match observations, we select as input one `standard' 1D Chandrasekhar-mass delayed-detonation hydrodynamical model, and then explore the sensitivity of radiation and gas properties of the ejecta on radiative-transfer modelling assumptions. The correct computation of SN Ia radiation is not exclusively a solution to an `opacity problem', characterized by the treatment of a large number of lines. We demonstrate that the key is to identify and treat important atomic processes consistently. This is not limited to treating line blanketing in non-LTE. We show that including forbidden-line transitions of metals, and in particular Co, is increasingly important for the temperature and ionization of the gas beyond maximum light. Non-thermal ionization and excitation are also critical since they affect the colour evolution and the ΔM15 decline rate of our model. While impacting little the bolometric luminosity, a more complete treatment of decay routes leads to enhanced line blanketing, e.g. associated with 48Ti in the U and B bands. Overall, we find that SN Ia radiation properties are influenced in a complicated way by the atomic data we employ, so that obtaining converged results is a real challenge. Nonetheless, with our fully fledged CMFGEN model, we obtain good agreement with the golden standard Type Ia SN 2005cf in the optical and near-IR, from 5 to 60 d after explosion, suggesting that assuming spherical symmetry is not detrimental to SN Ia radiative-transfer modelling at these times. Multi-D effects no doubt matter, but they are perhaps less important than accurately treating the non-LTE processes that are crucial to obtain reliable temperature and ionization structures.
Doped SiO2 telecommunication fibre as a 1-D detector for radiation therapy dosimetry
NASA Astrophysics Data System (ADS)
Abdul Rahman, A. T.; Abdul Sani, Siti Fairus; Bradley, D. A.
2012-02-01
Present studies concern Ge-doped SiO2 telecommunication fibre as a high spatial resolution 1-D thermoluminescence (TL) system for radiotherapeutic dosimetry. Using tube xray bremsstrahlung sources operating at kilovoltage energies, these fibres have been shown to offer linear response, from < 1Gy up to in excess of 30 Gy. Measurement of the photoelectron dose enhancement resulting from use of a moderately high atomic number medium (iodinated contrast media) demonstrates the fibres to have the local dose sensitivity required of interface dosimetry. In PMMA, the TL yield is ~60% greater in the presence of iodine than in its absence.
Radiative transfer of visible radiation in turbid atmosphere
NASA Technical Reports Server (NTRS)
Yamamoto, G.; Tanaka, M.
1974-01-01
Methods are presented for solving radiative transfer problems; they include the doubling method and the closely related matrix method, iterative method, Chandrasekhar's method of discrete ordinates, and Monte Carlo method. To consider radiation transport through turbid atmosphere, an atmospheric model was developed characterizing aerosols by parameters. Intensity and polarization of radiation in turbid atmospheres is discussed, as well as lower atmospheric heating due to solar radiation absorption by aerosols.
Anscher, Mitchell S. . E-mail: ansch001@notes.duke.edu; Thrasher, Bradley; Rabbani, Zahid; Teicher, Beverly; Vujaskovic, Zeljko
2006-07-01
Purpose The aim of this study was to determine whether a neutralizing transforming growth factor-{beta} (TGF{beta}) antibody can prevent radiation (RT) induced lung injury. Methods and Materials Fractionated and sham right lung irradiation in Fischer 344 rats was delivered to assess the radioprotective effect of the antibodies. Animals were divided into the following groups: (1) control (sham RT, control antibody 13C4); (2) RT (800cGy x 5)+13C4); (3) RT + 0.1 mg/kg 1D11 anti-TGF{beta} antibody; and (4) RT + 1 mg/kg 1D11 antibody. Antibodies were intraperitoneally administered immediately after the last fraction of RT. Animals were sacrificed at 6 and 26 weeks after irradiation. Lungs were assessed for histologic changes, activation of macrophages, expression/activation of TGF{beta} and its signal transduction pathway. Results At 6 weeks post-RT, there was a significant reduction in macrophage accumulation (p = 0.041), alveolar wall thickness (p = 0.0003), and TGF-{beta} activation (p = 0.032) in animals receiving 1.0 mg/kg 1D11 vs. in the control group. However, at 6 weeks, the low dose of 1D11 antibody (0.1 mg/kg) failed to produce any significant changes. At 6 months post-RT, radioprotection is apparent for the group receiving 1.0 mg/kg 1D11, with activated macrophages (p = 0.037), alveolar wall thickness (p = 0.0002), TGF{beta} activation (p = 0.002) and its signal transduction proteins (p < 0.05) compared with the control group. Conclusions Administration of a single dose of 1.0 mg/kg of the anti-TGF{beta} antibody 1D11 resulted in decreased morphologic changes, inflammatory response, and reduced expression and activation of TGF{beta} 6 weeks and 6 months after 40 Gy to the right hemithorax. Targeting the TGF{beta} pathway may be a useful strategy to prevent radiation-induced lung injury.
A fast method for computing 1-D wakefields due to coherent synchrotron radiation
NASA Astrophysics Data System (ADS)
Mitchell, Chad E.; Qiang, Ji; Ryne, Robert D.
2013-07-01
A method for computing the free-space longitudinal wakefield due to coherent synchrotron radiation (CSR) in a one-dimensional model is developed using a fast integrated Green function approach. This approach accurately captures the short-range behavior of the CSR interaction and does not require the numerical differentiation of a noisy longitudinal charge density. The transient wakefields that occur near bend entry and exit are included. This method can also be generalized to include the effect of upstream radiation that propagates through multiple lattice elements before interacting with the bunch.
Radiation heat transfer shapefactors for combustion systems
NASA Technical Reports Server (NTRS)
Emery, A. F.; Johansson, O.; Abrous, A.
1987-01-01
The computation of radiation heat transfer through absorbing media is commonly done through the zoning method which relies upon values of the geometric mean transmittance and absorptance. The computation of these values is difficult and expensive, particularly if many spectral bands are used. This paper describes the extension of a scan line algorithm, based upon surface-surface radiation, to the computation of surface-gas and gas-gas radiation transmittances.
Radiative Transfer In Gamma-Ray Bursts
NASA Astrophysics Data System (ADS)
Beloborodov, Andrei
We propose to develop state-of-the-art numerical tools for radiative transfer calculations in gamma-ray bursts (GRBs). We will investigate two problems: (1) Production and heating of photons at the early (opaque) stage of the explosion, which controls the brightness and spectral shape of the jet photospheric emission. (2) Transfer of GRB radiation through the external blast wave. Our recent results suggest that this transfer generates the GeV flash observed in GRBs, providing key information on the explosion and its progenitor. We will test our models against observations.
Discrete Space Theory of Radiative Transfer: Application
NASA Astrophysics Data System (ADS)
Rao, M. Srinivasa
2010-06-01
The method of obtaining the solution of radiative transfer equation using discrete space theory (DST) is described with (1) interaction principle for different geometries (2) star product (3) calculation of radiation field at internal points. Some of the important steps to obtain the solution of radiative transfer equation in spherical symmetry are also mentioned. Applications of DST are discussed with their results in two cases (a) study of reflection effect in close binary systems and (b) to compute KI 769.9 nm emission line profiles from N-type stars.
Discrete Space Theory of Radiative Transfer: Application
NASA Astrophysics Data System (ADS)
Rao, M. Srinivasa
The method of obtaining the solution of radiative transfer equation using discrete space theory (DST) is described with (1) interaction principle for different geometries (2) star product (3) calculation of radiation field at internal points. Some of the important steps to obtain the solution of radiative transfer equation in spherical symmetry are also mentioned. Applications of DST are discussed with their results in two cases (a) study of reflection effect in close binary systems and (b) to compute KI 769.9 nm emission line profiles from N-type stars.
Nonaxisymmetric radiative transfer in cylindrical enclosures
Moder, J.P.; Lee, H.S.; Chai, J.C.; Parthasarathy, G.; Patankar, S.V.
1996-12-31
A finite-volume method for radiative transfer in cylindrical enclosures is presented. Angular redistribution terms in the equation of transfer are avoided by defining radiation directions in terms of angular coordinates measured with respect to Cartesian base vectors; this definition of radiation directions can result in control angles which overlap control-volume faces, depending on the type of spatial and angular grids used in the azimuthal direction. A simple treatment for such control-angle overlaps is presented which is also applicable to nonorthogonal curvilinear spatial-coordinates. A comparison of the present procedure with other similar methods is given. Solutions are presented for axisymmetric transfer through a cylinder and nonaxisymmetric transfer through two- and three-dimensional annular sectors. Results show that the procedure produces reasonable solutions for transparent and participating media in axisymmetric and nonaxisymmetric cylindrical enclosures.
Transient radiative transfer through scattering absorbing media
Mitra, K.; Kumar, S.
1996-12-31
This paper outlines the formulation of the different methods for determining transient radiative transfer through scattering absorbing media. A boundary driven radiative problem is considered in a one-dimensional plane-parallel slab. The different methods of solving the transient radiative transfer equation include the P{sub 1}, P{sub 3}, and P{sub 5} approximations, two-flux method, and eight, twelve and sixteen discrete ordinates methods. In addition, the general transient radiative transfer equation is also solved by direct numerical integration without any simplifying assumptions. Different orders of approximation for the phase function are considered as is a parametric analysis of the different parameters such as the scattering albedo and optical depth is performed. The propagation speed obtained and the magnitude of the transmitted and back-scattered fluxes for different models obtained are a function of the approximation used to represent the intensity distribution.
Radiative transfer model: matrix operator method.
Liu, Q; Ruprecht, E
1996-07-20
A radiative transfer model, the matrix operator method, is discussed here. The matrix operator method is applied to a plane-parallel atmosphere within three spectral ranges: the visible, the infrared, and the microwave. For a homogeneous layer with spherical scattering, the radiative transfer equation can be solved analytically. The vertically inhomogeneous atmosphere can be subdivided into a set of homogeneous layers. The solution of the radiative transfer equation for the vertically inhomogeneous atmosphere is obtained recurrently from the analytical solutions for the subdivided layers. As an example for the application of the matrix operator method, the effects of the cirrus and the stratocumulus clouds on the net radiation at the surface and at the top of the atmosphere are investigated. The relationship between the polarization in the microwave range and the rain rates is also studied. Copies of the FORTRAN program and the documentation of the FORTRAN program on a diskette are available. PMID:21102832
NASA Technical Reports Server (NTRS)
Minow, Joseph I.; Coffey, Victoria N.; Parker, Linda N.; Blackwell, William C., Jr.; Jun, Insoo; Garrett, Henry B.
2007-01-01
The NUMIT 1-dimensional bulk charging model is used as a screening to ol for evaluating time-dependent bulk internal or deep dielectric) ch arging of dielectrics exposed to penetrating electron environments. T he code is modified to accept time dependent electron flux time serie s along satellite orbits for the electron environment inputs instead of using the static electron flux environment input originally used b y the code and widely adopted in bulk charging models. Application of the screening technique ts demonstrated for three cases of spacecraf t exposure within the Earth's radiation belts including a geostationa ry transfer orbit and an Earth-Moon transit trajectory for a range of orbit inclinations. Electric fields and charge densities are compute d for dielectric materials with varying electrical properties exposed to relativistic electron environments along the orbits. Our objectiv e is to demonstrate a preliminary application of the time-dependent e nvironments input to the NUMIT code for evaluating charging risks to exposed dielectrics used on spacecraft when exposed to the Earth's ra diation belts. The results demonstrate that the NUMIT electric field values in GTO orbits with multiple encounters with the Earth's radiat ion belts are consistent with previous studies of charging in GTO orb its and that potential threat conditions for electrostatic discharge exist on lunar transit trajectories depending on the electrical proper ties of the materials exposed to the radiation environment.
Session on modeling of radiative transfer processes
NASA Technical Reports Server (NTRS)
Flatau, Piotr
1993-01-01
The session on modeling of radiative transfer processes is reviewed. Six critical issues surfaced in the discussion concerning scale-interactive radiative processes relevent to the mesoscale convective systems (MCS's). These issues are the need to expand basic knowledge of how MCS's influence climate through extensive cloud shields and increased humidity in the upper troposphere; to improve radiation parameterizations used in mesoscale and General Circulation Model (GCM) models; to improve our basic understanding of the influence of radiation on MCS dynamics due to diabatic heating, production of condensate, and vertical and horizontal heat fluxes; to quantify our understanding of radiative impacts of MCS's on the surface and free atmosphere energy budgets; to quantify and identify radiative and microphysical processes important in the evolution of MCS's; and to improve the capability to remotely sense MCS radiative properties from space and ground-based systems.
NASA Astrophysics Data System (ADS)
Sen, K. K., Wilson, S. J.
The advancement of observational techniques over the years has led to the discovery of a large number of stars exhibiting complex spectral structures, thus necessitating the search for new techniques and methods to study radiative transfer in such stars with moving envelopes. This led to the introduction of the concept of "photon escape probability" and the wisdom of expressing the transfer equations in "comoving frames" (CMF). Radiative transfer problems in spherically moving media form a branch of mathematical physics which uses mathematics of a very distinctive kind. Radiative Transfer in Moving Media records the basic mathematical methodologies, both analytical and numerical, developed for solving radiation transfer problems in spherically symmetric moving media, in the consideration of macroscopic velocity fields only. Part I contains the basic notions of radiation-matter interaction in participating media and constructs the relevant transfer equations to be solved in the subsequent chapters. Part II considers the basic mathematical methods for solving the transfer problems in extensive moving atmospheres when it is observed in the lab frame. Part III introduces the analytical and numerical methods for solving radiative transfer problems in spherically symmetric moving atmospheres when expressed in the comoving frame. This book is addressed to graduate students and researchers in Astrophysics, in particular to those studying radiative transfer in stellar atmospheres.
A 1D radiative-convective model of H2O-CO2 atmospheres around young telluric planets: an update
NASA Astrophysics Data System (ADS)
Marcq, Emmanuel; Salvador, Arnaud; Massol, Hélène; Chassefière, Éric
2016-04-01
The study of the early phases of the evolution of terrestrial planets has recently known significant progress [1,2]. It appears that their cooling phase during the magma ocean stage is first dominated by a radiative cooling stage through its atmosphere. If the planet is able to reach radiative balance during this stage, then its further evolution is dominated by the escape flux, and no large scale condensation of water occurs (Hamano-type II planets). On the other hand, if the planet is far enough from the sun, then radiative equilibrium cannot be reached until the outgoing flux has fallen below the runaway greenhouse limit, implying the condensation of most atmospheric water vapor into a global water ocean, thus sheltering most water from atmospheric escape (Hamano-type I planet). In the solar system, Earth is clearly a type-I planet, whereas Venus was most likely a type-II planet from quite early on in its history [1,2]. In this presentation, we will deal with the atmospheric radiative model used by [2] and first described in [3]. After describing its recent improvements since [3] (pressure grid enabling an arbitrary total volatile amount, correction of the k-correlated radiative transfer in the thermal radiation, improvement of the numerical stability and integration scheme) and their consequences on the detectability of extrasolar type-I or type-II planets, we will deal with the possible improvements and extensions to such models, such as but not limited to: (1) adopting a 1D-spherical geometry suited for larger atmospheres around smaller planets, (2) improvement of the visible albedo parameterization based on recent 3D-modelling GCM [4]. [1] : K. Hamano et al., Nature (2013) [2] : T. Lebrun et al. JGR (2013) [3] : E. Marcq, JGR (2012) [4] : J. Leconte et al. (2015)
3D Thermal Infrared Radiative Transfer in Mountains
NASA Astrophysics Data System (ADS)
Lee, W.; Liou, K.; Hall, A.
2007-12-01
We developed a 3D Monte Carlo photon tracing program for radiative transfer in inhomogeneous and irregular terrain coupled with the correlated k-distribution method for gaseous absorption in the atmosphere for the calculation of broadband thermal infrared (IR) fluxes at mountain surfaces. The thermal IR radiative transfer program includes emission from the atmosphere to the surface and vice versa as well as emissions between mountain surfaces. Both the atmosphere and the land surface are discretized by using finite cubic cells characterized by the spectral optical properties of molecules and background aerosols (absorption coefficient, single-scattering albedo, and scattering phase function) and terrain configuration (albedo, elevation, slope, and orientation). The emissivity of gases is parameterized in terms of the vertical optical depth of cubic cell. We selected an area of 100×100 km2 in the Tibetan Plateau near Lhasa city with a horizontal resolution of 1 km2 and used the surface temperature and albedo available from MODIS/Terra dataset for this study. We show that surface temperature is the dominating factor in radiative transfer calculations and that subgrid variability of the net surface IR flux distribution relative to a flat surface (1D) with average elevation and temperature can be as large as 50 W/m2 at cold mountain surfaces.
NASA Astrophysics Data System (ADS)
Ashworth, K.; Chung, S. H.; Griffin, R. J.; Chen, J.; Forkel, R.; Bryan, A. M.; Steiner, A. L.
2015-11-01
Biosphere-atmosphere interactions play a critical role in governing atmospheric composition, mediating the concentrations of key species such as ozone and aerosol, thereby influencing air quality and climate. The exchange of reactive trace gases and their oxidation products (both gas and particle phase) is of particular importance in this process. The FORCAsT (FORest Canopy Atmosphere Transfer) 1-D model is developed to study the emission, deposition, chemistry and transport of volatile organic compounds (VOCs) and their oxidation products in the atmosphere within and above the forest canopy. We include an equilibrium partitioning scheme, making FORCAsT one of the few canopy models currently capable of simulating the formation of secondary organic aerosols (SOAs) from VOC oxidation in a forest environment. We evaluate the capability of FORCAsT to reproduce observed concentrations of key gas-phase species and report modeled SOA concentrations within and above a mixed forest at the University of Michigan Biological Station (UMBS) during the Community Atmosphere-Biosphere Interactions Experiment (CABINEX) field campaign in the summer of 2009. We examine the impact of two different gas-phase chemical mechanisms on modelled concentrations of short-lived primary emissions, such as isoprene and monoterpenes, and their oxidation products. While the two chemistry schemes perform similarly under high-NOx conditions, they diverge at the low levels of NOx at UMBS. We identify peroxy radical and alkyl nitrate chemistry as the key causes of the differences, highlighting the importance of this chemistry in understanding the fate of biogenic VOCs (bVOCs) for both the modelling and measurement communities.
Data Assimilation Using a Variational Method for a 1D Radiation Belt Diffusion Model
NASA Astrophysics Data System (ADS)
Marchand, R.; Degeling, A. W.; O'Donnell, S.; Rankin, R.; Kabin, K.
2009-12-01
A variational data assimilation algorithm has been developed to incorporate electron flux time-series data from satellites into a simple one dimensional diffusion model for the radial transport of radiation belt electrons. The model developed assumes a power law scaling for the radial diffusion coefficient with L shell. The effectiveness of this method is investigated by means of a series of identical twin numerical experiments. This involves using the diffusion model to produce synthetic observations along various satellite trajectories. These observations are in turn used to estimate time-dependent parameters input to the diffusion model, which are compared against the values initially used. The data assimilation algorithm considers the time dependent source located at the outer boundary as a function to be determined. Using synthetic satellite electron flux observations, the algorithm computes a source function that, when used as an input to the diffusion model, most closely reproduces the synthetic observations in a least-squares sense. Observational errors are taken into account, and an estimate of the uncertainty in the output source function is also produced. This uncertainty is found to consistently reflect the quality of the source function estimation during identical twin numerical experiments. Initial tests indicate that the quality of the outer boundary source estimation is strongly dependent on the satellite location, indicating that the outer boundary source estimation becomes poor as information relating to the outer boundary contained in the observations is reduced. The potential of using this data assimilation method to estimate one or more parameters that determine the radial diffusion coefficient, and the possibility of determining whether physical processes affecting the observations are missing in the dynamical model will be discussed.
Radiation-induced hydrogen transfer in metals
NASA Astrophysics Data System (ADS)
Tyurin, Yu I.; Vlasov, V. A.; Dolgov, A. S.
2015-11-01
The paper presents processes of hydrogen (deuterium) diffusion and release from hydrogen-saturated condensed matters in atomic, molecular and ionized states under the influence of the electron beam and X-ray radiation in the pre-threshold region. The dependence is described between the hydrogen isotope release intensity and the current density and the electron beam energy affecting sample, hydrogen concentration in the material volume and time of radiation exposure to the sample. The energy distribution of the emitted positive ions of hydrogen isotopes is investigated herein. Mechanisms of radiation-induced hydrogen transfer in condensed matters are suggested.
Fluss, M.J.; Berko, S.; Chakraborty, B.; Hoffmann, K.R.; Lippel, P.; Siegel, R.W.
1985-03-12
One- and two-dimensional angular correlation of positron-electron annihilation radiation (1D and 2D-ACAR) data have been obtained between 293 and 903 K for single crystals of aluminum. The peak counting rates vs temperature, which were measured using the 1D-ACAR technique, provide a model independent value for the temperature dependence of the positron trapping probability. Using these results it is possible to strip out the Bloch state contribution from the observed 2D-ACAR surfaces and then compare the resulting defect ACAR surfaces to calculated 2D-ACAR surfaces for positrons annihilating from the Bloch, monovacancy, and divacancy-trapped states. The result of this comparison is that the presence of an increasing equilibrium divacancy population is consistent with the observed temperature dependence of ACAR data at high temperature in Al and that the present results when compared to earlier studies on Al indicate that the ratio of the trapping rates at divacancies and monovacancies is of order two.
2-DUST: Dust radiative transfer code
NASA Astrophysics Data System (ADS)
Ueta, Toshiya; Meixner, Margaret
2016-04-01
2-DUST is a general-purpose dust radiative transfer code for an axisymmetric system that reveals the global energetics of dust grains in the shell and the 2-D projected morphologies of the shell that are strongly dependent on the mixed effects of the axisymmetric dust distribution and inclination angle. It can be used to model a variety of axisymmetric astronomical dust systems.
Three Dimensional Radiative Transfer In Tropical Deep Convective Clouds.
NASA Astrophysics Data System (ADS)
di Giuseppe, F.
In this study the focus is on the interaction between short-wave radiation with a field of tropical deep convective events generated using a 3D cloud resolving model (CRM) to assess the significance of 3D radiative transport (3DRT). It is not currently un- derstood what magnitude of error is involved when a two stream approximation is used to describe the radiative transfer through such a cloud field. It seems likely that deep convective clouds could be the most complex to represent, and that the error in neglecting horizontal transport could be relevant in these cases. The field here con- sidered has an extention of roughly 90x90 km, approximately equivalent to the grid box dimension of many global models. The 3DRT results are compared both with the calculations obtained by an Independent Pixel Approximation (IPA) approch and by the Plane Parallel radiative scheme (PP) implemented in ECMWF's Forecast model. The differences between the three calculations are used to assess both problems in current GCM's representation of radiative heating and inaccuracies in the dynamical response of CRM simulations due to the Independent Column Approximation (ICA). The understanding of the mechanisms involved in the main 3DRT/1D differences is the starting point for the future attempt to develop a parameterization procedure.
CRASH3: cosmological radiative transfer through metals
NASA Astrophysics Data System (ADS)
Graziani, L.; Maselli, A.; Ciardi, B.
2013-05-01
Here we introduce CRASH3, the latest release of the 3D radiative transfer code CRASH. In its current implementation, CRASH3 integrates into the reference algorithm the code CLOUDY to evaluate the ionization states of metals, self-consistently with the radiative transfer through H and He. The feedback of the heavy elements on the calculation of the gas temperature is also taken into account, making CRASH3 the first 3D code for cosmological applications which treats self-consistently the radiative transfer through an inhomogeneous distribution of metal-enriched gas with an arbitrary number of point sources and/or a background radiation. The code has been tested in idealized configurations, as well as in a more realistic case of multiple sources embedded in a polluted cosmic web. Through these validation tests, the new method has been proven to be numerically stable and convergent. We have studied the dependence of the results on a number of physical quantities such as the source characteristics (spectral range and shape, intensity), the metal composition, the gas number density and metallicity.
Infrared radiative energy transfer in gaseous systems
NASA Technical Reports Server (NTRS)
Tiwari, Surendra N.
1991-01-01
Analyses and numerical procedures are presented to investigate the radiative interactions in various energy transfer processes in gaseous systems. Both gray and non-gray radiative formulations for absorption and emission by molecular gases are presented. The gray gas formulations are based on the Planck mean absorption coefficient and the non-gray formulations are based on the wide band model correlations for molecular absorption. Various relations for the radiative flux and divergence of radiative flux are developed. These are useful for different flow conditions and physical problems. Specific plans for obtaining extensive results for different cases are presented. The procedure developed was applied to several realistic problems. Results of selected studies are presented.
Three Dimensional Atmospheric Radiative Transfer-Applications and Methods Comparison
NASA Technical Reports Server (NTRS)
Cahalan, Robert F.; Einaudi, Franco (Technical Monitor)
2001-01-01
We review applications of 3D radiative transfer in the atmosphere, emphasizing the wide spectrum of scales important to remote sensing and modeling of cloud fields, and the characteristic scales introduced into observed radiances and fluxes by the distribution of photon pathlengths at conservative and absorbing wavelengths. We define the "plane-parallel bias", which is a measure of the importance of 3D cloud structure in large-scale models, and the "independent pixel errors" that quantify the significance of 3D effects in remote sensing, and emphasize their relative magnitude and scale dependence. A variety of approaches in current use in 3D radiative transfer, and issues of speed, accuracy, and flexibility are summarized. We also describe a recently initiated "International Intercomparison of 3-Dimensional Radiation Codes", or I3RC. I3RC is a 3-phase effort that has as its goals to: (1) understand the errors and limits of 3D methods; (2) provide "baseline" cases for future 3D code development; (3) promote sharing of 3D tools; (4) derive guidelines for 3D tool selection; and (5) improve atmospheric science education in 3D radiative transfer. Selected results from Phases 1 and 2 of I3RC are discussed. These are taken from five cloud fields: a 1D field of bar clouds, a 2D radar-derived field, a 3D Landsat-derived field, a stratiform cloud from the model of C. Moeng, and a convective cloud from the model of B. Stevens. Computations have been carried out for three monochromatic wavelengths (one conservative, one absorptive, and one thermal) and two solar zenith angles (0, 60 degrees).
RRTM: A rapid radiative transfer model
Mlawer, E.J.; Taubman, S.J.; Clough, S.A.
1996-04-01
A rapid radiative transfer model (RRTM) for the calculation of longwave clear-sky fluxes and cooling rates has been developed. The model, which uses the correlated-k method, is both accurate and computationally fast. The foundation for RRTM is the line-by-line radiative transfer model (LBLRTM) from which the relevant k-distributions are obtained. LBLRTM, which has been extensively validated against spectral observations e.g., the high-resolution sounder and the Atmospheric Emitted Radiance Interferometer, is used to validate the flux and cooling rate results from RRTM. Validations of RRTM`s results have been performed for the tropical, midlatitude summer, and midlatitude winter atmospheres, as well as for the four Intercomparison of Radiation Codes in Climate Models (ICRCCM) cases from the Spectral Radiance Experiment (SPECTRE). Details of some of these validations are presented below. RRTM has the identical atmospheric input module as LBLRTM, facilitating intercomparisons with LBLRTM and application of the model at the Atmospheric Radiation Measurement Cloud and Radiation Testbed sites.
Modeling of Radiative Transfer in Protostellar Disks
NASA Technical Reports Server (NTRS)
VonAllmen, Paul; Turner, Neal
2007-01-01
This program implements a spectral line, radiative transfer tool for interpreting Spitzer Space Telescope observations by matching them with models of protostellar disks for improved understanding of planet and star formation. The Spitzer Space Telescope detects gas phase molecules in the infrared spectra of protostellar disks, with spectral lines carrying information on the chemical composition of the material from which planets form. Input to the software includes chemical models developed at JPL. The products are synthetic images and spectra for comparison with Spitzer measurements. Radiative transfer in a protostellar disk is primarily affected by absorption and emission processes in the dust and in molecular gases such as H2, CO, and HCO. The magnitude of the optical absorption and emission is determined by the population of the electronic, vibrational, and rotational energy levels. The population of the molecular level is in turn determined by the intensity of the radiation field. Therefore, the intensity of the radiation field and the population of the molecular levels are inter-dependent quantities. To meet the computational challenges of solving for the coupled radiation field and electronic level populations in disks having wide ranges of optical depths and spatial scales, the tool runs in parallel on the JPL Dell Cluster supercomputer with C++ and Fortran compiler with a Message Passing Interface. Because this software has been developed on a distributed computing platform, the modeling of systems previously beyond the reach of available computational resources is possible.
Radiative transfer in realistic planetary atmospheres. [bibliographies
NASA Technical Reports Server (NTRS)
Plass, G. N.; Kattawar, G. W.
1982-01-01
Some 40 publications that appeared in scientific journals from 1973 to 1981 as well as 45 scientific reports issued during the grant period are listed by title. Topics cover the development of a matrix operator theory of radiative transfer which made possible the exact model calculations of the radiance as a function of height in planetary atmospheres; calculation of the Mie phase matrix for various types of particles as well as for radiance and polarization in planetary atmospheres; analysis of high dispersion spectroscopic observations of Venus; calculation of curves of growth for Venus; the development of a theory for calculating radiative transfer in spherical shell atmospheres; investigations of zonal winds on Venus; and examination of Rayleigh scattering.
Introductory Tools for Radiative Transfer Models
NASA Astrophysics Data System (ADS)
Feldman, D.; Kuai, L.; Natraj, V.; Yung, Y.
2006-12-01
Satellite data are currently so voluminous that, despite their unprecedented quality and potential for scientific application, only a small fraction is analyzed due to two factors: researchers' computational constraints and a relatively small number of researchers actively utilizing the data. Ultimately it is hoped that the terabytes of unanalyzed data being archived can receive scientific scrutiny but this will require a popularization of the methods associated with the analysis. Since a large portion of complexity is associated with the proper implementation of the radiative transfer model, it is reasonable and appropriate to make the model as accessible as possible to general audiences. Unfortunately, the algorithmic and conceptual details that are necessary for state-of-the-art analysis also tend to frustrate the accessibility for those new to remote sensing. Several efforts have been made to have web- based radiative transfer calculations, and these are useful for limited calculations, but analysis of more than a few spectra requires the utilization of home- or server-based computing resources. We present a system that is designed to allow for easier access to radiative transfer models with implementation on a home computing platform in the hopes that this system can be utilized in and expanded upon in advanced high school and introductory college settings. This learning-by-doing process is aided through the use of several powerful tools. The first is a wikipedia-style introduction to the salient features of radiative transfer that references the seminal works in the field and refers to more complicated calculations and algorithms sparingly5. The second feature is a technical forum, commonly referred to as a tiki-wiki, that addresses technical and conceptual questions through public postings, private messages, and a ranked searching routine. Together, these tools may be able to facilitate greater interest in the field of remote sensing.
ASIMUT on line radiative transfer code
NASA Astrophysics Data System (ADS)
Vandaele, A. C.; Neary, L.; Robert, S.; Letocart, V.; Giuranna, M.; Kasaba, Y.
2015-10-01
The CROSS DRIVE project aims to develop an innovative collaborative workspace infrastructure for space missions that will allow distributed scientific and engineering teams to collectively analyse and interpret scientific data as well as execute operations of planetary spacecraft. ASIMUT will be one of the tools that will be made available to the users. Here we describe this radiative transfer code and how it will be integrated into the virtual environment developed within CROSS DRIVE.
Enhancing radiative energy transfer through thermal extraction
NASA Astrophysics Data System (ADS)
Tan, Yixuan; Liu, Baoan; Shen, Sheng; Yu, Zongfu
2016-06-01
Thermal radiation plays an increasingly important role in many emerging energy technologies, such as thermophotovoltaics, passive radiative cooling and wearable cooling clothes [1]. One of the fundamental constraints in thermal radiation is the Stefan-Boltzmann law, which limits the maximum power of far-field radiation to P0 = σT4S, where σ is the Boltzmann constant, S and T are the area and the temperature of the emitter, respectively (Fig. 1a). In order to overcome this limit, it has been shown that near-field radiations could have an energy density that is orders of magnitude greater than the Stefan-Boltzmann law [2-7]. Unfortunately, such near-field radiation transfer is spatially confined and cannot carry radiative heat to the far field. Recently, a new concept of thermal extraction was proposed [8] to enhance far-field thermal emission, which, conceptually, operates on a principle similar to oil immersion lenses and light extraction in light-emitting diodes using solid immersion lens to increase light output [62].Thermal extraction allows a blackbody to radiate more energy to the far field than the apparent limit of the Stefan-Boltzmann law without breaking the second law of thermodynamics. Thermal extraction works by using a specially designed thermal extractor to convert and guide the near-field energy to the far field, as shown in Fig. 1b. The same blackbody as shown in Fig. 1a is placed closely below the thermal extractor with a spacing smaller than the thermal wavelength. The near-field coupling transfers radiative energy with a density greater than σT4. The thermal extractor, made from transparent and high-index or structured materials, does not emit or absorb any radiation. It transforms the near-field energy and sends it toward the far field. As a result, the total amount of far-field radiative heat dissipated by the same blackbody is greatly enhanced above SσT4, where S is the area of the emitter. This paper will review the progress in thermal
Principles of Invariance in Radiative Transfer
NASA Astrophysics Data System (ADS)
Peraiah, A.
1999-09-01
We have reviewed the principle of invariance, its applications and its usefulness for obtaining the radiation field in semi-infinite and finite atmospheres. Various laws of scattering in dispersive media and the consequent radiation field are studied. The H-functions and X- and Y-functions in semi-infinite and finite media respectively are derived in a few cases. The Discrete Space Theory (DST) which is a general form of the Principle of Invariance is described. The method of addition of layers with general properties, is shown to describe all the properties of multiple scattering. A few examples of the application of DST such as polarization, line formation in expanding stellar atmospheres, etc., and a numerical analysis of DST are presented. Other developments in the theory of radiative transfer are briefly described.
Parameterization and Analysis of 3-D Solar Radiative Transfer in Clouds: Final Report
Jerry Y. Harrington
2012-09-21
This document reports on the research that we have done over the course of our two-year project. The report also covers the research done on this project during a 1 year no-cost extension of the grant. Our work has had two main, inter-related thrusts: The first thrust was to characterize the response of stratocumulus cloud structure and dynamics to systematic changes in cloud infrared radiative cooling and solar heating using one-dimensional radiative transfer models. The second was to couple a three-dimensional (3-D) solar radiative transfer model to the Large Eddy Simulation (LES) model that we use to simulate stratocumulus. The purpose of the studies with 3-D radiative transfer was to examine the possible influences of 3-D photon transport on the structure, evolution, and radiative properties of stratocumulus. While 3-D radiative transport has been examined in static cloud environments, few studies have attempted to examine whether the 3-D nature of radiative absorption and emission influence the structure and evolution of stratocumulus. We undertook this dual approach because only a small number of LES simulations with the 3-D radiative transfer model are possible due to the high computational costs. Consequently, LES simulations with a 1-D radiative transfer solver were used in order to examine the portions of stratocumulus parameter space that may be most sensitive to perturbations in the radiative fields. The goal was then to explore these sensitive regions with LES using full 3-D radiative transfer. Our overall goal was to discover whether 3-D radiative processes alter cloud structure and evolution, and whether this may have any indirect implications for cloud radiative properties. In addition, we collaborated with Dr. Tamas Varni, providing model output fields for his attempt at parameterizing 3-D radiative effects for cloud models.
Radiative transfer effects in primordial hydrogen recombination
Ali-Haiemoud, Yacine; Hirata, Christopher M.; Grin, Daniel
2010-12-15
The calculation of a highly accurate cosmological recombination history has been the object of particular attention recently, as it constitutes the major theoretical uncertainty when predicting the angular power spectrum of cosmic microwave background anisotropies. Lyman transitions, in particular the Lyman-{alpha} line, have long been recognized as one of the bottlenecks of recombination, due to their very low escape probabilities. The Sobolev approximation does not describe radiative transfer in the vicinity of Lyman lines to a sufficient degree of accuracy, and several corrections have already been computed in other works. In this paper, we compute the impact of some radiative transfer effects that were previously ignored, or for which previous treatments were incomplete. First, the effect of Thomson scattering in the vicinity of the Lyman-{alpha} line is evaluated, using a full redistribution kernel incorporated into a radiative transfer code. The effect of feedback of distortions generated by the optically thick deuterium Lyman-{alpha} line blueward of the hydrogen line is investigated with an analytic approximation. It is shown that both effects are negligible during cosmological hydrogen recombination. Second, the importance of high-lying, nonoverlapping Lyman transitions is assessed. It is shown that escape from lines above Ly{gamma} and frequency diffusion in Ly{beta} and higher lines can be neglected without loss of accuracy. Third, a formalism generalizing the Sobolev approximation is developed to account for the overlap of the high-lying Lyman lines, which is shown to lead to negligible changes to the recombination history. Finally, the possibility of a cosmological hydrogen recombination maser is investigated. It is shown that there is no such maser in the purely radiative treatment presented here.
Composite biasing in Monte Carlo radiative transfer
NASA Astrophysics Data System (ADS)
Baes, Maarten; Gordon, Karl D.; Lunttila, Tuomas; Bianchi, Simone; Camps, Peter; Juvela, Mika; Kuiper, Rolf
2016-05-01
Biasing or importance sampling is a powerful technique in Monte Carlo radiative transfer, and can be applied in different forms to increase the accuracy and efficiency of simulations. One of the drawbacks of the use of biasing is the potential introduction of large weight factors. We discuss a general strategy, composite biasing, to suppress the appearance of large weight factors. We use this composite biasing approach for two different problems faced by current state-of-the-art Monte Carlo radiative transfer codes: the generation of photon packages from multiple components, and the penetration of radiation through high optical depth barriers. In both cases, the implementation of the relevant algorithms is trivial and does not interfere with any other optimisation techniques. Through simple test models, we demonstrate the general applicability, accuracy and efficiency of the composite biasing approach. In particular, for the penetration of high optical depths, the gain in efficiency is spectacular for the specific problems that we consider: in simulations with composite path length stretching, high accuracy results are obtained even for simulations with modest numbers of photon packages, while simulations without biasing cannot reach convergence, even with a huge number of photon packages.
Accurate radiative transfer calculations for layered media.
Selden, Adrian C
2016-07-01
Simple yet accurate results for radiative transfer in layered media with discontinuous refractive index are obtained by the method of K-integrals. These are certain weighted integrals applied to the angular intensity distribution at the refracting boundaries. The radiative intensity is expressed as the sum of the asymptotic angular intensity distribution valid in the depth of the scattering medium and a transient term valid near the boundary. Integrated boundary equations are obtained, yielding simple linear equations for the intensity coefficients, enabling the angular emission intensity and the diffuse reflectance (albedo) and transmittance of the scattering layer to be calculated without solving the radiative transfer equation directly. Examples are given of half-space, slab, interface, and double-layer calculations, and extensions to multilayer systems are indicated. The K-integral method is orders of magnitude more accurate than diffusion theory and can be applied to layered scattering media with a wide range of scattering albedos, with potential applications to biomedical and ocean optics. PMID:27409700
Yashchuk, Valeriy V.; McKinney, Wayne R.; Takacs, Peter Z.
2008-05-19
We suggest and describe the use of a binary pseudo-random grating as a standard test surface for calibration of the modulation transfer function of microscopes. Results from calibration of a MicromapTM-570 interferometric microscope are presented.
Analysis for radiative heat transfer in a circulating fluidized bed
Steward, F.R.; Couturier, M.F.; Poolpol, S.
1995-12-31
The radiative heat transfer from the particles within a circulating fluidized bed has been determined for a number of different assumptions. Based on temperature profiles measured in an operating circulating fluidized bed burning coal, a procedure for predicting the radiative transfer from the solid particles to a cold wall is recommended. The radiative transfer from the solid particles to a cold wall makes up approximately 50% of the total heat transfer to the wall in a circulating fluidized bed combustor.
Studies of radiative transfer in planetary atmospheres
NASA Technical Reports Server (NTRS)
Irvine, W. M.; Schloerb, F. P.
1984-01-01
Progress is reported in modeling cometary emission in the 18-cm OH transition with specific application and predictions for Comet Halley. Radiative transfer is also being studied in rough and porous media. The kinematics of the cold, dark interstellar cloud Li34N were examined, and CO monitoring of Venus and Mars continues. Analysis of 3.4 mm maps of the lunar surface shows thermal anomalies associated with such surface features as the Crater Copernicus, Mare Imbrium, Mare Nubium, Mare Serenitatis, and Mare Tranquillatis.
Diffusion model for lightning radiative transfer
NASA Technical Reports Server (NTRS)
Koshak, William J.; Solakiewicz, Richard J.; Phanord, Dieudonne D.; Blakeslee, Richard J.
1994-01-01
A one-speed Boltzmann transport theory, with diffusion approximations, is applied to study the radiative transfer properties of lightning in optically thick thunderclouds. Near-infrared (lambda = 0.7774 micrometers) photons associated with a prominent oxygen emission triplet in the lightning spectrum are considered. Transient and spatially complex lightning radiation sources are placed inside a rectangular parallelepiped thundercloud geometry and the effects of multiple scattering are studied. The cloud is assumed to be composed of a homogeneous collection of identical spherical water droplets, each droplet a nearly conservative, anisotropic scatterer. Conceptually, we treat the thundercloud like a nuclear reactor, with photons replaced by neutrons, and utilize standard one-speed neutron diffusion techniques common in nuclear reactor analyses. Valid analytic results for the intensity distribution (expanded in spherical harmonics) are obtained for regions sufficiently far from sources. Model estimates of the arrival-time delay and pulse width broadening of lightning signals radiated from within the cloud are determined and the results are in good agreement with both experimental data and previous Monte Carlo estimates. Additional model studies of this kind will be used to study the general information content of cloud top lightning radiation signatures.
Radiative energy transfer in disordered photonic crystals.
Erementchouk, M V; Deych, L I; Noh, H; Cao, H; Lisyansky, A A
2009-04-29
The difficulty of description of the radiative transfer in disordered photonic crystals arises from the necessity to consider on an equal footing the wave scattering by periodic modulations of the dielectric function and by its random inhomogeneities. We resolve this difficulty by approaching this problem from the standpoint of the general multiple scattering theory in media with an arbitrary regular profile of the dielectric function. We use the general asymptotic solution of the Bethe-Salpeter equation in order to show that for a sufficiently weak disorder the diffusion limit in disordered photonic crystals is presented by incoherent superpositions of the modes of the ideal structure with weights inversely proportional to the respective group velocities. The radiative transfer and the diffusion equations are derived as a relaxation of long scale deviations from this limiting distribution. In particular, it is shown that in general the diffusion is anisotropic unless the crystal has sufficiently rich symmetry, say, the square lattice in 2D or the cubic lattice in 3D. In this case, the diffusion is isotropic and only in this case can the effect of the disorder be characterized by a single mean free path depending on frequency. PMID:21825416
Radiative Transfer and Retrievals in EOF Domain
NASA Technical Reports Server (NTRS)
Liu, Xu; Zhou, Daniel K.; Larar, Allen; Smith, William L.; Schluessel, Peter
2008-01-01
The Infrared Atmospheric Sounding Interferometer (IASI) is a hyperspectral sensor with 8461 spectral channels and a nominal spectral resolution of 0.25 cm(sup -1). It is computationally intensive to perform radiative transfer calculations and inversions using all these channels. We will present a Principal Component-based Radiative Transfer Model (PCRTM) and a retrieval algorithm which perform all the necessary calculations in EOF domain. Since the EOFs are orthogonal to each other, only about 100 principal components are needed to represent the information content of the 8461 channels. The PCRTM provides the EOF coefficients and associated derivatives with respect to atmospheric and surface parameters needed by the inversion algorithm. The inversion algorithm is based on a non-linear Levenberg-Marquardt method with climatology covariance and a priori information as constraints. The retrieved parameters include atmospheric temperature, moisture and ozone profiles, cloud parameters, surface skin temperature, and surface emissivities. To make the retrieval system even more compact and stable. The atmospheric vertical profiles are compressed into the EOF space as well. The surface emissivities are also compressed into EOF space.
Microwave radiative transfer studies of precipitation
NASA Technical Reports Server (NTRS)
Bringi, V. N.; Vivekanandan, J.; Turk, F. Joseph
1993-01-01
Since the deployment of the DMSP SSM/I microwave imagers in 1987, increased utilization of passive microwave radiometry throughout the 10 - 100 GHz spectrum has occurred for measurement of atmospheric constituents and terrestrial surfaces. Our efforts have focused on observations and analysis of the microwave radiative transfer behavior of precipitating clouds. We have focused particular attention on combining both aircraft and SSM/I radiometer imagery with ground-based multiparameter radar observations. As part of this and the past NASA contract, we have developed a multi-stream, polarized radiative transfer model which incorporates scattering. The model has the capability to be initialized with cloud model output or multiparameter radar products. This model provides the necessary 'link' between the passive microwave radiometer and active microwave radar observations. This unique arrangement has allowed the brightness temperatures (TB) to be compared against quantities such as rainfall, liquid/ice water paths, and the vertical structure of the cloud. Quantification of the amounts of ice and water in precipitating clouds is required for understanding of the global energy balance.
A rapid radiative transfer model for reflection of solar radiation
NASA Technical Reports Server (NTRS)
Xiang, X.; Smith, E. A.; Justus, C. G.
1994-01-01
A rapid analytical radiative transfer model for reflection of solar radiation in plane-parallel atmospheres is developed based on the Sobolev approach and the delta function transformation technique. A distinct advantage of this model over alternative two-stream solutions is that in addition to yielding the irradiance components, which turn out to be mathematically equivalent to the delta-Eddington approximation, the radiance field can also be expanded in a mathematically consistent fashion. Tests with the model against a more precise multistream discrete ordinate model over a wide range of input parameters demonstrate that the new approximate method typically produces average radiance differences of less than 5%, with worst average differences of approximately 10%-15%. By the same token, the computational speed of the new model is some tens to thousands times faster than that of the more precise model when its stream resolution is set to generate precise calculations.
New mechanism for mass gap scaling and transfer-matrix study for (1+1)D directed percolation
NASA Astrophysics Data System (ADS)
Henkel, Malte; Privman, Vladimir
1990-10-01
The standard scaling mechanism for matching asymptotic behaviors on the approach to criticality fails for directed percolation, as indicated by transfer-matrix calculations for sizes up to N=15. A formulation is proposed incorporating irrelevant-variable corrections in a new pattern of asymptotic scaling. Numerical data support these theoretical predictions and, among other results, yield a new estimate, λ=0.37+/-0.04, for the leading irrelevant-variable exponent.
Radiative energy transfer in molecular gases
NASA Technical Reports Server (NTRS)
Tiwari, Surendra N.
1992-01-01
Basic formulations, analyses, and numerical procedures are presented to study radiative interactions in gray as well as nongray gases under different physical and flow conditions. After preliminary fluid-dynamical considerations, essential governing equations for radiative transport are presented that are applicable under local and nonlocal thermodynamic equilibrium conditions. Auxiliary relations for relaxation times and spectral absorption models are also provided. For specific applications, several simple gaseous systems are analyzed. The first system considered consists of a gas bounded by two parallel plates having the same temperature. Within the gas there is a uniform heat source per unit volume. For this system, both vibrational nonequilibrium effects and radiation conduction interactions are studied. The second system consists of fully developed laminar flow and heat transfer in a parallel plate duct under the boundary condition of a uniform surface heat flux. For this system, effects of gray surface emittance are studied. With the single exception of a circular geometry, the third system is considered identical to the second system. Here, the influence of nongray walls is also studied.
Radiative Transfer on Mesoscopic Spatial Scales
NASA Astrophysics Data System (ADS)
Gardner, Adam Ronald
Accurate predictions of light transport produced by illumination of turbid media such as biological tissues, cloudy atmospheres, terrestrial surfaces, and soft matter is essential in many applications including remote sensing, functional optical imaging, realistic image synthesis, and materials characterization. The inability to model light transport on mesoscopic scales limits the spatial resolution and information content that can be extracted from optical measurements. While effective approaches exist to model light transport in singly- and diffusely-scattering regimes, modeling light propagation over the mesoscopic spatial scales remains an important challenge. Radiative transfer on these scales must account for the complete 5-dimensional spatial and angular distributions of the radiant field. Here, we present novel stochastic and analytic methods to analyze and predict light propagation in turbid media generated by collimated illumination on mesoscopic scales. We also consider coupled transport problems, resulting from illumination and detection, to facilitate measurement design and inverse problems. Specifically, we introduce a coupled Forward-Adjoint Monte Carlo (cFAMC) method that leverages generalized optical reciprocity to enable the computation of spatially-resolved distributions of light interrogation for specific source-detector pairs. cFAMC can aid the design of optical diagnostic measurements by tailoring the light field to interrogate specific sub-volumes of interest. We use cFAMC to examine the effects of angular resolution on the resulting interrogation distributions and analyze a diagnostically-relevant compact fiber probe design for the detection of epithelial precancer. While Monte Carlo simulation is considered a gold standard method to solve the equation of radiative transfer (ERT), it is computationally expensive. Thus, methods to obtain ERT solutions at lower computational cost are valuable. We introduce a general analytical framework to
A stochastic formation of radiative transfer in clouds
Stephens, G.L.; Gabriel, P.M.
1993-03-01
The research carried out under this award dealt with issues involving deterministic radiative transfer, remote sensing, Stochastic radiative transfer, and parameterization of cloud optical properties. A number of different forms of radiative transfer models in one, two, and three dimensions were developed in an attempt to build an understanding of the radiative transfer in clouds with realistic spatial structure and to determine the key geometrical parameter that influence this transfer. The research conducted also seeks to assess the relative importance of these geometrical effects in contrast to microphysical effects of clouds. The main conclusion of the work is that geometry has a profound influence on all aspects of radiative transfer and the interpretation of this transfer. We demonstrate how this geometry can influence estimate of particle effective radius to the 30-50% level and also how geometry can significantly bias the remote sensing of cloud optical depth.
Planetary Atmosphere Dynamics and Radiative Transfer
NASA Technical Reports Server (NTRS)
Atkinson, David H.
1996-01-01
This research program has dealt with two projects in the field of planetary atmosphere dynamics and radiative energy transfer, one theoretical and one experimental. The first project, in radiative energy transfer, incorporated the capability to isolate and quantify the contribution of individual atmospheric components to the Venus radiative balance and thermal structure to greatly improve the current understanding of the radiative processes occurring within the Venus atmosphere. This is possible by varying the mixing ratios of each gas species, and the location, number density and aerosol size distributions of the clouds. This project was a continuation of the work initiated under a 1992 University Consortium Agreement. Under the just completed grant, work has continued on the use of a convolution-based algorithm that provided the capability to calculate the k coefficients of a gas mixture at different temperatures, pressures and spectral intervals from the separate k-distributions of the individual gas species. The second primary goal of this research dealt with the Doppler wind retrieval for the Successful Galileo Jupiter probe mission in December, 1995. In anticipation of the arrival of Galileo at Jupiter, software development continued to read the radioscience and probe/orbiter trajectory data provided by the Galileo project and required for Jupiter zonal wind measurements. Sample experiment radioscience data records and probe/orbiter trajectory data files provided by the Galileo Radioscience and Navigation teams at the Jet Propulsion Laboratory, respectively, were used for the first phase of the software development. The software to read the necessary data records was completed in 1995. The procedure by which the wind retrieval takes place begins with initial consistency checks of the raw data, preliminary data reductions, wind recoveries, iterative reconstruction of the probe descent profile, and refined wind recoveries. At each stage of the wind recovery
NASA Astrophysics Data System (ADS)
Barlakas, Vasileios; Macke, Andreas; Wendisch, Manfred
2016-07-01
Non-spherical particles in the atmosphere absorb and scatter solar radiation. They change the polarization state of solar radiation depending on their shape, size, chemical composition and orientation. To quantify polarization effects, a new three-dimensional (3D) vector radiative transfer model, SPARTA (Solver for Polarized Atmospheric Radiative Transfer Applications) is introduced and validated against benchmark results. SPARTA employs the statistical forward Monte Carlo technique for efficient column-response pixel-based radiance calculations including polarization for 3D inhomogeneous cloudless and cloudy atmospheres. A sensitivity study has been carried out and exemplarily results are presented for two lidar-based mineral dust fields. The scattering and absorption properties of the dust particles have been computed for spheroids and irregular shaped particles. Polarized radiance fields in two-dimensional (2D) and one-dimensional (1D) inhomogeneous Saharan dust fields have been calculated at 532 nm wavelength. The domain-averaged results of the normalized reflected radiance are almost identical for the 1D and 2D modes. In the areas with large spatial gradient in optical thickness with expected significant horizontal photon transport, the radiance fields of the 2D mode differ by about ±12% for the first Stokes component (radiance, I) and ±8% for the second Stokes component (linear polarization, Q) from the fields of the 1D mode.
Plasma effects in high frequency radiative transfer
Alonso, C.T.
1981-02-08
This paper is intended as a survey of collective plasma processes which can affect the transfer of high frequency radiation in a hot dense plasma. We are rapidly approaching an era when this subject will become important in the laboratory. For pedagogical reasons we have chosen to examine plasma processes by relating them to a particular reference plasma which will consist of fully ionized carbon at a temperature kT=1 KeV (10/sup 70/K) and an electron density N = 3 x 10/sup 23/cm/sup -3/, (which corresponds to a mass density rho = 1 gm/cm/sup 3/ and an ion density N/sub i/ = 5 x 10/sup 22/ cm/sup -3/). We will consider the transport in such a plasma of photons ranging from 1 eV to 1 KeV in energy. Such photons will probably be frequently used as diagnostic probes of hot dense laboratory plasmas.
Radiative Transfer Simulations of Infrared Dark Clouds
NASA Astrophysics Data System (ADS)
Pavlyuchenkov, Yaroslav; Wiebe, Dmitry; Fateeva, Anna; Vasyunina, Tatiana
2011-04-01
The determination of prestellar core structure is often based on observations of (sub)millimeter dust continuum. However, recently the Spitzer Space Telescope provided us with IR images of many objects not only in emission but also in absorption. We developed a technique to reconstruct the density and temperature distributions of protostellar objects based on radiation transfer (RT) simulations both in mm and IR wavelengths. Best-fit model parameters are obtained with the genetic algorithm. We apply the method to two cores of Infrared Dark Clouds and show that their observations are better reproduced by a model with an embedded heating source despite the lack of 70 μm emission in one of these cores. Thus, the starless nature of massive cores can only be established with the careful case-by-case RT modeling.
Flare loop radiative hydrodynamics. III - Nonlocal radiative transfer effects
NASA Technical Reports Server (NTRS)
Canfield, R. C.; Fisher, G. H.; Mcclymont, A. N.
1983-01-01
The study has three goals. The first is to demonstrate that processes exist whose intrinsic nonlocal nature cannot be represented by local approximations. The second is to elucidate the physical nature and origins of these nonlocal processes. The third is to suggest that the methods and results described here may prove useful in constructing semiempirical models of the chromosphere by means more efficient than trial and error. Matrices are computed that describe the effect of a temperature perturbation at an arbitrary point in the loop on density, hydrogen ionized fraction, total radiative loss rate, and radiative loss rate of selected hydrogen lines and continua at all other points. It is found that the dominant nonlocal radiative transfer effects can be separated into flux divergence coefficient effects and upper level population effects. The former are most important when the perturbation takes place in a region of significant opacity. Upper level population effects arise in both optically thick and thin regions in response to nonlocal density, ionization, and interlocking effects.
Simulation of solar radiative transfer in cumulus clouds
Zuev, V.E.; Titov, G.A.
1996-04-01
This work presents a 3-D model of radiative transfer which is used to study the relationship between the spatial distribution of cumulus clouds and fluxes (albedo and transmittance) of visible solar radiation.
Comparison of 1D stagnation solutions to 3D wire-array Z pinch simulations in absence of radiation
NASA Astrophysics Data System (ADS)
Yu, Edmund; Velikovich, Alexander; Maron, Yitzhak
2013-10-01
In the idealized picture of a Z pinch, a cylindrically symmetric plasma shell implodes towards axis. In this 1D (radial) picture, the resulting stagnation is very efficient: all the kinetic energy of the shell converts to internal energy, as for instance in the Noh shock solution or the homogeneous stagnation flow. If we generalize the problem to 2D by deforming the shell from perfectly circular to oblate, the resulting stagnation will not be as efficient. As in the Hiemenz flow, in which a jet of fluid strikes a rigid flat boundary and squirts out to the sides, the more complicated flows allowed in 2D allow flow kinetic energy to redirect rather than stagnate. With this picture in mind, we might expect the stagnation of a wire-array Z pinch, which in actuality forms a highly distorted 3D imploding plasma, to dissipate its kinetic energy inefficiently due to the lack of symmetry, and be indescribable by means of the idealized 1D stagnation solutions. On the other hand, one might expect that if the imploding plasma is sufficiently messy, the non-uniformities might ``wash out,'' allowing a quasi-1D description of the averaged quantities of plasma. In this work we explore this idea, comparing predictions of 1D stagnation solutions with 3D simulation. Sandia is a multiprogram laboratory operated by Sandia Corporation, a Lockheed Martin Company, for the United States Department of Energy under contract DE-AC0 4-94AL85000.
Radiative Transfer Modeling and Retrievals for Advanced Hyperspectral Sensors
NASA Technical Reports Server (NTRS)
Liu, Xu; Zhou, Daniel K.; Larar, Allen M.; Smith, William L., Sr.; Mango, Stephen A.
2009-01-01
A novel radiative transfer model and a physical inversion algorithm based on principal component analysis will be presented. Instead of dealing with channel radiances, the new approach fits principal component scores of these quantities. Compared to channel-based radiative transfer models, the new approach compresses radiances into a much smaller dimension making both forward modeling and inversion algorithm more efficient.
Grant, K.E.; Taylor, K.E.; Ellis, J.S.; Wuebbles, D.J.
1987-07-01
The authors have implemented a series of state of the art radiation transport submodels in previously developed one dimensional and two dimensional chemical transport models of the troposphere and stratosphere. These submodels provide the capability of calculating accurate solar and infrared heating rates. They are a firm basis for further radiation submodel development as well as for studying interactions between radiation and model dynamics under varying conditions of clear sky, clouds, and aerosols. 37 refs., 3 figs.
Polar firn layering in radiative transfer models
NASA Astrophysics Data System (ADS)
Linow, Stefanie; Hoerhold, Maria
2016-04-01
For many applications in the geosciences, remote sensing is the only feasible method of obtaining data from large areas with limited accessibility. This is especially true for the cryosphere, where light conditions and cloud coverage additionally limit the use of optical sensors. Here, instruments operating at microwave frequencies become important, for instance in polar snow parameters / SWE (snow water equivalent) mapping. However, the interaction between snow and microwave radiation is a complex process and still not fully understood. RT (radiative transfer) models to simulate snow-microwave interaction are available, but they require a number of input parameters such as microstructure and density, which are partly ill-constrained. The layering of snow and firn introduces an additional degree of complexity, as all snow parameters show a strong variability with depth. Many studies on RT modeling of polar firn deal with layer variability by using statistical properties derived from previous measurements, such as the standard deviations of density and microstructure, to configure model input. Here, the variability of microstructure parameters, such as density and particle size, are usually assumed to be independent of each other. However, in the case of the firn pack of the polar ice sheets, we observe that microstructure evolution depends on environmental parameters, such as temperature and snow deposition. Accordingly, density and microstructure evolve together within the snow and firn. Based on CT (computer tomography) microstructure measurements of antarctic firn, we can show that: first, the variability of density and effective grain size are linked and can thus be implemented in the RT models as a coupled set of parameters. Second, the magnitude of layering is captured by the measured standard deviation. Based on high-resolution density measurements of an Antarctic firn core, we study the effect of firn layering at different microwave wavelengths. By means of
Radiative transfer during the reflooding step of a LOCA
NASA Astrophysics Data System (ADS)
Gérardin, J.; Seiler, N.; Ruyer, P.; Boulet, P.
2013-10-01
Within the evaluation of the heat transfer downstream a quench front during the reflood phase of a Loss of Coolant Accident (LOCA) in a nuclear power plant, a numerical study has been conducted on radiative transfer through a vapor-droplet medium. The non-grey behavior of the medium is obvious since it can be optically thin or thick depending on the wavelength. A six wide bands model has been tested, providing a satisfactory accuracy for the description of the radiative properties. Once the radiative properties of the medium computed, they have been introduced in a model solving the radiative heat transfer based on the Improved Differential Approximation. The fluxes and the flux divergence have been computed on a geometry characteristic of the reactor core showing that radiative transfer plays a relevant role, quite as important as convective heat transfer.
NASA Astrophysics Data System (ADS)
Tong, J.; Yang, J.; Hu, B.
2013-12-01
In this paper, a mixing layer theory is coupled with the Hydrus-1D to predict solute transfer from an initially saturated soil into surface runoff water under a controlled drainage condition with low infiltration rate, and various ponding waterwith the increasing depths of the ponding water on the soil surface before the surface runoff. A sand experiment is used as the reference experiment to identify the parameters for the water flow model and solute transfer model with increasing mixing layer depth after the surface runoff water. Based on these identified parameters andan increasing depth of the mixing layer, the model well predicts the experimental results, and much better than the a simulation results by thewith a constant mixing layer depth under the controlled drainage without drainage theory. The model is also applied to another initially saturated sand experiment under the controlled drainage water condition, and the simulation resultsdata agree well with the observed data. These results suggest that the model can efficiently simulate the solute transfer from the initially saturated soil into the surface runoff water. It is found from this study that with the increase of the surface runoff time or/and the decrease of the mixing layer, the solute loss in the mixing layer before the surface runoff will increase, the solute source in the mixing layer for the surface runoff water will decrease after the surface runoff and the increase rate of the mixing layer depth will increase. Simple sketch of the model Hydraulic parameters in the reference sand experiment
Studies of radiative transfer in planetary atmospheres
NASA Technical Reports Server (NTRS)
Irvine, W. M.; Schloerb, F. P.
1986-01-01
Schloerb and Claussen continued their analysis of the very high quality data set obtained on the 18 centimeter OH line from the Comet P/Halley with the NRAO 43 meter antenna. The high spectral resolution (0.22 km/sec) and high signal-to-noise of the OH spectra make them ideal for the study of kinematics in the coma. Synthetic profiles were initiated for comparison with the data. A vectorial model was developed using the Monte Carlo techniques originated by Combi and Delsemme. Analysis of the millimeter wavelength observations of HCN emission from P/Halley obtained throughout much of the recent apparition were continued using the University of Massachusetts 14 millimeter-wavelength (FCRAO) antenna. A detailed analysis of the HCN lineshpaes was performed over the last six months. The excitation of HCN in the coma was studied to obtain a detailed match to the observed spectra. The passive millimeter wave radiometer was used to probe the physical and chemical nature of comets from spacecraft. Work was continued on an improved theory of radiative transfer for rough and porous surfaces, such as the regoliths of satellites, asteroids, and comets.
Wavelets in the solution of nongray radiative heat transfer equation
Bayazitoglu, Y.; Wang, B.Y.
1996-12-31
The wavelet basis functions are introduced into the radiative transfer equation in the frequency domain. The intensity of radiation is expanded in terms of Daubechies` wrapped around wavelet functions. It is shown that the wavelet basis approach to modeling nongrayness can be incorporated into any solution method for the equation of transfer. In this paper the resulting system of equations is solved for the one-dimensional radiative equilibrium problem using the P-N approximation.
SHDOM: Spherical Harmonic Discrete Ordinate Method for atmospheric radiative transfer
NASA Astrophysics Data System (ADS)
Evans, K. Franklin
2015-08-01
The Spherical Harmonic Discrete Ordinate Method (SHDOM) radiative transfer model computes polarized monochromatic or spectral band radiative transfer in a one, two, or three-dimensional medium for either collimated solar and/or thermal emission sources of radiation. The model is written in a variant of Fortran 77 and in Fortran90 and requires a Fortran 90 compiler. Also included are programs for generating the optical property files input to SHDOM from physical properties of water cloud particles and aerosols.
Application of ray tracing in radiation heat transfer
NASA Technical Reports Server (NTRS)
Baumeister, Joseph F.
1993-01-01
This collection of presentation figures displays the capabilities of ray tracing for radiation propagation calculations as compared to an analytical approach. The goal is to introduce the terminology and solution process used in ray tracing, and provide insight into radiation heat transfer principles and analysis tools. A thermal analysis working environment is introduced that solves demanding radiation heat transfer problems based on ray tracing. This information may serve as a reference for designing and building ones own analysis environment.
Radiative heat transfer in the extreme near field.
Kim, Kyeongtae; Song, Bai; Fernández-Hurtado, Víctor; Lee, Woochul; Jeong, Wonho; Cui, Longji; Thompson, Dakotah; Feist, Johannes; Reid, M T Homer; García-Vidal, Francisco J; Cuevas, Juan Carlos; Meyhofer, Edgar; Reddy, Pramod
2015-12-17
Radiative transfer of energy at the nanometre length scale is of great importance to a variety of technologies including heat-assisted magnetic recording, near-field thermophotovoltaics and lithography. Although experimental advances have enabled elucidation of near-field radiative heat transfer in gaps as small as 20-30 nanometres (refs 4-6), quantitative analysis in the extreme near field (less than 10 nanometres) has been greatly limited by experimental challenges. Moreover, the results of pioneering measurements differed from theoretical predictions by orders of magnitude. Here we use custom-fabricated scanning probes with embedded thermocouples, in conjunction with new microdevices capable of periodic temperature modulation, to measure radiative heat transfer down to gaps as small as two nanometres. For our experiments we deposited suitably chosen metal or dielectric layers on the scanning probes and microdevices, enabling direct study of extreme near-field radiation between silica-silica, silicon nitride-silicon nitride and gold-gold surfaces to reveal marked, gap-size-dependent enhancements of radiative heat transfer. Furthermore, our state-of-the-art calculations of radiative heat transfer, performed within the theoretical framework of fluctuational electrodynamics, are in excellent agreement with our experimental results, providing unambiguous evidence that confirms the validity of this theory for modelling radiative heat transfer in gaps as small as a few nanometres. This work lays the foundations required for the rational design of novel technologies that leverage nanoscale radiative heat transfer. PMID:26641312
Study on radiation transfer in human skin for cosmetics
NASA Astrophysics Data System (ADS)
Yamada, Jun; Kawamura, Ayumu; Miura, Yoshimasa; Takata, Sadaki; Ogawa, Katsuki
2005-06-01
In order to design cosmetics producing the optical properties that are required for a beautiful skin, the radiation transfer in the skin has been numerically investigated by the Monte Carlo method and the effects of skin texture and cosmetics on the radiation transfer have been empirically investigated using an artificial skin. The numerical analysis showed that the total internal reflection suppresses large portion of radiation going out through the skin surface Additionally, the experimental study revealed that skin texture and cosmetics not only diffusely reflect the incoming radiation, but also lead the internally reflected radiation to the outside of the skin.
NASA Technical Reports Server (NTRS)
Egan, Michael P.; Leung, Chun Ming; Spagna, George F., Jr.
1988-01-01
The program solves the radiation transport problem in a dusty medium with one-dimensional planar, spherical or cylindrical geometry. It determines self-consistently the effects of multiple scattering, absorption, and re-emission of photons on the temperature of dust grains and the characteristics of the internal radiation field. The program can treat radiation field anisotropy, linear anisotropic scattering, and multi-grain components. The program output consists of the dust-temperature distribution, flux spectrum, surface brightness at each frequency and the observed intensities (involving a convolution with a telescope beam pattern).
Radiative transfer in atmosphere-sea ice-ocean system
Jin, Z.; Stamnes, K.; Weeks, W.F.; Tsay, S.C.
1996-04-01
Radiative energy is critical in controlling the heat and mass balance of sea ice, which significantly affects the polar climate. In the polar oceans, light transmission through the atmosphere and sea ice is essential to the growth of plankton and algae and, consequently, to the microbial community both in the ice and in the ocean. Therefore, the study of radiative transfer in the polar atmosphere, sea ice, and ocean system is of particular importance. Lacking a properly coupled radiative transfer model for the atmosphere-sea ice-ocean system, a consistent study of the radiative transfer in the polar atmosphere, snow, sea ice, and ocean system has not been undertaken before. The radiative transfer processes in the atmosphere and in the ice and ocean have been treated separately. Because the radiation processes in the atmosphere, sea ice, and ocean depend on each other, this separate treatment is inconsistent. To study the radiative interaction between the atmosphere, clouds, snow, sea ice, and ocean, a radiative transfer model with consistent treatment of radiation in the coupled system is needed and is under development.
A 3D radiative transfer framework. VI. PHOENIX/3D example applications
NASA Astrophysics Data System (ADS)
Hauschildt, P. H.; Baron, E.
2010-01-01
Aims: We demonstrate the application of our 3D radiative transfer framework in the model atmosphere code PHOENIX for a number of spectrum synthesis calculations for very different conditions. Methods: The 3DRT framework discussed in the previous papers of this series was added to our general-purpose model atmosphere code PHOENIX/1D and an extended 3D version PHOENIX/3D was created. The PHOENIX/3D code is parallelized via the MPI library using a hierarchical domain decomposition and displays very good strong scaling. Results: We present the results of several test cases for widely different atmosphere conditions and compare the 3D calculations with equivalent 1D models to assess the internal accuracy of the 3D modeling. In addition, we show the results for a number of parameterized 3D structures. Conclusions: With presently available computational resources it is possible to solve the full 3D radiative transfer (including scattering) problem with the same micro-physics as included in 1D modeling.
Radiation heat transfer in two-phase media
Adzerikho, K.S.
1988-05-01
The state of the art of approximate and numerical methods of the theory of radiation heat transfer is analyzed. The principles for producing engineering methods of computing the radiation heat-transfer characteristics in power plants are examined. These principles include: the integration of the transport equation, computing the radiation heat transfer in nonisothermal two-phase media bounded by emitting and reflecting surfaces, the thermal efficiency of screens as a function of the optical properties of the boundary surfaces and the furnace medium, the scattering processes, temperature distribution, and a program NOTAK in the FORTRAN-IV language.
NASA Astrophysics Data System (ADS)
Wan Ismail, Wan Zakiah; Goldys, Ewa M.; Dawes, Judith M.
2016-02-01
We demonstrate long-wavelength operation (>700 nm) of random dye lasers (using a methylene blue dye) with the addition of rhodamine 6G and titania, enabled by radiative and non-radiative energy transfer. The pump energy is efficiently absorbed and transferred to the acceptors, to support lasing in random dye lasers in the near infrared. The optimum random laser performance with the highest emission intensity and the lowest lasing threshold was achieved for a concentration of methylene blue as the acceptor equal to 6× the concentration of rhodamine 6G (donor). Excessive levels of methylene blue increased the lasing threshold and broadened the methylene blue emission linewidth due to dye quenching from re-absorption. This is due to competition between the donor emission and energy transfer and between absorption loss and fluorescence quenching. The radiative and non-radiative energy transfer is analyzed as a function of the acceptor concentration and pump energy density, with consideration of the spectral overlap. The dependence of the radiative and non-radiative transfer efficiency on the acceptor concentration is obtained, and the energy transfer parameters, including the radiative and non-radiative energy transfer rate constants ( K R and K NR), are investigated using Stern-Volmer analysis. The analysis indicates that radiative energy transfer is the dominant energy transfer mechanism in this system.
Shape-Independent Limits to Near-Field Radiative Heat Transfer.
Miller, Owen D; Johnson, Steven G; Rodriguez, Alejandro W
2015-11-13
We derive shape-independent limits to the spectral radiative heat transfer rate between two closely spaced bodies, generalizing the concept of a blackbody to the case of near-field energy transfer. Through conservation of energy and reciprocity, we show that each body of susceptibility χ can emit and absorb radiation at enhanced rates bounded by |χ|(2)/Im χ, optimally mediated by near-field photon transfer proportional to 1/d(2) across a separation distance d. Dipole-dipole and dipole-plate structures approach restricted versions of the limit, but common large-area structures do not exhibit the material enhancement factor and thus fall short of the general limit. By contrast, we find that particle arrays interacting in an idealized Born approximation (i.e., neglecting multiple scattering) exhibit both enhancement factors, suggesting the possibility of orders-of-magnitude improvement beyond previous designs and the potential for radiative heat transfer to be comparable to conductive heat transfer through air at room temperature, and significantly greater at higher temperatures. PMID:26613444
A Monte Carlo method for 3D thermal infrared radiative transfer
NASA Astrophysics Data System (ADS)
Chen, Y.; Liou, K. N.
2006-09-01
A 3D Monte Carlo model for specific application to the broadband thermal radiative transfer has been developed in which the emissivities for gases and cloud particles are parameterized by using a single cubic element as the building block in 3D space. For spectral integration in the thermal infrared, the correlated k-distribution method has been used for the sorting of gaseous absorption lines in multiple-scattering atmospheres involving 3D clouds. To check the Monte-Carlo simulation, we compare a variety of 1D broadband atmospheric fluxes and heating rates to those computed from the conventional plane-parallel (PP) model and demonstrate excellent agreement between the two. Comparisons of the Monte Carlo results for broadband thermal cooling rates in 3D clouds to those computed from the delta-diffusion approximation for 3D radiative transfer and the independent pixel-by-pixel approximation are subsequently carried out to understand the relative merits of these approaches.
Using 1D theory to understand 3D stagnation of a wire-array Z pinch in the absence of radiation
NASA Astrophysics Data System (ADS)
Yu, Edmund
2015-11-01
Many high-energy-density systems implode towards the axis of symmetry, where it collides on itself, forming a hot plasma. However, experiments show these imploding plasmas develop three-dimensional (3D) structures. As a result, the plasma cannot completely dissipate its kinetic energy at stagnation, instead retaining significant 3D flow. A useful tool for understanding the effects of this residual flow is 3D simulation, but the amount and complexity of information can be daunting. To address this problem, we explore the connection between 3D simulation and one-dimensional (1D) theory. Such a connection, if it exists, is mutually beneficial: 1D theory can provide a clear picture of the underlying dynamics of 3D stagnation. On the other hand, deviations between theory and simulation suggest how 1D theory must be modified to account for 3D effects. In this work, we focus on a 3D, magnetohydrodynamic simulation of a compact wire-array Z pinch. To provide a simpler background against which to test our ideas, we artificially turn off radiation during the stagnation phase. Examination of the initial accumulation of mass on axis reveals oblique collision between jets, shock accretion, and vortex formation. Despite evidence for shock-dominated stagnation, a 1D shockless stagnation solution is more appropriate for describing the global dynamics, in that it reproduces the increase of on-axis density with time. However, the 1D solution must be modified to account for 3D effects: the flows suggest enhanced thermal transport as well as centrifugal force. Upon reaching peak compression, the stagnation transitions to a second phase, in which the high-pressure core on axis expands outward into the remaining imploding plasma. During this phase, a 1D shock solution describes the growth of the shock accretion region, as well as the decrease of on-axis density with time. However, the effect of 3D flows is still present: the on-axis temperature does not cool during expansion, which
grtrans: Polarized general relativistic radiative transfer via ray tracing
NASA Astrophysics Data System (ADS)
Dexter, Jason
2016-05-01
grtrans calculates ray tracing radiative transfer in the Kerr metric, including the full treatment of polarised radiative transfer and parallel transport along geodesics, for comparing theoretical models of black hole accretion flows and jets with observations. The code is written in Fortran 90 and parallelizes with OpenMP; the full code and several components have Python interfaces. grtrans includes Geokerr (ascl:1011.015) and requires cfitsio (ascl:1010.001) and pyfits (ascl:1207.009).
Discrete diffusion Monte Carlo for frequency-dependent radiative transfer
Densmore, Jeffrey D; Kelly, Thompson G; Urbatish, Todd J
2010-11-17
Discrete Diffusion Monte Carlo (DDMC) is a technique for increasing the efficiency of Implicit Monte Carlo radiative-transfer simulations. In this paper, we develop an extension of DDMC for frequency-dependent radiative transfer. We base our new DDMC method on a frequency-integrated diffusion equation for frequencies below a specified threshold. Above this threshold we employ standard Monte Carlo. With a frequency-dependent test problem, we confirm the increased efficiency of our new DDMC technique.
grtrans: Polarized general relativistic radiative transfer via ray tracing
NASA Astrophysics Data System (ADS)
Dexter, Jason
2016-05-01
grtrans calculates ray tracing radiative transfer in the Kerr metric, including the full treatment of polarised radiative transfer and parallel transport along geodesics, for comparing theoretical models of black hole accretion flows and jets with observations. The code is written in Fortran 90 and parallelizes with OpenMP; the full code and several components have Python interfaces. grtrans requires Geokerr (ascl:1011.015), cfitsio (ascl:1010.001), and pyfits (ascl:1207.009).
Effect of radiation heat transfer on thermal diffusivity measurements
NASA Astrophysics Data System (ADS)
Araki, N.
1990-03-01
Experimental data on thermal conductivity and thermal diffusivity of a semitransparent material generally include an error due to the radiation heat transfer. This error varies in accordance with the experimental conditions such as the temperature level of the sample and the measuring method. In this paper, research on the influence of radiation heat transfer on thermal diffusivity are reviewed, and as an example, the method to correct the radiation component in the apparent thermal diffusivity measured by the stepwise heating technique is presented. The transient heat transfer by simultaneous thermal conduction and radiation in a semitransparent material is analyzed when the front surface is subjected to stepwise heating. The apparent thermal diffusivity, which includes the radiation component, is calculated for various parameters.
COMPARING THE EFFECT OF RADIATIVE TRANSFER SCHEMES ON CONVECTION SIMULATIONS
Tanner, Joel D.; Basu, Sarbani; Demarque, Pierre
2012-11-10
We examine the effect of different radiative transfer schemes on the properties of three-dimensional (3D) simulations of near-surface stellar convection in the superadiabatic layer, where energy transport transitions from fully convective to fully radiative. We employ two radiative transfer schemes that fundamentally differ in the way they cover the 3D domain. The first solver approximates domain coverage with moments, while the second solver samples the 3D domain with ray integrations. By comparing simulations that differ only in their respective radiative transfer methods, we are able to isolate the effect that radiative efficiency has on the structure of the superadiabatic layer. We find the simulations to be in good general agreement, but they show distinct differences in the thermal structure in the superadiabatic layer and atmosphere.
NASA Astrophysics Data System (ADS)
Marcq, E.
2012-01-01
In order to understand the early history of telluric interiors and atmospheres during the ocean magma stage, a coupled interior-atmosphere-escape model is being developed. This paper describes the atmospheric part and its first preliminary results. A unidimensional, radiative-convective, H2O-CO2 atmosphere is modeled following a vertical T(z) profile similar to Kasting (1988) and Abe and Matsui (1988). Opacities in the thermal IR are then computed using a k-correlated code (KSPECTRUM), tabulated continuum opacities for H2O-H2O and CO2-CO2 absorption, and water or sulphuric acid clouds in the moist convective zone (whenever present). The first results show the existence of two regimes depending on the relative value of the surface temperature Ts compared to a threshold temperature Tc depending on the total gaseous inventory. For Ts < Tc, efficient blanketing results in a cool upper atmosphere, a cloud cover, and a long lifetime for the underneath magma ocean with a net thermal IR flux between 160 and 200 Wm-2. For Ts > Tc, the blanketing is not efficient enough to prevent large radiative heat loss to space through a hot, cloudless atmosphere. Our current calculations may underestimate the thermal flux in the case of hot surfaces with little gaseous content in the atmosphere.
Verley, Jason C.; Axness, Carl L.; Hembree, Charles Edward; Keiter, Eric Richard; Kerr, Bert
2012-04-01
Photocurrent generated by ionizing radiation represents a threat to microelectronics in radiation environments. Circuit simulation tools such as SPICE [1] can be used to analyze these threats, and typically rely on compact models for individual electrical components such as transistors and diodes. Compact models consist of a handful of differential and/or algebraic equations, and are derived by making simplifying assumptions to any of the many semiconductor transport equations. Historically, many photocurrent compact models have suffered from accuracy issues due to the use of qualitative approximation, rather than mathematically correct solutions to the ambipolar diffusion equation. A practical consequence of this inaccuracy is that a given model calibration is trustworthy over only a narrow range of operating conditions. This report describes work to produce improved compact models for photocurrent. Specifically, an analytic model is developed for epitaxial diode structures that have a highly doped subcollector. The analytic model is compared with both numerical TCAD calculations, as well as the compact model described in reference [2]. The new analytic model compares well against TCAD over a wide range of operating conditions, and is shown to be superior to the compact model from reference [2].
NASA Technical Reports Server (NTRS)
Tripp, T. M.; Shemansky, D. E.; James, G. K.; Ajello, J. M.
1991-01-01
Laboratory measurements of EUV emission from electron-excited N2 have been obtained at medium resolution, providing N II EUV emission cross section measurements and allowing the confirmation of recent calculations by Fawcett (1987) indicating the presence of a strong intercombination line in N II at 748.37 A. The most recently available data are used to predict the basic collisional and radiative properties of N II, the plasma diagnostic properties are briefly explored, and radiative cooling coefficients are given. Some basic properties of electron-excited N II and N2 are examined in the EUV in order to diagnose emission spectra of the earth and Titan. The N II emissions in the earth dayglow, particularly at 916 A, are much brighter than current estimates of source rates. The N II 1085 A line in the dayglow contains a significant component from dissociative photoionization excitation. The N II 1085 A, 916 A, and 670 A lines in the Titan dayglow spectrum appear to be compatible with direct electron excitation of N2.
Tripp, T.M.; Shemansky, D.E.; James, G.K.; Ajello, J.M. JPL, Pasadena, CA )
1991-02-01
Laboratory measurements of EUV emission from electron-excited N2 have been obtained at medium resolution, providing N II EUV emission cross section measurements and allowing the confirmation of recent calculations by Fawcett (1987) indicating the presence of a strong intercombination line in N II at 748.37 A. The most recently available data are used to predict the basic collisional and radiative properties of N II, the plasma diagnostic properties are briefly explored, and radiative cooling coefficients are given. Some basic properties of electron-excited N II and N2 are examined in the EUV in order to diagnose emission spectra of the earth and Titan. The N II emissions in the earth dayglow, particularly at 916 A, are much brighter than current estimates of source rates. The N II 1085 A line in the dayglow contains a significant component from dissociative photoionization excitation. The N II 1085 A, 916 A, and 670 A lines in the Titan dayglow spectrum appear to be compatible with direct electron excitation of N2. 61 refs.
General Relativistic Radiative Transfer: Applications to Black-Hole Systems
NASA Technical Reports Server (NTRS)
Wu, Kinwah; Fuerst, Steven V.; Mizuno, Yosuke; Nishikawa, Ken-Ichi; Branduardi-Raymont, Graziella; Lee, Khee-Gan
2007-01-01
We present general relativistic radiation transfer formulations which include opacity effects due to absorption, emission and scattering explicitly. We consider a moment expansions for the transfer in the presence of scattering. The formulation is applied to calculation emissions from accretion and outflows in black-hole systems. Cases with thin accretion disks and accretion tori are considered. Effects, such as emission anisotropy, non-stationary flows and geometrical self-occultation are investigated. Polarisation transfer in curved space-time is discussed qualitatively.
Thermal radiation heat transfer (3rd revised and enlarged edition)
NASA Astrophysics Data System (ADS)
Siegel, Robert; Howell, John R.
This book first reviews the overall aspects and background information related to thermal radiation heat transfer and incorporates new general information, advances in analytical and computational techniques, and new reference material. Coverage focuses on radiation from opaque surfaces, radiation interchange between various types of surfaces enclosing a vacuum or transparent medium, and radiation including the effects of partially transmitting media, such as combustion gases, soot, or windows. Boundary conditions and multiple layers are discussed with information on radiation in materials with nonunity refractive indices.
Radiative Transfer Reconsidered as a Quantum Kinetic Theory Problem
NASA Astrophysics Data System (ADS)
Rosato, J.
2015-12-01
We revisit the radiative transfer theory from first principles approach, inspired from quantum kinetic theory. The radiation field is described within the second quantization formalism. A master equation for the radiation density operator is derived and transformed into a balance relation in the phase space, which involves nonlocal terms owing to radiation coherence. In a perturbative framework, we focus on the lowest order term in ℏ-expansion and show that the radiation coherence results in an alteration of the photon group velocity. An application to the formation of hydrogen lines in stellar atmospheres is performed as an illustration.
A study of Monte Carlo radiative transfer through fractal clouds
Gautier, C.; Lavallec, D.; O`Hirok, W.; Ricchiazzi, P.
1996-04-01
An understanding of radiation transport (RT) through clouds is fundamental to studies of the earth`s radiation budget and climate dynamics. The transmission through horizontally homogeneous clouds has been studied thoroughly using accurate, discreet ordinates radiative transfer models. However, the applicability of these results to general problems of global radiation budget is limited by the plane parallel assumption and the fact that real clouds fields show variability, both vertically and horizontally, on all size scales. To understand how radiation interacts with realistic clouds, we have used a Monte Carlo radiative transfer model to compute the details of the photon-cloud interaction on synthetic cloud fields. Synthetic cloud fields, generated by a cascade model, reproduce the scaling behavior, as well as the cloud variability observed and estimated from cloud satellite data.
Radiation Transfer in the Atmosphere: Scattering
NASA Technical Reports Server (NTRS)
Mishchenko, M.; Travis, L.; Lacis, Andrew A.
2014-01-01
Sunlight illuminating the Earth's atmosphere is scattered by gas molecules and suspended particles, giving rise to blue skies, white clouds, and optical displays such as rainbows and halos. By scattering and absorbing the shortwave solar radiation and the longwave radiation emitted by the underlying surface, cloud and aerosol particles strongly affect the radiation budget of the terrestrial climate system. As a consequence of the dependence of scattering characteristics on particle size, morphology, and composition, scattered light can be remarkably rich in information on particle properties and thus provides a sensitive tool for remote retrievals of macro- and microphysical parameters of clouds and aerosols.
NASA Astrophysics Data System (ADS)
Olson, Gordon L.
2012-04-01
When using polynomial expansions for the angular variables in the radiation transport equation, the usual procedure is to truncate the series by setting all higher order terms to zero. At low order, such simple closures may not give the optimum solution. This work tests alternate closures that scale either the time- or spatial-derivatives in the highest order equation. These scale factors can be chosen such that waves propagate at exactly the speed of light in optically thin media. Alternatively, they may be chosen to significantly improve the accuracy of low-order solutions with no additional computational cost. The same scaling procedure and scale factors work in one- and multi-dimensions. In multidimensions, reducing the order of a solution can save significant amounts of computer time.
Prediction of radiative heat transfer in rectangular enclosures
Jamaluddin, A.S.; Smith, P.J.
1987-01-01
Discrete ordinates solutions of the radiative transport equation have been obtained for two- and three-dimensional rectangular enclosures using the S/sub 2/ and S/sub 4/ approximations. Limited evaluations indicate that both S/sub 2/ and S/sub 4/ are suitable for predicting radiative transfer in two-dimensional enclosures. However, for the three-dimensional enclosures the S/sub 2/ approximation is found inadequate. It is inferred that S/sub 4/ or higher order approximations should be used to accurately predict radiative heat transfer in three-dimensional rectangular enclosures.
Spectrally-Invariant Approximation Within Atmospheric Radiative Transfer
NASA Technical Reports Server (NTRS)
Marshak, A.; Knyazikhin, Y.; Chiu, J. C.; Wiscombe, W. J.
2011-01-01
Certain algebraic combinations of single scattering albedo and solar radiation reflected from, or transmitted through, vegetation canopies do not vary with wavelength. These "spectrally invariant relationships" are the consequence of wavelength independence of the extinction coefficient and scattering phase function in vegetation. In general, this wavelength independence does not hold in the atmosphere, but in clouddominated atmospheres the total extinction and total scattering phase function vary only weakly with wavelength. This paper identifies the atmospheric conditions under which the spectrally invariant approximation can accurately describe the extinction. and scattering properties of cloudy atmospheres. The validity of the assumptions and the accuracy of the approximation are tested with ID radiative transfer calculations using publicly available radiative transfer models: Discrete Ordinate Radiative Transfer (DISORT) and Santa Barbara DISORT Atmospheric Radiative Transfer (SBDART). It is shown for cloudy atmospheres with cloud optical depth above 3, and for spectral intervals that exclude strong water vapor absorption, that the spectrally invariant relationships found in vegetation canopy radiative transfer are valid to better than 5%. The physics behind this phenomenon, its mathematical basis, and possible applications to remote sensing and climate are discussed.
Radiative heat transfer in low-dimensional systems -- microscopic mode
NASA Astrophysics Data System (ADS)
Woods, Lilia; Phan, Anh; Drosdoff, David
2013-03-01
Radiative heat transfer between objects can increase dramatically at sub-wavelength scales. Exploring ways to modulate such transport between nano-systems is a key issue from fundamental and applied points of view. We advance the theoretical understanding of radiative heat transfer between nano-objects by introducing a microscopic model, which takes into account the individual atoms and their atomic polarizabilities. This approach is especially useful to investigate nano-objects with various geometries and give a detailed description of the heat transfer distribution. We employ this model to study the heat exchange in graphene nanoribbon/substrate systems. Our results for the distance separations, substrates, and presence of extended or localized defects enable making predictions for tailoring the radiative heat transfer at the nanoscale. Financial support from the Department of Energy under Contract No. DE-FG02-06ER46297 is acknowledged.
Li, Xiangming; Kawamura, Akira; Andrews, Chasity D; Miller, Jessica L; Wu, Douglass; Tsao, Tiffany; Zhang, Min; Oren, Deena; Padte, Neal N; Porcelli, Steven A; Wong, Chi-Huey; Kappe, Stefan H I; Ho, David D; Tsuji, Moriya
2015-09-15
A CD1d-binding glycolipid, α-Galactosylceramide (αGalCer), activates invariant NK T cells and acts as an adjuvant. We previously identified a fluorinated phenyl ring-modified αGalCer analog, 7DW8-5, displaying nearly 100-fold stronger CD1d binding affinity. In the current study, 7DW8-5 was found to exert a more potent adjuvant effect than αGalCer for a vaccine based on radiation-attenuated sporozoites of a rodent malaria parasite, Plasmodium yoelii, also referred to as irradiated P. yoelii sporozoites (IrPySpz). 7DW8-5 had a superb adjuvant effect only when the glycolipid and IrPySpz were conjointly administered i.m. Therefore, we evaluated the effect of distinctly different biodistribution patterns of αGalCer and 7DW8-5 on their respective adjuvant activities. Although both glycolipids induce a similar cytokine response in sera of mice injected i.v., after i.m. injection, αGalCer induces a systemic cytokine response, whereas 7DW8-5 is locally trapped by CD1d expressed by dendritic cells (DCs) in draining lymph nodes (dLNs). Moreover, the i.m. coadministration of 7DW8-5 with IrPySpz results in the recruitment of DCs to dLNs and the activation and maturation of DCs. These events cause the potent adjuvant effect of 7DW8-5, resulting in the enhancement of the CD8(+) T cell response induced by IrPySpz and, ultimately, improved protection against malaria. Our study is the first to show that the colocalization of a CD1d-binding invariant NK T cell-stimulatory glycolipid and a vaccine, like radiation-attenuated sporozoites, in dLN-resident DCs upon i.m. conjoint administration governs the potency of the adjuvant effect of the glycolipid. PMID:26254338
NASA Astrophysics Data System (ADS)
Paletou, F.; Vial, Jean-Claude; Auer, L. H.
1992-11-01
The two dimensional (2D) PRD (Partial Redistribution) radiative transfer code of Auer and Paletou was used to compute the resonance lines of HI, MgII and CaII in quiescent prominences, which are modeled as isothermal freestanding slabs illuminated from the sides as well as from below. PRD and 2D effects are evidenced and compared to Complete Redistribution (CRD) computations for both 1D and 2D geometries. Important edge variations are found at the bottom and the top that should be observed with a spatial resolution of one arcsecond. As in 1D, PRD effects allow for greater penetration of the incident radiation into the layer. The 2D code computes both the radial emergent intensity and the amount of radiation backscattered into the chromosphere. It can accordingly, be used to estimate the visibility of filaments. It will be of special interest to build nonisothermal models and compare, for example the Ly alpha profiles with the SUMER/SOHO (Solar Ultraviolet Measurement of Emitted Radiation)/(Solar and Heliospheric Observatory) observations.
3D Radiative Transfer in Cloudy Atmospheres
NASA Astrophysics Data System (ADS)
Marshak, Alexander; Davis, Anthony
Developments in three-dimensional cloud radiation over the past few decades are assessed and distilled into this contributed volume. Chapters are authored by subject-matter experts who address a broad audience of graduate students, researchers, and anyone interested in cloud-radiation processes in the solar and infrared spectral regions. After two introductory chapters and a section on the fundamental physics and computational techniques, the volume extensively treats two main application areas: the impact of clouds on the Earth's radiation budget, which is an essential aspect of climate modeling; and remote observation of clouds, especially with the advanced sensors on current and future satellite missions. http://www.springeronline.com/alert/article?a=3D1_1fva7w_1j826l_41z_6
Light-Cone Effect of Radiation Fields in Cosmological Radiative Transfer Simulations
NASA Astrophysics Data System (ADS)
Ahn, Kyungjin
2015-02-01
We present a novel method to implement time-delayed propagation of radiation fields in cosmo-logical radiative transfer simulations. Time-delayed propagation of radiation fields requires construction of retarded-time fields by tracking the location and lifetime of radiation sources along the corresponding light-cones. Cosmological radiative transfer simulations have, until now, ignored this "light-cone effect" or implemented ray-tracing methods that are computationally demanding. We show that radiative trans-fer calculation of the time-delayed fields can be easily achieved in numerical simulations when periodic boundary conditions are used, by calculating the time-discretized retarded-time Green's function using the Fast Fourier Transform (FFT) method and convolving it with the source distribution. We also present a direct application of this method to the long-range radiation field of Lyman-Werner band photons, which is important in the high-redshift astrophysics with first stars.
Mishra, Subhash C. . E-mail: scm_iitg@yahoo.com; Roy, Hillol K.
2007-04-10
The lattice Boltzmann method (LBM) was used to solve the energy equation of a transient conduction-radiation heat transfer problem. The finite volume method (FVM) was used to compute the radiative information. To study the compatibility of the LBM for the energy equation and the FVM for the radiative transfer equation, transient conduction and radiation heat transfer problems in 1-D planar and 2-D rectangular geometries were considered. In order to establish the suitability of the LBM, the energy equations of the two problems were also solved using the FVM of the computational fluid dynamics. The FVM used in the radiative heat transfer was employed to compute the radiative information required for the solution of the energy equation using the LBM or the FVM (of the CFD). To study the compatibility and suitability of the LBM for the solution of energy equation and the FVM for the radiative information, results were analyzed for the effects of various parameters such as the scattering albedo, the conduction-radiation parameter and the boundary emissivity. The results of the LBM-FVM combination were found to be in excellent agreement with the FVM-FVM combination. The number of iterations and CPU times in both the combinations were found comparable.
User's Manual: Routines for Radiative Heat Transfer and Thermometry
NASA Technical Reports Server (NTRS)
Risch, Timothy K.
2016-01-01
Determining the intensity and spectral distribution of radiation emanating from a heated surface has applications in many areas of science and engineering. Areas of research in which the quantification of spectral radiation is used routinely include thermal radiation heat transfer, infrared signature analysis, and radiation thermometry. In the analysis of radiation, it is helpful to be able to predict the radiative intensity and the spectral distribution of the emitted energy. Presented in this report is a set of routines written in Microsoft Visual Basic for Applications (VBA) (Microsoft Corporation, Redmond, Washington) and incorporating functions specific to Microsoft Excel (Microsoft Corporation, Redmond, Washington) that are useful for predicting the radiative behavior of heated surfaces. These routines include functions for calculating quantities of primary importance to engineers and scientists. In addition, the routines also provide the capability to use such information to determine surface temperatures from spectral intensities and for calculating the sensitivity of the surface temperature measurements to unknowns in the input parameters.
Axness, Carl L.; Keiter, Eric Richard; Kerr, Bert
2011-04-01
Circuit simulation tools (e.g., SPICE) have become invaluable in the development and design of electronic circuits in radiation environments. These codes are often employed to study the effect of many thousands of devices under transient current conditions. Device-scale simulation tools (e.g., MEDICI) are commonly used in the design of individual semiconductor components, but require computing resources that make their incorporation into a circuit code impossible for large-scale circuits. Analytic solutions to the ambipolar diffusion equation, an approximation to the carrier transport equations, may be used to characterize the transient currents at nodes within a circuit simulator. We present new transient 1D excess carrier density and photocurrent density solutions to the ambipolar diffusion equation for low-level radiation pulses that take into account a finite device geometry, ohmic fields outside the depleted region, and an arbitrary change in the carrier lifetime due to neutron irradiation or other effects. The solutions are specifically evaluated for the case of an abrupt change in the carrier lifetime during or after, a step, square, or piecewise linear radiation pulse. Noting slow convergence of the raw Fourier series for certain parameter sets, we use closed-form formulas for some of the infinite sums to produce 'partial closed-form' solutions for the above three cases. These 'partial closed-form' solutions converge with only a few tens of terms, which enables efficient large-scale circuit simulations.
Axness, Carl L.; Keiter, Eric Richard; Kerr, Bert
2010-07-01
Circuit simulation codes, such as SPICE, are invaluable in the development and design of electronic circuits in radiation environments. These codes are often employed to study the effect of many thousands of devices under transient current conditions. Device-scale simulation codes are commonly used in the design of individual semiconductor components, but computational requirements limit their use to small-scale circuits. Analytic solutions to the ambipolar diffusion equation, an approximation to the carrier transport equations, may be used to characterize the transient currents at nodes within a circuit simulator. We present new analytic transient excess carrier density and photocurrent solutions to the ambipolar diffusion equation for 1-D abrupt-junction pn diodes. These solutions incorporate low-level radiation pulses and take into account a finite device geometry, ohmic fields outside the depleted region, and an arbitrary change in the carrier lifetime due to neutron irradiation or other effects. The solutions are specifically evaluated for the case of an abrupt change in the carrier lifetime during or after, a step, square, or piecewise linear radiation pulse. Noting slow convergence of the Fourier series solutions for some parameters sets, we evaluate portions of the solutions using closed-form formulas, which result in a two order of magnitude increase in computational efficiency.
Radiative interactions in transient energy transfer in gaseous systems
NASA Technical Reports Server (NTRS)
Tiwari, S. N.
1985-01-01
Analyses and numerical procedures are presented to investigate the radiative interactions in transient energy transfer processes in gaseous systems. The nongray radiative formulations are based on the wide-band model correlations for molecular absorption. Various relations for the radiative flux are developed; these are useful for different flow conditions and physical problems. Specific plans for obtaining extensive results for different cases are presented. The methods presented in this study can be extended easily to investigate the radiative interactions in realistic flows of hydrogen-air species in the scramjet engine.
Heat Transfer Analysis of a Closed Brayton Cycle Space Radiator
NASA Technical Reports Server (NTRS)
Juhasz, Albert J.
2007-01-01
This paper presents a mathematical analysis of the heat transfer processes taking place in a radiator for a closed cycle gas turbine (CCGT), also referred to as a Closed Brayton Cycle (CBC) space power system. The resulting equations and relationships have been incorporated into a radiator sub-routine of a numerical triple objective CCGT optimization program to determine operating conditions yielding maximum cycle efficiency, minimum radiator area and minimum overall systems mass. Study results should be of interest to numerical modeling of closed cycle Brayton space power systems and to the design of fluid cooled radiators in general.
Radiation Heat Transfer Procedures for Space-Related Applications
NASA Technical Reports Server (NTRS)
Chai, John C.
2000-01-01
Over the last contract year, a numerical procedure for combined conduction-radiation heat transfer using unstructured grids has been developed. As a result of this research, one paper has been published in the Numerical Heat Transfer Journal. One paper has been accepted for presentation at the International Center for Heat and Mass Transfer's International Symposium on Computational Heat Transfer to be held in Australia next year. A journal paper is under review by my NASA's contact. A conference paper for the ASME National Heat Transfer conference is under preparation. In summary, a total of four (4) papers (two journal and two conference) have been published, accepted or are under preparation. There are two (2) to three (3) more papers to be written for the project. In addition to the above publications, one book chapter, one journal paper and six conference papers have been published as a result of this project. Over the last contract year, the research project resulted in one Ph.D. thesis and partially supported another Ph.D. student. My NASA contact and myself have formulated radiation heat transfer procedures for materials with different indices of refraction and for combined conduction-radiation heat transfer. We are trying to find other applications for the procedures developed under this grant.
Partial moment entropy approximation to radiative heat transfer
Frank, Martin . E-mail: frank@mathematik.uni-kl.de; Dubroca, Bruno . E-mail: Bruno.Dubroca@math.u-bordeaux.fr; Klar, Axel . E-mail: klar@mathematik.uni-kl.de
2006-10-10
We extend the half moment entropy closure for the radiative heat transfer equations presented in Dubroca and Klar [B. Dubroca, A. Klar, Half moment closure for radiative transfer equations, J. Comput. Phys. 180 (2002) 584-596] and Turpault et al. [R. Turpault, M. Frank, B. Dubroca, A. Klar, Multigroup half space moment approximations to the radiative heat transfer equations, J. Comput. Phys. 198 (2004) 363-371] to multi-D. To that end, we consider a partial moment system with general partitions of the unit sphere closed by an entropy minimization principle. We give physical and mathematical reasons for this choice of model and study its properties. Several numerical examples in different physical regimes are presented.
Radiative transfer theory for polarimetric remote sensing of pine forest
NASA Technical Reports Server (NTRS)
Hsu, C. C.; Han, H. C.; Shin, Robert T.; Kong, Jin AU; Beaudoin, A.; Letoan, T.
1992-01-01
The radiative transfer theory is applied to interpret polarimetric radar backscatter from pine forest with clustered vegetation structures. To take into account the clustered structures with the radiative transfer theory, the scattering function of each cluster is calculated by incorporating the phase interference of scattered fields from each component. Subsequently, the resulting phase matrix is used in the radiative transfer equations to evaluate the polarimetric backscattering coefficients from random medium layers embedded with vegetation clusters. Upon including the multi-scale structures, namely, trunks, primary and secondary branches, as well as needles, we interpret and simulate the polarimetric radar responses from pine forest for different frequencies and looking angles. The preliminary results are shown to be in good agreement with the measured backscattering coefficients at the Landes maritime pine forest during the MAESTRO-1 experiment.
Estimation of radiative heat transfer using a geometric human model.
Kakuta, N; Yokoyama, S; Nakamura, M; Mabuchi, K
2001-03-01
In order to provide a detailed estimate of radiative heat transfer between a human body and its surrounding environment, we have developed a geometric model of a human form and an algorithm. The model closely resembles the actual shape of a human body and is composed of small quadrilateral surfaces. Dealing with an object or a space with an arbitrary shape, the developed algorithm can judge efficiently whether there is an obstruction between a pair of surfaces. As a result, the angle factors between a pair of surfaces that only occur during radiative heat transfer can be defined. The distribution of the radiative heat transfer rates shows the characteristics of body shape and variations in posture. PMID:11327500
Many-body radiative heat transfer theory.
Ben-Abdallah, Philippe; Biehs, Svend-Age; Joulain, Karl
2011-09-01
In this Letter, an N-body theory for the radiative heat exchange in thermally nonequilibrated discrete systems of finite size objects is presented. We report strong exaltation effects of heat flux which can be explained only by taking into account the presence of many-body interactions. Our theory extends the standard Polder and van Hove stochastic formalism used to evaluate heat exchanges between two objects isolated from their environment to a collection of objects in mutual interaction. It gives a natural theoretical framework to investigate the photon heat transport properties of complex systems at the mesoscopic scale. PMID:22026672
Radiative heat transfer in coal furnaces
Ahluwalia, R.K.; Im, K.H.
1992-01-01
A hybrid technique has been developed to solve three-dimensional spectral radiation transport equations for absorbing, emitting and anisotropically scattering media. An optimal mix of computational speed and accuracy is obtained by combining the discrete ordinate method (S{sub 4}), modified differential approximation (MDA) and P{sub 1} approximation for use in different range of optical thicknesses. The technique is used in conjunction with a char burnout model and spectroscopic data for H{sub 2}O, CO{sub 2}, CO, char, soot and ash to determine the influence of ash composition, ash content and coal preparation on furnace heat absorption.
Radiative heat transfer in coal furnaces
Ahluwalia, R.K.; Im, K.H.
1992-09-01
A hybrid technique has been developed to solve three-dimensional spectral radiation transport equations for absorbing, emitting and anisotropically scattering media. An optimal mix of computational speed and accuracy is obtained by combining the discrete ordinate method (S{sub 4}), modified differential approximation (MDA) and P{sub 1} approximation for use in different range of optical thicknesses. The technique is used in conjunction with a char burnout model and spectroscopic data for H{sub 2}O, CO{sub 2}, CO, char, soot and ash to determine the influence of ash composition, ash content and coal preparation on furnace heat absorption.
Combined conduction and radiation heat transfer in concentric cylindrical media
NASA Technical Reports Server (NTRS)
Pandey, D. K.
1987-01-01
The exact radiative transfer expressions for gray and nongray gases which are absorbing, emitting and nonscattering, contained between infinitely long concentric cylinders with black surfaces, are given in local thermodynamic equilibrium. Resulting energy equations due to the combination of conduction and radiation modes of heat transfer, under steady state conditions for gray and nongray media, are solved numerically using the undetermined parameters method. A single 4.3-micron band of CO2 is considered for the nongray problems. The present solutions for gray and nongray gases obtained in the plane-parallel limit (radius ratio approaches to one) are compared with the plane-parallel results reported in the literature.
Computation of Radiation Heat Transfer in Aeroengine Combustors
NASA Technical Reports Server (NTRS)
Patankar, S. V.
1996-01-01
In this report the highlights of the research completed for the NASA are summarized. This research has been completed in the form of two Ph.D. theses by Chai (1994) and Parthasarathy (1996). Readers are referred to these theses for a complete details of the work and lists of references. In the following sections, first objectives of this research are introduced, then the finite-volume method for radiation heat transfer is described, and finally computations of radiative heat transfer in non-gray participating media is presented.
A simplified scheme for computing radiation transfer in the troposphere
NASA Technical Reports Server (NTRS)
Katayama, A.
1973-01-01
A scheme is presented, for the heating of clear and cloudy air by solar and infrared radiation transfer, designed for use in tropospheric general circulation models with coarse vertical resolution. A bulk transmission function is defined for the infrared transfer. The interpolation factors, required for computing the bulk transmission function, are parameterized as functions of such physical parameters as the thickness of the layer, the pressure, and the mixing ratio at a reference level. The computation procedure for solar radiation is significantly simplified by the introduction of two basic concepts. The first is that the solar radiation spectrum can be divided into a scattered part, for which Rayleigh scattering is significant but absorption by water vapor is negligible, and an absorbed part for which absorption by water vapor is significant but Rayleigh scattering is negligible. The second concept is that of an equivalent cloud water vapor amount which absorbs the same amount of radiation as the cloud.
Introduction of Parallel GPGPU Acceleration Algorithms for the Solution of Radiative Transfer
NASA Technical Reports Server (NTRS)
Godoy, William F.; Liu, Xu
2011-01-01
General-purpose computing on graphics processing units (GPGPU) is a recent technique that allows the parallel graphics processing unit (GPU) to accelerate calculations performed sequentially by the central processing unit (CPU). To introduce GPGPU to radiative transfer, the Gauss-Seidel solution of the well-known expressions for 1-D and 3-D homogeneous, isotropic media is selected as a test case. Different algorithms are introduced to balance memory and GPU-CPU communication, critical aspects of GPGPU. Results show that speed-ups of one to two orders of magnitude are obtained when compared to sequential solutions. The underlying value of GPGPU is its potential extension in radiative solvers (e.g., Monte Carlo, discrete ordinates) at a minimal learning curve.
Theory of heat transfer and hydraulic resistance of oil radiators
NASA Technical Reports Server (NTRS)
Mariamov, N B
1942-01-01
In the present report the coefficients of heat transfer and hydraulic resistance are theoretically obtained for the case of laminar flow of a heated viscous liquid in a narrow rectangular channel. The results obtained are applied to the computation of oil radiators, which to a first approximation may be considered as made up of a system of such channels. In conclusion, a comparison is given of the theoretical with the experimental results obtained from tests on airplane oil radiators.
Modelling of Radiation Heat Transfer in Reacting Hot Gas Flows
NASA Astrophysics Data System (ADS)
Thellmann, A.; Mundt, C.
2009-01-01
In this work the interaction between a turbulent flow including chemical reactions and radiation transport is investigated. As a first step, the state-of-the art radiation models P1 based on the moment method and Discrete Transfer Model (DTM) based on the discrete ordinate method are used in conjunction with the CFD code ANSYS CFX. The absorbing and emitting medium (water vapor) is modeled by Weighted Sum of Gray Gases. For the chemical reactions the standard Eddy dissipation model combined with the two equation turbulence model k-epsilon is employed. A demonstration experiment is identified which delivers temperature distribution, species concentration and radiative intensity distribution in the investigated combustion enclosure. The simulation results are compared with the experiment and reveals that the P1 model predicts the location of the maximal radiation intensity unphysically. On the other hand the DTM model does better but over predicts the maximum value of the radiation intensity. This radiation sensitivity study is a first step on the way to identify a suitable radiation transport and spectral model in order to implement both in an existing 3D Navier-Stokes Code. Including radiation heat transfer we intend to investigate the influence on the overall energy balance in a hydrogen/oxygen rocket combustion chamber.
Super-Eddington radiation transfer in soft gamma repeaters
NASA Astrophysics Data System (ADS)
Ulmer, Andrew
1994-12-01
Bursts from soft gamma repeaters (SGRs) have been shown to be super-Eddington by a factor of 1000 and have been persuasively associated with compact objects. Super-Eddington radiation transfer on the surface of a strongly magnetic (greater than or equal to 1013 G) neutron star is studied and related to the observational constraints on SGRs. In strong magnetic fields, Thompson scattering is suppressed in one polarization state, so super-Eddington fluxes can be radiated while the plasma remains in hydrostatic equilibrium. We discuss a model which offers a somewhat natural explanation for the observation that the energy spectra of bursts with varying intensity are similar. The radiation produced is found to be linearly polarized to one part in 1000 in a direction determined by the local magnetic field, and intensity variations between bursts are understood as a change in the radiating area on the source. The net polarization is inversely correlated with burst intensity. Further, it is shown that for radiation transfer calculations in limit of superstrong magnetic fields, it is sufficient to solve the radiation transfer for the low opacity state rather than the coupled equations for both. With this approximation, standard stellar atmosphere techniques are utilized to calculate the model energy spectrum.
Gopinath, A.; Sadhal, S.S.; Jones, P.D.; Seyed-Yagoobi, J.; Woodbury, K.A.
1996-12-31
In the first section on heat transfer in microgravity, the papers cover phase-change phenomena and thermocapillary flows and surface effects. In the second section, several papers cover solution methods for radiative heat transfer while the rest cover heat transfer in low-temperature environments. The last section covers papers containing valuable information for thermal contact conductance of various materials plus papers on inverse problems in heat transfer. Separate abstracts were prepared for most papers in this volume.
Mesoscopic near-field radiative heat transfer at low temperatures
NASA Astrophysics Data System (ADS)
Maasilta, Ilari; Geng, Zhuoran; Chaudhuri, Saumyadip; Koppinen, Panu
2015-03-01
Near-field radiative heat transfer has mostly been discussed at room temperatures and/or macroscopic scale geometries. Here, we discuss our recent theoretical and experimental advances in understanding near-field transfer at ultra-low temperatures below 1K. As the thermal wavelengths increase with lowering temperature, we show that with sensitive tunnel junction bolometers it is possible to study near-field transfer up to distances ~ 10 μm currently, even though the power levels are low. In addition, these type of experiments correspond to the extreme near-field limit, as the near-field region starts at ~ mm distances at 0.1 K, and could have theoretical power enhancement factors of the order of 1010. Preliminary results on heat transfer between two parallel metallic wires are presented. We also comment on possible areas were such heat transfer might be relevant, such as densely packed arrays of low-temperature detectors.
Validation of the Poisson Stochastic Radiative Transfer Model
NASA Technical Reports Server (NTRS)
Zhuravleva, Tatiana; Marshak, Alexander
2004-01-01
A new approach to validation of the Poisson stochastic radiative transfer method is proposed. In contrast to other validations of stochastic models, the main parameter of the Poisson model responsible for cloud geometrical structure - cloud aspect ratio - is determined entirely by matching measurements and calculations of the direct solar radiation. If the measurements of the direct solar radiation is unavailable, it was shown that there is a range of the aspect ratios that allows the stochastic model to accurately approximate the average measurements of surface downward and cloud top upward fluxes. Realizations of the fractionally integrated cascade model are taken as a prototype of real measurements.
A modular radiative transfer program for gas filter correlation radiometry
NASA Technical Reports Server (NTRS)
Casas, J. C.; Campbell, S. A.
1977-01-01
The fundamentals of a computer program, simulated monochromatic atmospheric radiative transfer (SMART), which calculates atmospheric path transmission, solar radiation, and thermal radiation in the 4.6 micrometer spectral region, are described. A brief outline of atmospheric absorption properties and line by line transmission calculations is explained in conjunction with an outline of the SMART computational procedures. Program flexibility is demonstrated by simulating the response of a gas filter correlation radiometer as one example of an atmospheric infrared sensor. Program limitations, input data requirements, program listing, and comparison of SMART transmission calculations are presented.
Fractional integration and radiative transfer in a multifractal atmosphere
Naud, C.; Schertzer, D.; Lovejoy, S.
1996-04-01
Recently, Cess et al. (1995) and Ramathan et al. (1995) cited observations which exhibit an anomalous absorption of cloudy skies in comparison with the value predicted by usual models and which thus introduce large uncertainties for climatic change assessments. These observation raise questions concerning the way general circulation models have been tuned for decades, relying on classical methods, of both radiative transfer and dynamical modeling. The observations also tend to demonstrate that homogeneous models are simply not relevant in relating the highly variable properties of clouds and radiation fields. However smoothed, the intensity of cloud`s multi-scattered radiation fields reflect this extreme variability.
Debris disk radiative transfer simulation tool (DDS)
NASA Astrophysics Data System (ADS)
Wolf, S.; Hillenbrand, L. A.
2005-10-01
A WWW interface for the simulation of spectral energy distributions of optically thin dust configurations with an embedded radiative source is presented. The density distribution, radiative source, and dust parameters can be selected either from an internal database or defined by the user. This tool is optimized for studying circumstellar debris disks where large grains (a ≫1 μm) are expected to determine the far-infrared through millimeter dust reemission spectral energy distribution. The tool is available at http://aida28.mpia-hd.mpg.de/~swolf/dds. Catalogue identifier:ADVV Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADVV Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland Licensing provisions:none Computers:PC with Intel(R) XEON(TM) 2.80 GHz processor Operating systems or monitors under which the program has been tested:SUSE Linux 9.1 Programming language used:Fortran 90 (for the main program; furthermore Perl, CGI and HTML) Memory required to execute with typical data:108 words No. of bits in a word:8 No. of lines in distributed program, including test data, etc.:44 636 No. of bytes in distributed program, including test data, etc.: 4 806 280 Distribution format:tar.gz Nature of the physical problem:Simulation of scattered light and thermal reemission in arbitrary optically dust distributions with spherical, homogeneous grains where the dust parameters (optical properties, sublimation temperature, grain size) and SED of the illuminating/heating radiative source can be arbitrarily defined (example application: [S. Wolf, L.A. Hillenbrand, Astrophys. J. 596 (2003) 603]). The program is optimized for studying circumstellar debris disks where large grains (i.e. with large size parameters) are expected to determine the far-infrared through millimeter dust reemission spectral energy distribution. Method of solution:Calculation of the dust temperature distribution and dust reemission and scattering spectrum in the
Radiative transfer simulations of magnetar flare beaming
NASA Astrophysics Data System (ADS)
van Putten, T.; Watts, A. L.; Baring, M. G.; Wijers, R. A. M. J.
2016-05-01
Magnetar giant flares show oscillatory modulations in the tails of their light curves, which can only be explained via some form of beaming. The fireball model for magnetar bursts has been used successfully to fit the phase-averaged light curves of the tails of giant flares, but so far no attempts have been made to fit the pulsations. We present a relatively simple numerical model to simulate beaming of magnetar flare emission. In our simulations, radiation escapes from the base of a fireball trapped in a dipolar magnetic field, and is scattered through the optically thick magnetosphere of the magnetar until it escapes. Beaming is provided by the presence of a relativistic outflow, as well as by the geometry of the system. We find that a simple picture for the relativistic outflow is enough to create the pulse fraction and sharp peaks observed in pulse profiles of magnetar flares, while without a relativistic outflow the beaming is insufficient to explain giant flare rotational modulations.
Radiative transfer simulations of magnetar flare beaming
NASA Astrophysics Data System (ADS)
van Putten, T.; Watts, A. L.; Baring, M. G.; Wijers, R. A. M. J.
2016-09-01
Magnetar giant flares show oscillatory modulations in the tails of their light curves, which can only be explained via some form of beaming. The fireball model for magnetar bursts has been used successfully to fit the phase-averaged light curves of the tails of giant flares, but so far no attempts have been made to fit the pulsations. We present a relatively simple numerical model to simulate beaming of magnetar flare emission. In our simulations, radiation escapes from the base of a fireball trapped in a dipolar magnetic field, and is scattered through the optically thick magnetosphere of the magnetar until it escapes. Beaming is provided by the presence of a relativistic outflow, as well as by the geometry of the system. We find that a simple picture for the relativistic outflow is enough to create the pulse fraction and sharp peaks observed in pulse profiles of magnetar flares, while without a relativistic outflow the beaming is insufficient to explain giant flare rotational modulations.
One-dimensional transient radiative transfer by lattice Boltzmann method.
Zhang, Yong; Yi, Hongliang; Tan, Heping
2013-10-21
The lattice Boltzmann method (LBM) is extended to solve transient radiative transfer in one-dimensional slab containing scattering media subjected to a collimated short laser irradiation. By using a fully implicit backward differencing scheme to discretize the transient term in the radiative transfer equation, a new type of lattice structure is devised. The accuracy and computational efficiency of this algorithm are examined firstly. Afterwards, effects of the medium properties such as the extinction coefficient, the scattering albedo and the anisotropy factor, and the shapes of laser pulse on time-resolved signals of transmittance and reflectance are investigated. Results of the present method are found to compare very well with the data from the literature. For an oblique incidence, the LBM results in this paper are compared with those by Monte Carlo method generated by ourselves. In addition, transient radiative transfer in a two-Layer inhomogeneous media subjected to a short square pulse irradiation is investigated. At last, the LBM is further extended to study the transient radiative transfer in homogeneous medium with a refractive index discontinuity irradiated by the short pulse laser. Several trends on the time-resolved signals different from those for refractive index of 1 (i.e. refractive-index-matched boundary) are observed and analysed. PMID:24150298
RADMC: A 2-D Continuum Radiative Transfer Tool
NASA Astrophysics Data System (ADS)
Dullemond, C. P.
2011-08-01
RADMC is a 2-D Monte-Carlo code for dust continuum radiative transfer circumstellar disks and envelopes. It is based on the method of Bjorkman & Wood (ApJ 2001, 554, 615), but with several modifications to produce smoother results with fewer photon packages.
Radiation effects in low-thrust orbit transfers
Pollard, James E.
1998-01-15
A low-thrust orbit transfer vehicle (OTV) and its payload must be designed to survive in the near-Earth radiation environment for a much longer duration than a conventional upper stage. This paper examines the effects of natural radiation on OTV's using data that have become available since 1991 from the CRRES and APEX satellites. Dose rates for microelectronics in LEO-to-GEO missions are calculated for spiral orbit raising and for multi-impulse transfers. Semiconductor devices that are shielded by less than 2.5 mm of aluminum (0.69 g/cm{sup 2}) are inappropriate for spiral transfers, because they require hardness levels >100 krad (Si). Shield thicknesses of 6-12 mm reduce this requirement to about 10 krad (Si), which is still an order of magnitude higher than the radiation dose in a 10-year mission at GEO with similar shielding. The dose for a multi-impulse LEO-to-GEO transfer is about 10 times smaller than for a spiral transfer. Estimates of single event upset rates and photovoltaic array degradation are also provided.
Optical and radiative-transfer properties of mixed atmospheric aerosols
NASA Astrophysics Data System (ADS)
Degheidy, A. R.; Sallah, M.; Elgarayhi, A.; Shaaban, S. M.
2015-04-01
The optical and radiative-transfer properties of mixed atmospheric aerosols have been investigated. The aerosol medium is considered as a plane-parallel anisotropic scattering medium with diffusive reflecting boundaries and containing an internal radiation source. The basic components are defined by their complex refractive index, a lognormal size distribution and humidity dependence in hygroscopic particles. The aerosol particles are assumed to be spherical, so the scattering parameters in the form of single scattering albedo, asymmetry factor, scattering, absorption, extinction efficiencies and linear anisotropic coefficient are calculated using the Mie theory. The calculations have been performed for individual aerosol particles, internal and external mixing media. Radiation transfer problem through the considered aerosol medium has been solved in terms of the solution of the corresponding source-free problem with simple boundary conditions. For the solution of the source-free problem, the Variational Pomraning-Eddington technique has been employed. The variation of the radiative-transfer properties (partial radiative fluxes at the medium boundaries) have been calculated and represented graphically for the different aerosols with their different mixing states. A comparison of the obtained results versus available published data has been performed and a very good agreement was observed.
Cirrus microphysics and radiative transfer: A case study
NASA Technical Reports Server (NTRS)
Kinne, Stefan A.; Ackerman, Thomas P.; Heymsfield, Andrew J.
1990-01-01
During the Cirrus Intensive Field Operations of FIRE, data collected by the NCAR King Air in the vicinity of Wausau, WI on October 28 were selected to study the influence of cirrus cloud microphysics on radiative transfer and the role of microphysical approximations in radiative transfer models. The instrumentation of the King Air provided, aside from temperature and wind data, up-and downwelling broadband solar and infrared fluxes as well as detailed microphysical data. The aircraft data, supplied every second, are averaged over the 7 legs to represent the properties for that altitude. The resulting vertical profiles, however, suffer from the fact that each leg represents a different cloud column path. Based on the measured microphysical data particle size distributions of equivalent spheres for each cloud level are developed. Accurate radiative transfer calculations are performed, incorporating atmospheric and radiative data from the ground and the stratosphere. Comparing calculated to the measured up- and downwelling fluxes at the seven cloud levels for both the averaged and the three crossover data will help to assess the validity of particle size and shape approximation as they are frequently used to model cirrus clouds. Once agreement is achieved the model results may be applied to determine, in comparison to a cloudfree case, the influence of this particular cirrus on the radiation budget of the earth atmosphere system.
Relativistic radiative transfer and relativistic plane-parallel flows
NASA Astrophysics Data System (ADS)
Fukue, Jun
2015-04-01
Relativistic radiative transfer and relativistic plane-parallel flows accelerated from their base like accretion disk winds are numerically examined under the special relativistic treatment. We first solve the relativistic transfer equation iteratively, using a given velocity field, and obtain specific intensities as well as moment quantities. Using the obtained flux, we then solve the hydrodynamical equation, and obtain the new velocity field and the mass-loss rate as an eigen value. We repeat these double-iteration processes until both the intensity and velocity profiles converge. Under this double iteration, we solve the relativistic radiative transfer equation and relativistic flows in the vertical direction, simultaneously. The flows are gradually accelerated, as the optical depth decreases towards the surface. The mass-loss rate dot{J} is roughly expressed in terms of the optical depth τb and terminal speed βs of the flow as dot{J} ˜ 10 τ_b β _s^{-3/4}.
A New Look into the Effect of Large Drops on Radiative Transfer Process
NASA Technical Reports Server (NTRS)
Marshak, Alexander
2003-01-01
Recent studies indicate that a cloudy atmosphere absorbs more solar radiation than any current 1D or 3D radiation model can predict. The excess absorption is not large, perhaps 10-15 W/sq m or less, but any such systematic bias is of concern since radiative transfer models are assumed to be sufficiently accurate for remote sensing applications and climate modeling. The most natural explanation would be that models do not capture real 3D cloud structure and, as a consequence, their photon path lengths are too short. However, extensive calculations, using increasingly realistic 3D cloud structures, failed to produce photon paths long enough to explain the excess absorption. Other possible explanations have also been unsuccessful so, at this point, conventional models seem to offer no solution to this puzzle. The weakest link in conventional models is the way a size distribution of cloud particles is mathematically handled. Basically, real particles are replaced with a single average particle. This "ensemble assumption" assumes that all particle sizes are well represented in any given elementary volume. But the concentration of larger particles can be so low that this assumption is significantly violated. We show how a different mathematical route, using the concept of a cumulative distribution, avoids the ensemble assumption. The cumulative distribution has jumps, or steps, corresponding to the rarer sizes. These jumps result in an additional term, a kind of Green's function, in the solution of the radiative transfer equation. Solving the cloud radiative transfer equation with the measured particle distributions, described in a cumulative rather than an ensemble fashion, may lead to increased cloud absorption of the magnitude observed.
Radiative transfer in a polluted urban planetary boundary layer
NASA Technical Reports Server (NTRS)
Viskanta, R.; Johnson, R. O.; Bergstrom, R. W.
1977-01-01
Radiative transfer in a polluted urban atmosphere is studied using a dynamic model. The diurnal nature of radiative transfer for summer conditions is simulated for an urban area 40 km in extent and the effects of various parameters arising in the problem are investigated. The results of numerical computations show that air pollution has the potential of playing a major role in the radiative regime of the urban area. Absorption of solar energy by aerosols in realistic models of urban atmosphere are of the same order of magnitude as that due to water vapor. The predicted effect of the air pollution aerosol in the city is to warm the earth-atmosphere system, and the net effect of gaseous pollutant is to warm the surface and cool the planetary boundary layer, particularly near the top.
A Thermokinetic Approach to Radiative Heat Transfer at the Nanoscale
Pérez-Madrid, Agustín; Lapas, Luciano C.; Rubí, J. Miguel
2013-01-01
Radiative heat exchange at the nanoscale presents a challenge for several areas due to its scope and nature. Here, we provide a thermokinetic description of microscale radiative energy transfer including phonon-photon coupling manifested through a non-Debye relaxation behavior. We show that a lognormal-like distribution of modes of relaxation accounts for this non-Debye relaxation behavior leading to the thermal conductance. We also discuss the validity of the fluctuation-dissipation theorem. The general expression for the thermal conductance we obtain fits existing experimental results with remarkable accuracy. Accordingly, our approach offers an overall explanation of radiative energy transfer through micrometric gaps regardless of geometrical configurations and distances. PMID:23527019
Realistic three-dimensional radiative transfer simulations of observed precipitation
NASA Astrophysics Data System (ADS)
Adams, I. S.; Bettenhausen, M. H.
2013-12-01
Remote sensing observations of precipitation typically utilize a number of instruments on various platforms. Ground validation campaigns incorporate ground-based and airborne measurements to characterize and study precipitating clouds, while the precipitation measurement constellation envisioned by the Global Precipitation Measurement (GPM) mission includes measurements from differing space-borne instruments. In addition to disparities such as frequency channel selection and bandwidth, measurement geometry and resolution differences between observing platforms result in inherent inconsistencies between data products. In order to harmonize measurements from multiple passive radiometers, a framework is required that addresses these differences. To accomplish this, we have implemented a flexible three-dimensional radiative transfer model. As its core, the radiative transfer model uses the Atmospheric Radiative Transfer Simulator (ARTS) version 2 to solve the radiative transfer equation in three dimensions using Monte Carlo integration. Gaseous absorption is computed with MonoRTM and formatted into look-up tables for rapid processing. Likewise, scattering properties are pre-computed using a number of publicly available codes, such as T-Matrix and DDSCAT. If necessary, a melting layer model can be applied to the input profiles. Gaussian antenna beams estimate the spatial resolutions of the passive measurements, and realistic bandpass characteristics can be included to properly account for the spectral response of the simulated instrument. This work presents three-dimensional simulations of WindSat brightness temperatures for an oceanic rain event sampled by the Tropical Rainfall Measuring Mission (TRMM) satellite. The 2B-31 combined Precipitation Radar / TRMM Microwave Imager (TMI) retrievals provide profiles that are the input to the radiative transfer model. TMI brightness temperatures are also simulated. Comparisons between monochromatic, pencil beam simulations and
A Fast Infrared Radiative Transfer Model for Overlapping Clouds
NASA Technical Reports Server (NTRS)
Niu, Jianguo; Yang, Ping; Huang, Huang-Lung; Davies, James E.; Li, Jun; Baum, Bryan A.; Hu, Yong X.
2006-01-01
A fast infrared radiative transfer model (FIRTM2) appropriate for application to both single-layered and overlapping cloud situations is developed for simulating the outgoing infrared spectral radiance at the top of the atmosphere (TOA). In FIRTM2 a pre-computed library of cloud reflectance and transmittance values is employed to account for one or two cloud layers, whereas the background atmospheric optical thickness due to gaseous absorption can be computed from a clear-sky radiative transfer model. FIRTM2 is applicable to three atmospheric conditions: 1) clear-sky, 2) single-layered ice or water cloud, and 3) two simultaneous cloud layers in a column (e.g., ice cloud overlying water cloud). Moreover, FIRTM2 outputs the derivatives (i.e., Jacobians) of the TOA brightness temperature with respect to cloud optical thickness and effective particle size. Sensitivity analyses have been carried out to assess the performance of FIRTM2 for two spectral regions, namely the longwave (LW) band (587.3 - 1179.5/cm) and the short-to-medium wave (SMW) band (1180.1 - 2228.9/cm). The assessment is carried out in terms of brightness temperature differences (BTD) between FIRTM2 and the well-known discrete ordinates radiative transfer model (DISORT), henceforth referred to as BTD (F-D). The BTD (F-D) values for single-layered clouds are generally less than 0.8 K. For the case of two cloud layers (specifically ice cloud over water cloud), the BTD(F-D) values are also generally less than 0.8 K except for the SMW band for the case of a very high altitude (>15 km) cloud comprised of small ice particles. Note that for clear-sky atmospheres, FIRTM2 reduces to the clear-sky radiative transfer model that is incorporated into FIRTM2, and the errors in this case are essentially those of the clear-sky radiative transfer model.
Coupling radiative heat transfer in participating media with other heat transfer modes
Tencer, John; Howell, John R.
2015-09-28
The common methods for finding the local radiative flux divergence in participating media through solution of the radiative transfer equation are outlined. The pros and cons of each method are discussed in terms of their speed, ability to handle spectral properties and scattering phenomena, as well as their accuracy in different ranges of media transport properties. The suitability of each method for inclusion in the energy equation to efficiently solve multi-mode thermal transfer problems is discussed. Lastly, remaining topics needing research are outlined.
3D-radiative transfer in terrestrial atmosphere: An efficient parallel numerical procedure
NASA Astrophysics Data System (ADS)
Bass, L. P.; Germogenova, T. A.; Nikolaeva, O. V.; Kokhanovsky, A. A.; Kuznetsov, V. S.
2003-04-01
Light propagation and scattering in terrestrial atmosphere is usually studied in the framework of the 1D radiative transfer theory [1]. However, in reality particles (e.g., ice crystals, solid and liquid aerosols, cloud droplets) are randomly distributed in 3D space. In particular, their concentrations vary both in vertical and horizontal directions. Therefore, 3D effects influence modern cloud and aerosol retrieval procedures, which are currently based on the 1D radiative transfer theory. It should be pointed out that the standard radiative transfer equation allows to study these more complex situations as well [2]. In recent year the parallel version of the 2D and 3D RADUGA code has been developed. This version is successfully used in gammas and neutrons transport problems [3]. Applications of this code to radiative transfer in atmosphere problems are contained in [4]. Possibilities of code RADUGA are presented in [5]. The RADUGA code system is an universal solver of radiative transfer problems for complicated models, including 2D and 3D aerosol and cloud fields with arbitrary scattering anisotropy, light absorption, inhomogeneous underlying surface and topography. Both delta type and distributed light sources can be accounted for in the framework of the algorithm developed. The accurate numerical procedure is based on the new discrete ordinate SWDD scheme [6]. The algorithm is specifically designed for parallel supercomputers. The version RADUGA 5.1(P) can run on MBC1000M [7] (768 processors with 10 Gb of hard disc memory for each processor). The peak productivity is equal 1 Tfl. Corresponding scalar version RADUGA 5.1 is working on PC. As a first example of application of the algorithm developed, we have studied the shadowing effects of clouds on neighboring cloudless atmosphere, depending on the cloud optical thickness, surface albedo, and illumination conditions. This is of importance for modern satellite aerosol retrieval algorithms development. [1] Sobolev
Radiation transfer in plant canopies - Scattering of solar radiation and canopy reflectance
NASA Technical Reports Server (NTRS)
Verstraete, Michel M.
1988-01-01
The one-dimensional vertical model of radiation transfer in a plant canopy described by Verstraete (1987) is extended to account for the transfer of diffuse radiation. This improved model computes the absorption and scattering of both visible and near-infrared radiation in a multilayer canopy as a function of solar position and leaf orientation distribution. Multiple scattering is allowed, and the spectral reflectance of the vegetation stand is predicted. The results of the model are compared to those of other models and actual observations.
Cloud Property Retrieval and 3D Radiative Transfer
NASA Technical Reports Server (NTRS)
Cahalan, Robert F.
2003-01-01
Cloud thickness and photon mean-free-path together determine the scale of "radiative smoothing" of cloud fluxes and radiances. This scale is observed as a change in the spatial spectrum of cloud radiances, and also as the "halo size" seen by off beam lidar such as THOR and WAIL. Such of beam lidar returns are now being used to retrieve cloud layer thickness and vertical scattering extinction profile. We illustrate with recent measurements taken at the Oklahoma ARM site, comparing these to the-dependent 3D simulations. These and other measurements sensitive to 3D transfer in clouds, coupled with Monte Carlo and other 3D transfer methods, are providing a better understanding of the dependence of radiation on cloud inhomogeneity, and to suggest new retrieval algorithms appropriate for inhomogeneous clouds. The international "Intercomparison of 3D Radiation Codes" or I3RC, program is coordinating and evaluating the variety of 3D radiative transfer methods now available, and to make them more widely available. Information is on the Web at: http://i3rc.gsfc.nasa.gov/. Input consists of selected cloud fields derived from data sources such as radar, microwave and satellite, and from models involved in the GEWEX Cloud Systems Studies. Output is selected radiative quantities that characterize the large-scale properties of the fields of radiative fluxes and heating. Several example cloud fields will be used to illustrate. I3RC is currently implementing an "open source" 3d code capable of solving the baseline cases. Maintenance of this effort is one of the goals of a new 3DRT Working Group under the International Radiation Commission. It is hoped that the 3DRT WG will include active participation by land and ocean modelers as well, such as 3D vegetation modelers participating in RAMI.
Dana E. Veron
2012-04-09
This project had two primary goals: (1) development of stochastic radiative transfer as a parameterization that could be employed in an AGCM environment, and (2) exploration of the stochastic approach as a means for representing shortwave radiative transfer through mixed-phase layer clouds. To achieve these goals, climatology of cloud properties was developed at the ARM CART sites, an analysis of the performance of the stochastic approach was performed, a simple stochastic cloud-radiation parameterization for an AGCM was developed and tested, a statistical description of Arctic mixed phase clouds was developed and the appropriateness of stochastic approach for representing radiative transfer through mixed-phase clouds was assessed. Significant progress has been made in all of these areas and is detailed in the final report.
Veron, Dana E
2009-03-12
This project had two primary goals: 1) development of stochastic radiative transfer as a parameterization that could be employed in an AGCM environment, and 2) exploration of the stochastic approach as a means for representing shortwave radiative transfer through mixed-phase layer clouds. To achieve these goals, an analysis of the performance of the stochastic approach was performed, a simple stochastic cloud-radiation parameterization for an AGCM was developed and tested, a statistical description of Arctic mixed phase clouds was developed and the appropriateness of stochastic approach for representing radiative transfer through mixed-phase clouds was assessed. Significant progress has been made in all of these areas and is detailed below.
NASA Astrophysics Data System (ADS)
Zhao, J. M.; Tan, J. Y.; Liu, L. H.
2012-02-01
Light transport in graded index media follows a curved trajectory determined by Fermat's principle. Besides the effect of variation of the refractive index on the transport of radiative intensity, the curved ray trajectory will induce geometrical effects on the transport of polarization ellipse. This paper presents a complete derivation of vector radiative transfer equation for polarized radiation transport in absorption, emission and scattering graded index media. The derivation is based on the analysis of the conserved quantities for polarized light transport along curved trajectory and a novel approach. The obtained transfer equation can be considered as a generalization of the classic vector radiative transfer equation that is only valid for uniform refractive index media. Several variant forms of the transport equation are also presented, which include the form for Stokes parameters defined with a fixed reference and the Eulerian forms in the ray coordinate and in several common orthogonal coordinate systems.
NASA Astrophysics Data System (ADS)
Widlowski, J.-L.; Taberner, M.; Pinty, B.; Bruniquel-Pinel, V.; Disney, M.; Fernandes, R.; Gastellu-Etchegorry, J.-P.; Gobron, N.; Kuusk, A.; Lavergne, T.; Leblanc, S.; Lewis, P. E.; Martin, E.; Mõttus, M.; North, P. R. J.; Qin, W.; Robustelli, M.; Rochdi, N.; Ruiloba, R.; Soler, C.; Thompson, R.; Verhoef, W.; Verstraete, M. M.; Xie, D.
2007-05-01
The Radiation Transfer Model Intercomparison (RAMI) initiative benchmarks canopy reflectance models under well-controlled experimental conditions. Launched for the first time in 1999, this triennial community exercise encourages the systematic evaluation of canopy reflectance models on a voluntary basis. The first phase of RAMI focused on documenting the spread among radiative transfer (RT) simulations over a small set of primarily 1-D canopies. The second phase expanded the scope to include structurally complex 3-D plant architectures with and without background topography. Here sometimes significant discrepancies were noted which effectively prevented the definition of a reliable "surrogate truth," over heterogeneous vegetation canopies, against which other RT models could then be compared. The present paper documents the outcome of the third phase of RAMI, highlighting both the significant progress that has been made in terms of model agreement since RAMI-2 and the capability of/need for RT models to accurately reproduce local estimates of radiative quantities under conditions that are reminiscent of in situ measurements. Our assessment of the self-consistency and the relative and absolute performance of 3-D Monte Carlo models in RAMI-3 supports their usage in the generation of a "surrogate truth" for all RAMI test cases. This development then leads (1) to the presentation of the "RAMI Online Model Checker" (ROMC), an open-access web-based interface to evaluate RT models automatically, and (2) to a reassessment of the role, scope, and opportunities of the RAMI project in the future.
SGPGET: AN SBDART Module for Aerosol Radiative Transfer
McComiskey, A.; Ricchiazzi, P.; Ogren, J.A.; Dutton, E.
2005-03-18
Quantification of the aerosol direct effect and climate sensitivity requires accurate estimates of optical properties as inputs to a radiative transfer model. Long-term measurements of aerosol properties at the Southern Great Plains (SGP) site can be used as an improvement over a best guess or global average for optical properties (e.g., asymmetry factor of 0.7) used in Atmospheric Radiation Measurement (ARM) products such as the Broadband Heating Rate Profile VAP. To make this information readily available to the ARM community and others, an add-on module for a commonly used radiative transfer model, SBDART (Ricchiazzi et al. 1998), is being developed. A look up table and algorithm will provide aerosol related model inputs including aerosol optical and atmospheric state properties at high temporal resolution. These inputs can be used in conjunction with any mode of operation and with any other information, for example, cloud properties, in SBDART or any other radiative transfer model. Aerosol properties measured at three visible wavelengths are extrapolated so that flux calculations can be made in any desired wavelength across the shortwave spectrum. Several sources of uncertainty contribute to degraded accuracy of the aerosol property estimation. The effect of these uncertainties is shown through error analysis and comparisons of modeled and observed surface irradiance. A module is also being developed for the North Slope of Alaska site.
Radiative heat transfer between two dielectric-filled metal gratings
NASA Astrophysics Data System (ADS)
Dai, J.; Dyakov, S. A.; Yan, M.
2016-04-01
Nanoscale surface corrugation is known to be able to drastically enhance radiative heat transfer between two metal plates. Here we numerically calculate the radiative heat transfer between two dielectric-filled metal gratings at dissimilar temperatures based on a scattering approach. It is demonstrated that, compared to unfilled metal gratings, the heat flux for a fixed geometry can be further enhanced, by up to 650% for the geometry separated by a vacuum gap of g =1 μ m and temperature values concerned in our study. The enhancement in radiative heat transfer is found to depend on refractive index of the filling dielectric, the specific grating temperatures, and naturally the gap size between the two gratings. The enhancement can be understood through examining the transmission factor spectra, especially the spectral locations of the spoof surface plasmon polariton modes. Of more practical importance, it's shown that the radiative heat flux can exceed that between two planar SiC plates with same thickness, separation, and temperature settings over a wide temperature range. This reaffirms that one can harness rich electromagnetic modal properties in nanostructured materials for efficient thermal management at nanoscale.
Radiative heat transfer in 2D Dirac materials
Rodriguez-López, Pablo; Tse, Wang -Kong; Dalvit, Diego A. R.
2015-05-12
We compute the radiative heat transfer between two sheets of 2D Dirac materials, including topological Chern insulators and graphene, within the framework of the local approximation for the optical response of these materials. In this approximation, which neglects spatial dispersion, we derive both numerically and analytically the short-distance asymptotic of the near-field heat transfer in these systems, and show that it scales as the inverse of the distance between the two sheets. In conclusion, we discuss the limitations to the validity of this scaling law imposed by spatial dispersion in 2D Dirac materials.
SPHRAY: A Smoothed Particle Hydrodynamics Ray Tracer for Radiative Transfer
NASA Astrophysics Data System (ADS)
Altay, Gabriel; Croft, Rupert A. C.; Pelupessy, Inti
2011-03-01
SPHRAY, a Smoothed Particle Hydrodynamics (SPH) ray tracer, is designed to solve the 3D, time dependent, radiative transfer (RT) equations for arbitrary density fields. The SPH nature of SPHRAY makes the incorporation of separate hydrodynamics and gravity solvers very natural. SPHRAY relies on a Monte Carlo (MC) ray tracing scheme that does not interpolate the SPH particles onto a grid but instead integrates directly through the SPH kernels. Given initial conditions and a description of the sources of ionizing radiation, the code will calculate the non-equilibrium ionization state (HI, HII, HeI, HeII, HeIII, e) and temperature (internal energy/entropy) of each SPH particle. The sources of radiation can include point like objects, diffuse recombination radiation, and a background field from outside the computational volume. The MC ray tracing implementation allows for the quick introduction of new physics and is parallelization friendly. A quick Axis Aligned Bounding Box (AABB) test taken from computer graphics applications allows for the acceleration of the raytracing component. We present the algorithms used in SPHRAY and verify the code by performing all the test problems detailed in the recent Radiative Transfer Comparison Project of Iliev et. al. The Fortran 90 source code for SPHRAY and example SPH density fields are made available online.
An Improved Radiative Transfer Model for Climate Calculations
NASA Technical Reports Server (NTRS)
Bergstrom, Robert W.; Mlawer, Eli J.; Sokolik, Irina N.; Clough, Shepard A.; Toon, Owen B.
1998-01-01
This paper presents a radiative transfer model that has been developed to accurately predict the atmospheric radiant flux in both the infrared and the solar spectrum with a minimum of computational effort. The model is designed to be included in numerical climate models To assess the accuracy of the model, the results are compared to other more detailed models for several standard cases in the solar and thermal spectrum. As the thermal spectrum has been treated in other publications, we focus here on the solar part of the spectrum. We perform several example calculations focussing on the question of absorption of solar radiation by gases and aerosols.
Fire Intensity Data for Validation of the Radiative Transfer Equation
Blanchat, Thomas K.; Jernigan, Dann A.
2016-01-01
A set of experiments and test data are outlined in this report that provides radiation intensity data for the validation of models for the radiative transfer equation. The experiments were performed with lightly-sooting liquid hydrocarbon fuels that yielded fully turbulent fires 2 m diameter). In addition, supplemental measurements of air flow and temperature, fuel temperature and burn rate, and flame surface emissive power, wall heat, and flame height and width provide a complete set of boundary condition data needed for validation of models used in fire simulations.
Preliminary results of a three-dimensional radiative transfer model
O`Hirok, W.
1995-09-01
Clouds act as the primary modulator of the Earth`s radiation at the top of the atmosphere, within the atmospheric column, and at the Earth`s surface. They interact with both shortwave and longwave radiation, but it is primarily in the case of shortwave where most of the uncertainty lies because of the difficulties in treating scattered solar radiation. To understand cloud-radiative interactions, radiative transfer models portray clouds as plane-parallel homogeneous entities to ease the computational physics. Unfortunately, clouds are far from being homogeneous, and large differences between measurement and theory point to a stronger need to understand and model cloud macrophysical properties. In an attempt to better comprehend the role of cloud morphology on the 3-dimensional radiation field, a Monte Carlo model has been developed. This model can simulate broadband shortwave radiation fluxes while incorporating all of the major atmospheric constituents. The model is used to investigate the cloud absorption anomaly where cloud absorption measurements exceed theoretical estimates and to examine the efficacy of ERBE measurements and cloud field experiments. 3 figs.
Backward and Forward Monte Carlo Method in Polarized Radiative Transfer
NASA Astrophysics Data System (ADS)
Yong, Huang; Guo-Dong, Shi; Ke-Yong, Zhu
2016-03-01
In general, the Stocks vector cannot be calculated in reverse in the vector radiative transfer. This paper presents a novel backward and forward Monte Carlo simulation strategy to study the vector radiative transfer in the participated medium. A backward Monte Carlo process is used to calculate the ray trajectory and the endpoint of the ray. The Stocks vector is carried out by a forward Monte Carlo process. A one-dimensional graded index semi-transparent medium was presented as the physical model and the thermal emission consideration of polarization was studied in the medium. The solution process to non-scattering, isotropic scattering, and the anisotropic scattering medium, respectively, is discussed. The influence of the optical thickness and albedo on the Stocks vector are studied. The results show that the U, V-components of the apparent Stocks vector are very small, but the Q-component of the apparent Stocks vector is relatively larger, which cannot be ignored.
Spatial-multiblock procedure for radiation heat transfer
Chai, J.C.; Moder, J.P.
1996-12-31
A spatial-multiblock procedure for radiation heat transfer is presented in this article. The proposed procedure is applicable to isothermal or nonisothermal, absorbing, emitting and scattering of transparent media with black or reflecting walls. Although not shown in this article, the procedure is also applicable to nongray conditions. The proposed procedure can be used with the discrete ordinates method and the finite volume method. The heat transfer rate, net radiation power and other full-range and half-range moments are conserved between spatial blocks by the proposed procedure. The utilities of the proposed procedure are shown using four sample problems. The solutions indicate that the multiblock procedure can reproduce the results of a single-block procedure even when very coarse spatial grids are used in the multiblock procedure.
Application of nonlinear Krylov acceleration to radiative transfer problems
Till, A. T.; Adams, M. L.; Morel, J. E.
2013-07-01
The iterative solution technique used for radiative transfer is normally nested, with outer thermal iterations and inner transport iterations. We implement a nonlinear Krylov acceleration (NKA) method in the PDT code for radiative transfer problems that breaks nesting, resulting in more thermal iterations but significantly fewer total inner transport iterations. Using the metric of total inner transport iterations, we investigate a crooked-pipe-like problem and a pseudo-shock-tube problem. Using only sweep preconditioning, we compare NKA against a typical inner / outer method employing GMRES / Newton and find NKA to be comparable or superior. Finally, we demonstrate the efficacy of applying diffusion-based preconditioning to grey problems in conjunction with NKA. (authors)
Radiative transfer theory for polarimetric remote sensing of pine forest
NASA Technical Reports Server (NTRS)
Hsu, C. C.; Han, H. C.; Shin, R. T.; Kong, J. A.; Beaudoin, A.; Le Toan, T.
1992-01-01
The radiative transfer theory is applied to interpret polarimetric radar backscatter from pine forest with clustered vegetation structures. The scattering function of each cluster is calculated by incorporating the phase interference of scattered fields from each component. The resulting phase matrix is used in the radiative transfer equations to evaluate the polarimetric backscattering coefficients from random medium layers embedded with vegetation clusters. Upon including multiscale structures (trunks, primary and secondary branches, and needles), polarimetric radar responses from pine forest for different frequencies and looking angles are interpreted and simulated. Preliminary results are shown to be in good agreement with the measured backscattering coefficients at the Landes maritime pine forest during the MAESTRO-1 experiment.
A Radiation Transfer Solver for Athena Using Short Characteristics
NASA Astrophysics Data System (ADS)
Davis, Shane W.; Stone, James M.; Jiang, Yan-Fei
2012-03-01
We describe the implementation of a module for the Athena magnetohydrodynamics (MHD) code that solves the time-independent, multi-frequency radiative transfer (RT) equation on multidimensional Cartesian simulation domains, including scattering and non-local thermodynamic equilibrium (LTE) effects. The module is based on well known and well tested algorithms developed for modeling stellar atmospheres, including the method of short characteristics to solve the RT equation, accelerated Lambda iteration to handle scattering and non-LTE effects, and parallelization via domain decomposition. The module serves several purposes: it can be used to generate spectra and images, to compute a variable Eddington tensor (VET) for full radiation MHD simulations, and to calculate the heating and cooling source terms in the MHD equations in flows where radiation pressure is small compared with gas pressure. For the latter case, the module is combined with the standard MHD integrators using operator splitting: we describe this approach in detail, including a new constraint on the time step for stability due to radiation diffusion modes. Implementation of the VET method for radiation pressure dominated flows is described in a companion paper. We present results from a suite of test problems for both the RT solver itself and for dynamical problems that include radiative heating and cooling. These tests demonstrate that the radiative transfer solution is accurate and confirm that the operator split method is stable, convergent, and efficient for problems of interest. We demonstrate there is no need to adopt ad hoc assumptions of questionable accuracy to solve RT problems in concert with MHD: the computational cost for our general-purpose module for simple (e.g., LTE gray) problems can be comparable to or less than a single time step of Athena's MHD integrators, and only few times more expensive than that for more general (non-LTE) problems.
A RADIATION TRANSFER SOLVER FOR ATHENA USING SHORT CHARACTERISTICS
Davis, Shane W.; Stone, James M.; Jiang Yanfei
2012-03-01
We describe the implementation of a module for the Athena magnetohydrodynamics (MHD) code that solves the time-independent, multi-frequency radiative transfer (RT) equation on multidimensional Cartesian simulation domains, including scattering and non-local thermodynamic equilibrium (LTE) effects. The module is based on well known and well tested algorithms developed for modeling stellar atmospheres, including the method of short characteristics to solve the RT equation, accelerated Lambda iteration to handle scattering and non-LTE effects, and parallelization via domain decomposition. The module serves several purposes: it can be used to generate spectra and images, to compute a variable Eddington tensor (VET) for full radiation MHD simulations, and to calculate the heating and cooling source terms in the MHD equations in flows where radiation pressure is small compared with gas pressure. For the latter case, the module is combined with the standard MHD integrators using operator splitting: we describe this approach in detail, including a new constraint on the time step for stability due to radiation diffusion modes. Implementation of the VET method for radiation pressure dominated flows is described in a companion paper. We present results from a suite of test problems for both the RT solver itself and for dynamical problems that include radiative heating and cooling. These tests demonstrate that the radiative transfer solution is accurate and confirm that the operator split method is stable, convergent, and efficient for problems of interest. We demonstrate there is no need to adopt ad hoc assumptions of questionable accuracy to solve RT problems in concert with MHD: the computational cost for our general-purpose module for simple (e.g., LTE gray) problems can be comparable to or less than a single time step of Athena's MHD integrators, and only few times more expensive than that for more general (non-LTE) problems.
TWILIGHT: A Cellular Framework for Three-Dimensional Radiative Transfer
NASA Astrophysics Data System (ADS)
Khatami, David; Madore, Barry
2015-01-01
We describe a new framework for solving three-dimensional radiative transfer of arbitrary geometries, including a full characterisation of the wavelength-dependent anisotropic scattering, absorption, and thermal reemission of light by dust. By adopting a cellular approach to discretising the light and dust, the problem can be efficiently solved through a fully deterministic iterative process. As a proof of concept we present TWILIGHT, our implementation of the cellular approach, in order to demonstrate and benchmark the new method. TWILIGHT simultaneously renders over one hundred unique images of a given environment with no additional slowdown, enabling a close study of inclination effects of three-dimensional dust geometries. In addition to qualitative rendering tests, TWILIGHT is successfully tested against two Monte-Carlo radiative transfer benchmarks, producing similar brightness profiles at varying inclinations. With the proof-of-concept established, we describe the improvements and current developments underway using the cellular framework, including a technique to resolve the subgrid physics of dust radiative transfer from micron-scale grain models to kiloparsec-sized dust environments.
Interpreting snowpack radiometry using currently existing microwave radiative transfer models
NASA Astrophysics Data System (ADS)
Kang, Do-Hyuk; Tang, Shurun; Kim, Edward J.
2015-10-01
A radiative transfer model (RTM) to calculate the snow brightness temperatures (Tb) is a critical element in terrestrial snow parameter retrieval from microwave remote sensing observations. The RTM simulates the Tb based on a layered snow by solving a set of microwave radiative transfer equations. Even with the same snow physical inputs to drive the RTM, currently existing models such as Microwave Emission Model of Layered Snowpacks (MEMLS), Dense Media Radiative Transfer (DMRT-QMS), and Helsinki University of Technology (HUT) models produce different Tb responses. To backwardly invert snow physical properties from the Tb, differences from RTMs are first to be quantitatively explained. To this end, this initial investigation evaluates the sources of perturbations in these RTMs, and reveals the equations where the variations are made among the three models. Modelling experiments are conducted by providing the same but gradual changes in snow physical inputs such as snow grain size, and snow density to the 3 RTMs. Simulations are conducted with the frequencies consistent with the Advanced Microwave Scanning Radiometer- E (AMSR-E) at 6.9, 10.7, 18.7, 23.8, 36.5, and 89.0 GHz. For realistic simulations, the 3 RTMs are simultaneously driven by the same snow physics model with the meteorological forcing datasets and are validated against the snow insitu samplings from the CLPX (Cold Land Processes Field Experiment) 2002-2003, and NoSREx (Nordic Snow Radar Experiment) 2009-2010.
Interpreting snowpack radiometry using currently existing microwave radiative transfer models
NASA Astrophysics Data System (ADS)
Kang, D. H.; Tan, S.; Kim, E. J.
2015-12-01
A radiative transfer model (RTM) to calculate a snow brightness temperature (Tb) is a critical element to retrieve terrestrial snow from microwave remote sensing observations. The RTM simulates the Tb based on a layered snow by solving a set of microwave radiative transfer formulas. Even with the same snow physical inputs used for the RTM, currently existing models such as Microwave Emission Model of Layered Snowpacks (MEMLS), Dense Media Radiative Transfer (DMRT-Tsang), and Helsinki University of Technology (HUT) models produce different Tb responses. To backwardly invert snow physical properties from the Tb, the differences from the RTMs are to be quantitatively explained. To this end, the paper evaluates the sources of perturbations in the RTMs, and reveals the equations where the variations are made among three models. Investigations are conducted by providing the same but gradual changes in snow physical inputs such as snow grain size, and snow density to the 3 RTMs. Simulations are done with the frequencies consistent with the Advanced Microwave Scanning Radiometer-E (AMSR-E) at 6.9, 10.7, 18.7, 23.8, 36.5, and 89.0 GHz. For realistic simulations, the 3 RTMs are simultaneously driven by the same snow physics model with the meteorological forcing datasets and are validated from the snow core samplings from the CLPX (Cold Land Processes Field Experiment) 2002-2003, and NoSREx (Nordic Snow Radar Experiment) 2009-2010.
Millimeter wave radiative transfer studies for precipitation measurements
NASA Technical Reports Server (NTRS)
Vivekanandan, J.; Evans, Frank
1989-01-01
Scattering calculations using the discrete dipole approximation and vector radiative transfer calculations were performed to model multiparameter radar return and passive microwave emission for a simple model of a winter storm. The issue of dendrite riming was addressed by computing scattering properties of thin ice disks with varying bulk density. It was shown that C-band multiparameter radar contains information about particle density and the number concentration of the ice particles. The radiative transfer modeling indicated that polarized multifrequency passive microwave emission may be used to infer some properties of ice hydrometers. Detailed radar modeling and vector radiative transfer modeling is in progress to enhance the understanding of simultaneous radar and radiometer measurements, as in the case of the proposed TRMM field program. A one-dimensional cloud model will be used to simulate the storm structure in detail and study the microphysics, such as size and density. Multifrequency polarized radiometer measurements from the SSMI satellite instrument will be analyzed in relation to dual-frequency and dual-polarization radar measurements.
Lattice Boltzmann method for one-dimensional vector radiative transfer.
Zhang, Yong; Yi, Hongliang; Tan, Heping
2016-02-01
A one-dimensional vector radiative transfer (VRT) model based on lattice Boltzmann method (LBM) that considers polarization using four Stokes parameters is developed. The angular space is discretized by the discrete-ordinates approach, and the spatial discretization is conducted by LBM. LBM has such attractive properties as simple calculation procedure, straightforward and efficient handing of boundary conditions, and capability of stable and accurate simulation. To validate the performance of LBM for vector radiative transfer, four various test problems are examined. The first case investigates the non-scattering thermal-emitting atmosphere with no external collimated solar. For the other three cases, the external collimated solar and three different scattering types are considered. Particularly, the LBM is extended to solve VRT in the atmospheric aerosol system where the scattering function contains singularities and the hemisphere space distributions for the Stokes vector are presented and discussed. The accuracy and computational efficiency of this algorithm are discussed. Numerical results show that the LBM is accurate, flexible and effective to solve one-dimensional polarized radiative transfer problems. PMID:26906779
Radiation heat transfer within an optical fiber draw tower furnace
Issa, J.; Jaluria, Y.; Polymeropoulos, C.E.; Yin, Z.
1995-12-31
Study of the thermal transport and material flow processes associated with the drawing of optical fiber in a graphite draw furnace requires modeling of the heat transfer from the furnace wall. Previous work has shown that accurate knowledge of the furnace heater element axial temperature distribution is essential for proper modeling of the radiative transfer process. The present work is aimed at providing this information, as well as generating a set of data for the study of radiation exchange in the furnace cavity. The experimental procedure involved measuring the centerline temperature distribution in graphite and fused silica rods inserted into an optical fiber draw tower furnace. The temperature measurements were then used along with a model for radiative-convective heat transfer in the furnace in order to obtain the furnace temperature profile. This is an inverse problem since the centerline temperature in the rod is known whereas the furnace thermal conditions are not. The results obtained showed that the furnace temperature distribution was independent of rod material and size. The shape of the computed temperature distributions suggest that they can be well represented by a Gaussian function.
A conjugate gradient method for solving the non-LTE line radiation transfer problem
NASA Astrophysics Data System (ADS)
Paletou, F.; Anterrieu, E.
2009-12-01
This study concerns the fast and accurate solution of the line radiation transfer problem, under non-LTE conditions. We propose and evaluate an alternative iterative scheme to the classical ALI-Jacobi method, and to the more recently proposed Gauss-Seidel and successive over-relaxation (GS/SOR) schemes. Our study is indeed based on applying a preconditioned bi-conjugate gradient method (BiCG-P). Standard tests, in 1D plane parallel geometry and in the frame of the two-level atom model with monochromatic scattering are discussed. Rates of convergence between the previously mentioned iterative schemes are compared, as are their respective timing properties. The smoothing capability of the BiCG-P method is also demonstrated.
Radiative Heat Transfer in a Hydrous Transition Zone
NASA Astrophysics Data System (ADS)
Thomas, S.; Bina, C. R.; Jacobsen, S. D.; Goncharov, A. F.
2012-12-01
The structure and dynamics of Earth's interior depend crucially upon heat flow and thus upon the thermal conductivity of its constituents. The bulk thermal conductivity has two components: lattice conductivity (klat) and radiative conductivity (krad) [1,2]. Whereas lattice conductivity is governed by phonon propagation, radiative conductivity arises from heat transport by emission and absorption of photons. The latter, therefore, can be indirectly measured by analyzing the visible and infrared (VIS-IR) regions of a material's optical absorption spectrum. Thermal conductivity in the mantle is controlled by temperature, pressure, the electronic structure and concentration of transition metal ions (such as iron), and the water content of the material [1,3]. The radiative component has generally been assumed to be negligible, as most ferromagnesian minerals become opaque in the VIS-IR range at high pressures due to intensification and red-shift of electronic charge-transfer bands [4, 5]. However, more recent studies have suggested that mantle minerals may, in fact, remain relatively transparent at high pressures, thereby allowing for a potentially significant contribution to thermal conductivity from the radiative component [6]. We measured optical absorbance spectra of hydrous wadsleyite and hydrous ringwoodite at simultaneous high-pressure and high-temperature conditions up to 26 GPa and 823 K in order to determine their radiative conductivities and to study the potential influence of hydration in the transition zone on thermal conductivity of the mantle. We report radiative thermal conductivities of 1.5 ± 0.2 Wm-1K-1 for hydrous wadsleyite and 1.2 ± 0.1 Wm-1K-1 for hydrous ringwoodite at transition zone conditions. The analytically derived radiative thermal conductivities of anhydrous wadsleyite and ringwoodite are 2.1 ± 0.2 Wm-1K-1 and 1.6 ± 0.2 Wm-1K-1, respectively. Our results imply that a water content of ~1 wt% H2O lowers the thermal radiative conductivity
Radiative heat transfer in rocket thrust chambers and nozzles
NASA Technical Reports Server (NTRS)
Hammad, K. J.; Naraghi, M. H. N.
1989-01-01
Numerical models based on the discrete exchange factor (DEF) and the zonal methods for radiative analysis of rocket engines containing a radiatively participating medium have been developed. These models implement a new technique for calculating the direct exchange factors to account for possible blockage by the nozzle throat. Given the gas and surface temperature distributions, engine geometry, and radiative properties, the models compute the wall radiative heat fluxes at different axial positions. The results of sample calculations for a typical rocket engine (engine 700 at NASA), which uses RP-1 (a kerosene-type propellant), are presented for a wide range of surface and gas properties. It is found that the heat transfer by radiation can reach up to 50 percent of that due to convection. The maximum radiative heat flux is at the inner side of the engine, where the gas temperature is the highest. While the results of both models are in excellent agreement, the computation time of the DEF method is found to be much smaller.
IPRT polarized radiative transfer model intercomparison project - Phase A
NASA Astrophysics Data System (ADS)
Emde, Claudia; Barlakas, Vasileios; Cornet, Céline; Evans, Frank; Korkin, Sergey; Ota, Yoshifumi; Labonnote, Laurent C.; Lyapustin, Alexei; Macke, Andreas; Mayer, Bernhard; Wendisch, Manfred
2015-10-01
The polarization state of electromagnetic radiation scattered by atmospheric particles such as aerosols, cloud droplets, or ice crystals contains much more information about the optical and microphysical properties than the total intensity alone. For this reason an increasing number of polarimetric observations are performed from space, from the ground and from aircraft. Polarized radiative transfer models are required to interpret and analyse these measurements and to develop retrieval algorithms exploiting polarimetric observations. In the last years a large number of new codes have been developed, mostly for specific applications. Benchmark results are available for specific cases, but not for more sophisticated scenarios including polarized surface reflection and multi-layer atmospheres. The International Polarized Radiative Transfer (IPRT) working group of the International Radiation Commission (IRC) has initiated a model intercomparison project in order to fill this gap. This paper presents the results of the first phase A of the IPRT project which includes ten test cases, from simple setups with only one layer and Rayleigh scattering to rather sophisticated setups with a cloud embedded in a standard atmosphere above an ocean surface. All scenarios in the first phase A of the intercomparison project are for a one-dimensional plane-parallel model geometry. The commonly established benchmark results are available at the IPRT website.
Three-Dimensional Radiative Transfer on a Massively Parallel Computer.
NASA Astrophysics Data System (ADS)
Vath, Horst Michael
1994-01-01
We perform three-dimensional radiative transfer calculations on the MasPar MP-1, which contains 8192 processors and is a single instruction multiple data (SIMD) machine, an example of the new generation of massively parallel computers. To make radiative transfer calculations efficient, we must re-consider the numerical methods and methods of storage of data that have been used with serial machines. We developed a numerical code which efficiently calculates images and spectra of astrophysical systems as seen from different viewing directions and at different wavelengths. We use this code to examine a number of different astrophysical systems. First we image the HI distribution of model galaxies. Then we investigate the galaxy NGC 5055, which displays a radial asymmetry in its optical appearance. This can be explained by the presence of dust in the outer HI disk far beyond the optical disk. As the formation of dust is connected to the presence of stars, the existence of dust in outer regions of this galaxy could have consequences for star formation at a time when this galaxy was just forming. Next we use the code for polarized radiative transfer. We first discuss the numerical computation of the required cyclotron opacities and use them to calculate spectra of AM Her systems, binaries containing accreting magnetic white dwarfs. Then we obtain spectra of an extended polar cap. Previous calculations did not consider the three -dimensional extension of the shock. We find that this results in a significant underestimate of the radiation emitted in the shock. Next we calculate the spectrum of the intermediate polar RE 0751+14. For this system we obtain a magnetic field of ~10 MG, which has consequences for the evolution of intermediate polars. Finally we perform 3D radiative transfer in NLTE in the two-level atom approximation. To solve the transfer equation in this case, we adapt the short characteristic method and examine different acceleration methods to obtain the
Radiative transfer in layered media: Application to the radar sounding of Martian polar ices. II
NASA Astrophysics Data System (ADS)
Ilyushin, Ya. A.
2007-01-01
The paper addresses the problem of the retrieval of the physical parameters of the Martian polar ices from the ground-penetrating radar (GPR) measurements. The results of the previous paper [Ilyushin, Y.A., Seu, R., Phillips, R.J., 2005. Subsurface radar sounding of the Martian polar cap: radiative transfer approach. Planet. Space Sci. 53(14-15), 1427-1436, doi: 10.1016/j.pss.2005.08.002] are further developed. Several electrical models of the Martian polar caps are considered, and several radiative transfer models, appropriate for layered media, are reviewed and compared to each other. Both 1D case, applicable to the orbital GPR instruments, and 2D and 3D geometry, applicable to the landed radars, are studied. A technique for retrieval of the average loss in the medium, which can be attributed to the impurities content in the ice, is proposed. The retrieval technique is numerically tested with the exact solutions of electromagnetic equations for the waves in the layered media. Numerical tests show that the loss in the medium can be determined from the observed radar signal in the whole range of values, expected to be the case in Martian polar caps.
A linear stability analysis for nonlinear, grey, thermal radiative transfer problems
Wollaber, Allan B.; Larsen, Edward W.
2011-02-20
We present a new linear stability analysis of three time discretizations and Monte Carlo interpretations of the nonlinear, grey thermal radiative transfer (TRT) equations: the widely used 'Implicit Monte Carlo' (IMC) equations, the Carter Forest (CF) equations, and the Ahrens-Larsen or 'Semi-Analog Monte Carlo' (SMC) equations. Using a spatial Fourier analysis of the 1-D Implicit Monte Carlo (IMC) equations that are linearized about an equilibrium solution, we show that the IMC equations are unconditionally stable (undamped perturbations do not exist) if {alpha}, the IMC time-discretization parameter, satisfies 0.5 < {alpha} {<=} 1. This is consistent with conventional wisdom. However, we also show that for sufficiently large time steps, unphysical damped oscillations can exist that correspond to the lowest-frequency Fourier modes. After numerically confirming this result, we develop a method to assess the stability of any time discretization of the 0-D, nonlinear, grey, thermal radiative transfer problem. Subsequent analyses of the CF and SMC methods then demonstrate that the CF method is unconditionally stable and monotonic, but the SMC method is conditionally stable and permits unphysical oscillatory solutions that can prevent it from reaching equilibrium. This stability theory provides new conditions on the time step to guarantee monotonicity of the IMC solution, although they are likely too conservative to be used in practice. Theoretical predictions are tested and confirmed with numerical experiments.
NASA Astrophysics Data System (ADS)
Havemann, Stephan; Thelen, Jean-Claude; Taylor, Jonathan P.; Keil, Andreas
2009-03-01
The Havemann-Taylor Fast Radiative Transfer Code (HT-FRTC) has been developed for the simulation of highly spectrally resolved measurements from satellite based (i.e. Infrared Atmospheric Sounding Interferometer (IASI), Atmospheric Infrared Sounder (AIRS)) and airborne (i.e. Atmospheric Research Interferometer Evaluation System (ARIES)) instruments. The use of principle components enables the calculation of a complete spectrum in less than a second. The principal compoents are derived from a diverse training set of atmospheres and surfaces and contain their spectral characteristics in a highly compressed form. For any given atmosphere/surface, the HT-FRTC calculates the weightings (also called scores) of a few hundred principal components based on selected monochromatic radiative transfer calculations, which is far cheaper than thousands of channel radiance calculations. By intercomparison with line-by-line and other fast models the HT-FRTC has been shown to be accurate. The HT-FRTC has been successfully applied to simultaneous variational retrievals of atmospheric temperature and humidity profiles, surface temperature and surface emissivity over land. This is the subject of another presentation at this conference. The HT-FRTC has now also been extended to include an exact treatment of scattering by aerosols/clouds. The radiative transfer problem is solved using a discrete ordinate method (DISORT). Modelling results at high-spectral resolution for non-clear sky atmospheres obtained with the HT-FRTC are presented.
Direct transfer of solar radiation to high temperature applications
NASA Astrophysics Data System (ADS)
Rahou, Maryam; Andrews, John; Rosengarten, Gary
2013-12-01
This paper reviews the different methods of directly transferring solar radiation from concentrated solar collectors to medium to high temperature thermal absorbers, at temperatures ranging from 100 to 400°. These methods are divided into four main categories associated with the radiation transfer medium: optical fibres, photonic crystal fibres, metal waveguides and light guides. The reviewed methods are novel compared to most rooftop solar concentrators that have a receiver and a thermal storage unit coupled by heat transfer fluids. Bundled optical fibres have the capability of transferring concentrated solar energy across the full wavelength spectrum with the maximum optical efficiency. In this study two different types of optical bundle, including hard polymer cladding silica (HPCS) and polymer clad silica (PCS) fibres are introduced which offer a broad spectrum transmission range from 300 to 1700 nm, low levels of losses through attenuation and the best resistance to heating. These fibres are able to transmit about 94% of the solar radiation over a distance of 10 m. The main parameters that determine the overall efficiency of the system are the concentration ratio, the acceptance angle of the fibres, and the matching of the diameter of the focus spot of the concentrator and the internal diameter of the fibre. In order to maximize the coupling efficiency of the system, higher levels of concentration are required which can be achieved through lenses or other non-imaging concentrators. However, these additional components add to the cost and complexity of the system. To avoid this problem we use tapered bundles of optical fibres that enhance the coupling efficiency by increasing the acceptance angle and consequently the coupling efficiency of the system.
Clouds Radiative Transfer Study at Microwave Region-RTM
NASA Astrophysics Data System (ADS)
Heredia, S. D.; Masuelli, S.; Caranti, G. M.; Jones, L.
2011-12-01
The objective of the recently launched SAC-D/Aquarius satellite mission is to globally and indirectly measure certain geophysical parameters such as: sea surface salinity (Sal), column water vapor (CWV), column liquid water (CLW), rain rate (RR), wind speed (WS), wind direction (WD), ice concentration (SIC) and others. On board the satellite there are several instruments designed for specific purposes like the passive microwave sensor MWR (Fig. 1) whose specifications are shown in Table 1. The aim of the latter is to determine the following parameters: CWV, CLW, RR, WS, WD and SIC. The MWR sensor measures brightness temperatures at two frequencies: 23.8 and 36.5GHz. In the case of 36.5GHz, it measures both polarizations (vertical and horizontal) while for 23.8GHz it only measures the horizontal component. Since this sensor measures brightness temperatures and not geophysical variables, it is necessary to establish a relationship that links both. These relationships are determined by radiative transfer models (RTM). In remote sensing there are two types of models, namely: Forward and Inverse Model. The radiative transfer model in the forward direction obtains brightness temperatures for a given configuration within the pixel (geophysical variables). The most important applications of these models are: * Simulator Development: spectral bands selection to meet the high-level requirements within the expected error. * Intercalibration: in the calculation of corrections due to differences in incidence angles and frequencies between sensors involved in this process. * Inverse Radiative Transfer Models to obtain geophysical variables from brightness temperatures. In this paper, we developed a module that simulates the interaction of radiation with cloud droplets and raindrops. These modules were incorporated into a radiative transfer model from CFRSL (Central Florida Remote Sensing Lab) to calculate the brightness temperatures that would measure a passive microwave sensor
Three-dimensional radiative transfer on a massively parallel computer
NASA Astrophysics Data System (ADS)
Vath, H. M.
1994-04-01
We perform 3D radiative transfer calculations in non-local thermodynamic equilibrium (NLTE) in the simple two-level atom approximation on the Mas-Par MP-1, which contains 8192 processors and is a single instruction multiple data (SIMD) machine, an example of the new generation of massively parallel computers. On such a machine, all processors execute the same command at a given time, but on different data. To make radiative transfer calculations efficient, we must re-consider the numerical methods and storage of data. To solve the transfer equation, we adopt the short characteristic method and examine different acceleration methods to obtain the source function. We use the ALI method and test local and non-local operators. Furthermore, we compare the Ng and the orthomin methods of acceleration. We also investigate the use of multi-grid methods to get fast solutions for the NLTE case. In order to test these numerical methods, we apply them to two problems with and without periodic boundary conditions.
Computing Radiative Transfer in a 3D Medium
NASA Technical Reports Server (NTRS)
Von Allmen, Paul; Lee, Seungwon
2012-01-01
A package of software computes the time-dependent propagation of a narrow laser beam in an arbitrary three- dimensional (3D) medium with absorption and scattering, using the transient-discrete-ordinates method and a direct integration method. Unlike prior software that utilizes a Monte Carlo method, this software enables simulation at very small signal-to-noise ratios. The ability to simulate propagation of a narrow laser beam in a 3D medium is an improvement over other discrete-ordinate software. Unlike other direct-integration software, this software is not limited to simulation of propagation of thermal radiation with broad angular spread in three dimensions or of a laser pulse with narrow angular spread in two dimensions. Uses for this software include (1) computing scattering of a pulsed laser beam on a material having given elastic scattering and absorption profiles, and (2) evaluating concepts for laser-based instruments for sensing oceanic turbulence and related measurements of oceanic mixed-layer depths. With suitable augmentation, this software could be used to compute radiative transfer in ultrasound imaging in biological tissues, radiative transfer in the upper Earth crust for oil exploration, and propagation of laser pulses in telecommunication applications.
Conjugate conductive, convective, and radiative heat transfer in rocket engines
Naraghi, M.H.N.; DeLise, J.C.
1995-12-31
A comprehensive conductive, convective and radiative model for thermal analysis of rocket thrust chambers and nozzles is presented. In this model, the rocket thrust chamber and nozzle are subdivided into a number of stations along the longitudinal direction. At each station a finite element scheme is used to evaluate wall temperature distribution. The hot-gas-side convective heat transport is evaluated by numerically solving the compressible boundary layer equations and the radiative fluxes are evaluated by implementing an exchange factor scheme. The convective heat flux in the cooling channel is modeled based on the existing closed form correlations for rocket cooling channels. The conductive, convective and radiative processes are conjugated through an iterative procedure. The hot-gas-side heat transfer coefficients evaluated based on this model are compared to the experimental results reported in the literature. The computed convective heat transfer coefficients agree very well with experimental data for most of the engine except the throat where a discrepancy of approximately 20% exists. The model is applied to a typical regeneratively cooled rocket engine and the resulting wall temperature and heat flux distribution are presented.
Radiative transfer in cylindrical threads with incident radiation. II. 2D azimuth-dependent case
NASA Astrophysics Data System (ADS)
Gouttebroze, P.
2005-05-01
A method is proposed for the solution of NLTE radiative transfer equations in long cylinders with an external incident radiation that varies with direction. This method is designed principally for the modelling of elongated structures imbedded in the solar corona (loops, prominence threads). The radiative transfer problem under consideration is a 2D one, since the source functions and absorption coefficients vary with both distance to axis and azimuth. The method is based on the general principles of finite-differences and accelerated Λ-iteration. A Fourier series is used for interpolation in azimuth. The method is applied to a line emitted by a two-level atom with complete frequency redistribution. Convergence properties of the method and influence of the inclination angle on the source function are discussed.
Radiative transfer of X-rays in the solar corona
NASA Technical Reports Server (NTRS)
Acton, L. W.
1978-01-01
The problem of resonance scattering of X-ray emission lines in the solar corona is investigated. For the resonance lines of some helium-like ions, significant optical depths are reached over distances small compared with the size of typical coronal features. A general integral equation for the transfer of resonance-line radiation under solar coronal conditions is derived. This expression is in a form useful for modeling the complex three-dimensional temperature and density structure of coronal active regions. The transfer equation is then cast in a form illustrating the terms which give rise to the attenuation or enhancement of the resonance-line intensity. The source function for helium-like oxygen (O VII) under coronal conditions is computed and discussed in terms of the relative importance of scattering.
Comparison of vibrational conductivity and radiative energy transfer methods
NASA Astrophysics Data System (ADS)
Le Bot, A.
2005-05-01
This paper is concerned with the comparison of two methods well suited for the prediction of the wideband response of built-up structures subjected to high-frequency vibrational excitation. The first method is sometimes called the vibrational conductivity method and the second one is rather known as the radiosity method in the field of acoustics, or the radiative energy transfer method. Both are based on quite similar physical assumptions i.e. uncorrelated sources, mean response and high-frequency excitation. Both are based on analogies with some equations encountered in the field of heat transfer. However these models do not lead to similar results. This paper compares the two methods. Some numerical simulations on a pair of plates joined along one edge are provided to illustrate the discussion.
NASA Astrophysics Data System (ADS)
Braghiere, Renato; Quaife, Tristan; Black, Emily
2016-04-01
Incoming shortwave radiation is the primary source of energy driving the majority of the Earth's climate system. The partitioning of shortwave radiation by vegetation into absorbed, reflected, and transmitted terms is important for most of biogeophysical processes, including leaf temperature changes and photosynthesis, and it is currently calculated by most of land surface schemes (LSS) of climate and/or numerical weather prediction models. The most commonly used radiative transfer scheme in LSS is the two-stream approximation, however it does not explicitly account for vegetation architectural effects on shortwave radiation partitioning. Detailed three-dimensional (3D) canopy radiative transfer schemes have been developed, but they are too computationally expensive to address large-scale related studies over long time periods. Using a straightforward one-dimensional (1D) parameterisation proposed by Pinty et al. (2006), we modified a two-stream radiative transfer scheme by including a simple function of Sun zenith angle, so-called "structure factor", which does not require an explicit description and understanding of the complex phenomena arising from the presence of vegetation heterogeneous architecture, and it guarantees accurate simulations of the radiative balance consistently with 3D representations. In order to evaluate the ability of the proposed parameterisation in accurately represent the radiative balance of more complex 3D schemes, a comparison between the modified two-stream approximation with the "structure factor" parameterisation and state-of-art 3D radiative transfer schemes was conducted, following a set of virtual scenarios described in the RAMI4PILPS experiment. These experiments have been evaluating the radiative balance of several models under perfectly controlled conditions in order to eliminate uncertainties arising from an incomplete or erroneous knowledge of the structural, spectral and illumination related canopy characteristics typical
Global sensitivity analysis of the radiative transfer model
NASA Astrophysics Data System (ADS)
Neelam, Maheshwari; Mohanty, Binayak P.
2015-04-01
With the recently launched Soil Moisture Active Passive (SMAP) mission, it is very important to have a complete understanding of the radiative transfer model for better soil moisture retrievals and to direct future research and field campaigns in areas of necessity. Because natural systems show great variability and complexity with respect to soil, land cover, topography, precipitation, there exist large uncertainties and heterogeneities in model input factors. In this paper, we explore the possibility of using global sensitivity analysis (GSA) technique to study the influence of heterogeneity and uncertainties in model inputs on zero order radiative transfer (ZRT) model and to quantify interactions between parameters. GSA technique is based on decomposition of variance and can handle nonlinear and nonmonotonic functions. We direct our analyses toward growing agricultural fields of corn and soybean in two different regions, Iowa, USA (SMEX02) and Winnipeg, Canada (SMAPVEX12). We noticed that, there exists a spatio-temporal variation in parameter interactions under different soil moisture and vegetation conditions. Radiative Transfer Model (RTM) behaves more non-linearly in SMEX02 and linearly in SMAPVEX12, with average parameter interactions of 14% in SMEX02 and 5% in SMAPVEX12. Also, parameter interactions increased with vegetation water content (VWC) and roughness conditions. Interestingly, soil moisture shows an exponentially decreasing sensitivity function whereas parameters such as root mean square height (RMS height) and vegetation water content show increasing sensitivity with 0.05 v/v increase in soil moisture range. Overall, considering the SMAPVEX12 fields to be water rich environment (due to higher observed SM) and SMEX02 fields to be energy rich environment (due to lower SM and wide ranges of TSURF), our results indicate that first order as well as interactions between the parameters change with water and energy rich environments.
NASA Astrophysics Data System (ADS)
Chevallier, L.
2010-11-01
Tests are presented of the 1D Accelerated Lambda Iteration method, which is widely used for solving the radiative transfer equation for a stellar atmosphere. We use our ARTY code as a reference solution and tables for these tests are provided. We model a static idealized stellar atmosphere, which is illuminated on its inner face and where internal sources are distributed with weak or strong gradients. This is an extension of published tests for a slab without incident radiation and gradients. Typical physical conditions for the continuum radiation and spectral lines are used, as well as typical values for the numerical parameters in order to reach a 1% accuracy. It is shown that the method is able to reach such an accuracy for most cases but the spatial discretization has to be refined for strong gradients and spectral lines, beyond the scope of realistic stellar atmospheres models. Discussion is provided on faster methods.
A multilevel method for conductive-radiative heat transfer
Banoczi, J.M.; Kelley, C.T.
1996-12-31
We present a fast multilevel algorithm for the solution of a system of nonlinear integro-differential equations that model steady-state combined radiative-conductive heat transfer. The equations can be formulated as a compact fixed point problem with a fixed point map that requires both a solution of the linear transport equation and the linear heat equation for its evaluation. We use fast transport solvers developed by the second author, to construct an efficient evaluation of the fixed point map and then apply the Atkinson-Brakhage, method, with Newton-GMRES as the coarse mesh solver, to the full nonlinear system.
Peregrinations through topics in light scattering and radiative transfer
NASA Astrophysics Data System (ADS)
Kattawar, George W.
2016-07-01
In this van de Hulst essay, I have taken the liberty to present a journey through some topics in light scattering and radiative transfer which I feel were major contributions to the field but the number of topics I would like to cover is far more numerous than I have the time or the space to present. I also wanted to share with the reader some heartwarming memories I have of my wonderful friend and truly distinguished colleague Hendrik Christoffel van de Hulst (affectionately known to his colleagues as "Henk") whom I consider to be one of the preeminent scientists of his era.
Radiative Transfer and Absorbing Structures in the Transition Region
NASA Astrophysics Data System (ADS)
Plovanic, Jacob; Kankelborg, C. C.
2012-05-01
A fully satisfactory explanation for the anomalous He II 304 Å intensity in the solar transition region has yet to be offered. As an extension of previous work, we use a full radiative transfer code to build a more consistent model of the transition region that allows the He II line to form with low filling factor and low opacity. Our results are constrained by the quiet sun center-to-limb profile of He II 304 Å obtained from the MOSES sounding rocket mission and by AIA full-disk data.
3D Monte Carlo radiation transfer modelling of photodynamic therapy
NASA Astrophysics Data System (ADS)
Campbell, C. Louise; Christison, Craig; Brown, C. Tom A.; Wood, Kenneth; Valentine, Ronan M.; Moseley, Harry
2015-06-01
The effects of ageing and skin type on Photodynamic Therapy (PDT) for different treatment methods have been theoretically investigated. A multilayered Monte Carlo Radiation Transfer model is presented where both daylight activated PDT and conventional PDT are compared. It was found that light penetrates deeper through older skin with a lighter complexion, which translates into a deeper effective treatment depth. The effect of ageing was found to be larger for darker skin types. The investigation further strengthens the usage of daylight as a potential light source for PDT where effective treatment depths of about 2 mm can be achieved.
The diffusion approximation. An application to radiative transfer in clouds
NASA Technical Reports Server (NTRS)
Arduini, R. F.; Barkstrom, B. R.
1976-01-01
It is shown how the radiative transfer equation reduces to the diffusion equation. To keep the mathematics as simple as possible, the approximation is applied to a cylindrical cloud of radius R and height h. The diffusion equation separates in cylindrical coordinates and, in a sample calculation, the solution is evaluated for a range of cloud radii with cloud heights of 0.5 km and 1.0 km. The simplicity of the method and the speed with which solutions are obtained give it potential as a tool with which to study the effects of finite-sized clouds on the albedo of the earth-atmosphere system.
Odyssey: Ray tracing and radiative transfer in Kerr spacetime
NASA Astrophysics Data System (ADS)
Pu, Hung-Yi; Yun, Kiyun; Younsi, Ziri; Yoon, Suk-Jin
2016-01-01
Odyssey is a GPU-based General Relativistic Radiative Transfer (GRRT) code for computing images and/or spectra in Kerr metric describing the spacetime around a rotating black hole. Odyssey is implemented in CUDA C/C++. For flexibility, the namespace structure in C++ is used for different tasks; the two default tasks presented in the source code are the redshift of a Keplerian disk and the image of a Keplerian rotating shell at 340GHz. Odyssey_Edu, an educational software package for visualizing the ray trajectories in the Kerr spacetime that uses Odyssey, is also available.
A field test of a simple stochastic radiative transfer model
Byrne, N.
1995-09-01
The problem of determining the effect of clouds on the radiative energy balance of the globe is of well-recognized importance. One can in principle solve the problem for any given configuration of clouds using numerical techniques. This knowledge is not useful however, because of the amount of input data and computer resources required. Besides, we need only the average of the resulting solution over the grid scale of a general circulation model (GCM). Therefore, we are interested in estimating the average of the solutions of such fine-grained problems using only coarse grained data, a science or art called stochastic radiation transfer. Results of the described field test indicate that the stochastic description is a somewhat better fit to the data than is a fractional cloud cover model, but more data are needed. 1 ref., 3 figs.
Radiative transfer for a three-dimensional raining cloud
NASA Technical Reports Server (NTRS)
Haferman, J. L.; Krajewski, W. F.; Smith, T. F.; Sanchez, A.
1993-01-01
Satellite-sensor-based microwave brightness temperatures for a three-dimensional raining cloud over a reflecting surface are computed by using a radiative transfer model based on the discrete-ordinates solution procedure. The three-dimensional model applied to a plane layer is validated by comparison with results from a one-dimensional model that is available in the literature. Results examining the effects of cloud height, rainfall rate, surface reflectance, rainfall footprint area, and satellite viewing position on one- and three-dimensional brightness temperature calculations are reported. The numerical experiments indicate that, under certain conditions, three-dimensional effects are significant in the analysis of satellite-sensor-based rainfall retrieval algorithms. The results point to the need to consider carefully three-dimensional effects as well as surface reflectance effects when interpreting satellite-measured radiation data.
Analytical properties of the radiance in atmospheric radiative transfer theory
NASA Astrophysics Data System (ADS)
Otto, Sebastian
2014-01-01
It is demonstrated mathematically strictly that state density functions, as the radiance (specific intensity), exist to describe certain state properties of transported photons on microscopic and the state of the radiation field on macroscopic scale, which have independent physical meanings. Analytical properties as boundedness, continuity, differentiability and integrability of these functions to describe the photon transport are discussed. It is shown that the density functions may be derived based on the assumption of photons as real particles of non-zero and finite size, independently of usual electrodynamics, and certain historically postulated functional relationships between them were proved, that is, these functions can be derived mathematically strictly and consistently within the framework of the theory of the phenomenological radiative transfer if one takes the theory seriously by really assuming photons as particles. In this sense these functions may be treated as fundamental physical quantities within the scope of this theory, if one considers the possibility of the existence of photons.
The Chandrasekhar method and its applications to atmospheric radiative transfer
Stamnes, K.
1994-12-31
Problems involving radiation and particle transport in a host medium require solution of the linear (or linearized) Boltzmann equation. A convenient strategy for solving such problems is to apply a multigroup procedure in which the problem is reformulated as a series of one-group problems in such a way that each one-group problem may be cast into a form identical to the monochromatic radiative transfer equation. In essence, Chandrasekhar`s method consists of converting the integro-differential equation for the resulting one-group problem into a system of coupled differential equations for which eigensolutions are sought. The basic method is well described in Chandrasekhar`s classic text in which applications to simple problems were used to demonstrate the potential power of the method before the advent of the modern computer.
NASA Astrophysics Data System (ADS)
Espy, P. J.; Daae, M.; Shprits, Y.
2010-12-01
The correlation between the inner edge of the outer radiation belt phase space density (PSD) and the plasmapause location (Lpp) using reanalysis is investigated. A large data set is applied for the statistical analysis, using data from 1990-1991 from the CRRES satellite, GEO 1989, GPS-ns18 and Akebono. These data are incorporated into reanalysis by means of a Kalman filter with the UCLA 1-D VERB code. The result is a continuous radial and temporal distribution of the PSD from L*=3 to L*=7. The innovation vector of the reconstructed PSD can give us information about regions where local loss or source processes are dominating. We analyze both the PSD and the innovation vector by binning them into slots of Dst and Kp values. This has been done by finding the time for when the Dst (Kp) is within each bin-size of 20 nT (1) from 10 nT to -130 nT (1 to 8). The PSD and innovation vector was then averaged over each of those times. The result shows a good correlation between the location of the inner edge of the outer radiation belt in the PSD and the location of the plasmapause, which is consistent with previous observations. The boundary between the inner edge of the radiation belt and the Lpp becomes sharper, and the radiation belt becomes thinner, during times of high geomagnetic activity. The innovation vector shows that the inner edge of the source region also lines up well with the Lpp, and further showing a battle between losses and sources during active times. This study also illustrates how data assimilation in the radiation belts can be used to understand the underlining processes of acceleration and loss in the inner magnetosphere.
II. The Second Law in Relation to Thermal Radiative Transfer
NASA Astrophysics Data System (ADS)
Jesudason, Christopher G.
2011-12-01
Planck introduced the quantum hypothesis from his Blackbody radiation studies, where he and subsequent workers opined that classical mechanics and electrodynamical theories could not account for the phenomenon. Hence a statistical mechanics with an appropriate Second law entropy was invented and coupled to the First law to account for quantum effects. Here, as an academic exercise we derive the quantum of energy by considering two structures, that of the dipole oscillators on a 2-D surface and the scattering of radiation into the 3-D cavity. Previous derivations are briefly cited and reviewed where none followed this approach. One prediction from this first order Brownian motion development is that a 2-D sheet of oscillators should emit radiation largely with energy density factor T1 of the Kelvin temperature T, rather than that deduced as T4 from detailed balance. Preliminary measurements conducted here seemed to verify the the T1 density. The first order theory also admits a possibility of nonlinear quanta and the consequences are explored briefly. It was noticed in passing during the experimentation that certain bodies suspended in a vacuum exhibited small persistent temperature differentials. A Second law statement is presented for such cases and consequences explored for processes that are not coupled by Newtonian momentum energy transfer mechanisms, such as for the radiation field as deduced by Planck. The different forms of heat transfer due to different laws (e.g. gravity waves and electromagnetic waves) are strictly separable and cannot be confused or forced to an equivalence. We generalize on the Zeroth law, the Kirchoff law and postulate an appropriate entropy form due to these generalizations.
Radiative Transfer, Black Hole Growth, AGN Feedback in Galaxies
NASA Astrophysics Data System (ADS)
Novak, Gregory
2013-01-01
We have performed 3D hydrodynamic simulations of black hole fueling and AGN feedback using a novel method for treating the radial forces on interstellar gas due to absorption of photons by dust grains. The method provides a solution to the radiative transfer equation and hence computes forces on the gas self-consistently by first solving for the radiation field taking into account radiation sources, absorption, and scattering. The algorithm gives the correct behavior in all of the relevant limits (dominated by the central point source; dominated by the distributed isotropic source; optically thin; optically thick to UV/optical; optically thick to IR) and reasonably interpolates between the limits when necessary. The simulations allow us to study gas flows and feedback processes over length scales from ~1 pc to ~100 kpc. We find that the dynamics and final state of simulations are measurably but only moderately affected by radiative forces on dust, even when assumptions about the dust-to-gas ratio are varied from zero to a value appropriate for the Milky Way. In simulations with high gas densities designed to mimic ULIRGs with a star formation rate of several hundred solar masses per year, dust makes a more substantial contribution to the dynamics and outcome of the simulation.
Radiative properties and heat transfer analysis of fibrous insulations
Yeh, H.Y.
1986-01-01
The spectral radiative properties of pink fiberglass insulation were determined from monochromatic directional-hemispherical reflectance or monochromatic specular reflectance data coupled with an analytical model in the wavelength range of 3 to 80 microns. The fiberglass samples were 3.18 cm square with thicknesses varying from 0.10 to 1.28 cm. Three types of spectral data were recorded for the two sets of samples. First, normal transmittance measurements in the wavelength range of 2.5-40 microns were made with a Perkin-Elmer M521 infrared grating spectrophotometer. Second, the Willey 318S Fourier transform spectrometer (FTS) was utilized to measure the directional-hemispherical reflectance in the wavelength range of 3-20 microns. The nonlinear least squares approach coupled with Chandrasekhar's discrete ordinate method for isotropic and anisotropic radiative transport was used in determining the radiative properties. Third, the specular reflectance at an incidence angle of 16 degrees was measured with a Digilab FTS-20 in the wavelength range of 3-80 microns. A method to invert the radiation properties from the experimental data was developed. Typical residential attic fiberglass insulation was chosen for the heat transfer analysis; results are compared to experimental data.
Numerical Radiative Transfer and the Hydrogen Reionization of the Universe
NASA Astrophysics Data System (ADS)
Petkova, M.
2011-03-01
One of the most interesting questions in cosmology is to understand how the Universe evolved from its nearly uniform and simple state briefly after the Big Bang to the complex state we see around us today. In particular, we would like to explain how galaxies have formed, and why they have the properties that we observe in the local Universe. Computer simulations play a highly important role in studying these questions, because they allow one to follow the dynamical equations of gravity and hydrodynamics well into the non-linear regime of the growth of cosmic structures. The current generation of simulation codes for cosmological structure formation calculates the self-gravity of dark matter and cosmic gas, and the fluid dynamics of the cosmic gas, but radiation processes are typically not taken into account, or only at the level of a spatially uniform, externally imposed background field. However, we know that the radiation field has been highly inhomogeneous during certain phases of the growth of structure, and may have in fact provided important feedback effects for galaxy formation. In particular, it is well established that the diffuse gas in the universe was nearly fully neutral after recombination at very high redshift, but today this gas is highly ionized. Sometime during the evolution, a transition to the ionized state must have occurred, a process we refer to as reionization. The UV radiation responsible for this reionization is now permeating the universe and may in part explain why small dwarf galaxies have so low luminosities. It is therefore clear that accurate and self-consistent studies of galaxy formation and of the dynamics of the reionization process should ideally be done with simulation codes that directly include a treatment of radiative transfer, and that account for all relevant source and sink terms of the radiation. We present a novel numerical implementation of radiative transfer in the cosmological smoothed particle hydrodynamics (SPH
Energy Transfer Based Nanocomposite Scintillator for Radiation Detection
NASA Astrophysics Data System (ADS)
Aslam, Soha; Sahi, Sunil; Chen, Wei; Ma, Lun; Kenarangui, Rasool
2014-09-01
Scintillators are the materials that emit light upon irradiation with high energy radiation like X-ray or gamma-ray. Inorganic single crystal and organic (plastic and liquid) are the two most used scintillator types. Both of these scintillator kinds have advantages and disadvantages. Inorganic single crystals are expensive and difficult to grow in desire shape and size. Also, single crystal scintillator such as NaI and CsI are very hygroscopic. On the other hand, organic scintillators have low density which limits their applications in gamma spectroscopy. Due to high quantum yield and size dependent emission, nanoparticles have attracted interested in various field of research. Here, we have studies the nanoparticles for radiation detection. We have synthesized nanoparticles of Cerium fluoride (CeF3), Zinc Oxide (ZnO), Cadmium Telluride (CdTe), Copper complex and Zinc sulfide (ZnS). We have used Fluorescence Resonance Energy Transfer (FRET) principle to enhance the luminescence properties of nanocomposite scintillator. Nanocomposites scintillators are structurally characterized with X-ray diffraction (XRD) and Transmission Electron Microscopy (TEM). Optical properties are studied using Photoluminescence, UV-Visible and X-ray. Enhancements in the luminescence are observed under UV and X-ray excitation. Preliminary studies shows nanocomposite scintillators are promising for radiation detection. Scintillators are the materials that emit light upon irradiation with high energy radiation like X-ray or gamma-ray. Inorganic single crystal and organic (plastic and liquid) are the two most used scintillator types. Both of these scintillator kinds have advantages and disadvantages. Inorganic single crystals are expensive and difficult to grow in desire shape and size. Also, single crystal scintillator such as NaI and CsI are very hygroscopic. On the other hand, organic scintillators have low density which limits their applications in gamma spectroscopy. Due to high quantum
THE RADIATIVE TRANSFER OF SYNCHROTRON RADIATION THROUGH A COMPRESSED RANDOM MAGNETIC FIELD
Cawthorne, T. V.; Hughes, P. A.
2013-07-01
This paper examines the radiative transfer of synchrotron radiation in the presence of a magnetic field configuration resulting from the compression of a highly disordered magnetic field. It is shown that, provided Faraday rotation and circular polarization can be neglected, the radiative transfer equations for synchrotron radiation separate for this configuration, and the intensities and polarization values for sources that are uniform on large scales can be found straightforwardly in the case where opacity is significant. Although the emission and absorption coefficients must, in general, be obtained numerically, the process is much simpler than a full numerical solution to the transfer equations. Some illustrative results are given and an interesting effect, whereby the polarization increases while the magnetic field distribution becomes less strongly confined to the plane of compression, is discussed. The results are of importance for the interpretation of polarization near the edges of lobes in radio galaxies and of bright features in the parsec-scale jets of active galactic nuclei, where such magnetic field configurations are believed to exist.
NASA Astrophysics Data System (ADS)
Guerlet, Sandrine; Spiga, A.; Sylvestre, M.; Fouchet, T.; Millour, E.; Wordsworth, R.; Leconte, J.; Forget, F.
2013-10-01
Recent observations of Saturn’s stratospheric thermal structure and composition revealed new phenomena: an equatorial oscillation in temperature, reminiscent of the Earth's Quasi-Biennal Oscillation ; strong meridional contrasts of hydrocarbons ; a warm “beacon” associated with the powerful 2010 storm. Those signatures cannot be reproduced by 1D photochemical and radiative models and suggest that atmospheric dynamics plays a key role. This motivated us to develop a complete 3D General Circulation Model (GCM) for Saturn, based on the LMDz hydrodynamical core, to explore the circulation, seasonal variability, and wave activity in Saturn's atmosphere. In order to closely reproduce Saturn's radiative forcing, a particular emphasis was put in obtaining fast and accurate radiative transfer calculations. Our radiative model uses correlated-k distributions and spectral discretization tailored for Saturn's atmosphere. We include internal heat flux, ring shadowing and aerosols. We will report on the sensitivity of the model to spectral discretization, spectroscopic databases, and aerosol scenarios (varying particle sizes, opacities and vertical structures). We will also discuss the radiative effect of the ring shadowing on Saturn's atmosphere. We will present a comparison of temperature fields obtained with this new radiative equilibrium model to that inferred from Cassini/CIRS observations. In the troposphere, our model reproduces the observed temperature knee caused by heating at the top of the tropospheric aerosol layer. In the lower stratosphere (20mbar
radiative heating/cooling by trace
Suomi NPP VIIRS Striping Analysis using Radiative Transfer Model Calculations
NASA Astrophysics Data System (ADS)
Wang, Z.; Cao, C.
2015-12-01
Modern satellite radiometers such as VIIRS have many detectors with slightly different relative spectral response (RSR). These differences can introduce artifacts such as striping in the imagery. In recent studies we have analyzed the striping pattern related to the detector level RSR difference in VIIRS Thermal Emissive Bands (TEB) M15 and M16, which includes line-by-line radiative transfer model (LBLRTM) detector level response study and onboard detector stability evaluation using the solar diffuser. Now we extend these analysis to the Reflective Solar Bands (RSB) using MODTRAN atmospheric radiative transfer model (RTM) for detector level radiance simulation. Previous studies analyzed the striping pattern in the images of VIIRS ocean color and reflectance in RSB, further studies about the root cause for striping are still needed. In this study, we will use the MODTRAN model at spectral resolution of 1 cm^-1 under different atmospheric conditions for VIIRS RSB, for example band M1 centered at 410nm which is used for Ocean Color product retrieval. The impact of detector level RSR difference, atmospheric dependency, and solar geometry on the striping in VIIRS SDR imagery will be investigated. The cumulative histogram method used successfully for the TEB striping analysis will be used to quantify the striping. These analysis help S-NPP and J1 to better understand the root cause for VIIRS image artifacts and reduce the uncertainties in geophysical retrievals to meet the user needs.
Radiative transfer code SHARM for atmospheric and terrestrial applications
NASA Astrophysics Data System (ADS)
Lyapustin, A. I.
2005-12-01
An overview of the publicly available radiative transfer Spherical Harmonics code (SHARM) is presented. SHARM is a rigorous code, as accurate as the Discrete Ordinate Radiative Transfer (DISORT) code, yet faster. It performs simultaneous calculations for different solar zenith angles, view zenith angles, and view azimuths and allows the user to make multiwavelength calculations in one run. The Δ-M method is implemented for calculations with highly anisotropic phase functions. Rayleigh scattering is automatically included as a function of wavelength, surface elevation, and the selected vertical profile of one of the standard atmospheric models. The current version of the SHARM code does not explicitly include atmospheric gaseous absorption, which should be provided by the user. The SHARM code has several built-in models of the bidirectional reflectance of land and wind-ruffled water surfaces that are most widely used in research and satellite data processing. A modification of the SHARM code with the built-in Mie algorithm designed for calculations with spherical aerosols is also described.
Lattice Boltzmann model for a steady radiative transfer equation.
Yi, Hong-Liang; Yao, Feng-Ju; Tan, He-Ping
2016-08-01
A complete lattice Boltzmann model (LBM) is proposed for the steady radiative transfer equation (RTE). The RTE can be regarded as a pure convection equation with a source term. To derive the expressions for the equilibrium distribution function and the relaxation time, an artificial isotropic diffusion term is introduced to form a convection-diffusion equation. When the dimensionless relaxation time has a value of 0.5, the lattice Boltzmann equation (LBE) is exactly applicable to the original steady RTE. We also perform a multiscale analysis based on the Chapman-Enskog expansion to recover the macroscopic RTE from the mesoscopic LBE. The D2Q9 model is used to solve the LBE, and the numerical results obtained by the LBM are comparable to the results obtained by other methods or analytical solutions, which demonstrates that the proposed model is highly accurate and stable in simulating multidimensional radiative transfer. In addition, we find that the convergence rate of the LBM depends on the transport properties of RTE: for diffusion-dominated RTE with a large optical thickness, the LBM shows a second-order convergence rate in space, while for convection-dominated RTE with a small optical thickness, a lower convergence rate is observed. PMID:27627417
Infrared radiative transfer through a regular array of cuboidal clouds
NASA Technical Reports Server (NTRS)
HARSHVARDHAN; Weinman, J. A.
1981-01-01
Infrared radiative transfer through a regular array of cuboidal clouds is studied and the interaction of the sides of the clouds with each other and the ground is considered. The theory is developed for black clouds and is extended to scattering clouds using a variable azimuth two-stream approximation. It is shown that geometrical considerations often dominate over the microphysical aspects of radiative transfer through the clouds. For example, the difference in simulated 10 micron brightness temperature between black isothermal cubic clouds and cubic clouds of optical depth 10, is less than 2 deg for zenith angles less than 50 deg for all cloud fractions when viewed parallel to the array. The results show that serious errors are made in flux and cooling rate computations if broken clouds are modeled as planiform. Radiances computed by the usual practice of area-weighting cloudy and clear sky radiances are in error by 2 to 8 K in brightness temperature for cubic clouds over a wide range of cloud fractions and zenith angles. It is also shown that the lapse rate does not markedly affect the exiting radiances for cuboidal clouds of unit aspect ratio and optical depth 10.
Rabacus: A Python package for analytic cosmological radiative transfer calculations
NASA Astrophysics Data System (ADS)
Altay, G.; Wise, J. H.
2015-04-01
We describe RABACUS, a Python package for calculating the transfer of hydrogen ionizing radiation in simplified geometries relevant to astronomy and cosmology. We present example solutions for three specific cases: (1) a semi-infinite slab gas distribution in a homogeneous isotropic background, (2) a spherically symmetric gas distribution with a point source at the center, and (3) a spherically symmetric gas distribution in a homogeneous isotropic background. All problems can accommodate arbitrary spectra and density profiles as input. The solutions include a treatment of both hydrogen and helium, a self-consistent calculation of equilibrium temperatures, and the transfer of recombination radiation. The core routines are written in Fortran 90 and then wrapped in Python leading to execution speeds thousands of times faster than equivalent routines written in pure Python. In addition, all variables have associated units for ease of analysis. The software is part of the Python Package Index and the source code is available on Bitbucket at https://bitbucket.org/galtay/rabacus. In addition, installation instructions and a detailed users guide are available at http://pythonhosted.org//rabacus.
NASA Technical Reports Server (NTRS)
Weston, K. C.; Reynolds, A. C., Jr.; Alikhan, A.; Drago, D. W.
1974-01-01
Numerical solutions for radiative transport in a class of anisotropically scattering materials are presented. Conditions for convergence and divergence of the iterative method are given and supported by computed results. The relation of two flux theories to the equation of radiative transfer for isotropic scattering is discussed. The adequacy of the two flux approach for the reflectance, radiative flux and radiative flux divergence of highly scattering media is evaluated with respect to solutions of the radiative transfer equation.
Preliminary design for Arctic atmospheric radiative transfer experiments
NASA Technical Reports Server (NTRS)
Zak, B. D.; Church, H. W.; Stamnes, K.; Shaw, G.; Filyushkin, V.; Jin, Z.; Ellingson, R. G.; Tsay, S. C.
1995-01-01
If current plans are realized, within the next few years, an extraordinary set of coordinated research efforts focusing on energy flows in the Arctic will be implemented. All are motivated by the prospect of global climate change. SHEBA (Surface Energy Budget of the Arctic Ocean), led by the National Science Foundation (NSF) and the Office of Naval Research (ONR), involves instrumenting an ice camp in the perennial Arctic ice pack, and taking data for 12-18 months. The ARM (Atmospheric Radiation Measurement) North Slope of Alaska and Adjacent Arctic Ocean (NSA/AAO) Cloud and Radiation Testbed (CART) focuses on atmospheric radiative transport, especially in the presence of clouds. The NSA/AAO CART involves instrumenting a sizeable area on the North Slope of Alaska and adjacent waters in the vicinity of Barrow, and acquiring data over a period of about 10 years. FIRE (First ISCCP (International Satellite Cloud Climatology Program) Regional Experiment) Phase 3 is a program led by the National Aeronautics and Space Administration (NASA) which focuses on Arctic clouds, and which is coordinated with SHEBA and ARM. FIRE has historically emphasized data from airborne and satellite platforms. All three program anticipate initiating Arctic data acquisition during spring, 1997. In light of his historic opportunity, the authors discuss a strawman atmospheric radiative transfer experimental plan that identifies which features of the radiative transport models they think should be tested, what experimental data are required for each type of test, the platforms and instrumentation necessary to acquire those data, and in general terms, how the experiments could be conducted. Aspects of the plan are applicable to all three programs.
Preliminary design for Arctic atmospheric radiative transfer experiments
Zak, B.D.; Church, H.W.; Stamnes, K.; Shaw, G.; Filyushkin, V.; Jin, Z.; Ellingson, R.G.; Tsay, S.C.
1995-04-01
If current plans are realized, within the next few years, an extraordinary set of coordinated research efforts focusing on energy flows in the Arctic will be implemented. All are motivated by the prospect of global climate change. SHEBA (Surface Energy Budget of the Arctic Ocean), led by the National Science Foundation (NSF) and the Office of Naval Research (ONR), involves instrumenting an ice camp in the perennial Arctic ice pack, and taking data for 12--18 months. The ARM (Atmospheric Radiation Measurement) North Slope of Alaska and Adjacent Arctic Ocean (NSA/AAO) Cloud and Radiation Testbed (CART) focuses on atmospheric radiative transport, especially in the presence of clouds. The NSA/AAO CART involves instrumenting a sizeable area on the North Slope of Alaska and adjacent waters in the vicinity of Barrow, and acquiring data over a period of about 10 years. FIRE (First ISCCP [International Satellite Cloud Climatology Program] Regional Experiment) Phase 3 is a program led by the National Aeronautics and Space Administration (NASA) which focuses on Arctic clouds, and which is coordinated with SHEBA and ARM. FIRE has historically emphasized data from airborne and satellite platforms. All three program anticipate initiating Arctic data acquisition during spring, 1997. In light of his historic opportunity, the authors discuss a strawman atmospheric radiative transfer experimental plan that identifies which features of the radiative transport models they think should be tested, what experimental data are required for each type of test, the platforms and instrumentation necessary to acquire those data, and in general terms, how the experiments could be conducted. Aspects of the plan are applicable to all three programs.
Chetia, Tridip Ranjan; Ansari, Mohammad Shaad; Qureshi, Mohammad
2016-02-21
Mesoporous and hollow zinc oxide (ZnO) hierarchical superstructures assembled with compact 1D building blocks that provide an efficient and faster transport pathway for photo-generated charge carriers have been synthesized using a biomass derived polysaccharide "alginic acid". To understand the interactions between the organic bio-template and inorganic growth units of ZnO in aqueous medium, the effects of additives such as the alginate ion (ALGI) and ammonium hydroxide (NH4OH), along with the controlled reaction conditions, are investigated using Field Emission Scanning Electron Microscopy (FESEM) and powder X-ray diffraction. Dynamic and steady-state photoluminescence measurements are carried out to understand the charge transfer processes in the compact 1D superstructures. Experimental analyses reveal that the alginate ions, under hydrothermal reaction conditions, act as a structure directing agent and assemble 1D ZnO nanorods (NRs) hierarchically while NH4OH assists the formation of ZnO growth units. A plausible mechanism for ZnO cage formation is proposed based on the experimental observations. Morphology dependent photovoltaic properties of ZnO heterostructures, i.e., for ZnO cages, ZnO NRs and ZnO PNPs, have been studied along with electrochemical impedance spectroscopy (EIS). Enhancement of ∼ 60% and ∼ 35% in power conversion efficiency (PCE) is observed in ZnO cage based devices as compared to ZnO NR- and ZnO PNP-based devices, respectively. PMID:26818181
Detectivity of gas leakage based on electromagnetic radiation transfer
NASA Astrophysics Data System (ADS)
Long, Yunting; Wang, Lingxue; Li, Jiakun; Zhang, Changxing; Zhang, Bei
2011-05-01
Standoff detection of gas leakage is a fundamental need in petrochemical and power industries. The passive gas imaging system using thermal imager has been proven to be efficient to visualize leaking gas which is not visible to the naked eye. The detection probability of gas leakage is the basis for designing a gas imaging system. Supposing the performance parameters of the thermal imager are known, the detectivity based on electromagnetic radiation transfer model to image gas leakage is analyzed. This model takes into consideration a physical analysis of the gas plume spread in the atmosphere-the interaction processes between the gas and its surrounding environment, the temperature of the gas and the background, the background surface emissivity, and also gas concentration, etc. Under a certain environmental conditions, through calculating the radiation reaching to the detector from the camera's optical field of view, we obtain an entity "Gas Equivalent Blackbody Temperature Difference (GEBTD)" which is the radiation difference between the on-plume and off-plume regions. Comparing the GEBTD with the Noise Equivalent Temperature Difference (NETD) of the thermal imager, we can know whether the system can image the gas leakage. At last, an example of detecting CO2 gas by JADE MWIR thermal imager with a narrow band-pass filter is presented.
Wu, Xinyuan; Huang, Tao; Lekich, Travis T; Sommer, Roger D; Weare, Walter W
2015-06-01
Heterobimetallic complexes composed only of first-row transition metals [(TMTAA)V(IV)═O→M(II)Py5Me2](OTf)2 (TMTAA = 7,16-dihydro-6,8,15,17-tetramethyldibenzo[b,i][1,4,8,11]tetraazacyclotetradecine; Py5Me2 = 2,6-bis(1,1-bis(2-pyridyl)ethyl)pyridine; M = Mn(II), Fe(II), Co(II), Ni(II), Cu(II); OTf = trifluoromethanesulfonate) have been synthesized through a dative interaction between a terminal oxido and M(II) metal centers. This is the first series of V(IV)═O→M(II) heterobimetallic complexes containing an unsupported oxido bridge. Among these five complexes, only V(IV)═O→Fe(II) (3b) has a clear new absorption band upon formation of the dinuclear species (502 nm, ε = 1700 M(-1) cm(-1)). This feature is assigned to a metal-to-metal charge transfer (MMCT) transition from V(IV) to Fe(II), which forms a V(V)-O-Fe(I) excited state. This assignment is supported by electrochemical data, electronic absorption profiles, and resonance Raman spectroscopy and represents the first report of visible-light induced MMCT in a heterobimetallic oxido-bridged molecule where the electron originates on a d(1) metal center. PMID:25970134
NASA Technical Reports Server (NTRS)
Verstraete, Michel M.
1987-01-01
Understanding the details of the interaction between the radiation field and plant structures is important climatically because of the influence of vegetation on the surface water and energy balance, but also biologically, since solar radiation provides the energy necessary for photosynthesis. The problem is complex because of the extreme variety of vegetation forms in space and time, as well as within and across plant species. This one-dimensional vertical multilayer model describes the transfer of direct solar radiation through a leaf canopy, accounting explicitly for the vertical inhomogeneities of a plant stand and leaf orientation, as well as heliotropic plant behavior. This model reproduces observational results on homogeneous canopies, but it is also well adapted to describe vertically inhomogeneous canopies. Some of the implications of leaf orientation and plant structure as far as light collection is concerned are briefly reviewed.
Radiative transfer theory applied to ocean bottom modeling.
Quijano, Jorge E; Zurk, Lisa M
2009-10-01
Research on the propagation of acoustic waves in the ocean bottom sediment is of interest for active sonar applications such as target detection and remote sensing. The interaction of acoustic energy with the sea floor sublayers is usually modeled with techniques based on the full solution of the wave equation, which sometimes leads to mathematically intractable problems. An alternative way to model wave propagation in layered media containing random scatterers is the radiative transfer (RT) formulation, which is a well established technique in the electromagnetics community and is based on the principle of conservation of energy. In this paper, the RT equation is used to model the backscattering of acoustic energy from a layered elastic bottom sediment containing distributions of independent scatterers due to a constant single frequency excitation in the water column. It is shown that the RT formulation provides insight into the physical phenomena of scattering and conversion of energy between waves of different polarizations. PMID:19813787