Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 325 records · Page 18

Geometrical optics design of a compact range Gregorian subreflector system by the principle of the central ray

In recent years, the compact range has become very popular for measuring Radar Cross Section (RCS) and antenna patterns. The compact range, in fact, offers several advantages due to reduced size, a controlled environment, and privacy. On the other hand, it has some problems of its own, which must be solved properly in order to achieve high quality measurement results. For example, diffraction from the edges of the main reflector corrupts the plane wave in the target zone and creates spurious scattering centers in RCS measurements. While diffraction can be minimized by using rolled edges, the field of an offset single reflector compact range is corrupted by three other errors: the taper of the reflected field, the cross polarization introduced by the tilt of the feed and the aperture blockage introduced by the feed itself. These three errors can be eliminated by the use of a subreflector system. A properly designed subreflector system offers very little aperture blockage, no cross-polarization introduced and a minimization of the taper of the reflected field. A Gregorian configuration has been adopted in order to enclose the feed and the ellipsoidal subreflector in a lower chamber, which is isolated by absorbers from the upper chamber, where the main parabolic reflector and the target zone are enclosed. The coupling between the two rooms is performed through a coupling aperture. The first cut design for such a subreflector system is performed through Geometrical Optics ray tracing techniques (GO), and is greatly simplified by the use of the concept of the central ray introduced by Dragone. The purpose of the GO design is to establish the basic dimensions of the main reflector and subreflector, the size of the primary and secondary illuminating surfaces, the tilt angles of the subreflector and feed, and estimate the feed beamwidth. At the same time, the shape of the coupling aperture is initially determined.

Clerici, Giancarlo↗

IFSAR Simulation Using the Shooting and Bouncing Ray Technique

Interferometric Synthetic Aperture Radar (IFSAR) is a technique that allows an automated way to carry out terrain mapping. IFSAR is carried out by first generating a SAR image pair from two antennas that are spatially separated. The phase difference between the SAR image pair is proportional to the topography. After registering the SAR images, the difference in phase in each pixel is extracted to generate an interferogram. Since the phase can only be measured within 2pi radians, phase unwrapping is carried out to extract the absolute phase for each pixel that will be proportional to the local height. While IFSAR algorithm is typically applied to measurement data, it is useful to develop an IFSAR simulator to develop a better understanding of the IFSAR technique. The IFSAR simulator can be used in choosing system parameters, experimenting with processing procedures and mission planning. In this paper we will present an IFSAR simulation methodology to simulate the interferogram based on the shooting and bouncing ray (SBR) technique. SBR is a standard ray-tracing technique used to simulate scattering from large, complex targets. SBR is carried out by shooting rays at the target or scene. At the exit point of each ray, a ray-tube integration is done to find its contribution to the total field. A fast algorithm has been developed for the SBR for simulating SAR images of complex targets. In the IFSAR simulation, we build upon the fast SAR simulation technique. Given the antenna pair configuration, radar system parameters and the geometrical description of the scene, we first simulate two SAR images from each antenna. After post processing the two SAR images, we generate an interferogram. Phase unwrapping is then performed on the interferogram to arrive at the desired terrain map. We will present results from the SBR-based IFSAR simulator. The results will include terrain map reconstruction of urban environments. The reconstruction will be compared to the ground truth to examine the fidelity of the simulation. We will also investigate the effect of multi-bounce scattering in urban environments on phase unwrapping and reconstruction.

Houshmand, Bijan↗

An Analytical Model for the Prediction of a Micro-Dosimeter Response Function

