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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Radiation transfer in gamma-ray bursts

We have calculated gamma-ray radiative transport in regions of high-energy density, such as gamma-ray burst source regions, using a discrete ordinate, discrete energy group method. The calculations include two-photon pair production and annihilation, as well as three-photon annihilation. The radiation fields itself acts as an absorbing medium, and the optical depth depends on its intensity, so the problem is intrinsically nonlinear. Spherical divergence produces effective collimation of the flux. At high optical depth the high energy (E is greater than 1 MeV) portion of the emergent spectrum assumes a nearly universal form. An approximate limit is derived for the high-energy flux from a gamma-ray burst source region of given size, and the implications of this limit for the distance to the 1979 March 5 event are briefly discussed. We discuss more generally the problem of very luminous bursts, and implications of Galactic halo distances for flare models.

Carrigan, B. J.↗
NASA NTRS1972

Discrete ordinates-Monte Carlo coupling: A comparison of techniques in NERVA radiation analysis

In the radiation analysis of the NERVA nuclear rocket system, two-dimensional discrete ordinates calculations are sufficient to provide detail in the pressure vessel and reactor assembly. Other parts of the system, however, require three-dimensional Monte Carlo analyses. To use these two methods in a single analysis, a means of coupling was developed whereby the results of a discrete ordinates calculation can be used to produce source data for a Monte Carlo calculation. Several techniques for producing source detail were investigated. Results of calculations on the NERVA system are compared and limitations and advantages of the coupling techniques discussed.

Lindstrom, D. G.↗
NASA NTRS2004

A Numerical Investigation of the Extinction of Low Strain Rate Diffusion Flames by an Agent in Microgravity

Our goal has been to investigate the influence of both dilution and radiation on the extinction process of nonpremixed flames at low strain rates. Simulations have been performed by using a counterflow code and three radiation models have been included in it, namely, the optically thin, the narrowband, and discrete ordinate models. The counterflow flame code OPPDIFF was modified to account for heat transfer losses by radiation from the hot gases. The discrete ordinate method (DOM) approximation was first suggested by Chandrasekhar for solving problems in interstellar atmospheres. Carlson and Lathrop developed the method for solving multi-dimensional problem in neutron transport. Only recently has the method received attention in the field of heat transfer. Due to the applicability of the discrete ordinate method for thermal radiation problems involving flames, the narrowband code RADCAL was modified to calculate the radiative properties of the gases. A non-premixed counterflow flame was simulated with the discrete ordinate method for radiative emissions. In comparison with two other models, it was found that the heat losses were comparable with the optically thin and simple narrowband model. The optically thin model had the highest heat losses followed by the DOM model and the narrow-band model.

Puri, Ishwar K.↗
DOE OSTI2025

Single Grid Error Estimation for Neutron Transport Solvers

The method of nearby problems (MNP) is a solution verification technique that does not require the use of multiple spatial grids. To estimate spatial discretization error without requiring a high-fidelity spatial grid, an analytical curve fit is interpolated from the numerical solution. The residual between the curve fit solution and numerical solution is calculated and added as an additional source term to the governing equation. The nearby solution is estimated using the updated source term and boundary conditions to remain consistent with the curve fit interpolation. The nearby solution can be compared to the curve fit solution as a discretization error estimation while using a single spatial grid. Without the use of higher fidelity spatial grids, the MNP is able to approximate the spatial discretization error, a facet of solution verification. The application of the method of nearby problems is presented for one- and two-dimensional neutron transport problems for both fixed source and criticality problems on the spatial variable. The fixed source results demonstrate the effectiveness of nearby problems for spatial error identification using the discrete ordinates method. Criticality results are shown to identify area of high spatial error for the C5G7 problem as well as for the discrete ordinates solver. A novel approach of combining the capabilities of Monte Carlo with the discrete ordinates nearby problems is presented for one- and two-dimensional fixed source problems. In conclusion, the MNP demonstrates its effectiveness at identifying spatial error on a single structured grid with a wide variety of neutron transport problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2016

Fast and Accurate Hybrid Stream PCRTMSOLAR Radiative Transfer Model for Reflected Solar Spectrum Simulation in the Cloudy Atmosphere

