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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Refraction effects on the Galileo probe telemetry carrier frequency

As the Galileo probe relay radio link (RRL) signal propagates outward through the Jovian atmosphere, the atmosphere will manifest itself in two ways. First, the geometric path length of the signal is increased, resulting in a small change of the RRL signal departure angle from the proble (transmitter). Secondly, the velocity of the signal is decreased. For a spherical, static atmosphere with a known profile of refractivity versus altitude the effects of refraction on the RRL frequency can be found using a variation of standard ray-tracing techniques, whereby the ray departure angle is found by an iterative process. From the dispersive characteristics of a mixture of hydrogen and helium with trace amounts of methane and ammonia a simple model of the Jovian atmosphere is constructed assuming spherical symmetry and uniform mixing. The contribution to the RRL Doppler frequency arising from refraction is calculated, and its effect on the Doppler wind measurements is discussed.

Atkinson, D. H.↗
NASA NTRS2014

The Prediction and Analysis of Jet Flows and Scattered Turbulent Mixing Noise About Flight Vehicle Airframes

Jet flows interacting with nearby surfaces exhibit a complex behavior in which acoustic and aerodynamic characteristics are altered. The physical understanding and prediction of these characteristics are essential to designing future low noise aircraft. A new approach is created for predicting scattered jet mixing noise that utilizes an acoustic analogy and steady Reynolds-averaged Navier-Stokes solutions. A tailored Green's function accounts for the propagation of mixing noise about the air-frame and is calculated numerically using a newly developed ray tracing method. The steady aerodynamic statistics, associated unsteady sound source, and acoustic intensity are examined as jet conditions are varied about a large at plate. A non-dimensional number is proposed to estimate the effect of the aerodynamic noise source relative to jet operating condition and airframe position. The steady Reynolds-averaged Navier-Stokes solutions, acoustic analogy, tailored Green's function, non- dimensional number, and predicted noise are validated with a wide variety of measurements. The combination of the developed theory, ray tracing method, and careful implementation in a stand-alone computer program result in an approach that is more first principles oriented than alternatives, computationally efficient, and captures the relevant physics of fluid-structure interaction.

Miller, Steven A.↗
NASA NTRS2014

The Prediction and Analysis of Jet Flows and Scattered Turbulent Mixing Noise about Flight Vehicle Airframes

Jet flows interacting with nearby surfaces exhibit a complex behavior in which acoustic and aerodynamic characteristics are altered. The physical understanding and prediction of these characteristics are essential to designing future low noise aircraft. A new approach is created for predicting scattered jet mixing noise that utilizes an acoustic analogy and steady Reynolds-averaged Navier-Stokes solutions. A tailored Green's function accounts for the propagation of mixing noise about the airframe and is calculated numerically using a newly developed ray tracing method. The steady aerodynamic statistics, associated unsteady sound source, and acoustic intensity are examined as jet conditions are varied about a large flat plate. A non-dimensional number is proposed to estimate the effect of the aerodynamic noise source relative to jet operating condition and airframe position.The steady Reynolds-averaged Navier-Stokes solutions, acoustic analogy, tailored Green's function, non-dimensional number, and predicted noise are validated with a wide variety of measurements. The combination of the developed theory, ray tracing method, and careful implementation in a stand-alone computer program result in an approach that is more first principles oriented than alternatives, computationally efficient, and captures the relevant physics of fluid-structure interaction.

Miller, Steven A. E.↗
NASA NTRS1992

Funnel-shaped, low-frequency equatorial waves

