Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ray tracing”

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

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

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1978

Winston Solar Concentrators and Evaluation Support. Phase 2: Non-imaging Concentrators for Space Applications

A 4.67X, plus or minus 5 deg. compound parabolic concentrator (CPC) for a large photovoltaic array in space was analyzed. The design was demonstrated to be effective in achieving a net power gain which can be varied from more than a factor of 3 down to approximately unity. A method for reducing nonuniformities in illumination to a given desired level was found. The effectiveness of this method, which involves the introduction of a degree of non-specularity in the reflector surface, was confirmed by direct measurements with prepared foil reflectors in a CPC in terrestrial sunshine as well as by computer ray tracing. Further ray tracing confirms that the CPC design is extremely tolerant to pointing and alignment errors, minor distortions, etc. A two stage non-imaging design was shown, by preliminary measurements and analysis, to provide both the desired angular tolerance and the required degree of intensity uniformity if higher concentrations (4X-10X) are necessary.

Winston, R.↗
NASA NTRS1965

Computer programs simplify optical system analysis

The optical ray-trace computer program performs geometrical ray tracing. The energy-trace program calculates the relative monochromatic flux density on a specific target area. This program uses the ray-trace program as a subroutine to generate a representation of the optical system.

Source record↗
NASA NTRS1999

Wave Dynamics and Transport in the Stratosphere

The report discusses: (1) Gravity waves generated by tropical convection: A study in which a two-dimensional cloud-resolving model was used to examine the possible role of gravity waves generated by a simulated tropical squall line in forcing the quasi-biennial oscillation was completed. (2) Gravity wave ray tracing studies:It was developed a linear ray tracing model of gravity wave propagation to extend the nonlinear storm model results into the mesosphere and thermosphere. (3) tracer filamentation: Vertical soundings of stratospheric ozone often exhibit laminated tracer structures characterized by strong vertical tracer gradients. (4) Mesospheric gravity wave modeling studies: Although our emphasis in numerical simulation of gravity waves generated by convection has shifted from simulation of idealized two-dimensional squall lines to the most realistic (and complex) study of wave generation by three-dimensional storms. (5) Gravity wave climatology studies: Mr. Alexander applied a linear gravity wave propagation model together with observations of the background wind and stability fields to compute climatologies of gravity wave activity for comparison to observations. (6) Convective forcing of gravity waves: Theoretical study of gravity wave forcing by convective heat sources has completed. (7) Gravity waves observation from UARS: The objective of this work is to apply ray tracing, and other model technique, in order to determine to what extend the horizontal and vertical variation in satellite observed distribution of small-scale temperature variance can be attributed to gravity waves from particular sources. (8) The annual and interannual variations in temperature and mass flux near the tropical tropopause. and (9) Three dimensional cloud model.

Holton, James R.↗
NASA NTRS1982

Tolerance analysis of optical telescopes using coherent addition of wavefront errors

A near diffraction-limited telescope requires that tolerance analysis be done on the basis of system wavefront error. One method of analyzing the wavefront error is to represent the wavefront error function in terms of its Zernike polynomial expansion. A Ramsey-Korsch ray trace package, a computer program that simulates the tracing of rays through an optical telescope system, was expanded to include the Zernike polynomial expansion up through the fifth-order spherical term. An option to determine a 3 dimensional plot of the wavefront error function was also included in the Ramsey-Korsch package. Several assimulation runs were analyzed to determine the particular set of coefficients in the Zernike expansion that are effected by various errors such as tilt, decenter and despace. A 3 dimensional plot of each error up through the fifth-order spherical term was also included in the study. Tolerance analysis data are presented.

Davenport, J. W.↗
NASA NTRS1994

Refraction of high frequency noise in an arbitrary jet flow

Refraction of high frequency noise by mean flow gradients in a jet is studied using the ray-tracing methods of geometrical acoustics. Both the two-dimensional (2D) and three-dimensional (3D) formulations are considered. In the former case, the mean flow is assumed parallel and the governing propagation equations are described by a system of four first order ordinary differential equations. The 3D formulation, on the other hand, accounts for the jet spreading as well as the axial flow development. In this case, a system of six first order differential equations are solved to trace a ray from its source location to an observer in the far field. For subsonic jets with a small spreading angle both methods lead to similar results outside the zone of silence. However, with increasing jet speed the two prediction models diverge to the point where the parallel flow assumption is no longer justified. The Doppler factor of supersonic jets as influenced by the refraction effects is discussed and compared with the conventional modified Doppler factor.