A hybrid stream PCRTM-SOLAR model has been proposed for fast and accurate radiative transfer simulation. It calculates the reflected solar (RS) radiances with a fast coarse way and then, with the help of a pre-saved matrix, transforms the results to obtain the desired high accurate RS spectrum. The methodology has been demonstrated with the hybrid stream discrete ordinate (HSDO) radiative transfer (RT) model. The HSDO method calculates the monochromatic radiances using a 4-stream discrete ordinate method, where only a small number of monochromatic radiances are simulated with both 4-stream and a larger N-stream (N ≥ 16) discrete ordinate RT algorithm. The accuracy of the obtained channel radiance is comparable to the result from N-stream moderate resolution atmospheric transmission version 5 (MODTRAN5). The root-mean-square errors are usually less than 5x10(exp -4) mW/sq cm/sr/cm. The computational speed is three to four-orders of magnitude faster than the medium speed correlated-k option MODTRAN5. This method is very efficient to simulate thousands of RS spectra under multi-layer clouds/aerosols and solar radiation conditions for climate change study and numerical weather prediction applications.

Yang, Qiguang↗
DOE OSTI2026

The Physics Imposed on a Streaming Operator by Spherical Transport Problems

The streaming operator, which generates a displacement of a particle on a straight line at a constant speed in transport theory, is derived algebraically from a spherical coordinate formulation of Newton’s second law. This derivation leads to an operator that has more partial derivatives than a Cartesian coordinate formulation of the operator. The additional partial derivatives, which are with respect to the normalized velocity variables of a particle, take into account the intrinsic curvature of a ball. Moreover, these partial derivatives mitigate ray effects, which arise when a finite number of normalized velocities (also called directions or discrete ordinates) are used to simulate a continuous S 2 sphere of directions, by rotating the polar axis of the S 2 sphere into the radial direction of the coordinate system. As a consequence of this rotation, the number of actual discrete ordinates is greatly amplified to an enormous number of effective discrete ordinates by a multiplier that is equal to the number of patches that partitions a spherical surface. In addition to the derivation of the streaming operator, we provide in closed form a solution to the system of characteristic equations that is equivalent to the streaming operator. Furthermore, the solution to the system of characteristic equations enables the construction of an integral operator that is the inverse to the streaming operator. Examples in which ray effects are immensely mitigated by spherical coordinates are presented.

integral operator↗
NASA NTRS1978

The phase-integral method for radiative transfer problems with highly-peaked phase functions

Complete solutions to the radiative transfer equation, including both azimuth and depth dependence, are provided by the discrete-ordinate method of Chandrasekhar, but these solutions are often limited because of large computer requirements. This paper presents a 'phase-integral' method which greatly reduces the number of discrete ordinates needed in the solution for highly peaked phase functions. A composite quadrature method is shown to be effective in further reducing the number of discrete ordinates required for highly anisotropic phase functions. Examples are given to indicate convergence requirements and expected accuracy in the complete solution for Henyey-Greenstein and cloud-type phase functions.

Fricke, C. L.↗
NASA NTRSDate not supplied

Data Simulation Using VLIDORT and PCRTM and its Use for Radiometric Calibration of Hyperspectral UV-NIR-SWIR Sensors

Using the vector linearized discrete ordinate radiative transfer model (VLIDORT), a multiple-scattering multi-layer discrete ordinate scattering codes, coupled with a fast principal component (PC) radiative transfer model (PCRTM), we can simulate the radiation fields for hyperspectral UV-NIR-SWIR sensors from 250nm to 2500nm with a very high speed. Since the accurate gas absorptions, molecules and aerosol/cloud particle scatterings, and the surface reflectance have been considered in these simulations, these model data can be used for satellite or airborne sensor design, characterization, or end-to-end system level uncertainty assessment. Comparison of the modelled data with measured results has been a general practice for operational satellite sensor calibrations, such as the Ozone Mapping and Profiler Suite (OMPS) onboard on SNPP and NOAA-20. This presentation will give more details of this simulation system and models. Some examples of the generation of simulations data for CLARREO PathFinder (CPF) and the Surface Biology and Geology (SBG), as well as some examples to use these data for CPF and OMPS radiometric calibration will be presented.

Xiaozhen Xiong↗
NASA NTRS2011

Analysis of the Radiative Transfer Equation with Highly Asymmetric Phase Function