A rapid analytical procedure for the prediction of a micro-dosimeter response function in low Earth orbit (LEO), correlated with the Space Transportation System (STS, shuttle) Tissue Equivalent Proportional Counter (TEPC) measurements is presented. The analytical model takes into consideration the energy loss straggling and chord length distribution of the detector, and is capable of predicting energy deposition fluctuations in a cylindrical micro-volume of arbitrary aspect ratio (height/diameter) by incoming ions through both direct and indirect (ray) events. At any designated (ray traced) target point within the vehicle, the model accepts the differential flux spectrum of Galactic Cosmic Rays (GCR) and/or trapped protons at LEO as input. On a desktop PC, the response function of TEPC for each ion in the GCR/trapped field is computed at the average rate of 30 seconds/ion. The ionizing radiation environment at LEO is represented by O'Neill fs GCR model (2004), covering charged particles in the 1 less than or equal to Z less than or equal to 28. O'Neill's free space GCR model is coupled with the Langley Research Center (LaRC) angular dependent geomagnetic cutoff model to compute the transmission coefficient in LEO. The trapped proton environment is represented by a LaRC developed time dependent procedure which couples the AP8MIN/AP8MAX, Deep River Neutron Monitor (DRNM) and F10.7 solar radio frequency measurements. The albedo neutron environment is represented by the extrapolation of the Atmospheric Ionizing Radiation (AIR) measurements. The charged particle transport calculations correlated with STS 51 and 114 flights are accomplished by using the most recent version (2005) of the LaRC deterministic High charge (Z) and Energy TRaNsport (HZETRN) code. We present the correlations between the TEPC model predictions (response function) and TEPC measured differential/integral spectra in the lineal energy (y) domain for both GCR and trapped protons, with the conclusion that the model correctly accounts for the increase in flux at low y values where energetic ions are the primary contributor. We further discuss that, even with the incorporation of angular dependency in the cutoffs, comparison of the GCR differential/integral flux between STS 51 and 114 TEPC measured data and current calculations indicates that there still exists an underestimation by the simulations at low to mid range y values. This underestimation is partly related the exclusion of the secondary pion particle production from the current version of HZETRN.

Badavi, Francis F.↗

Theoretical design and analysis of the layered synthetic microstructure optic for the dual path X-ray telescope

A ray tracing analysis was performed for several configurations for the inner channel of the dual path X-ray telescope, which is proposed to use the second mirror of the Stanford/MSFC Wolter-Schwarzchild telescope and a normal incident layered synthetic microstructure (LSM) mirror to form a secondary image near the front of the telescope. The LSM mirror shapes considered were spherical, ellipsoid, hyperboloid, and constant optical path length (OPL) aspheric. Only the constant OPL case gave good axial resolution. All cases had poor off axis resolution as judged by the RMS blur circle radius.

Shealy, D. L.↗

Computing material volume fractions on a superimposed mesh as applied to Monte Carlo particle transport simulations

Here, we present a newly implemented ray tracing algorithm in OpenMC for efficiently computing material volume fractions on superimposed meshes in complex geometries. By firing rays along each coordinate direction through the geometry, the approach accumulates track-length data in each mesh element, thereby determining the fractional composition of each material. Scaling studies on three different models—a random tetrahedra configuration, the Frascati Neutron Generator ITER dose rate benchmark, and a stellarator design—show excellent parallel performance, with nearly linear speedup on modern multi-threaded and distributed-memory systems. An analysis of the residual error relative to high-resolution reference solutions demonstrated that under optimal conditions it decreases as 1/R, where R is the number of rays fired, making it straightforward to achieve user-prescribed accuracy. This new functionality enables practical, mesh-based approaches for detailed nuclear analyses in production Monte Carlo workflows without resorting to expensive, fully conformal or unstructured meshing.

Monte Carlo↗

Observation of synergy between lower hybrid waves at two frequencies in EAST

Synergistic effects between two frequencies of lower hybrid (LH) waves—operating at 2.45 and 4.6 GHz—were observed in experiment on EAST for the first time. At low density (ne,lin≈2.0×1019 m−3), simultaneous injection of a 65/35 mix of 2.45/4.6 GHz power achieved an lower hybrid current drive efficiency that was 25% higher than what should be expected from the linear combination of the two sources. The experiment was interpreted with time-dependent simulations, using the equilibrium and transport solver TRANSP, coupled with the ray-tracing code GENRAY and the Fokker–Planck solver CQL3D. For each discharge, profiles of current and hard x-rays from simulation and measurement agree within uncertainties. An examination of the electron distribution function indicates that the LH synergy is supported by the increased width of the LH resonance plateau in the simultaneous injection case compared to independent injection.