Khavaran, Abbas↗
NASA NTRS2019

Understanding the Impact of Meteorological Spatiotemporal Variability on Distant Focusing Overpressure Risk

Since the 1950s, the amplification of blast energy from explosions has been understood to be a significant hazard to public safety at launch ranges. Historically, the risk assessment of the Distant Focusing Overpressure (DFO) hazard started with a single temperature and wind profile (e.g., a radiosonde launch) as input to acoustic 1-D ray tracing models. By analyzing rays propagation and performing ray tracing, population centers under high DFO risk can be identified. Although this method is useful, less is known about how the blast waves can be focused when the spatiotemporal variability of the input profiles are considered. In summary, this work aims to consider how realistic atmospheric boundary layer variability (e.g., turbulence, land-surface contrasts) may affect blast waves propagation and focusing and, as a result, DFO risk assessment results.

Blumberg, W. Greg↗
NASA NTRS1972

Moment expansion for ionospheric range error

On a plane earth, the ionospheric or tropospheric range error depends only on the total refractivity content or zeroth moment of the refracting layer and the elevation angle. On a spherical earth, however, the dependence is more complex; so for more accurate results it has been necessary to resort to complex ray-tracing calculations. A simple, high-accuracy alternative to the ray-tracing calculation is presented. By appropriate expansion of the angular dependence in the ray-tracing integral in a power series in height, an expression is obtained for the range error in terms of a simple function of elevation angle, E, at the expansion height and of the mth moment of the refractivity, N, distribution about the expansion height. The rapidity of convergence is heavily dependent on the choice of expansion height. For expansion heights in the neighborhood of the centroid of the layer (300-490 km), the expansion to N = 2 (three terms) gives results accurate to about 0.4% at E = 10 deg. As an analytic tool, the expansion affords some insight on the influence of layer shape on range errors in special problems.

Mallinckrodt, A.↗
NASA NTRS1998

Polarization Considerations

As light passes through a optical system the reflections and refractions will in general change the polarization state of the light. If we assume that all of the materials in the thin film coatings and substrate are isotropic and homogeneous then calculating the amount of "instrumental" polarization is a relatively straight forward task. In the following sections we will present a of the steps required to perform a 'polarization ray trace' calculation for a single ray and monochromatic and hence polarized light. The thin film portion of the calculation is also shown. The reason for explicitly showing the thin film equations is that there are sign conventions imposed on the boundary value equations by the orientation and handedness of the various coordinate frames which are attached to the geometric rays. The attenuation of light through a optical system, is relatively simple, and requires at the very least a lens (average) reflectivity or transmissivity. Determining the polarization sensitivity of a optical system is still relatively straight forward requiring at least a knowledge of the behavior of the "s" and "p" components at each interface for the chief ray. Determining the thin film induced aberrations of a optical system are somewhat more demanding. Questions about the arithmetic sign of the phase factors and how this relates to the overall "OPD" of a ray are ubiquitous. Many rays are required to construct a wavefront. Thin film codes which modify the OPD's of rays are a requirement for this last mentioned computation. This requires a consistent scheme of coordinate frames and sign conventions and is probably the most demanding task of a polarization ray trace. Only the electric field will used in the discussion. This is not a restriction as the Stokes parameters are functions of the electric field. The following does not attempt to explain, but only to present all of the required concepts and formulas.

Waluschka, Eugene↗
NASA NTRS1993

TransFit: Finite element analysis data fitting software

The Advanced X-Ray Astrophysics Facility (AXAF) mission support team has made extensive use of geometric ray tracing to analyze the performance of AXAF developmental and flight optics. One important aspect of this performance modeling is the incorporation of finite element analysis (FEA) data into the surface deformations of the optical elements. TransFit is software designed for the fitting of FEA data of Wolter I optical surface distortions with a continuous surface description which can then be used by SAO's analytic ray tracing software, currently OSAC (Optical Surface Analysis Code). The improved capabilities of Transfit over previous methods include bicubic spline fitting of FEA data to accommodate higher spatial frequency distortions, fitted data visualization for assessing the quality of fit, the ability to accommodate input data from three FEA codes plus other standard formats, and options for alignment of the model coordinate system with the ray trace coordinate system. TransFit uses the AnswerGarden graphical user interface (GUI) to edit input parameters and then access routines written in PV-WAVE, C, and FORTRAN to allow the user to interactively create, evaluate, and modify the fit. The topics covered include an introduction to TransFit: requirements, designs philosophy, and implementation; design specifics: modules, parameters, fitting algorithms, and data displays; a procedural example; verification of performance; future work; and appendices on online help and ray trace results of the verification section.