This paper considers a scalar radiative transfer problem with high scattering anisotropy, Two computational methods are presented based on decomposition of the diffuse light field into a regular and anisotropic part. The first algorithm (DOMAS) singles out the anisotropic radiance in the forward scattering peak using the Small-Angle Modification of RTE. The second algorithm (DOM2+) separates the single scattering radiance as an anisotropic part, which largely defines the fine detail of the total radiance in the backscattering directions. In both cases, the anisotropic part is represented analytically. With anisotropy subtraction, the regular part of the signal. which requires a numerical solution, is essentially smoothed as a function of angles. Further, the transport equation is obtained for the regular part that contains an additional source function from the anisotropic part of the signal. This equation is solved with the discrete ordinates method. A conducted numerical analysis of this work showed that algorithm DOMAS has a strong advantage as compared to the standard discrete ordinates method for simulation of the radiance transmission, and DOM2 + is the best of the three for the reflection computations. Both algorithms offer at least a factor of three acceleration of convergence of the azimuthal series for highly anisotropic phase functions.

Korkin, Sergey V.↗
NASA NTRS1972

The design of asymmetric 4 pi shields for space reactors

A one dimensional shield optimization program based on the method of discrete ordinates has been developed and is used to determine material thicknesses used in asymmetric 4 pion shields for space power reactors. The two dimensional discrete ordinates program DOT is used to check the design, and the information generated in the DOT calculation is used as a guide in shaping the shield which may be considered a first step in two dimensional shield optimization.

Engle, W. W., Jr.↗
NASA NTRS2018

Modified Geometric Truncation of the Scattering Phase Function

Phase function of light scattering on large atmospheric particles has very strong peak in forward direction constituting a challenge for accurate numerical calculations of radiance required in remote sensing problems. Scaling transformation replaces original phase function with a sum of the delta function and a new regular smooth phase function. Geometric truncation is one of the ways to construct such a smooth function. The replacement phase function coincides with the original one outside the forward cone and preserves the asymmetry parameter. It has discontinuity at the cone. Another simple functional form of the replacement phase function within the cone is suggested. It enables continuity and allows for a number of modifications. Three of them are considered in this study: preserving asymmetry parameter, providing continuity of the 1st derivative of the phase function, and preserving mean scattering angle. Yet another problem addressed in this study is objective selection of the width of the forward cone. That angle affects truncation fraction and values of the phase function within the cone. A heuristic approach providing unambiguous criterion of selection of the truncation angle is proposed. The approach has easy numerical implementation. Suggested modifications were tested on cloud phase function using discrete ordinates and Monte Carlo methods. It was shown that the modifications provide better accuracy of the radiance computation compare to the original geometric truncation with discrete ordinates while continuous derivative approach provides significant gain in computer time with Monte Carlo simulations.

Alexander Radkevich↗
NASA NTRS1993

On the accuracy of the Eddington approximation for radiative transfer in the microwave frequencies

The paper examines how well an Eddington approximation can reproduce brightness temperatures obtained from a more complete, N-stream discrete ordinate solution in the microwave regime. Radiation propagation through a plane parallel medium is considered. Although model discrepancies are complicated functions of the cloud constituents, the differences between an eight-stream discrete ordinate solution and an analytical Eddington solution were found to be generally small, ranging from 0 to 6 K when only one uniform layer of hydrometeors was considered. When realistic multilayered cloud hydrometeor profiles were used, the differences between these two models never exceeded 3 K over the entire range of microwave frequencies considered (6.6-183 GHz). The models agreed to within 0.2 K in the absence of scattering constituents.

Kummerow, Christian↗
NASA NTRS2001

Satellite Remote Sensing of Tropical Precipitation and Ice Clouds for GCM Verification

