Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ray trace”

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 523 records · Page 29

MPACT 4.4 Theory Manual

MPACT is a three-dimensional (3D) full-core neutron transport code capable of calculating subpin power distributions. Calculations are based on the Boltzmann transport equation for neutron fluxes for problems in which the detailed geometrical configuration of fuel components such as the pellet and cladding are explicitly retained. The cross-section data needed for the neutron transport calculation are obtained directly from a multigroup cross section library, which has traditionally been used by lattice physics codes to generate few-group homogenized cross sections for nodal core simulators. Hence, MPACT assumes neither a priori homogenization nor group condensation for the full core spatial solution. The 3D MPACT transport solution can be obtained using the method of characteristics (MOC), which employs discrete ray tracing within each fuel pin. However, for practical reactor applications, the direct application of MOC to 3D core configurations requires an excessive amount of memory and computing time due to the very large number of rays. For practical 3D full-core calculations, MPACT commonly uses an approximate “2D/1D” method that treats the radial (x and y) variables differently from the axial (z) variable. In particular, the radial dependence of the solution is calculated using transport theory, and the axial dependence is calculated using diffusion or P 3 theory. The 2D/1D method requires the core to be divided into a vertical stack of axial slices with a thickness of Δ z ≈ 5–10 cm. Each axial slice is divided radially into coarse spatial cells with boundaries that usually constitute the pin cell boundaries, for which Δ x = Δ y ≈ 1.5 cm. Then, each coarse radial cell (pin cell) is divided into 50–100 fine radial cells, which resolve the angular flux in the fuel, cladding, and moderator regions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Modelling of capillary optics as a focussing hard X-ray concentrator

The development of glass microcapillary arrays by Kumakhov and co-workers in the Soviet Union suggests that it is possible to construct novel devices for concentrating and focussing X-ray beams. Hypothetical focussing concentrators made up of capillary arrays have been configured to operate over a broad range of X-ray energy and its characteristics have been studied with a Monte Carlo ray tracing code. These modeling results indicate that the concentrators would have significant effective area up to 100 keV while their dimensions are similar to conventional grazing incidence telescopes that operate in the soft X-ray band.

Gorenstein, P.↗

A GENERAL RELATIVISTIC MODEL OF ACCRETION DISKS WITH CORONAE SURROUNDING KERR BLACK HOLES

We calculate the structure of a standard accretion disk with a corona surrounding a massive Kerr black hole in the general relativistic frame, in which the corona is assumed to be heated by the reconnection of the strongly buoyant magnetic fields generated in the cold accretion disk. The emergent spectra of accretion disk-corona systems are calculated by using the relativistic ray-tracing method. We propose a new method to calculate the emergent Comptonized spectra from the coronae. The spectra of disk-corona systems with a modified {alpha}-magnetic stress show that both the hard X-ray spectral index and the hard X-ray bolometric correction factor L{sub bol}/L{sub X,2-10keV} increase with the dimensionless mass accretion rate, which is qualitatively consistent with the observations of active galactic nuclei. The fraction of the power dissipated in the corona decreases with increasing black hole spin parameter a, which leads to lower electron temperatures of the coronae for rapidly spinning black holes. The X-ray emission from the coronae surrounding rapidly spinning black holes becomes weak and soft. The ratio of the X-ray luminosity to the optical/UV luminosity increases with the viewing angle, while the spectral shape in the X-ray band is insensitive to the viewing angle. We find that the spectral index in the infrared waveband depends on the mass accretion rate and the black hole spin a, which deviates from the f{sub {nu}}{proportional_to}{nu}{sup 1/3} relation expected by the standard thin disk model.

79 ASTRONOMY AND ASTROPHYSICS↗

A model investigation of the impact of lower hybrid wave scattering angle on current drive profile in EAST and Alcator C-Mod

