Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “discrete ordinates”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

Coupled atmosphere/canopy model for remote sensing of plant reflectance features

Solar radiative transfer through a coupled system of atmosphere and plant canopy is modeled as a multiple-scattering problem through a layered medium of random scatterers. The radiative transfer equation is solved by the discrete-ordinates finite-element method. Analytic expressions are derived that allow the calculation of scattering and absorption cross sections for any plant canopy layer form measurable biophysical parameters such as the leaf area index, leaf angle distribution, and individual leaf reflectance and transmittance data. An expression for a canopy scattering phase function is also given. Computational results are in good agreement with spectral reflectance measurements directly above a soybean canopy, and the concept of greenness- and brightness-transforms of Landsat MSS data is reconfirmed with the computed results. A sensitivity analysis with the coupled atmosphere/canopy model quantifies how satellite-sensed spectral radiances are affected by increased atmospheric aerosols, by varying leaf area index, by anisotropic leaf scattering, and by non-Lambertian soil boundary conditions. Possible extensions to a 2-D model are also discussed.

Gerstl, S. A.↗

Multidimensional Modeling of Atmospheric Effects and Surface Heterogeneities on Remote Sensing

The overall goal of this project is to establish a modeling capability that allows a quantitative determination of atmospheric effects on remote sensing including the effects of surface heterogeneities. This includes an improved understanding of aerosol and haze effects in connection with structural, angular, and spatial surface heterogeneities. One important objective of the research is the possible identification of intrinsic surface or canopy characteristics that might be invariant to atmospheric perturbations so that they could be used for scene identification. Conversely, an equally important objective is to find a correction algorithm for atmospheric effects in satellite-sensed surface reflectances. The technical approach is centered around a systematic model and code development effort based on existing, highly advanced computer codes that were originally developed for nuclear radiation shielding applications. Computational techniques for the numerical solution of the radiative transfer equation are adapted on the basis of the discrete-ordinates finite-element method which proved highly successful for one and two-dimensional radiative transfer problems with fully resolved angular representation of the radiation field.

Gerstl, S. A. W.↗

Comparative accuracy of the Albedo, transmission and absorption for selected radiative transfer approximations

Illustrations of both the relative and absolute accuracy of eight different radiative transfer approximations as a function of optical thickness, solar zenith angle and single scattering albedo are given. Computational results for the plane albedo, total transmission and fractional absorption were obtained for plane-parallel atmospheres composed of cloud particles. These computations, which were obtained using the doubling method, are compared with comparable results obtained using selected radiative transfer approximations. Comparisons were made between asymptotic theory for thick layers and the following widely used two stream approximations: Coakley-Chylek's models 1 and 2, Meador-Weaver, Eddington, delta-Eddington, PIFM and delta-discrete ordinates.

King, M. D.↗

Introduction to the Theory of Atmospheric Radiative Transfer

The fundamental physical and mathematical principles governing the transmission of radiation through the atmosphere are presented, with emphasis on the scattering of visible and near-IR radiation. The classical two-stream, thin-atmosphere, and Eddington approximations, along with some of their offspring, are developed in detail, along with the discrete ordinates method of Chandrasekhar. The adding and doubling methods are discussed from basic principles, and references for further reading are suggested.

Buglia, J. J.↗

Angular dependent transport of auroral electrons in the upper atmosphere

The transport of auroral electrons through the upper atmosphere is analyzed. The transport equation is solved using a discrete-ordinate method, including elastic and inelastic scattering of electrons (resulting in changes of pitch angle) and degradation in energy as the electrons penetrate into the atmosphere. The transport equation is solved numerically for the electron intensity as a function of altitude, pitch angle, and energy. In situ measurements of the pitch angle and energy distribution of precipitating electrons over an auroral arc provide boundary conditions for the calculation. Model calculations were carried out with various different phase functions for elastic and inelastic collisions to attempt changing the angular scattering, but the observed pitch angle distributions remain unexplained. It is suggested that mechanisms other than collisional scattering influence the angular distribution of auroral electrons at or below 300 km altitude in the low-energy domain.

Lummerzheim, D.↗

Flow of rarefied gases over two-dimensional bodies

A kinetic-theory analysis is made of the flow of rarefied gases over two-dimensional bodies of arbitrary curvature. The Boltzmann equation simplified by a model collision integral is written in an arbitrary orthogonal curvilinear coordinate system, and solved by means of finite-difference approximation with the discrete ordinate method. A numerical code is developed which can be applied to any two-dimensional submerged body of arbitrary curvature for the flow regimes from free-molecular to slip at transonic Mach numbers. Predictions are made for the case of a right circular cylinder.

Jeng, Duen-Ren↗

Rarefied gas flow through two-dimensional nozzles