Funnel-shaped, low-frequency radiation, as observed in frequency time spectrograms, is frequently found at the earth's magnetic equator which extends from the proton-cyclotron frequency up to the lower hybrid frequency. Ray-tracing calculations can qualitatively reproduce the observed frequency-time characteristics of these emissions if the waves are propagating in the fast magnetosonic mode starting with wave normal angles of about 88 deg at the magnetic equator. The funnel-shaped emissions are consistent with generation by protons with a ring-type velocity space distribution. A ring-shaped region of positive slope in the velocity space density distribution of protons is observed near the Alfven velocity, indicating that the ring protons strongly interact with the waves. Ray-tracing calculations show that for similar equatorial wave normal angles lower-frequency fast magnetosonic waves are more closely confined to the magnetic equator than higher-frequency fast magnetosonic waves. For waves refracted back toward the equator at similar magnetic latitudes, the lower-frequency waves experience stronger damping in the vicinity of the equator than higher-frequency waves. Also, wave growth is restricted to higher frequencies at larger magnetic latitudes. Wave damping at the equator and wave growth off the equator favors equatorial wave normal angle distributions which lead to the funnel-shaped frequency time characteristic.

Boardsen, S. A.↗
DOE OSTI2024

Simple strategy for the simulation of axially symmetric large-area metasurfaces

Metalenses are composed of nanostructures for focusing light and have been widely explored in many exciting applications. However, their expanding dimensions pose simulation challenges. We propose a method to simulate metalenses in a timely manner using vectorial wave and ray tracing models. We sample the metalens’s radial phase gradient and locally approximate the phase profile by a linear phase response. Each sampling point is modeled as a binary blazed grating, employing the chosen nanostructure, to build a transfer function set. The metalens transmission or reflection is then obtained by applying the corresponding transfer function to the incoming field on the regions surrounding each sampling point. Fourier optics is used to calculate the scattered fields under arbitrary illumination for the vectorial wave method, and a Monte Carlo algorithm is used in the ray tracing formalism. We validated our method against finite-difference time domain simulations at 632 nm, and we were able to simulate metalenses larger than 3000 wavelengths in diameter on a personal computer.

Martins, Augusto (ORCID:0000000295546481)↗
NASA NTRS1993

Spectral radiation analyses of the GOES solar illuminated hexagonal cell scan mirror back

A ray tracing analytical tool has been developed for the simulation of spectral radiation exchange in complex systems. Algorithms are used to account for heat source spectral energy, surface directional radiation properties, and surface spectral absorptivity properties. This tool has been used to calculate the effective solar absorptivity of the geostationary operational environmental satellites (GOES) scan mirror in the calibration position. The development and design of Sounder and Imager instruments on board GOES is reviewed and the problem of calculating the effective solar absorptivity associated with the GOES hexagonal cell configuration is presented. The analytical methodology based on the Monte Carlo ray tracing technique is described and results are presented and verified by experimental measurements for selected solar incidence angles.

Fantano, Louis G.↗
DOE OSTI2014

Application of the TITAN Hybrid Deterministic Transport Code to Medical Physics - Paper 118

The TITAN code is a 3-dimensional parallel deterministic transport code that solves the linear Boltzmann equation. The TITAN code has hybrid algorithms that are ideal for medical physics applications: i) SN and characteristics method and ii) SN with a fictitious quadrature set and ray-tracing. These hybrid algorithms have been used in the simulation of computed tomography (CT) and single photon emission computed tomography (SPECT) In this paper, the hybrid algorithms developed for medical physics applications in the TITAN code are discussed. First, the TITAN code's block-oriented approach for utilizing SN and characteristic solvers for CT is presented and results for a simple test case are given. Second, the paper describes a hybrid formulation for SPECT with an SN solver in the phantom, a fictitious quadrature set for directions from the phantom surface passing through the collimator, and a simplified ray-tracing algorithm for particle transport to the detector. Finally, ongoing work is described and concluding remarks are given. The TITAN code's hybrid algorithms allow it to efficiently solve problems with low-scatter regions, which are frequently found in medical physics applications. The TITAN code has been used to solve a simplified CT problem and generate projection images for SPECT simulations. The WCOS technique was developed to improve the TITAN code's collimator representation. For these applications, the TITAN code was shown to require significantly less computation time than MCNP5. Work is currently underway to implement an iterative image reconstruction algorithm in the TITAN code. A script that calls the TITAN code to complete the forward projection step in SPECT has already been used to reconstruction a 2-dimensional phantom. The ability of TITAN to fully model particle transport in the patient in short computation times is expected to be a great advantage to iterative reconstruction. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗
NASA NTRS2016