Freeman, Mark↗
NASA NTRS1976

The effects of atmospheric refraction on the accuracy of laser ranging systems

Correction formulas derived by Saastamoinen and Marini, and the ray traces through the refractivity profiles all assume a spherically symmetric refractivity profile. The errors introduced by this assumption were investigated by ray tracing through three-dimensional profiles. The results of this investigation indicate that the difference between ray traces through the spherically symmetric and three-dimensional profiles is approximately three centimeters at 10 deg and decreases to less than one half of a centimeter at 80 deg. If the accuracy desired in future laser ranging systems is less than a few centimeters, Saastamoinen and Marini's formulas must be altered to account for the fact that the refractivity profile is not spherically symmetric.

Zanter, D. L.↗
NASA NTRS1996

Evaluation of Two Computational Techniques of Calculating Multipath Using Global Positioning System Carrier Phase Measurements

Two computational techniques are used to calculate differential phase errors on Global Positioning System (GPS) carrier war phase measurements due to certain multipath-producing objects. The two computational techniques are a rigorous computati electromagnetics technique called Geometric Theory of Diffraction (GTD) and the other is a simple ray tracing method. The GTD technique has been used successfully to predict microwave propagation characteristics by taking into account the dominant multipath components due to reflections and diffractions from scattering structures. The ray tracing technique only solves for reflected signals. The results from the two techniques are compared to GPS differential carrier phase ns taken on the ground using a GPS receiver in the presence of typical International Space Station (ISS) interference structures. The calculations produced using the GTD code compared to the measured results better than the ray tracing technique. The agreement was good, demonstrating that the phase errors due to multipath can be modeled and characterized using the GTD technique and characterized to a lesser fidelity using the DECAT technique. However, some discrepancies were observed. Most of the discrepancies occurred at lower devations and were either due to phase center deviations of the antenna, the background multipath environment, or the receiver itself. Selected measured and predicted differential carrier phase error results are presented and compared. Results indicate that reflections and diffractions caused by the multipath producers, located near the GPS antennas, can produce phase shifts of greater than 10 mm, and as high as 95 mm. It should be noted tl the field test configuration was meant to simulate typical ISS structures, but the two environments are not identical. The GZ and DECAT techniques have been used to calculate phase errors due to multipath o the ISS configuration to quantify the expected attitude determination errors.

Gomez, Susan F.↗
NASA NTRS2006

Self-Consistent Model of Magnetospheric Ring Current and Propagating Electromagnetic Ion Cyclotron Waves: Waves in Multi Ion Magnetosphere - 1

The further development of a self-consistent theoretical model of interacting ring current ions and electromagnetic ion cyclotron waves [Khazanov et al., 2003] is presented. In order to adequately take into account the wave propagation and refraction in a multi-ion plasmasphere, we explicitly include the ray tracing equations in our previous self-consistent model and use the general form of the wave kinetic equation. This is a major new feature of the present model and, to the best of our knowledge, the ray tracing equations for the first time are explicitly employed on a global magnetospheric scale in order to self-consistently simulate spatial, temporal, and spectral evolutions of the ring current and electromagnetic ion cyclotron waves. To demonstrate the effects of EMIC wave propagation and refraction on the EMIC wave energy distributions and evolution we simulate the May 1998 storm. The main findings of our simulation can be summarized as follows. First, due to the density gradient at the plasmapause, the net wave refraction is suppressed, and He(+)-mode grows preferably at plasmapause. This result is in a total agreement with the previous ray tracing studies, and very clear observed in presented B-field spectrograms. Second, comparison the global wave distributions with the results from other ring current model [Kozyra et al., 1997] reveals that our model provides more intense and higher plasmapause organized distributions during the May, 1998 storm period. Finally, the found He(+)-mode energy distributions are not Gaussian distributions, and most important that wave energy can occupy not only the region of generation, i. e. the region of small wave normal angles, but the entire wave normal angle region and even only the region near 90 degrees. The latter is extremely crucial for energy transfer to thermal plasmaspheric electrons by resonant Landau damping, and subsequent downward heat transport and excitation of stable auroral red arcs.

