Results for “Discrete ordinates method” — Page 7

NASA NTRS1986

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.↗

NASA NTRS2016

Accuracy of RT Code SORD for Realistic Atmospheric Profiles

We discuss accuracy of our recently developed RT code SORD using 2 benchmark scenarios published by the IPRT group in 2015. These scenarios define atmospheres with a complicate dependence of scattering and absorption properties over height (profile). Equal step, dh=1km, is assumed in the profiles. We developed subroutines that split such atmospheres into layers of the same optical thickness, d tau. We provide full text of the subroutines with comments in Appendix. The d tau is a step for vertical integration in the method of successive orders. Modification of the input profiles from "equal step over h" to "equal step over tau" changes input for RT simulations. This may cause errors at or above the acceptable level of the measurement uncertainty. We show errors of the RT code SORD for both intensity and polarization. In addition to that, using our discrete ordinates RT code IPOL, we discuss one more IPRT scenario, in which changes in height profile indeed cause unacceptable errors. Clear understanding of source and magnitude of these errors is important, e.g. for the AERONET retrieval algorithm. SORD is available from ftp://maiac.gsfc.nasa.gov/pub/skorkin/SORD_IP_16B/ or by email request from the first author.

polarization↗

NASA NTRS1978

Photoelectrons in the upper atmosphere - A formulation incorporating effects of transport

A technique for solving the problem of photoelectron spatial transport and spectral degradation is presented. The developed model, which allows for a relatively complete description of the relevant transport and collisional physics, can be applied to study processes at either low or high altitude, and may be used to investigate interhemispheric transport. The technique involves dividing the spatial dimension into a set of discrete ordinates and calculating the photoelectron intensity for a specific energy at each point in space and at a set of angle points. The selected range of the energy spectrum was covered completely by a set of energy bins, and the photoelectron intensity in the highest energy bin was calculated first by assuming that the only source term was derived from photoionization. The method is similar to a method used in radiative transfer theory. Spatial and angular redistribution of electrons play a dominant role at heights above 400 km.

Prather, M. J.↗

NASA NTRS2001

Gamma-Weighted Discrete Ordinate Two-Stream Approximation for Computation of Domain Averaged Solar Irradiance

An algorithm is developed for the gamma-weighted discrete ordinate two-stream approximation that computes profiles of domain-averaged shortwave irradiances for horizontally inhomogeneous cloudy atmospheres. The algorithm assumes that frequency distributions of cloud optical depth at unresolved scales can be represented by a gamma distribution though it neglects net horizontal transport of radiation. This algorithm is an alternative to the one used in earlier studies that adopted the adding method. At present, only overcast cloudy layers are permitted.

Kato, S.↗

NASA NTRS1987

Optimal application of Morrison's iterative noise removal for deconvolution. Appendices

Morrison's iterative method of noise removal, or Morrison's smoothing, is applied in a simulation to noise-added data sets of various noise levels to determine its optimum use. Morrison's smoothing is applied for noise removal alone, and for noise removal prior to deconvolution. For the latter, an accurate method is analyzed to provide confidence in the optimization. The method consists of convolving the data with an inverse filter calculated by taking the inverse discrete Fourier transform of the reciprocal of the transform of the response of the system. Various length filters are calculated for the narrow and wide Gaussian response functions used. Deconvolution of non-noisy data is performed, and the error in each deconvolution calculated. Plots are produced of error versus filter length; and from these plots the most accurate length filters determined. The statistical methodologies employed in the optimizations of Morrison's method are similar. A typical peak-type input is selected and convolved with the two response functions to produce the data sets to be analyzed. Both constant and ordinate-dependent Gaussian distributed noise is added to the data, where the noise levels of the data are characterized by their signal-to-noise ratios. The error measures employed in the optimizations are the L1 and L2 norms. Results of the optimizations for both Gaussians, both noise types, and both norms include figures of optimum iteration number and error improvement versus signal-to-noise ratio, and tables of results. The statistical variation of all quantities considered is also given.

Ioup, George E.↗

NASA NTRS1972

Common radiation analysis model for 75,000 pound thrust NERVA engine (1137400E)

The mathematical model and sources of radiation used for the radiation analysis and shielding activities in support of the design of the 1137400E version of the 75,000 lbs thrust NERVA engine are presented. The nuclear subsystem (NSS) and non-nuclear components are discussed. The geometrical model for the NSS is two dimensional as required for the DOT discrete ordinates computer code or for an azimuthally symetrical three dimensional Point Kernel or Monte Carlo code. The geometrical model for the non-nuclear components is three dimensional in the FASTER geometry format. This geometry routine is inherent in the ANSC versions of the QAD and GGG Point Kernal programs and the COHORT Monte Carlo program. Data are included pertaining to a pressure vessel surface radiation source data tape which has been used as the basis for starting ANSC analyses with the DASH code to bridge into the COHORT Monte Carlo code using the WANL supplied DOT angular flux leakage data. In addition to the model descriptions and sources of radiation, the methods of analyses are briefly described.

Warman, E. A.↗

NASA NTRSDate not supplied

Using Polarized PCRTM-Solar Model to Validate EMIT Surface Reflectance L2A Product

