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 271 records · Page 15

Recent Trends in the Analysis of Quasioptical Systems

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

Hoppe, D. J.↗

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

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

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Fabrication and testing of a newly designed slit system for depth-resolved X-ray diffraction measurements

A new system of slits called `spiderweb slits' have been developed for depth-resolved powder or polycrystalline X-ray diffraction measurements. The slits act on diffracted X-rays to select a particular gauge volume of sample, while absorbing diffracted X-rays from outside of this volume. Although the slit geometry is to some extent similar to that of previously developed conical slits or spiral slits, this new design has advantages over the previous ones in use for complex heterogeneous materials and in situ and operando diffraction measurements. For example, the slits can measure a majority of any diffraction cone for any polycrystalline material, over a continuous range of diffraction angles, and work for X-ray energies of tens to hundreds of kiloelectronvolts. The design is generated and optimized using ray-tracing simulations, and fabricated through laser micromachining. The first prototype was successfully tested at the X17A beamline at the National Synchrotron Light Source, and shows similar performance to simulations, demonstrating gauge volume selection for standard powders, for all diffraction peaks over angles of 2–10°. A similar, but improved, design will be implemented at the X-ray Powder Diffraction beamline at the National Synchrotron Light Source II.

depth-resolved↗

Ray methods in random media

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

Hornstein, J.↗

Integrated Modeling Activities for the James Webb Space Telescope (JWST): Structural-Thermal-Optical Analysis

This is a continuation of a series of papers on modeling activities for JWST. The structural-thermal- optical, often referred to as "STOP", analysis process is used to predict the effect of thermal distortion on optical performance. The benchmark STOP analysis for JWST assesses the effect of an observatory slew on wavefront error. The paper begins an overview of multi-disciplinary engineering analysis, or integrated modeling, which is a critical element of the JWST mission. The STOP analysis process is then described. This process consists of the following steps: thermal analysis, structural analysis, and optical analysis. Temperatures predicted using geometric and thermal math models are mapped to the structural finite element model in order to predict thermally-induced deformations. Motions and deformations at optical surfaces are input to optical models and optical performance is predicted using either an optical ray trace or WFE estimation techniques based on prior ray traces or first order optics. Following the discussion of the analysis process, results based on models representing the design at the time of the System Requirements Review. In addition to baseline performance predictions, sensitivity studies are performed to assess modeling uncertainties. Of particular interest is the sensitivity of optical performance to uncertainties in temperature predictions and variations in metal properties. The paper concludes with a discussion of modeling uncertainty as it pertains to STOP analysis.

Johnston, John D.↗

Resolving broadening of ECRH deposition due to edge turbulence in DIII-D by heat deposition measurement

Interaction between microwave power, used for local heating and mode control, and density fluctuations can produce a broadening of the injected beam, as confirmed in experiment and simulation. Increased power deposition width could impact suppression of tearing mode structures on ITER[1]. This work discusses the experimental portion of an effort to understand scattering of injected microwaves by turbulence on the DIII-D tokamak. The corresponding theoretical modeling work can be found in M.B. Thomas et. al.: Submitted to Nuclear Fusion (2017)[Author Note - this paper to be published in same journal]. In a set of perturbative heat transport experiments, tokamak edge millimeter-scale fluctuation levels and microwave heat deposition are measured simultaneously. Beam broadening is separated from heat transport through fitting of modulated fluxes[2]. Electron temperature measurements from a 500 kHz, 48-channel ECE radiometer are Fourier analyzed and used to calculate a deposition-dependent flux[3]. Consistency of this flux with a transport model is evaluated. A diffusive(∝ ∇˜ T e ) and convective(∝ ˜ T e ) transport solution is linearized and compared with energy conservation-derived fluxes. Comparison between these two forms of heat flux is used to evaluate the quality of ECRF deposition profiles, and a Χ 2 minimization finds a significant broadening of 1D equilibrium ray tracing calculations from the benchmarked TORAY-GA ray tracing code[4] is needed. The physical basis, cross-validation, and application of the heat flux method is presented. The method is applied to a range of DIII-D discharges and finds a broadening factor of the deposition profile width which scales linearly with edge density fluctuation level. Furthermore, these experimental results are found to be consistent with the full-wave beam broadening measured by the 3D full wave simulations in the same discharges[5].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SU-F-T-575: Verification of a Monte-Carlo Small Field SRS/SBRT Dose Calculation System

Purpose: To develop a methodology for validation of a Monte-Carlo dose calculation model for robotic small field SRS/SBRT deliveries. Methods: In a robotic treatment planning system, a Monte-Carlo model was iteratively optimized to match with beam data. A two-part analysis was developed to verify this model. 1) The Monte-Carlo model was validated in a simulated water phantom versus a Ray-Tracing calculation on a single beam collimator-by-collimator calculation. 2) The Monte-Carlo model was validated to be accurate in the most challenging situation, lung, by acquiring in-phantom measurements. A plan was created and delivered in a CIRS lung phantom with film insert. Separately, plans were delivered in an in-house created lung phantom with a PinPoint chamber insert within a lung simulating material. For medium to large collimator sizes, a single beam was delivered to the phantom. For small size collimators (10, 12.5, and 15mm), a robotically delivered plan was created to generate a uniform dose field of irradiation over a 2×2cm{sup 2} area. Results: Dose differences in simulated water between Ray-Tracing and Monte-Carlo were all within 1% at dmax and deeper. Maximum dose differences occurred prior to dmax but were all within 3%. Film measurements in a lung phantom show high correspondence of over 95% gamma at the 2%/2mm level for Monte-Carlo. Ion chamber measurements for collimator sizes of 12.5mm and above were within 3% of Monte-Carlo calculated values. Uniform irradiation involving the 10mm collimator resulted in a dose difference of ∼8% for both Monte-Carlo and Ray-Tracing indicating that there may be limitations with the dose calculation. Conclusion: We have developed a methodology to validate a Monte-Carlo model by verifying that it matches in water and, separately, that it corresponds well in lung simulating materials. The Monte-Carlo model and algorithm tested may have more limited accuracy for 10mm fields and smaller.

60 APPLIED LIFE SCIENCES↗

A pencil beam dose calculation model for CyberKnife system

Purpose: CyberKnife system is initially equipped with fixed circular cones for stereotactic radiosurgery. Two dose calculation algorithms, Ray-Tracing and Monte Carlo, are available in the supplied treatment planning system. A multileaf collimator system was recently introduced in the latest generation of system, capable of arbitrarily shaped treatment field. The purpose of this study is to develop a model based dose calculation algorithm to better handle the lateral scatter in an irregularly shaped small field for the CyberKnife system. Methods: A pencil beam dose calculation algorithm widely used in linac based treatment planning system was modified. The kernel parameters and intensity profile were systematically determined by fitting to the commissioning data. The model was tuned using only a subset of measured data (4 out of 12 cones) and applied to all fixed circular cones for evaluation. The root mean square (RMS) of the difference between the measured and calculated tissue-phantom-ratios (TPRs) and off-center-ratio (OCR) was compared. Three cone size correction techniques were developed to better fit the OCRs at the penumbra region, which are further evaluated by the output factors (OFs). The pencil beam model was further validated against measurement data on the variable dodecagon-shaped Iris collimators and a half-beam blocked field. Comparison with Ray-Tracing and Monte Carlo methods was also performed on a lung SBRT case. Results: The RMS between the measured and calculated TPRs is 0.7% averaged for all cones, with the descending region at 0.5%. The RMSs of OCR at infield and outfield regions are both at 0.5%. The distance to agreement (DTA) at the OCR penumbra region is 0.2 mm. All three cone size correction models achieve the same improvement in OCR agreement, with the effective source shift model (SSM) preferred, due to their ability to predict more accurately the OF variations with the source to axis distance (SAD). In noncircular field validation, the pencil beam calculated results agreed well with the film measurement of both Iris collimators and the half-beam blocked field, fared much better than the Ray-Tracing calculation. Conclusions: The authors have developed a pencil beam dose calculation model for the CyberKnife system. The dose calculation accuracy is better than the standard linac based system because the model parameters were specifically tuned to the CyberKnife system and geometry correction factors. The model handles better the lateral scatter and has the potential to be used for the irregularly shaped fields. Comprehensive validations on MLC equipped system are necessary for its clinical implementation. It is reasonably fast enough to be used during plan optimization.

60 APPLIED LIFE SCIENCES↗

Spacecraft observations of man-made whistler-mode signals near the electron gyrofrequency

The reported investigation extends the range of whistler-mode wave observations to a wave frequency/electron gyrofrequency ratio of about 0.9, where an abrupt cutoff is observed. This cutoff can be explained entirely in terms of accessibility and hence, if there is damping, it must be limited to normalized frequencies above 0.9. In connection with a study of the behavior of the signal intensity, ray tracings were carried out at 80 kHz. The ray-tracing calculations were carried out with the aid of a computer program written by Walter (1969) and modified by Angerami (1970).

Dunckel, N.↗

Effects of refractive index mismatch in optical CT imaging of polymer gel dosimeters

Purpose: Proposing an image reconstruction technique, algebraic reconstruction technique-refraction correction (ART-rc). The proposed method takes care of refractive index mismatches present in gel dosimeter scanner at the boundary, and also corrects for the interior ray refraction. Polymer gel dosimeters with high dose regions have higher refractive index and optical density compared to the background medium, these changes in refractive index at high dose results in interior ray bending. Methods: The inclusion of the effects of refraction is an important step in reconstruction of optical density in gel dosimeters. The proposed ray tracing algorithm models the interior multiple refraction at the inhomogeneities. Jacob’s ray tracing algorithm has been modified to calculate the pathlengths of the ray that traverses through the higher dose regions. The algorithm computes the length of the ray in each pixel along its path and is used as the weight matrix. Algebraic reconstruction technique and pixel based reconstruction algorithms are used for solving the reconstruction problem. The proposed method is tested with numerical phantoms for various noise levels. The experimental dosimetric results are also presented. Results: The results show that the proposed scheme ART-rc is able to reconstruct optical density inside the dosimeter better than the results obtained using filtered backprojection and conventional algebraic reconstruction approaches. The quantitative improvement using ART-rc is evaluated using gamma-index. The refraction errors due to regions of different refractive indices are discussed. The effects of modeling of interior refraction in the dose region are presented. Conclusions: The errors propagated due to multiple refraction effects have been modeled and the improvements in reconstruction using proposed model is presented. The refractive index of the dosimeter has a mismatch with the surrounding medium (for dry air or water scanning). The algorithm reconstructs the dose profiles by estimating refractive indices of multiple inhomogeneities having different refractive indices and optical densities embedded in the dosimeter. This is achieved by tracking the path of the ray that traverses through the dosimeter. Extensive simulation studies have been carried out and results are found to be matching that of experimental results.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Reflectance model of a plant leaf

A light ray, incident at 5 deg to the normal, is geometrically plotted through the drawing of the cross section of a soybean leaf using Fresnel's Equations and Snell's Law. The optical mediums of the leaf considered for ray tracing are: air, cell sap, chloroplast, and cell wall. The above ray is also drawn through the same leaf cross section considering cell wall and air as the only optical mediums. The values of the reflection and transmission found from ray tracing agree closely with the experimental results obtained using a Beckman DK-2A Spectroreflectometer. Similarly a light ray, incident at about 60 deg to the normal, is drawn through the palisade cells of a soybean leaf to illustrate the pathway of light, incident at an oblique angle, through the palisade cells.

Kumar, R.↗

Quasioptical modeling of wave beams with and without mode conversion. III. Numerical simulations of mode-converting beams

This work continues a series of papers where we propose an algorithm for quasi-optical modeling of electromagnetic beams with and without mode conversion. The general theory was reported in the first paper of this series, where a parabolic partial differential equation was derived for the field envelope that may contain one or multiple modes with close group velocities. In the second paper, we presented a corresponding code PARADE (PAraxial RAy DEscription) and its test applications to single-mode beams. Here, we report quasi-optical simulations of mode-converting beams for the first time. Here, we also demonstrate that PARADE can model splitting of two-mode beams. The numerical results produced by PARADE show good agreement with those of one-dimensional full-wave simulations and also with conventional ray tracing (to the extent that one-dimensional and ray-tracing simulations are applicable).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electron cyclotron emission measurements and modeling in optically semi-transparent plasma on the HSX stellarator

On the Helically Symmetric Experiment (HSX), electron cyclotron emission (ECE) at the second harmonic of the electron gyrofrequency is measured with a 16-channel heterodyne receiver. An ECE radiometer diagnostic has been calibrated with a broadband frequency source, allowing independent measurements of the electron temperature. The antenna pattern has a significant impact on the measurements, and modeling its effect is necessary to produce an ECE radiation temperature spectrum. Measurements, calibration results, and modeling of the ECE data are presented in this paper. At a relatively high plasma density, HSX plasmas are thermal, and the radiation temperature measured by the radiometer is in good agreement with the electron temperature measured by a Thomson scattering diagnostic. At a relatively low plasma density, the deviation from a Maxwellian distribution is large, and asymmetry in the radiation temperature profile is measured. These measurements are supported by single-pass and multi-pass ray tracing calculations of the electron cyclotron emission using the TRAVIS ray tracing code. Thermal and nonthermal electron populations exist in HSX, and the effect of these populations on the asymmetry in the measured emission is presented.

Cyclotron radiation↗

Prediction of regional infrasound produced by supersonic sources using a ray-based Mach cone source

The geometry of the Mach cone produced by a supersonic source is analyzed and mapped into initial conditions used in acoustic ray tracing. The resulting source model is combined with spherical geometry ray tracing methods to enable propagation simulations for infrasonic signals produced by bolides, space debris, rockets, aircraft, and other fast-than-sound sources out to typical infrasonic observation distances of hundreds or thousands of kilometers. Idealized linear and parabolic trajectories typical of bolides and rockets, respectively, are used to demonstrate the calculation of regional infrasonic signals produced by such sources and characteristics of the radiated infrasonic waves are found to vary strongly with the geometry of the trajectory and atmospheric structure. Predicted regional infrasonic signals are compared with those observed from a November 2020 bolide that passed over Scandinavia using a combination of institutionally maintained infrasound stations and “citizen scientist” data from the Raspberry Shake data repository.

47 OTHER INSTRUMENTATION↗

Quasioptical modeling of wave beams with and without mode conversion: III. Numerical simulations of mode-converting beams

This work continues a series of papers where we propose an algorithm for quasioptical modeling of electromagnetic beams with and without mode conversion. The general theory was reported in the first paper of this series, where a parabolic partial differential equation was derived for the field envelope that may contain one or multiple modes with close group velocities. In the second paper, we presented a corresponding code PARADE (PAraxial RAy DEscription) and its test applications to single-mode beams. Here, we report quasioptical simulations of mode-converting beams for the first time. We also demonstrate that PARADE can model splitting of two-mode beams. The numerical results produced by PARADE show good agreement with those of one-dimensional full-wave simulations and also with conventional ray tracing (to the extent that one-dimensional and ray-tracing simulations are applicable).

ECRH↗

Prediction of regional infrasound produced by supersonic sources using a ray-based Mach cone source

The geometry of the Mach cone produced by a supersonic source is analyzed and mapped into initial conditions used in acoustic ray tracing. The resulting source model is combined with spherical geometry ray tracing methods to enable propagation simulations for infrasonic signals produced by bolides, space debris, rockets, aircraft, and other fast-than-sound sources out to typical infrasonic observation distances of hundreds or thousands of kilometers. Idealized linear and parabolic trajectories typical of bolides and rockets, respectively, are used to demonstrate the calculation of regional infrasonic signals produced by such sources and characteristics of the radiated infrasonic waves are found to vary strongly with the geometry of the trajectory and atmospheric structure. Predicted regional infrasonic signals are compared with those observed from a November 2020 bolide that passed over Scandinavia using a combination of institutionally maintained infrasound stations and “citizen scientist” data from the Raspberry Shake data repository.

42 ENGINEERING↗

Concentrator optical characterization using computer mathematical modelling and point source testing

The optical characteristics of a paraboloidal solar concentrator are analyzed using the intercept factor curve (a format for image data) to describe the results of a mathematical model and to represent reduced data from experimental testing. This procedure makes it possible not only to test an assembled concentrator, but also to evaluate single optical panels or to conduct non-solar tests of an assembled concentrator. The use of three-dimensional ray tracing computer programs to calculate the mathematical model is described. These ray tracing programs can include any type of optical configuration from simple paraboloids to array of spherical facets and can be adapted to microcomputers or larger computers, which can graphically display real-time comparison of calculated and measured data.

Dennison, E. W.↗

Software to model AXAF-I image quality

A modular user-friendly computer program for the modeling of grazing-incidence type x-ray optical systems has been developed. This comprehensive computer software GRAZTRACE covers the manipulation of input data, ray tracing with reflectivity and surface deformation effects, convolution with x-ray source shape, and x-ray scattering. The program also includes the capabilities for image analysis, detector scan modeling, and graphical presentation of the results. A number of utilities have been developed to interface the predicted Advanced X-ray Astrophysics Facility-Imaging (AXAF-I) mirror structural and thermal distortions with the ray-trace. This software is written in FORTRAN 77 and runs on a SUN/SPARC station. An interactive command mode version and a batch mode version of the software have been developed.

Ahmad, Anees↗