Khazanov, G. V.↗
NASA NTRS2006

Self-Consistent Model of Magnetospheric Ring Current and Propagating Electromagnetic Ion Cyclotron Waves: Waves in Multi-Ion Magnetosphere

The further development of a self-consistent theoretical model of interacting ring current ions and electromagnetic ion cyclotron waves (Khazanov et al., 2003) is presented In order to adequately take into account wave propagation and refraction in a multi-ion magnetosphere, we explicitly include the ray tracing equations in our previous self-consistent model and use the general form of the wave kinetic equation. This is a major new feature of the present model and, to the best of our knowledge, the ray tracing equations for the first time are explicitly employed on a global magnetospheric scale in order to self-consistently simulate the spatial, temporal, and spectral evolution of the ring current and of electromagnetic ion cyclotron waves To demonstrate the effects of EMIC wave propagation and refraction on the wave energy distribution and evolution, we simulate the May 1998 storm. The main findings of our simulation can be summarized as follows. First, owing to the density gradient at the plasmapause, the net wave refraction is suppressed, and He+-mode grows preferably at the plasmapause. This result is in total agreement with previous ray tracing studies and is very clearly found in presented B field spectrograms. Second, comparison of global wave distributions with the results from another ring current model (Kozyra et al., 1997) reveals that this new model provides more intense and more highly plasmapause-organized wave distributions during the May 1998 storm period Finally, it is found that He(+)-mode energy distributions are not Gaussian distributions and most important that wave energy can occupy not only the region of generation, i.e., the region of small wave normal angles, but all wave normal angles, including those to near 90 . The latter is extremely crucial for energy transfer to thermal plasmaspheric electrons by resonant Landau damping and subsequent downward heat transport and excitation of stable auroral red arcs.

Khazanov, G. V.↗
NASA NTRS1996

A Simulated Spectrum of Convectively Generated Gravity Waves: Propagation from the Tropopause to the Mesopause and Effects on the Middle Atmosphere

This work evaluates the interaction of a simulated spectrum of convectively generated gravity waves with realistic middle atmosphere mean winds. The wave spectrum is derived from the nonlinear convection model described by Alexander et al. that simulated a two-dimensional midlatitude squall line. This spectrum becomes input to a linear ray tracing model for evaluation of wave propagation as a function of height through climatological background wind and buoyancy frequency profiles. The energy defined by the spectrum as a function of wavenumber and frequency is distributed spatially and temporally into wave packets for the purpose of estimating wave amplitudes at the lower boundary of the ray tracing model. A wavelet analysis provides an estimate of these wave packet widths in space and time. Without this redistribution of energies into wave packets the Fourier analysis alone inaccurately assumes the energy is evenly distributed throughout the storm model domain. The growth with height of wave amplitudes is derived from wave action flux conservation coupled to a convective instability saturation condition. Mean flow accelerations and wave energy dissipation profiles are derived from this analysis and compared to parameterized estimates of gravity wave forcing, providing a measure of the importance of the storm source to global gravity wave forcing. The results suggest that a single large convective storm system like the simulated squall line could provide a significant fraction of the zonal mean gravity wave forcing at some levels, particularly in the mesosphere. The vertical distributions of mean flow acceleration and energy dissipation do not much resemble the parameterized profiles in form because of the peculiarities of the spectral properties of the waves from the storm source. The ray tracing model developed herein provides a tool for examining the role of convectively generated waves in middle atmosphere physics.

Alexander, Joan↗
NASA NTRS1996

A Simulated Spectrum of Convectively Generated Gravity Waves: Propagation from the Tropopause to the Mesopause and Effects on the Middle Atmosphere