Lower hybrid current drive (LHCD) is beneficial for developing a steady-state operation scenario in a tokamak. This paper conducts a modeling investigation to identify an optimum rotation angle of the initial lower hybrid perpendicular (to the background magnetic field $\overrightarrow{B}$) wavevector for best matching the experimental RF current profile. It is hypothesized that central RF power deposition widely observed in the present-day LHCD experiments arises from wave scattering by turbulence. In a standard model without considering such interactions, the predicted power deposition profile is generally broad with off-axis peaking, not in agreement with experimental observations. A heuristic approach is adopted by introducing a spectral broadening mechanism by modifying the initial orientation of the perpendicular wavevector. The ray-tracing/Fokker–Planck solver GENRAY/CQL3D is utilized within the python-based π-scope framework. A focus is given to identify the perpendicular wavenumber orientation angle with respect to the magnetic surface normal vector at the initial ray location. Our modeling study shows that rotating the perpendicular wavevector in such a way as to increase the initial poloidal component is effective in reproducing the centrally peaked current profile observed in normal shear plasmas on both EAST and C-Mod. These waves can readily be absorbed to the central plasma, which reduces the sensitivity of the power deposition profile to a slight change of the plasma condition. The same approach is also found to help broaden the off-axis power deposition profile in a reverse-shear EAST plasma, leading to a better agreement with the experiment. The results presented here suggest that spectral modification arising from edge density fluctuations in a tokamak may need to be considered in understanding wave propagation and absorption. A further experimental and theoretical/modeling study is vital as a reverse approach is adopted in this study. Finally, our work suggests that mitigation or control measures are critical for parasitic effects occurring on the first pass in a reactor regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impact of Spacecraft Shielding on Direct Ionization Soft Error Rates for sub-130 nm Technologies

We use ray tracing software to model various levels of spacecraft shielding complexity and energy deposition pulse height analysis to study how it affects the direct ionization soft error rate of microelectronic components in space. The analysis incorporates the galactic cosmic ray background, trapped proton, and solar heavy ion environments as well as the October 1989 and July 2000 solar particle events.

Pellish, Jonathan A.↗

Impact of Spacecraft Shielding on Direct Ionization Soft Error Rates for Sub-130 nm Technologies

We use ray tracing software to model various levels of spacecraft shielding complexity and energy deposition pulse height analysis to study how it affects the direct ionization soft error rate of microelectronic components in space. The analysis incorporates the galactic cosmic ray background, trapped proton, and solar heavy ion environments as well as the October 1989 and July 2000 solar particle events.

Pellish, Jonathan A.↗

The global atmospheric response to low-frequency tropical forcing: Zonally averaged basic states

The extratropical response to localized, low-frequency tropical forcing is examined using a linearized, non-divergent barotropic model on a sphere. Zonal-mean basic states characterized by solid-body rotation or critical latitudes are considered. An analytical analysis based on WKB and ray tracing methods shows that, in contrast to stationary Rossby waves, westward moving, low-frequency Rossby waves can propagate through the tropical easterlies into the extratropics. It is shown analytically that the difference between the stationary and low-frequency ray paths is proportional to the forcing frequency and inversely proportional to the zonal wavenumber cubed. An expression for the disturbance amplitude is derived that shows the ability of the forced waves to maintain their strength well into middle latitudes depends on their meridional wave scale and northward group velocity, both of which are functions of the slowly varying background flow. A local energetics analysis shows that the combination of energy dispersion from the forcing region and energy extraction from the equatorward flank of the midlatitude jet produces disturbances that have the greatest impact on the extratropical circulation. Under the assumption that the forcing amplitude is independent of frequency, this impact is largest when the tropical forcing period is in the range 10-20 days.

Li, Long↗

Wave-current interaction near the Gulf Stream during the surface wave dynamics experiment

This paper presents a case study on the wave-current interaction near the local curvature of a Gulf Stream meander. The wave data were obtained from in situ measurements by a pitch-roll discus buoy during the Surface Wave Dynamics Experiment (SWADE) conducted off Wallops Island, Virginia, from October 1990 to March 1991. Owing to the advection of the Gulf Stream by the semidiurnal tide, the discus buoy was alternately located outside and inside the Gulf Stream. The directional wave measurements from the buoy show the changes in wave direction, wave energy, and directional spreading when waves encountered the current in the Gulf Stream meanders. A wave refraction model, using the ray-tracing method with an estimated Gulf Stream velocity field and meandering condition, was used to simulate wave refraction patterns and to estimate wave parameters at relative locations corresponding to buoy measurements. The numerical simulation shows that a focusing zone of wave rays was formed near the boundary and behind the crest of a simulated Gulf Stream meander. The focusing of wave rays causes changes in wave direction, increases in wave energy, and decreases in wave directional spreading, which are in good agreement with the results from the buoy measurements.