Using the vector linearized discrete ordinate radiative transfer model (VLIDORT), a multiple-scattering multi-layer scattering code, coupled with a fast principal component (PC) radiative transfer model (PCRTM), we can simulate hyperspectral radiation fields from visible to shortwave infrared (VSWIR), i.e., 250 nm to 2500 nm. Since only limited accurate line-by-line calculation using VLIDORT is needed, this coupled PCRTM-VLIDORT approach is very efficient for simulating hyperspectral VSWIR spectra. These simulations incorporate accurately defined optical properties (gas absorptions, molecular and aerosol/cloud particle scattering) and detailed surface reflectance inputs. In this study we used the current L2A products of the Earth Surface Mineral Dust Source Investigation (EMIT) imaging spectrometer as the inputs to PCRTM-VLIDORT to simulate EMIT spectra. Comparing the simulated spectra with EMIT measurements provides an efficient method to validate the current EMIT product.

Xiaozhen Xiong↗

NASA NTRS2017

Validation of a Solid Rocket Motor Internal Environment Model

In a prior effort, a thermal/fluid model of the interior of Penn State University's laboratory-scale Insulation Test Motor (ITM) was constructed to predict both the convective and radiative heat transfer to the interior walls of the ITM with a minimum of empiricism. These predictions were then compared to values of total and radiative heat flux measured in a previous series of ITM test firings to assess the capabilities and shortcomings of the chosen modeling approach. Though the calculated fluxes reasonably agreed with those measured during testing, this exercise revealed means of improving the fidelity of the model to, in the case of the thermal radiation, enable direct comparison of the measured and calculated fluxes and, for the total heat flux, compute a value indicative of the average measured condition. By replacing the P1-Approximation with the discrete ordinates (DO) model for the solution of the gray radiative transfer equation, the radiation intensity field in the optically thin region near the radiometer is accurately estimated, allowing the thermal radiation flux to be calculated on the heat-flux sensor itself, which was then compared directly to the measured values. Though the fully coupling the wall thermal response with the flow model was not attempted due to the excessive computational time required, a separate wall thermal response model was used to better estimate the average temperature of the graphite surfaces upstream of the heat flux gauges and improve the accuracy of both the total and radiative heat flux computations. The success of this modeling approach increases confidence in the ability of state-of-the-art thermal and fluid modeling to accurately predict SRM internal environments, offers corrections to older methods, and supplies a tool for further studies of the dynamics of SRM interiors.

Martin, Heath T.↗

NASA NTRS2017

Validation of a Solid Rocket Motor Internal Environment Model

In a prior effort, a thermal/fluid model of the interior of Penn State University's laboratory-scale Insulation Test Motor (ITM) was constructed to predict both the convective and radiative heat transfer to the interior walls of the ITM with a minimum of empiricism. These predictions were then compared to values of total and radiative heat flux measured in a previous series of ITM test firings to assess the capabilities and shortcomings of the chosen modeling approach. Though the calculated fluxes reasonably agreed with those measured during testing, this exercise revealed means of improving the fidelity of the model to, in the case of the thermal radiation, enable direct comparison of the measured and calculated fluxes and, for the total heat flux, compute a value indicative of the average measured condition. By replacing the P1-Approximation with the discrete ordinates (DO) model for the solution of the gray radiative transfer equation, the radiation intensity field in the optically thin region near the radiometer is accurately estimated, allowing the thermal radiation flux to be calculated on the heat-flux sensor itself, which was then compared directly to the measured values. Though the fully coupling the wall thermal response with the flow model was not attempted due to the excessive computational time required, a separate wall thermal response model was used to better estimate the average temperature of the graphite surfaces upstream of the heat flux gauges and improve the accuracy of both the total and radiative heat flux computations. The success of this modeling approach increases confidence in the ability of state-of-the-art thermal and fluid modeling to accurately predict SRM internal environments, offers corrections to older methods, and supplies a tool for further studies of the dynamics of SRM interiors.

Martin, Heath T.↗

NASA NTRS2016

An Optimal-Estimation-Based Aerosol Retrieval Algorithm Using OMI Near-UV Observations

An optimal-estimation(OE)-based aerosol retrieval algorithm using the OMI (Ozone Monitoring Instrument) near-ultraviolet observation was developed in this study. The OE-based algorithm has the merit of providing useful estimates of errors simultaneously with the inversion products. Furthermore, instead of using the traditional lookup tables for inversion, it performs online radiative transfer calculations with the VLIDORT (linearized pseudo-spherical vector discrete ordinate radiative transfer code) to eliminate interpolation errors and improve stability. The measurements and inversion products of the Distributed Regional Aerosol Gridded Observation Network campaign in northeast Asia (DRAGON NE-Asia 2012) were used to validate the retrieved aerosol optical thickness (AOT) and single scattering albedo (SSA). The retrieved AOT and SSA at 388 nm have a correlation with the Aerosol Robotic Network (AERONET) products that is comparable to or better than the correlation with the operational product during the campaign. The OEbased estimated error represented the variance of actual biases of AOT at 388 nm between the retrieval and AERONET measurements better than the operational error estimates. The forward model parameter errors were analyzed separately for both AOT and SSA retrievals. The surface reflectance at 388 nm, the imaginary part of the refractive index at 354 nm, and the number fine-mode fraction (FMF) were found to be the most important parameters affecting the retrieval accuracy of AOT, while FMF was the most important parameter for the SSA retrieval. The additional information provided with the retrievals, including the estimated error and degrees of freedom, is expected to be valuable for relevant studies. Detailed advantages of using the OE method were described and discussed in this paper.

Air Pollution↗