A kinetic theory analysis is made of the flow of a rarefied gas from one reservoir to another through two-dimensional nozzles with arbitrary curvature. The Boltzmann equation simplified by a model collision integral is solved by means of finite-difference approximations with the discrete ordinate method. The physical space is transformed by a general grid generation technique and the velocity space is transformed to a polar coordinate system. A numerical code is developed which can be applied to any two-dimensional passage of complicated geometry for the flow regimes from free-molecular to slip. Numerical values of flow quantities can be calculated for the entire physical space including both inside the nozzle and in the outside plume. Predictions are made for the case of parallel slots and compared with existing literature data. Also, results for the cases of convergent or divergent slots and two-dimensional nozzles with arbitrary curvature at arbitrary knudsen number are presented.

De Witt, Kenneth J.↗

Formation of cyclotron lines in gamma-ray burst spectra

A transmission model of gamma-ray burst sources is studied using the relativistic QED magnetic-resonant opacities including multiple photon scattering, incorporated into a discrete-ordinate radiative-transport scheme. The physics of the cyclotron line-producing region is discussed in general, and the expected line profiles, relative harmonic strengths, and polarizations are indicated under various conditions. The calculated spectra for these models show good agreement with the spectra reported from Ginga for GB 880205 and GB 870303.

Alexander, S. G.↗

A new spherical model for computing the radiation field available for photolysis and heating at twilight

Accurate computation of atmospheric photodissociation and heating rates is needed in photochemical models. These quantities are proportional to the mean intensity of the solar radiation penetrating to various levels in the atmosphere. For large solar zenith angles a solution of the radiative transfer equation valid for a spherical atmosphere is required in order to obtain accurate values of the mean intensity. Such a solution based on a perturbation technique combined with the discrete ordinate method is presented. Mean intensity calculations are carried out for various solar zenith angles. These results are compared with calculations from a plane parallel radiative transfer model in order to assess the importance of using correct geometry around sunrise and sunset. This comparison shows, in agreement with previous investigations, that for solar zenith angles less than 90 deg adequate solutions are obtained for plane parallel geometry as long as spherical geometry is used to compute the direct beam attenuation; but for solar zenith angles greater than 90 deg this pseudospherical plane parallel approximation overstimates the mean intensity.

Dahlback, Arne↗

Numerical solution of a three-dimensional cubic cavity flow by using the Boltzmann equation

A three-dimensional cubic cavity flow has been analyzed for diatomic gases by using the Boltzmann equation with the Bhatnagar-Gross-Krook (B-G-K) model. The method of discrete ordinate was applied, and the diffuse reflection boundary condition was assumed. The results, which show a consistent trend toward the Navier-Stokes solution as the Knudson number is reduced, give us confidence to apply the method to a three-dimensional geometry for practical predictions of rarefied-flow characteristics. The CPU time and the main memory required for a three-dimensional geometry using this method seem reasonable.

Hwang, Danny P.↗

How well is Martian ozone inferred with reflectance spectroscopy?

The Mariner 9 UV spectrometer scanned from 2100 to 3500 Angstroms in one of its two spectral channels every 3 seconds with a spectral resolution of 15 Angstroms and an effective field-of-view of approximately 300 km(exp 2). The only gaseous absorption in the 2000 to 3000 Angstrom region was assumed to come from the Hartley band system of ozone, and therefore the amount of ozone was inferred by fitting this absorption feature with laboratory data of ozone absorption. Mars O3 as inferred from these spectra shows strong seasonal and latitudinal variation, with column abundances ranging from 0.2 microns at equatorial latitudes to 60 microns over the northern winter polar latitudes. The detectability limit of the spectrometer was approximately 3 microns. I use a radiative transfer model based on the discrete ordinate method to calculate synthetic radiance spectra. When typical amounts of dust and cloud are present, significant underestimation of O3 occurs. A factor of 3 times as much O3 is needed to generate the same spectrum for cloudy, dusty atmospheres as for a clear atmosphere. If the scattering properties of Martian clouds and dust were well known, then their appearance would not be a problem, as a model would be capable of retrieving the O3 abundance. However, these properties are not well known, which raises doubts about the effectiveness of the current UV spectroscopy technique used to measure O3. Spatial and temporal variability in temperature and water vapor have been claimed to account for the scatter of the data points. However, water vapor is a small source of odd hydrogen in the winter polar atmosphere, and may not account for most of the variability. Masking by clouds and dust may also account for some of the observed O3 variability, because the nature and opacity of the clouds and dust in the polar hood change dramatically in latitude and even on a day-to-day basis. As the maximum O3 abundance resides near the surface, spacecraft must be able to observe through the entire cloud and dust abundance in order to actually see the total O3 column abundance. If reflectance spectroscopy is used, as on Mariner 9, then the cloud and the airborne dust must be traversed twice; first by the incoming solar flux down to the surface, and then once again upon reflection from the surface out to the spacecraft.