Wang, David W.↗

Solid, structured composite neutron detectors with high dynamic range capability

Neutron detectors are essential across disciplines such as fundamental science, nuclear security, safeguards, and civilian applications. While 3 He-filled gas proportional counters have long been revered for their efficacy in detecting thermal neutrons and praised for their efficiency, neutron/gamma discrimination, and stability, the scarcity of 3He has spurred a search for alternatives. Here, we explore a solid structured scintillating particle composite (SPC) consisting of 6 Li-containing scintillating glass particles within an acrylic matrix as a neutron detector for high dynamic range applications. We show for the first time that an SPC neutron detector can boast an intrinsic detection efficiency of 0.261% for pure 252 Cf fission neutrons and an overall neutron detection efficiency of (0.546 ± 0.003)% at the Neutron Free-in-Air facility while being able to function in an intense gamma-ray environment. We also show that the SPC neutron detector supports fast neutron capture times and enables a dual-readout scheme that extends the detector dynamic range to high incident neutron fluxes. A scalable fabrication process allows for tailoring the SPC detector properties to the requirements of specific applications. Good agreement is found between the experimental results taken with a National Institute of Standards and Technology traceable 252 Cf source and the coupled MCNP6 and optical-ray-tracing simulations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Angular distribution of polarized spontaneous emissions and its effect on light extraction behavior in InGaN-based light emitting diodes

Angular intensity distributions of differently polarized light sources in multiple quantum wells (MQWs) and their effects on extraction behavior of spontaneous emission from light emitting diode (LED) chips have been studied. Theoretical calculation based on k·p approximation, ray tracing simulation and angular electroluminescence measurement were applied in this work. It is found that the electron-hole recombination in the InGaN MQWs produces a spherical distribution of an s-polarized source and a dumbbell-shaped p-polarized source. Light rays from different polarized sources experience different extraction processes, determining the polarization degree of electro-luminescence and extraction efficiency of LEDs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Scattering of light by crystals - A modified Kirchhoff approximation

A modified Kirchhoff approximation (MKA) is developed for the scattering of light by randomly oriented crystals. The reflected and transmitted near fields are calculated from ray tracing; the corresponding far fields are then obtained via the vector Kirchhoff integral. On the shadow side of the particle, an additional near field exactly cancels the incident field and causes the forward diffraction. MKA contains a particle size dependence, which is not included in ray optics treatments, and satisfactory results can be obtained for size parameters larger than ten. The scattering phase functions and degrees of linear polarization are calculated for some hexagonal and cubic water ice crystals using MKA. The Kirchhoff approximation for particles other than crystals is discussed, and attention is paid to the backscattering enhancement due to the cyclic passage of internally or multiply externally reflected electromagnetic waves.

Muinonen, Karri↗

Design and development of conical X-ray imaging mirrors

The basic principles for the design of conical optics, and the development of an arc-minute quality conical mirror are described. Computer ray tracing data for over 100 conical mirror designs are analyzed in order to study the intrinsic capabilities of conical mirrors; the intrinsic spatial resolution and approximate filling factor are considered. The procedures for the fabrication of a conical mirror for the Broad-Band X-ray Telescope experiment to be flown on the Shuttle are examined. A f/10 conical mirror with 3.8 m focal length and a filling factor of about 60 percent, consisting of 18 nested reflector pairs has been developed. The difference between the intrinsic and measured half-power radii is investigated. Methods for improving the spatial resolution of conical mirrors are being researched. The advantages and disadvantages of conical mirrors are discussed. It is noted that conical mirrors are applicable for studies in which a broad bandpass, moderate spatial resolution, and a small number of mirrors are required.

Petre, R.↗

Low-Cost Heliostat for High-Flux Small-Area Receivers (Final Technical Report)