Monochromatic Measurements of the JPSS-1 VIIRS Polarization Sensitivity

Polarization sensitivity is a critical property that must be characterized for spaceborne remote sensing instruments designed to measure reflected solar radiation. Broadband testing of the first Joint Polar-orbiting Satellite System (JPSS-1) Visible Infrared Imaging Radiometer Suite (VIIRS) showed unexpectedly large polarization sensitivities for the bluest bands on VIIRS (centered between 400 and 600 nm). Subsequent ray trace modeling indicated that large diattenuation on the edges of the bandpass for these spectral bands was the driver behind these large sensitivities. Additional testing using the National Institute of Standards and Technologies Traveling Spectral Irradiance and Radiance Responsivity Calibrations Using Uniform Sources was added to the test program to verify and enhance the model. The testing was limited in scope to two spectral bands at two scan angles; nonetheless, this additional testing provided valuable insight into the polarization sensitivity. Analysis has shown that the derived diattenuation agreed with the broadband measurements to within an absolute difference of about0.4 and that the ray trace model reproduced the general features of the measured data. Additionally, by deriving the spectral responsivity, the linear diattenuation is shown to be explicitly dependent on the changes in bandwidth with polarization state.

Monochromatic↗
DOE OSTI2016

Sensitivity analysis of high resolution gamma-ray detection for safeguards monitoring at natural uranium conversion facilities

Under the policies proposed by recent International Atomic Energy Agency (IAEA) circulars and policy papers, implementation of safeguards exists when any purified aqueous uranium solution or uranium oxides suitable for isotopic enrichment or fuel fabrication exists. Under IAEA Policy Paper 18, the starting point for nuclear material under safeguards was reinterpreted, suggesting that purified uranium compounds should be subject to safeguards procedures no later than the first point in the conversion process. In response to this technical need, a combination of simulation models and experimental measurements were employed in previous work to develop and validate gamma-ray nondestructive assay monitoring systems in a natural uranium conversion plant (NUCP). In particular, uranyl nitrate (UO 2 (NO 3 ) 2 ) solution exiting solvent extraction was identified as a key measurement point (KMP). Passive nondestructive assay techniques using high resolution gamma-ray spectroscopy were evaluated to determine their viability as a technical means for drawing safeguards conclusions at NUCPs, and if the IAEA detection requirements of 1 significant quantity (SQ) can be met in a timely manner. Building upon the aforementioned previous validation work on detector sensitivity to varying concentrations of uranyl nitrate via a series of dilution measurements, this work investigates detector response parameter sensitivities to gamma-ray signatures of uranyl nitrate. The full energy peak efficiency of a detection system is dependent upon the sample, geometry, absorption, and intrinsic efficiency parameters. Perturbation of these parameters translates into corresponding variations of the 185.7 keV peak area of the 235 U in uranyl nitrate. Such perturbations in the assayed signature impact the quality or versatility of the safeguards conclusions drawn. Given the potentially high throughput of uranyl nitrate in NUCPs, the ability to assay 1 SQ of material requires uncertainty «1%. Accounting for material self-shielding properties, pipe thickness, and source-detector orientation is instrumental in determining the robustness of gamma-ray detection in the process monitoring of uranyl nitrate in NUCPs. Monte Carlo models and ray-tracing models were employed to determine the sensitivity of the detected 185.7 keV photon to self-shielding properties, pipe thickness, and source-detector geometry. Considering the implementation of the detection of 1 SQ, diversion of 1 SQ becomes essentially undetectable given the systematic uncertainty, in addition to considerations such as propagating uncertainties due to pipe offset/position, as well as minor variations in pipe thickness. Consequently, pipe thickness was the most sensitive variable in affecting full energy efficiency of the 185.7 keV signature peak with up to 8% variation in efficiency for ±0.5 mm changes in Schedule 40 304L stainless steel piping. Furthermore, computation of the attenuation correction factor of the uranyl nitrate solution [CF(AT) (i.e. εsample)] using Parker's method using with the approximation for the geometrical factor κ≈π/4 was validated through experimental, Monte Carlo and ray-tracing calculations for a uranyl nitrate filled transfer pipe segment. Furthermore, quantifying sensitivity in detector position, as well as voiding effects due to bubbly flow or laminar flow with an air gap in the uranyl nitrate becomes increasingly important as considerations from (static) design-scale measurements translate into (dynamic) field operations tests.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1993