Lindner, Bernhard Lee↗

The effect of polar clouds and dust on the radiative budget of the Martian polar cap

How ozone, clouds, and airborne dust affect CO2 ice formation and sublimation is examined to determine if they help explain the asymmetry in the residual polar caps. The radiative flux striking the surface is calculated for 57 solar and 10 IR wavelength intervals from 0 to 100 microns including the scatering, absorption, and emission by O3, CO2, clouds, and dust. UV absorption by H2O, O2, HO2, and H2O2 is also included, as is Rayleigh scattering. The discrete ordinate method of Stamnes et al. (1988) treats the emission, scattering, and absorption of monochromatic radiation.

Lindner, Bernhard L.↗

Radiative transfer solutions for cylindrical coordinates with emitting, absorbing, and anisotropic scattering medium

The radiative transfer equation is solved by the S-N discrete ordinates method in one-dimensional radial and two-dimensional r-z coordinates systems. The walls of the enclosures are diffuse, and the participating medium absorbs, emits and anisotropically scatters the radiative energy. Diffuse wall incidence, isothermal medium emission and collimated incidence problems are considered. Effects of the scattering phase functions on average incident radiation and net radiative heat fluxes are studied. In addition, the effects of scattering albedo, optical thickness and the wall emissivity are briefly discussed.

Jendoubi, Slah↗

A solar radiation model for use in climate studies

A solar radiation routine is developed for use in climate studies that includes absorption and scattering due to ozone, water vapor, oxygen, carbon dioxide, clouds, and aerosols. Rayleigh scattering is also included. Broadband parameterization is used to compute the absorption by water vapor in a clear atmosphere, and the k-distribution method is applied to compute fluxes in a scattering atmosphere. The reflectivity and transmissivity of a scattering layer are computed analytically using the delta-four-stream discrete-ordinate approximation. The two-stream adding method is then applied to compute fluxes for a composite of clear and scattering layers. Compared to the results of high spectral resolution and detailed multiple-scattering calculations, fluxes and heating rate are accurately computed to within a few percent. The high accuracy of the flux and heating-rate calculations is achieved with a reasonable amount of computing time. With the UV and visible region grouped into four bands, this solar radiation routine is useful not only for climate studies but also for studies on photolysis in the upper atmosphere and photosynthesis in the biosphere.

Chou, Ming-Dah↗

Estimating snow grain size using AVIRIS data

Estimates of snow grain size for the near-surface snow layer were calculated for the Tioga Pass region and Mammoth Mountain in the Sierra Nevada, California, using an inversion technique and data collected by the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS). The Tioga Pass and Mammoth Mountain single-band AVIRIS radiance images were atmospherically corrected to obtain surface reflectance. A discrete-ordinate model was used to calculate directional reflectance as a function of snowpack grain size for a wide range of snow grain radii. The resulting radius vs. reflectance curves were each fit using a nonlinear least squares technique which provided a means of transforming surface reflectance in each AVIRIS image to optically equivalent grain size on a per-pixel basis. The model results and grain size estimates derived from the AVIRIS data show that, for solar incidence angles between 0 and 30, the technique provides good estimates of grain size. This work provides the first quantitative estimates for grain size using data acquired from an airborne remote sensing instrument and is an important step in improving our ability to retrieve snow physical properties independent of field measurements.

Nolin, Anne W.↗

Radiative transfer for a three-dimensional raining cloud

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.

Haferman, J. L.↗

Radiative transfer in three-dimensional atmosphere-vegetation media

The problem of radiation transfer in three-dimensional atmosphere-vegetation media is formulated in this paper. The atmosphere is treated as a horizontally homogeneous medium but with spatially varying anisotropic lower boundary. The boundary condition is defined in terms of vegetation canopy bidirectional reflectance factors to be evaluated from a numerical solution of the three-dimensional canopy radiative transfer equation. Numerical results obtained with the discrete ordinates method are in good agreement with the Monte Carlo results of Pearce and field soybean reflectance data of Ranson et al.

Myneni, Ranga B.↗

Combustion chamber analysis code

A three-dimensional, time dependent, Favre averaged, finite volume Navier-Stokes code has been developed to model compressible and incompressible flows (with and without chemical reactions) in liquid rocket engines. The code has a non-staggered formulation with generalized body-fitted-coordinates (BFC) capability. Higher order differencing methodologies such as MUSCL and Osher-Chakravarthy schemes are available. Turbulent flows can be modeled using any of the five turbulent models present in the code. A two-phase, two-liquid, Lagrangian spray model has been incorporated into the code. Chemical equilibrium and finite rate reaction models are available to model chemically reacting flows. The discrete ordinate method is used to model effects of thermal radiation. The code has been validated extensively against benchmark experimental data and has been applied to model flows in several propulsion system components of the SSME and the STME.

Przekwas, A. J.↗