This work evaluates the interaction of a simulated spectrum of convectively generated gravity waves with realistic middle atmosphere mean winds. The wave spectrum is derived from the nonlinear convection model described by Alexander et al. [1995] that simulated a two-dimensional midlatitude squall line. This spectrum becomes input to a linear ray tracing model for evaluation of wave propagation as a function of height through climatological background wind and buoyancy frequency profiles. The energy defined by the spectrum as a function of wavenumber and frequency is distributed spatially and temporally into wave packets for the purpose of estimating wave amplitudes at the lower boundary of the ray tracing model. A wavelet analysis provides an estimate of these wave packet widths in space and time. Without this redistribution of energies into wave packets the Fourier analysis alone inaccurately assumes the energy is evenly distributed throughout the storm model domain. The growth with height of wave amplitudes is derived from wave action flux conservation coupled to a convective instability saturation condition. Mean flow accelerations and wave energy dissipation profiles are derived from this analysis and compared to parameterized estimates of gravity wave forcing, providing a measure of the importance of the storm source to global gravity wave forcing. The results suggest that a single large convective storm system like the simulated squall line could provide a significant fraction of the zonal mean gravity wave forcing at some levels, particularly in the mesosphere. The vertical distributions of mean flow acceleration and energy dissipation do not much resemble the parameterized profiles in form because of the peculiarities of the spectral properties of the waves from the storm source. The ray tracing model developed herein provides a tool for examining the role of convectively generated waves in middle atmosphere physics.

Alexander, M. Joan↗
NASA NTRS2000

Recent Trends in the Analysis of Quasioptical Systems

The recent trend in microwave instruments is the use of multiple millimeter and submillimeter wavelength bands. These systems are typically analyzed by using physical optics, Gaussian beams or ray tracing techniques. Physical optics offers high accuracy at the expense of computation time. This trade-off becomes particularly apparent in the analysis of multiple reflector antennas, such as beam waveguide antennas, where physical optics is used to compute the current on each reflector from the current on the previous reflector. At the other end of the spectrum is ray tracing approaches that ignore diffraction effects entirely. These methods are fast but sacrifice the ability to predict some effects accurately. An intermediate approach is to use an appropriate set of expansion functions to model the field between the reflectors. If the set is chosen wisely only a few coefficients need to be determined from each reflector current. The field is then computed at the next reflector through the use of the expansion functions and their coefficients rather than by using the previous reflector current. For a beam waveguide system with no enclosing tubes an excellent set of expansion functions is the Gaussian beam mode set. In many cases a preliminary design which includes the effects on diffraction may be obtained by considering only the fundamental mode and a thin lens model for the reflectors. Higher-order modes are included to model the effects of the curved reflector, which include asymmetric distortion of the beam, cross polarization, and beam truncation. This paper describes a computer code implementing higher-order Gaussian beam scattering by multiple reflector systems. There are four essential steps in the algorithm. (1) Compute the current on the first reflector using physical optics using either a feed model or by an incident set of Gaussian beam modes. (2) Find the direction of propagation for the reflected Gaussian beam-set using ray tracing. (3) Determine the waist size and location for the output beam set by examining the amplitude and phase distribution of the current on the reflector. (4) Compute the amplitudes of the individual modes in the output mode set. These steps are then repeated for each addition reflector in the chain. In each of these cases the previous Gaussian beam set provides the input field for the current calculation. Details of the four steps discussed above will be discussed. Examples will compare results from the Gaussian beam approach to pure physical optics, illustrating both its merits and limitations. Hybrid approaches capable of eliminating some of the limitations will also be discussed.

Hoppe, D. J.↗
DOE OSTI2024

Qualitative Conditions for Reasonable Estimates of Kerma Using Total and Kerma Cross Section Attenuation of Photons

Performing a full adjoint simulation in the Integrated Tiger Series (ITS) can be time consuming to obtain statistically significant results. The ray-trace capability in ITS allows for rapid scoping calculations of the uncollided kerma from photon sources without needing to run a full adjoint simulation. However, the uncollided estimate will always underestimate the full-physics kerma since it neglects scattered radiation, and under certain conditions, the result from the capability may not provide a sufficiently accurate estimate of the full-physics kerma. To exemplify the conditions in which the capability provides reasonable estimates, two problem geometries with different materials are simulated using the full adjoint capability as well as the ray-trace capability with a total cross section treatment and a new kerma-attenuation cross-section treatment. The results are then compared to show under which conditions the estimates are accurate. Reasonable estimates are provided from the ray-trace feature when there is minimal scattering into a region of interest occurring, such as with low-Z material and when the photons travel through small amounts of material.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1981

Ray methods in random media

We review ray-optical methods of analyzing short-wavelength propagation in random media. The advantages and limitations of ray methods are discussed, and results of the statistical theory of ray segment fluctuations pertinent to ray tracing are summarized. The standard method of Monte Carlo ray tracing is compared to a new method which takes into account recent results on the statistics of ray segment fluctuations.

Hornstein, J.↗