Choi, W. (ORCID:0000000268912300)↗

Localization of infrasonic sources via Bayesian back projection

SUMMARY A Bayesian framework is investigated for event-specific localization of infrasonic sources using back projection ray tracing. Direction-of-arrival information from array-based detection analysis is used to initialize a back projection ray path originating from the detecting array location and quantifying propagation characteristics from hypothetical source locations. The Fisher statistic, computed from the array’s beam coherence, is mapped into uncertainty in the launch angles of the ray path. Auxiliary parameters previously introduced for solving the Transport equation to compute geometric spreading along ray paths are used to map uncertainty in the ray launch angles into spatial and temporal uncertainties in the ray path. An atmospheric ensemble approach is applied to account for atmospheric uncertainty, and the relation between uncertainties in the atmospheric state and confidence in estimated localization are evaluated using several ensembles with specified variances. The method is evaluated using a synthetic event in the western United States constructed via forward propagation simulations as well as a single-station, multi-arrival detection from a surface explosion in the western United States. Localization results using this event-specific approach are more accurate and exhibit improved precision than existing Bayesian localization methods that leverage generalized, pre-computed propagation statistics.

58 GEOSCIENCES↗

Development of high resolution imaging detectors for x ray astronomy

This final report summarizes our past activities and discusses the work performed over the period of 1 April 1990 through 1 April 1991 on x-ray optics, soft x-ray (0.1 - 10 KeV) imaging detectors, and hard x-ray (10 - 300 KeV) imaging detectors. If microchannel plates (MCPs) can be used to focus x-rays with a high efficiency and good angular resolution, they will revolutionize the field of x-ray optics. An x-ray image of a point source through an array of square MCP pores compared favorably with our ray tracing model for the MCP. Initial analysis of this image demonstrates the feasibility of MCPs for soft x-rays. Our work continues with optimizing the performance of our soft x-ray MCP imaging detectors. This work involves readout technology that should provide improved MCP readout devices (thin film crossed grid, curved, and resistive sheets), defect removal in MCPs, and photocathode optimization. In the area of hard x-ray detector development we have developed two different techniques for producing a CsI photocathode thickness of 10 to 100 microns, such that it is thick enough to absorb the high energy x-rays and still allow the photoelectrons to escape to the top MCP of a modified soft x-ray imaging detector. The methods involve vacuum depositing a thick film of CsI on a strong back, and producing a converter device that takes the place of the photocathode.

Murray, S. S.↗

Study of turbulence-induced refraction of lower hybrid waves using synthetic scrape-off layer filaments

Here, turbulence-induced refraction effects to lower hybrid (LH) wave propagation and current drive are studied using synthetic scrape-off layer (SOL) blob/filament fields. A synthetic 3D, field-following, blob turbulence model is implemented in the ray-tracing/Fokker-Planck (RTFP) codes GENRAY/CQL3D. In Alcator C-Mod, the blob field is shown to significantly affect LH ray-trajectories, leading to increased on-axis damping and smoother current profiles. This effect depends on the average blob size and amplitude. In addition, the diffusion of ray-trajectories in phase-space caused by turbulence increases the robustness of the RTFP model. A modified $N_{||}$ launch spectrum, acting as a proxy for parametric decay instability (PDI) effects, is included in simulations with the blob model. A synergy between the modified launch spectrum and turbulence-induced refraction results in synthetic hard x-ray profiles that agree with experiment. Lastly, the blob model is used to predict the effect of SOL turbulence on DIII-D high-field side (HFS) LH launch. Assuming low turbulence amplitude in the HFS SOL (~5%), turbulence-induced refraction is predicted to have little effect on current drive efficiency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Observation of synergy between lower hybrid waves at two frequencies in EAST

