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 307 records · Page 17

Optimizing X-Ray Optical Prescriptions for Wide-Field Applications

X-ray telescopes with spatial resolution optimized over the field of view (FOV) are of special interest for missions, such as WFXT, focused on moderately deep and deep surveys of the x-ray sky, and for solar x-ray observations. Here we report on the present status of an on-going study of the properties of Wolter I and polynominal grazing incidence designs with a view to gain a deeper insight into their properties and simply the design process. With these goals in mind, we present some results in the complementary topics of (1) properties of Wolter I x-ray optics and polynominal x-ray optic ray tracing. Of crucial importance for the design of wide-field x-ray optics is the optimization criteria. Here we have adopted the minimization of a merit function, M, which measures the spatial resolution averaged over the FOV: M= ((integral of d phi) between the limits of 0 and 2 pi) (integral of d theta theta w(theta) sigma square (theta,phi) between the limits of 0 and theta(sub FOV)) (integral of d phi between the limits of 0 and phi/4) (Integral of d theta theta w(theta) between the limits of 0 and theta(sub FOV) where w(theta(sub 1) is a weighting function and Merit function: sigma-square (theta, phi) = summation of (x,y,z) [ -<(x,y,z)> (exp 2)] is the spatial variance for a point source on the sky at polar and azimuthal off-axis angles (theta,phi).

Elsner, R. F.↗

Tree-based solvers for adaptive mesh refinement code $\scriptsize{FLASH}$ – IV. An X-ray radiation scheme to couple discrete and diffuse X-ray emission sources to the thermochemistry of the interstellar medium

X-ray radiation, in particular radiation between 0.1 and 10 keV, is evident from both point-like sources, such as compact objects and T-Tauri young stellar objects, and extended emission from hot, cooling gas, such as in supernova remnants. The X-ray radiation is absorbed by nearby gas, providing a source of both heating and ionization. While protoplanetary chemistry models now often include X-ray emission from the central young stellar object, simulations of star-forming regions have yet to include X-ray emission coupled to the chemo-dynamical evolution of the gas. We present an extension of the $\scriptsize{TREERAY}$ reverse ray trace algorithm implemented in the flash magnetohydrodynamic code which enables the inclusion of X-ray radiation from 0.1 keV < E γ < 100 keV, dubbed $\scriptsize{XRAYTHESPOT}$. $\scriptsize{XRAYTHESPOT}$ allows for the use of an arbitrary number of bins, minimum and maximum energies, and both temperature-independent and temperature-dependent user-defined cross-sections, along with the ability to include both point and extended diffuse emission and is coupled to the thermochemical evolution. We demonstrate the method with several multibin benchmarks testing the radiation transfer solution and coupling to the thermochemistry. Finally, we show two example star formation science cases for this module: X-ray emission from protostellar accretion irradiating an accretion disc and simulations of molecular clouds with active chemistry, radiation pressure, and protostellar radiation feedback from infrared to X-ray radiation.

79 ASTRONOMY AND ASTROPHYSICS↗

Progress Report on Optimizing X-ray Optical Prescriptions for Wide-Field Applications

We report on the present status of our continuing efforts to develop a method for optimizing wide-field nested x-ray telescope mirror prescriptions. Utilizing extensive Monte-Carlo ray trace simulations, we find an analytic form for the root-mean-square dispersion of rays from a Wolter I optic on the surface of a flat focal plane detector as a function of detector tilt away from the nominal focal plane and detector displacement along the optical axis. The configuration minimizing the ray dispersion from a nested array of Wolter I telescopes is found by solving a linear system of equations for tilt and individual mirror pair displacement. Finally we outline our initial efforts at expanding this method to include higher order polynomial terms in the mirror prescriptions.

Elsner, R. F.↗

Advancing X-ray quantum imaging through Monte-Carlo simulations

Imaging with X-rays poses fundamental limits due to radiation damage of the highly energetic photons. This becomes problematic for sensitive biological systems such as subcellular structures. Lowering the radiation dose, without sacrificing the signal-to-noise ratio, would be desirable for any kind of imaging modalities involving X-rays. To achieve this goal, quantum imaging with entangled X-ray photons constitutes a promising route. Production of biphotons have been demonstrated in the X-ray regime by the process of Spontaneous Parametric Down-Conversion (SPDC). However, compared to SPDC in the regime of visible light, the production rate for X-ray biphotons is extremely low. With the introduction of new high average brightness X-ray sources, such as 4th generation synchrotrons and high repetition rate Free-Electron X-ray Lasers (XFEL), quantum imaging may become practical. We introduce a ray tracing approach using Monte-Carlo sampling, specifically designed for quantum imaging with entangled X-ray photons generated by SPDC. By simulation, the superior image quality of quantum over classical imaging methods is demonstrated using realistic experimental conditions available at high repetition rate XFELs. With these simulations, we can efficiently assist the design of future experiments at beam lines, which can substantially accelerate the advancement of X-ray quantum imaging and reduce costs.

Entangled Photons↗

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↗