This project, supported by the NASA New Investigator Program, has primarily been funding a graduate student, Darren McKague. Since August 1999 Darren has been working part time at Raytheon, while continuing his PhD research. Darren is planning to finish his thesis work in May 2001, thus some of the work described here is ongoing. The proposed research was to use GOES visible and infrared imager data and SSM/I microwave data to obtain joint distributions of cirrus cloud ice mass and precipitation for a study region in the Eastern Tropical Pacific. These joint distributions of cirrus cloud and rainfall were to be compared to those from the CSU general circulation model to evaluate the cloud microphysical amd cumulus parameterizations in the GCM. Existing algorithms were to be used for the retrieval of cloud ice water path from GOES (Minnis) and rainfall from SSM/I (Wilheit). A theoretical study using radiative transfer models and realistic variations in cloud and precipitation profiles was to be used to estimate the retrieval errors. Due to the unavailability of the GOES satellite cloud retrieval algorithm from Dr. Minnis (a co-PI), there was a change in the approach and emphasis of the project. The new approach was to develop a completely new type of remote sensing algorithm - one to directly retrieve joint probability density functions (pdf's) of cloud properties from multi-dimensional histograms of satellite radiances. The usual approach is to retrieve individual pixels of variables (i.e. cloud optical depth), and then aggregate the information. Only statistical information is actually needed, however, and so a more direct method is desirable. We developed forward radiative transfer models for the SSM/I and GOES channels, originally for testing the retrieval algorithms. The visible and near infrared ice scattering information is obtained from geometric ray tracing of fractal ice crystals (Andreas Macke), while the mid-infrared and microwave scattering is computed with Mie scattering. The radiative transfer is performed with the Spherical Harmonic Discrete Ordinate Method (developed by the PI), and infrared molecular absorption is included with the correlated k-distribution method. The SHDOM radiances have been validated by comparison to version 2 of DISORT (the community "standard" discrete-ordinates radiative transfer model), however we use SHDOM since it is computationally more efficient.

Evans, K. Franklin↗
NASA NTRS1972

Application of DOT-MORSE coupling to the analysis of three-dimensional SNAP shielding problems

The use of discrete ordinates and Monte Carlo techniques to solve radiation transport problems is discussed. A general discussion of two possible coupling schemes is given for the two methods. The calculation of the reactor radiation scattered from a docked service and command module is used as an example of coupling discrete ordinates (DOT) and Monte Carlo (MORSE) calculations.

Straker, E. A.↗
NASA NTRS1992

FDDO and DSMC analyses of rarefied gas flow through 2D nozzles

Two different approaches, the finite-difference method coupled with the discrete-ordinate method (FDDO), and the direct-simulation Monte Carlo (DSMC) method, are used in the analysis of the flow of a rarefied gas expanding through a two-dimensional nozzle and into a surrounding low-density environment. In the FDDO analysis, by employing the discrete-ordinate method, the Boltzmann equation simplified by a model collision integral is transformed to a set of partial differential equations which are continuous in physical space but are point functions in molecular velocity space. The set of partial differential equations are solved by means of a finite-difference approximation. In the DSMC analysis, the variable hard sphere model is used as a molecular model and the no time counter method is employed as a collision sampling technique. The results of both the FDDO and the DSMC methods show good agreement. The FDDO method requires less computational effort than the DSMC method by factors of 10 to 40 in CPU time, depending on the degree of rarefaction.

Chung, Chan-Hong↗
NASA NTRS2012

Computing Radiative Transfer in a 3D Medium

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.

Von Allmen, Paul↗
NASA NTRS1990

Numerical simulation of rarefied gas flow through a slit

Two different approaches, the finite-difference method coupled with the discrete-ordinate method (FDDO), and the direct-simulation Monte Carlo (DSMC) method, are used in the analysis of the flow of a rarefied gas from one reservoir to another through a two-dimensional slit. The cases considered are for hard vacuum downstream pressure, finite pressure ratios, and isobaric pressure with thermal diffusion, which are not well established in spite of the simplicity of the flow field. In the FDDO analysis, by employing the discrete-ordinate method, the Boltzmann equation simplified by a model collision integral is transformed to a set of partial differential equations which are continuous in physical space but are point functions in molecular velocity space. The set of partial differential equations are solved by means of a finite-difference approximation. In the DSMC analysis, three kinds of collision sampling techniques, the time counter (TC) method, the null collision (NC) method, and the no time counter (NTC) method, are used.

Keith, Theo G., Jr.↗
NASA NTRS1993

Multilevel methods for transport equations in diffusive regimes

We consider the numerical solution of the single-group, steady state, isotropic transport equation. An analysis by means of the moment equations shows that a discrete ordinate S(sub N) discretization in direction (angle) with a least squares finite element discretization in space does not behave properly in the diffusion limit. A scaling of the S(sub N) equations is introduced so that the least squares discretization has the correct diffusion limit. For the resulting discrete system a full multigrid algorithm was developed.

Manteuffel, Thomas A.↗