Propagation characteristics of Pc 3 compressional waves generated at the dayside magnetopause

New, 3D ray tracing of Pc 3 compressional waves from the magnetosheath reveals that the magnetosphere can present a major propagation barrier to the penetration of these waves to the plasmasphere. This barrier is the ion-ion cutoff between the He(+) and O(+) gyroresonances. As a result of the frequency-dependent location of this cutoff, the magnetosphere behaves like a filter for Pc 3 compressional waves, and only low-frequency components of Pc 3 compressional waves can penetrate to inner magnetosphere. Results are in agreement with previous satellite observations. This 'filter action' strongly depends on the relative concentration of He(+) and O(+) and is therefore sensitive to solar and magnetic activity. Ray-tracing results are based on a cold plasma dispersion relation, a semiempirical model of plasma density, and the Mead-Fairfield (1975) magnetic field model.

Zhang, X.↗
DOE OSTI2017

Full-sky Gravitational Lensing Simulation for Large-area Galaxy Surveys and Cosmic Microwave Background Experiments

We present 108 full-sky gravitational lensing simulation data sets generated by performing multiple-lens plane ray-tracing through high-resolution cosmological N-body simulations. The data sets include full-sky convergence and shear maps from redshifts z = 0.05 to 5.3 at intervals of $150\,{h}^{-1}\mathrm{Mpc}$ comoving radial distance (corresponding to a redshift interval of ${\rm{\Delta }}z\simeq 0.05$ at the nearby universe), enabling the construction of a mock shear catalog for an arbitrary source distribution up to z = 5.3. The dark matter halos are identified from the same N-body simulations with enough mass resolution to resolve the host halos of the Sloan Digital Sky Survey (SDSS) CMASS and luminous red galaxies (LRGs). Angular positions and redshifts of the halos are provided by a ray-tracing calculation, enabling the creation of a mock halo catalog to be used for galaxy–galaxy and cluster–galaxy lensing. The simulation also yields maps of gravitational lensing deflections for a source redshift at the last scattering surface, and we provide 108 realizations of lensed cosmic microwave background (CMB) maps in which the post-Born corrections caused by multiple light scattering are included. We present basic statistics of the simulation data, including the angular power spectra of cosmic shear, CMB temperature and polarization anisotropies, galaxy–galaxy lensing signals for halos, and their covariances. The angular power spectra of the cosmic shear and CMB anisotropies agree with theoretical predictions within 5% up to ${\ell }=3000$ (or at an angular scale $\theta \gt 0.5$ arcmin). The simulation data sets are generated primarily for the ongoing Subaru Hyper Suprime-Cam survey, but are freely available for download at http://cosmo.phys.hirosaki-u.ac.jp/takahasi/allsky_raytracing/.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Metaplectic geometrical optics for ray-based modeling of caustics: Theory and algorithms

The optimization of radio frequency-wave (RF) systems for fusion experiments is often performed using ray-tracing codes, which rely on the geometrical-optics (GO) approximation. However, GO fails at caustics such as cutoffs and focal points, erroneously predicting the wave intensity to be infinite. This is a critical shortcoming of GO, since the caustic wave intensity is often the quantity of interest, e.g., RF heating. Full-wave modeling can be used instead, but the computational cost limits the speed at which such optimizations can be performed. Here, we have developed a less expensive alternative called metaplectic geometrical optics (MGO). Instead of evolving waves in the usual x (coordinate) or k (spectral) representation, MGO uses a mixed X$\equiv$Ax+Bk representation. By continuously adjusting the matrix coefficients A and B along the rays, one can ensure that GO remains valid in the X coordinates without caustic singularities. The caustic-free result is then mapped back onto the original x space using metaplectic transforms. Here, we overview the MGO theory and review algorithms that will aid the development of an MGO-based ray-tracing code. We show how using orthosymplectic transformations leads to considerable simplifications compared to previously published MGO formulas. We also prove explicitly that MGO exactly reproduces standard GO when evaluated far from caustics (an important property that until now has only been inferred from numerical simulations), and we relate MGO to other semiclassical caustic-removal schemes published in the literature. Finally this discussion is then augmented by an explicit comparison of the computed spectrum for a wave bounded between two cutoffs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Quantifying electron cyclotron power deposition broadening in DIII-D and the potential consequences for the ITER EC system