This project analyzed a two-stage heliostat concept consisting of a tracking stage and a concentrating stage. The tracking stage uses mirrors mounted on a common drive that move to track the sun. The concentrating stage consists of stationary mirrors that each have a unique angle to direct rays towards a small-area, high-flux, point-focused receiver. By splitting the collection and concentrating process into two stages, multiple small, inexpensive mirrors can share a structure and be controlled by a single drive in the tracking stage. The project effort developed modeling techniques that were specifically relevant to this two-stage heliostat concept. Both field-level and unit-level models were developed. The field-level model does not explicitly consider unit-level losses which are predicted by the unit-level model and then integrated into the field-level model through a correlation referred to as an efficiency modifier. This approach is referred to as the two-model approach; the development and demonstration of this two-model approach for a multi-stage heliostat technology is a key outcome of this work. The field-level model is used to design a field that hits a specific design day power given a set of heliostat design parameters. An oversized field is simulated and then heliostat units are removed based on their annual energy production in order to generate the highest performing field. The field reduction procedure fits a smooth curve fit to annual energy production as a function of position in the field which has the effect of reducing the noise that is otherwise caused by the Monte Carlo ray tracing technique. This approach is referred to as the annual energy fit method and substantially reduces computational run time for a given field level modeling accuracy. The annual energy fit approach enables the selection of a properly sized, high-performing field using orders of magnitude fewer rays than would otherwise be possible and the development of this approach is a second key outcome of this work. These models are used within a genetic optimization algorithm in order to optimize the geometric parameters associated with a heliostat in order to achieve the lowest cost per unit of collected design day power. The cost modeling that underlies the optimization is a simple, scaling type analysis backed up by a much more detailed Design for Manufacture and Assembly (DFMA) analysis. Although the figure of merit used for optimization was not cost per mirror area, this metric is reasonable to use as a means of comparison. The optimally designed 500 kW design has a tracking mirror specific cost of $181.85/m 2 , which is significantly larger than the target value and also larger than the current state of the art. The cost of the torque-tube type linkages contributed substantially to the overall cost. Based on this observation, potentially attractive alternative design configuration utilizing a capstan type actuation system should be investigated. Finally, NREL compared the performance of the two-stage heliostat to the performance of a focused and different sized flat conventional heliostats and showed that, as expected, additional losses versus the convention heliostat caused by a worse cosine efficiency, two stages of reflection, and interstage interactions. The two-stage heliostat requires around 75% more reflective area than a flat 1x1 meter conventional heliostat (similar to a focused heliostat) and 40% more than a flat 2x2 meter conventional heliostat.

14 SOLAR ENERGY↗

Designing Aberration-Corrected Solid Unstable Resonators

In improved method of designing solid unstable resonator of laser diode, shapes of mirrors calculated to yield specified mode. Ray tracing used to compute shape of initially unspecified end mirror, given shape of initially specified end mirror and specified output mode. No need to accept aberrations or suboptimal circular shapes, to make iterative design computations in effort to converge on desired mode, or to assume paraxiality of rays: Angles between rays and optical axis large, cross sections of surfaces noncircular, and computed shape of end mirror exact. End mirror corrects for all aberrations.

Lang, Robert J.↗

Performance simulation of the soft gamma-ray concentrator

The soft gamma-ray concentrator is a telescope mission concept utilizing a suitable arrangement of bent multilayer structures of alternating low- and high-density materials. This lens is able to channel gamma-ray photons via total external reflection and concentrate the incident radiation to a point. The channeling technique offers the potential for concentrating gamma rays with focal lengths <10 m and energies >100 keV, beyond the reach of current grazing-incidence hard x-ray mirrors. For the performance estimation of such an instrument, we have developed a flexible set of computer modeling tools to compute the optical properties of multilayer structures, predict the channeling efficiency for a given multilayer configuration, and aid in the optimization of potential gamma-ray concentrator-based telescope designs. This modeling includes the multilayer optical properties calculated by the IMD software, the ray tracing using an IDL code, and the focal plane detector simulation by MEGAlib. We illustrate the potential of this approach by presenting simulated astronomical observations from a balloon-borne platform. The final result, including simulated effective area, instrument sensitivity, and polarization performance, shows that the gamma-ray concentrator will provide greatly increased sensitivity for next-generation soft gamma-ray missions with modest cost and complexity.

47 OTHER INSTRUMENTATION↗

Parallel 3-D method of characteristics in MPACT