Synergistic effects between two frequencies of lower hybrid (LH) waves—operating at 2.45 and 4.6 GHz—were observed in experiment on EAST for the first time. At low density (n_e,lin ≈ 2.0 × 10^19m^−3), simultaneous injection of a 65/35 mix of 2.45 GHz/4.6 GHz power achieved an LHCD efficiency that was 25% higher than what should be expected from the linear combination of the two sources. The experiment was interpreted with time-dependent simulations, using the equilibrium and transport solver TRANSP, coupled with the ray-tracing code GENRAY and the Fokker-Planck solver CQL3D. For each discharge, profiles of current and hard x-ray from simulation and measurement agree within uncertainties. An examination of the electron distribution function indicates that the LH synergy is supported by the increased width of the LH resonance plateau in the simultaneous injection case compared to independent injection.

Lower hybrid synergy↗

Characterization of NiC Multilayer Coating Stress for X-Ray Optic Applications

NiC multilayers have been identified as a promising coating design for hard X-ray astrophysical imaging applications enabling bandpass extension beyond the Pt K-edge of approximately 78 keV. However, these coatings are difficult to deposit with low interfacial roughness below a bilayer thickness of about 35 ˚A. Utilizing a DC magnetron sputtering system, NiC multilayer of varying d-spacings are deposited on flat Si wafers and characterized using 8.048 keV X-ray reflectometry measurements. The residual coating stress is also measured using interferometry. We investigate how deposition parameters affect both the coating quality (i.e. surface/interfacial roughness, density, etc.) and the residual coating stress. From these experimental results conducted on flat substrates, we employ FEM and ray trace simulations to determine how NiC multilayer stress could impact the figure, and therefore performance, of full shell X-ray optics.

mulitlayer↗

Surface Method Raytracer

Matlab based experiments that implement multiple ray tracing methods for deterministic radiation transport, including the Surface Method as described in the PhD Thesis, "CONSERVATIVE FIRSTCOLLISION SOURCE TREATMENT FOR RAY EFFECT MITIGATION IN DISCRETE-ORDINATE RADIATION TRANSPORT SOLUTIONS" by Alex Christensen.

Christensen, AlexB↗

Modeling of thermal lensing in side and end-pumped finite solid-state laser rods

An analytical expression for approximating the time-dependent thermal focal length in finite solid state laser rods was derived. The analysis is based on the temperature variation of the material refractive index caused by optical pumping of these rods. Several quantities were found to be relevant to this analysis. These quantities were the specific thermal profiles of the rods, type of optical pumping employed, type of cooling scheme employed (side and end-cooling parameters), and the specific material characteristics of the rods. The Thermal Lensing Model was formulated using the geometric ray tracing approach. The focal lengths are then approximated, by calculating the phase shift in the index of refraction, as the different rays of an incident plane wave are tracked through a lens-like crystal medium. The approach also applies in the case of Gaussian or parabolic pump beams. It is shown that the prediction of thermal focal length is in good quantitative agreement with experimentally obtained data.

Brackett, Vincent G.↗

Quasi-optical beam tracing module development for millimeter-wave high-wavenumber collective scattering on the NSTX-U and EAST tokamaks

A Python3-based beam tracing code utilizing Quasi-Optics has been developed to track both incident and receiving beams in high-k collective millimeter wave scattering systems within magnetic fusion plasmas. In contrast to existing ray tracing codes that solely consider refraction, this beam tracing code incorporates diffraction phenomena, providing a more comprehensive calculation. Here, this enhanced capability allows for a more accurate calculation of the scattering volume and spatial resolution in high-k collective scattering systems, crucial for evaluating system performance and facilitating data analysis. Unlike Geometrical Optics, Quasi-Optics employs the complex eikonal method, representing a Gaussian beam as a collection of coupled rays to accurately preserve diffraction characteristics. The developed code is intended for application in NSTX-Upgrade and EAST high-k beam tracing analyses, targeting frequencies of 693 GHz and 270 GHz, respectively. The high-k system's primary objective is the observation of electron-scale instabilities. Employing a symplectic integrator, the code ensures numerical accuracy, assessed through the conservation of the Hamiltonian. With its precision and efficiency, the code facilitates rapid inter-shot analyses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electrostatic wave propagation and trapping near the magnetic equator