Abstract The injection of electron cyclotron (EC) waves fulfills a number of important tasks in nuclear fusion devices for which detailed knowledge of the spatial power deposition profile is critical. This deposition profile is commonly determined using forward models such as beam or ray tracing. Recent numerical and experimental studies have shown that small-angle scattering of the EC beam as it passes through the turbulent plasma edge can cause significant broadening of the effective deposition profile, leading to considerable underestimation of the deposition width by forward methods. However, traditional inverse methods to determine the deposition profile from measurements overestimate the deposition profile width due to transport broadening. In this work, we implement three novel methods to resolve the EC power deposition profile from measurements that counteract transport broadening by simultaneously resolving transport and power deposition. We validate their assumptions and compare the results from these methods to the traditional break-in-slope method as well as to the TORAY ray-tracing code in a set of DIII-D discharges spanning five different confinement modes. We show that the four different inverse methods, novel and established, paint a consistent picture of deposition broadening. Specifically, we show that the measured power deposition profile is between 1.6 and 3.6 times wider than the TORAY profiles. Moreover, we show the considerable consequences that this level of broadening can have for ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Interpreting ~1 Hz magnetic compressional waves in Mercury's inner magnetosphere in terms of propagating ion‐Bernstein waves

Abstract We show that ~1 Hz magnetic compressional waves observed in Mercury's inner magnetosphere could be interpreted as ion‐Bernstein waves in a moderate proton beta ~0.1 plasma. An observation of a proton distribution with a large planetary loss cone is presented, and we show that this type of distribution is highly unstable to the generation of ion‐Bernstein waves with low magnetic compression. Ray tracing shows that as these waves propagate back and forth about the magnetic equator; they cycle between a state of low and high magnetic compression. The group velocity decreases during the high‐compression state leading to a pileup of compressional wave energy, which could explain the observed dominance of the highly compressional waves. This bimodal nature is due to the complexity of the index of refraction surface in a warm plasma whose upper branch has high growth rate with low compression, and its lower branch has low growth/damping rate with strong compression. Two different cycles are found: one where the compression maximum occurs at the magnetic equator and one where the compression maximum straddles the magnetic equator. The later cycle could explain observations where the maximum in compression straddles the equator. Ray tracing shows that this mode is confined within ±12° magnetic latitude which can account for the bulk of the observations. We show that the Doppler shift can account for the difference between the observed and model wave frequency, if the wave vector direction is in opposition to the plasma flow direction. We note that the Wentzel‐Kramers‐Brillouin approximation breaks down during the pileup of compressional energy and that a study involving full wave solutions is required.

Boardsen, S. A.↗
DOE OSTI2012

Zgoubi Users' Guide