A new parallel 3-D MOC kernel has been developed and implemented in MPACT which makes use of the modular ray tracing technique to reduce computational requirements and to facilitate parallel decomposition. The parallel model makes use of both distributed and shared memory parallelism which are implemented with the MPI and OpenMP standards, respectively. The kernel is capable of parallel decomposition of problems in space, angle, and by characteristic rays up to 0(104) processors. Initial verification of the parallel 3-D MOC kernel was performed using the Takeda 3-D transport benchmark problems. The eigenvalues computed by MPACT are within the statistical uncertainty of the benchmark reference and agree well with the averages of other participants. The MPACT k{sub eff} differs from the benchmark results for rodded and un-rodded cases by 11 and -40 pcm, respectively. The calculations were performed for various numbers of processors and parallel decompositions up to 15625 processors; all producing the same result at convergence. The parallel efficiency of the worst case was 60%, while very good efficiency (>95%) was observed for cases using 500 processors. The overall run time for the 500 processor case was 231 seconds and 19 seconds for the case with 15625 processors. Ongoing work is focused on developing theoretical performance models and the implementation of acceleration techniques to minimize the number of iterations to converge. (authors)

97 MATHEMATICS AND COMPUTING↗

SU-E-I-01: A Fast, Analytical Pencil Beam Based Method for First Order X-Ray Scatter Estimation of Kilovoltage Cone Beam X-Rays

Purpose: To analytically estimate first-order x-ray scatter for kV cone beam x-ray imaging with high computational efficiency. Methods: In calculating first-order scatter using the Klein-Nishina formula, we found that by integrating the point-to-point scatter along an interaction line, a “pencil-beam” scatter kernel (BSK) can be approximated to a quartic expression when the imaging field is small. This BSK model for monoenergetic, 100keV x-rays has been verified on homogeneous cube and cylinder water phantoms by comparing with the exact implementation of KN formula. For heterogeneous medium, the water-equivalent length of a BSK was acquired with an improved Siddon's ray-tracing algorithm, which was also used in calculating pre- and post- scattering attenuation. To include the electron binding effect for scattering of low-kV photons, the mean corresponding scattering angle is determined from the effective point of scattered photons of a BSK. The behavior of polyenergetic x-rays was also investigated for 120kV x-rays incident to a sandwiched infinite heterogeneous slab phantom, with the electron binding effect incorporated. Exact computation and Monte Carlo simulations were performed for comparisons, using the EGSnrc code package. Results: By reducing the 3D volumetric target (o(n{sup 3})) to 2D pencil-beams (o(n{sup 2})), the computation expense can be generally lowered by n times, which our experience verifies. The scatter distribution on a flat detector shows high agreement between the analytic BSK model and exact calculations. The pixel-to-pixel differences are within (-2%, 2%) for the homogeneous cube and cylinder phantoms and within (0, 6%) for the heterogeneous slab phantom. However, the Monte Carlo simulation shows increased deviation of the BSK model toward detector periphery. Conclusion: The proposed BSK model, accommodating polyenergetic x-rays and electron binding effect at low kV, shows great potential in efficiently estimating the first-order scatter from small imaging fields. We are investigating more thoroughly to improve performance and explore applications.

07 ISOTOPE AND RADIATION SOURCES↗

Predicted performance of a tangential viewing hard x-ray camera for the DIII-D high field side lower hybrid current drive experiment

High field side launch of lower hybrid current drive (LHCD) has improved accessibility and penetration over low field side launch on DIII-D. Simulations predict single pass absorption under a wide range of plasma conditions. Hard x-ray (HXR) measurement of LHCD generated fast electron bremsstrahlung (50–250 keV) will validate wave propagation and absorption. Emissivity profiles are recovered from one-dimensional inversion of HXR brightness to determine LH damping location, fast electron slowing down time, and some indication of the fast electron energy. The camera will be implemented by populating 32 tangential sightlines of the existing Gamma Ray Imager with Kromek SPEARTM Cadmium Zinc Telluride (CZT) detectors sensitive to 10–1000 keV photons with 10 keV energy resolution. Expected count rates allow for <0.5 ms time resolution. Pulses are processed using 50 ns shaping time Cremat CR-200 Gaussian shaping modules and are digitized by 25 MHz D-TACQ ACQ216 digitizers. The performance of the HXR camera is evaluated by comparing predicted fast electron density profiles and inverted synthetic brightnesses obtained from the ray-tracing/Fokker–Planck codes GENRAY/CQL3D. Inversions closely matched predicted fast electron profiles for a range of experimental parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