Results of a two-dimensional ray tracing computer code, based on Snell's law, for electrostatic wave propagation in a dipole magnetic field are discussed. A survey of possible ray paths varying a wide range of parameters is conducted for low-harmonic Bernstein modes in a high-density plasma. It is shown that the ray paths exhibit similarity with radial distance and that there exists the possibility of two classes of wave statistics of the equator: a broad emission region extending to about + or - 4 deg and a class of events restricted to the smaller region of 1-2 deg about the magnetic equator. The regulating parameter between these two types of events is the transition energy from the isotropic background electrons to the unstable distribution of superthermals. Ray paths for propagation in the magnetic equatorial plane are considered and an explanation is given for ray focusing in the equatorial plane based on electron gyroradius considerations.

Barbosa, D. D.↗

Angular Broadening: Effects of Nonzero, Spatially Varying Plasma Frequency Between the Source and Observer

Angular broadening of radiation due to scattering by density irregularities is usually described using geometric optics (GO) or the parabolic wave equation (PWE) with the assumptions that the radiation frequency f greatly exceeds the local plasma frequency f(sub p0) or that f(sub p0)/f is constant along the path. These assumptions are inappropriate for many solar system radio phenomena. Here the PWE and GO formalisms are extended to treat angular broadening in plasmas with nonzero, spatially varying ratios, f(sub p0)(z)/f < 1. The new PWE results show that the correlation function, scattered angular spectrum, and other quantities are modified by inclusion of a denominator factor [1 - f(sup 2, sub p0)(z prime)/f prime] inside the path integral over z prime, while the mean-square scattering angle depends on both f(sub p0)(z)/f at the observer and the foregoing factor inside the path integral. The PWE and GO predictions for are identical and involve equivalent assumptions. Previous GO and PWE results are recovered in the limits that f(sub p0)(z prime)/f is constant or zero. The new PWE and GO results will permit more accurate calculation of angular broadening for solar system and astrophysical sources. Moreover and importantly, due to the PWE and GO results for being identical, previous GO analyses of in solar system contexts are essentially correct, except for the neglect of or minor deficiencies in the treatment of nonzero, spatially varying f(sub p0)/f effects. The identical GO and PWE results for and the form of the PWE equation for the correlation function raise questions as to whether diffraction is unimportant for angular broadening (under the usual PWE conditions). Future direct comparisons of the PWE predictions with angular spectra calculated using existing GO ray-tracing codes should answer these questions. Diffraction effects are probably important when the medium and turbulence are not sufficiently homogeneous transverse to the central ray path. The implications are also discussed for studies of the 2-3.5 kHz radiation observed in the outer heliosphere.

Cairns, Iver H.↗

Computing Performance Of An Optical System

Collection of programs performs complicated ray-tracing and diffraction calculations. Optical Surface Analysis Code (OSAC) designed to provide comprehensive analysis of performance of optical system. Analyzes both conventional and x-ray systems. Composed of nine compatible programs executed in a sequence defined by user. Written in FORTRAN IV.

Glenn, P.↗

Grazing incidence relay optics

The necessity to work in the focal plane of the primary mirrors has been one of the factors limiting the utility of grazing incidence telescopes in X-ray astronomy. In connection with the reported investigation, computer ray tracing programs have been used to study the performance of several grazing incidence relay optics (GIRO) systems used together with a large nested solar X-ray telescope. It was found that GIRO magnifiers are useful to map appropriate sized regions of the sun onto available CCD detectors. GIRO collimators can be used together with an X-ray spectrometer to study the X-ray spectrum from very small regions on the sun. Attention is given to the stationary mode, the tracking mode, and the size of GIRO elements. It is found that for a given GIRO size and magnification a use of the diverging system has the advantage of reducing the overall length of the main telescope-GIRO combination. However, the resolution provided by the diverging GIRO may not be as good as that obtained with the corresponding converging GIRO.

Chase, R. C.↗