The initial version of zgoubi, dedicated to ray-tracing in magnetic fields, was developed by D. Garreta and J.C. Faivre at CEN-Saclay in the early 1970’s. It was perfected for the purpose of studying the four spectrometers SPES I, II, III, IV at the Laboratoire National SATURNE (CEA-Saclay, France), and SPEG at GANIL (Caen, France). It is being used since long in several national and foreign laboratories. The first manual was in French. Accounting for many developments and improvements, and in or der to facilitate access to the program an English version of the manual was written at TRIUMF with the assistance of J. Doornbos. P. Stewart prepared the manuscript for publication. An updating of the latter was necessary for accompanying the third version of the code which included developments regarding spin tracking and ray-tracing in combined electric and magnetic fields; this was done with the help of D. Bunel (SATURNE Laboratory, Saclay) for the preparation of the document and lead to the third release. In the mid-1990s, the computation of synchrotron radiation electromagnetic impulse and spectra was introduced for the purpose of studying interference effects at the LEP synchrotron radiation based diagnostic mini-wiggler. In the mean time, several new optical elements were added, such as electro-magnetic and other electrostatic lenses. Used since several years for special studies in periodic machines (e.g., SATURNE at Saclay, COSY at Julich, LEP and LHC at CERN, Neutrino Factory rings), zgoubi has also undergone extensive developments regarding storage ring related features. These developments of zgoubi have strongly benefited of the environment of the Groupe Théorie, Laboratoire National SATURNE, CEA/DSM-Saclay, in the years 1985-1995. The graphic interface to zgoubi (addressed in Part D) has also undergone concomitant extensive devel opments, which make it a performing tool for the post-processing of zgoubi outputs. The Users’ Guide is intended to describe the contents of the most recent version of zgoubi, which is far from being a “finished product”. Françios Méot

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Pre-commissioning of the wide-dynamic range halo monitor to be installed in the Fermilab MI-8 line

The beam halo can contribute to beam losses in accelerators and is very difficult to measure. With an increase in beam intensity following the PIP-II upgrade at Fermilab, the beam losses are expected to be higher with some coming from beam halo. Therefore, it is important to measure the sources of beam halo to minimize the beam losses. A modified Halo Monitor developed by J-PARC will be installed in Fermilab MI-8 transfer line to measure the beam halo. In this paper, an update on the beam profile monitor fabrication is covered. The updates include the location selection for the Halo Monitor in the MI-8 transfer line, shielding options for instrumentation, initial testing of equipment, and ray tracing simulations for the Offner optics and the targets used in the monitor.

Babacan, Betiay [IIT, Chicago; Fermilab]↗
DOE OSTI2013

EFFECT OF MASKED REGIONS ON WEAK-LENSING STATISTICS

Sky masking is unavoidable in wide-field weak-lensing observations. We study how masks affect the measurement of statistics of matter distribution probed by weak gravitational lensing. We first use 1000 cosmological ray-tracing simulations to examine in detail the impact of masked regions on the weak-lensing Minkowski Functionals (MFs). We consider actual sky masks used for a Subaru Suprime-Cam imaging survey. The masks increase the variance of the convergence field and the expected values of the MFs are biased. The bias then compromises the non-Gaussian signals induced by the gravitational growth of structure. We then explore how masks affect cosmological parameter estimation. We calculate the cumulative signal-to-noise ratio (S/N) for masked maps to study the information content of lensing MFs. We show that the degradation of S/N for masked maps is mainly determined by the effective survey area. We also perform simple {chi}{sup 2} analysis to show the impact of lensing MF bias due to masked regions. Finally, we compare ray-tracing simulations with data from a Subaru 2 deg{sup 2} survey in order to address if the observed lensing MFs are consistent with those of the standard cosmology. The resulting {chi}{sup 2}/n{sub dof} = 29.6/30 for three combined MFs, obtained with the mask effects taken into account, suggests that the observational data are indeed consistent with the standard {Lambda}CDM model. We conclude that the lensing MFs are a powerful probe of cosmology only if mask effects are correctly taken into account.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Alpha channeling with high-field launch of lower hybrid waves

Although lower hybrid waves are effective at driving currents in present-day tokamaks, they are expected to interact strongly with high-energy particles in extrapolating to reactors. In the presence of a radial alpha particle birth gradient, this interaction can take the form of wave amplification rather than damping. While it is known that this amplification more easily occurs when launching from the tokamak high-field side, the extent of this amplification has not been made quantitative. Here, by tracing rays launched from the high-field-side of a tokamak, the required radial gradients to achieve amplification are calculated for a temperature and density regime consistent with a hot-ion-mode fusion reactor. These simulations, while valid only in the linear regime of wave amplification, nonetheless illustrate the possibilities for wave amplification using high-field launch of the lower hybrid wave.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