Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linearized 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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Initial conditions for simulations of beam physics in linear induction accelerators

Flash radiography of hydrodynamic experiments driven by high explosives is a well-known diagnostic technique in use at many laboratories. At Los Alamos, the DualAxis Radiographic Hydrodynamic Test (DARHT) facility two linear induction electron accelerators (LIAs) make the bremsstrahlung radiographic source spots for point projection radiographs from orthogonal views. A new LIA, called Scorpius, is presently under development to advance this technology. To better understand electron-beam physics in these LIAs, numerical simulations are frequently performed with the objective of improving the radiography. At Los Alamos we frequently use the TRAK ray-trace and LSP particle-in-cell (PIC) codes to simulate the injector, and the XTR and LAMDA envelope/centroid codes along with LSP to simulate transport of the accelerated beam through the LIAs. The LIA simulations need the injected beam parameters as initial conditions for calculating beam transport and stability. The determination of these initial conditions is the topic of this note.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Correct Initial Conditions for Simulations of Beam Physics in Linear Induction Accelerators

Flash radiography of hydrodynamic experiments driven by high explosives is a well-known diagnostic technique in use at many laboratories. At Los Alamos, the Dual Axis Radiographic Hydrodynamic Test (DARHT) facility two linear induction electron accelerators (LIAs) make the bremsstrahlung radiographic source spots for point projection radiographs from orthogonal views. A new LIA, called Scorpius, is presently under development to advance this technology. To better understand electron-beam physics in these LIAs, numerical simulations are frequently performed with the objective of improving the radiography. At Los Alamos we frequently use the TRAK ray-trace and LSP particle-in-cell (PIC) codes to simulate the injector, and the XTR and LAMDA envelope/centroid codes along with LSP to simulate transport of the accelerated beam through the LIAs. The LIA simulations need the injected beam parameters as initial conditions for calculating beam transport and stability. The determination of these initial conditions is the topic of this note.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Single-Crystal X-Ray Spectropolarimeter

X-ray spectropolarimetry is a diagnostics technique suitable for measuring intense electron beams and strong magnetic fields in hot, dense plasmas. To increase the accuracy achievable, we investigated the possibility of measuring both orthogonal components of linear polarization simultaneously using a single crystal. We demonstrated experimentally for the first time that crystals with cubic symmetry can be used as polarizing beam splitters. We also found that the polarization-splitting planes in cubic crystals can be accessed in two crystal orientations, a property that can be used to decrease the clutter in polarized spectra. In parallel, we developed a ray-tracing code and are working to add the crystal reflectivity to instrument modeling codes, to help us interpret the spectra and to design spectroscopic and spectropolarimetric instruments. This work is motivated by the need to perform measurements of important quantities in difficult experimental conditions, such as the axial magnetic field at stagnation in magnetized liner inertial fusion loads.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Stochastic Averaging of Radiative Transfer Coefficients for Relativistic Electrons

Abstract Synchrotron emissivities, absorptivities, and Faraday rotation and conversion coefficients are needed in modeling a variety of astrophysical sources, including Event Horizon Telescope (EHT) sources. We develop a method for estimating transfer coefficients that exploits their linear dependence on the electron distribution function, decomposing the distribution function into a sum of parts each of whose emissivity can be calculated easily. We refer to this procedure as stochastic averaging and apply it in two contexts. First, we use it to estimate the emissivity of an isotropicκdistribution function with a high-energy cutoff. The resulting coefficients can be evaluated efficiently enough to be used directly in ray-tracing calculations, and we provide an example calculation. Second, we use stochastic averaging to assess the effect of subgrid turbulence on the volume-averaged emissivity and along the way provide a prescription for a turbulent emissivity. We find that for parameters appropriate to EHT sources turbulence reduces the emissivity slightly. In the infrared, turbulence can dramatically increase the emissivity.

Astronomy & Astrophysics↗
NASA NTRS1996

Modeling of Jovian Hectometric Radiation Source Locations: Ulysses Observations

The Unified Radio and Plasma Wave (URAP) experiment on Ulysses has provided unique high latitude measurements of Jovian hectometric radiation (HOM) during its encounter with Jupiter in February 1992. URAP was the first radio instrument in the Jovian environment with radio direction-finding capability, which was previously used to determine the HOM source locations in the Jovian magnetosphere. These initial source location determinations were based on several assumptions, including the neglect of refractive effects, which may be tested. We have, for the first time, combined the measured incident ray-direction at the spacecraft with a model magnetosphere to directly trace the rays back to the HOM source. We concentrate on the observations of HOM from high northern latitudes when Ulysses was at distances less than 15 R(sub j). The three- dimensional ray-tracing calculations presented here indicate that the HOM sources probably lie on L shells in the range 3 less than or approximately equal to L less than 7 (tilted dipole magnetic field model) consistent with previous determinations that ignored the effects of refraction. The ray-tracing results, however, indicate that wave refraction due to the Io torus and the magnetic field can significantly influence the precise source location. We show that constraints on the locations imposed by the gyroemission mechanism suggest that the lo torus density may have experienced temporal and/or spatial fluctuations during the Ulysses observations of HOM. Finally, in the cold plasma approximation we demonstrate that even if the emission were nearly linearly polarized near the source region, almost circular polarization will be observed at Ulysses, in agreement with observations.

Menietti, J. D.↗
DOE OSTI2015

Simulation of radiation damping in rings, using stepwise ray-tracing methods

The ray-tracing code Zgoubi computes particle trajectories in arbitrary magnetic and/or electric field maps or analytical field models. It includes a built-in fitting procedure, spin tracking many Monte Carlo processes. The accuracy of the integration method makes it an efficient tool for multi-turn tracking in periodic machines. Energy loss by synchrotron radiation, based on Monte Carlo techniques, had been introduced in Zgoubi in the early 2000s for studies regarding the linear collider beam delivery system. However, only recently has this Monte Carlo tool been used for systematic beam dynamics and spin diffusion studies in rings, including eRHIC electron-ion collider project at the Brookhaven National Laboratory. Some beam dynamics aspects of this recent use of Zgoubi capabilities, including considerations of accuracy as well as further benchmarking in the presence of synchrotron radiation in rings, are reported here.

43 PARTICLE ACCELERATORS↗
NASA NTRS1990

Optical modeling for dynamics and control analysis

This paper presents a coordinate-free ray-trace analysis of optical beam trains consisting of mirrors, lenses and reference surfaces. The analysis leads to optical models that can be directly integrated with standard structure and control models for integrated instrument design, analysis and simulation. This capability is required for a coming generation of spaceborn optical instruments which use controlled optical elements supported by flexible structures. New results include analytic formulas for optical sensitivities as functions of structural and control geometric parameters. An example problem calculates a linear (small-motion) optical model, derives optical gains for small- and large-angle controlled mirrors, and shows the propagation of a wavefront.

Redding, David C.↗
DOE OSTI2017

Linear feature detection algorithm for astronomical surveys – II. Defocusing effects on meteor tracks

Given the current limited knowledge of meteor plasma micro-physics and its interaction with the surrounding atmosphere and ionosphere, meteors are a highly interesting observational target for high-resolution wide-field astronomical surveys. Such surveys are capable of resolving the physical size of meteor plasma heads, but they produce large volumes of images that need to be automatically inspected for possible existence of long linear features produced by meteors. Here in this paper, we show how big aperture sky survey telescopes detect meteors as defocused tracks with a central brightness depression. We derive an analytic expression for a defocused point source meteor track and use it to calculate brightness profiles of meteors modelled as uniform brightness discs. We apply our modelling to meteor images as seen by the Sloan Digital Sky Survey and Large Synoptic Survey Telescope telescopes. The expression is validated by Monte Carlo ray-tracing simulations of photons travelling through the atmosphere and the Large Synoptic Survey Telescope telescope optics. We show that estimates of the meteor distance and size can be extracted from the measured full width at half-maximum and the strength of the central dip in the observed brightness profile. However, this extraction becomes difficult when the defocused meteor track is distorted by the atmospheric seeing or contaminated by a long-lasting glowing meteor trail. The full width at half-maximum of satellite tracks is distinctly narrower than meteor values, which enables removal of a possible confusion between satellites and meteors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A Scalable Gaussian Process Approach to Shear Mapping with MuyGPs

Analysis of cosmic shear is an integral part of understanding structure growth across cosmic time, which in turn provides us with information about the nature of dark energy. Conventional methods generate shear maps from which we can infer the matter distribution in the universe. Current methods (e.g., Kaiser–Squires inversion) for generating these maps, however, are tricky to implement and can introduce bias. Recent alternatives construct a spatial process prior for the lensing potential, which allows for inference of the convergence and shear parameters given lensing shear measurements. Realizing these spatial processes, however, scales cubically in the number of observations—an unacceptable expense as near-term surveys expect billions of correlated measurements. Therefore, we present a linearly scaling shear map construction alternative using a scalable Gaussian process prior called MuyGPs. MuyGPs avoids cubic scaling by conditioning interpolation on only nearest neighbors and fits hyperparameters using batched leave-one-out cross-validation. This work is the first step toward a full, scalable mass mapping method. We work in a simplified regime where we validate our method by interpolating and analyzing maps given noisy point-estimate data from all three shear fields, taken from a suite of N -body ray-tracing simulations. We also show that we can perform these operations at the scale of billions of galaxies on high-performance computing platforms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

An electron Bernstein wave heating scheme for the Proto-MPEX linear device

Linear plasma devices, such as the Prototype Material Plasma Exposure eXperiment (Proto-MPEX), are an economic method to study plasma-material interactions under high heat and particle fluxes. The Proto-MPEX device at Oak Ridge National Laboratory is a high-density helicon plasma generator with additional resonant electron heating to study plasma-material interactions in ITER-like conditions. Here, we present modeling results of an electron Bernstein wave (EBW) heating scheme for this device using a microwave-based system. A 28 GHz injection system is available to provide power levels up to 100 kW injected into low electron temperature and high electron density plasmas. A modified version of the GENRAY ray-tracing code, GENRAY-C, has been used to determine the EBW and electron cyclotron heating wave accessibility for these overdense plasmas. Modeling has shown that greater than 80% of mode conversion from the ordinary mode to the EBW is possible for 28 GHz injection in Proto-MPEX. Calculations show that by utilizing the effects of Doppler-shifted resonance absorption, power deposition near the core is possible. Finally, significant collisional damping, leading to edge absorption, is expected to occur for neutral pressures >1.0 mTorr; however, experiments can be designed to minimize these effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

The kinetic Sunyaev-Zel'dovich effect from reionization: simulated full-sky maps at arcminute resolution

The kinetic Sunyaev–Zel’dovich (kSZ) effect results from Thomson scattering by coherent flows in the reionized intergalactic medium. We present new results based on ray-tracing an 8 Gpc/h realization of reionization with resolution elements 2 Mpc/h (subtending ∼1′ at z = 6) on a side to create a full-sky kSZ map. The realization includes, self-consistently, the effects of reionization on scales corresponding to multipoles 10≲ℓ≲5000. We separate the kSZ map into Doppler (v), Ostriker–Vishniac (δv), patchy (xv), and third-order (xδv) components, and compute explicitly all the auto- and cross-correlations (e.g., 〈vv〉, 〈δvxv〉, etc.) that contribute to the total power. We find a complex and nonmonotonic dependence on the duration of reionization at ℓ∼300 and evidence for a non-negligible (10%–30%) contribution from connected four-point correlations, 〈xvxv〉{sub c}, usually neglected in analytical models. We also investigate the cross-correlation of linear matter and large-scale kSZ temperature fluctuations, focusing on (1) cross-power spectra with biased tracers of the matter density and (2) cold spots from infall onto large, rare H ii regions centered on peaks in the matter distribution at redshifts z>10 that are a generic non-Gaussian feature of patchy reionization. Finally, we show that the reionization history can be reconstructed at 5σ–10σ significance by correlating full-sky 21 cm maps stacked in bins with Δν = 10 MHz with existing cosmic microwave background (CMB) temperature maps at ℓ<500, raising the prospects for probing reionization by correlating CMB and LSS measurements. The resulting kSZ maps have been made publicly available at www.cita.utoronto.ca/∼malvarez/research/ksz-data/.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Power transport efficiency during O-X-B 2nd harmonic electron cyclotron heating in a helicon linear plasma device

The principal objective of this work is to report on the power coupled to a tungsten target in the prototype-material plasma exposure experiment device during oblique injection of a microwave beam (<70 kW at 28 GHz) into a high-power (~100 kW at 13.56 MHz) over-dense (${n_e} > 1 \times {10^{19}}{\text{ }}{{\text{m}}^{ - 3}}$) deuterium helicon plasma column. The experimental setup, electron heating system, electron heating scheme, and IR thermographic diagnostic for quantifying the power transport is described in detail. It is demonstrated that the power transported to the target can be effectively controlled by adjusting the magnetic field profile. Using this method, heat fluxes up to 22 MW m -2 and power transport efficiencies in the range of 17%–20% have been achieved using 70 kW of microwave power. It is observed that most of the heat flux is confined to a narrow region at the plasma periphery. Ray-tracing calculations are presented which indicate that the power is coupled to the plasma electrons via an O-X-B mode conversion process. Calculations indicate that the microwave power is absorbed in a single pass at the plasma periphery via collisions and in the over-dense region via 2nd harmonic cyclotron resonance of the electron Bernstein wave. The impact of our results is discussed in the context of MPEX.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

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↗
DOE OSTI2013

Maps of CMB lensing deflection from N-body simulations in Coupled Dark Energy Cosmologies

We produce lensing potential and deflection-angle maps in order to simulate the weak gravitational lensing of the Cosmic Microwave Background (CMB) via ray-tracing through the COupled Dark Energy Cosmological Simulations (CoDECS), the largest suite of N-body simulations to date for interacting Dark Energy cosmologies. The constructed maps faithfully reflect the N-body cosmic structures on a range of scales going from the arcminute to the degree scale, limited only by the resolution and extension of the simulations. We investigate the variation of the lensing pattern due to the underlying Dark Energy (DE) dynamics, characterised by different background and perturbation behaviours as a consequence of the interaction between the DE field and Cold Dark Matter (CDM). In particular, we study in detail the results from three cosmological models differing in the background and perturbations evolution at the epoch in which the lensing cross section is most effective, corresponding to a redshift of ∼ 1, with the purpose to isolate their imprints in the lensing observables, regardless of the compatibility of these models with present constraints. The scenarios investigated here include a reference ΛCDM cosmology, a standard coupled DE (cDE) scenario, and a ''bouncing'' cDE scenario. For the standard cDE scenario, we find that typical differences in the lensing potential result from two effects: the enhanced growth of linear CDM density fluctuations with respect to the ΛCDM case, and the modified nonlinear dynamics of collapsed structures induced by the DE-CDM interaction. As a consequence, CMB lensing highlights the DE impact in the cosmological expansion, even in the degenerate case where the amplitude of the linear matter density perturbations, parametrised through σ{sub 8}, is the same in both the standard cDE and ΛCDM cosmologies. For the ''bouncing'' scenario, we find that the two opposite behaviours of the lens density contrast and of the matter abundance lead to a counter-intuitive effect, making the power of the lensing signal in this model lower by 10% than in the ΛCDM scenario. Moreover, we compare the behaviour of CDM and baryons in CoDECS separately, in order to isolate effects coming from the coupling with the DE component. We find that, in the bouncing scenario, baryons show an opposite trend with respect to CDM, due to the coupling of the latter with the DE component. These results confirm the relevance of CMB lensing as a probe for DE at the early stages of cosmic acceleration, and demonstrate the reliability of N-body based large scale CMB lensing simulations in the context of DE studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Modelling of the EAST lower-hybrid current drive experiment using GENRAY/CQL3D and TORLH/CQL3D

The coupled GENRAY-CQL3D code has been used to do systematic ray-tracing and Fokker–Planck analysis for EAST Lower Hybrid wave Current Drive (LHCD) experiments. Despite being in the weak absorption regime, the experimental level of LH current drive is successfully simulated, by taking into account the variations in the parallel wavenumber due to the toroidal effect. The effect of radial transport of the fast LH electrons in EAST has also been studied, which shows that a modest amount of radial transport diffusion can redistribute the fast LH current significantly. Taking advantage of the new capability in GENRAY, the actual Scrape Off Layer (SOL) model with magnetic field, density, temperature, and geometry is included in the simulation for both the lower and the higher density cases, so that the collisional losses of Lower Hybrid Wave (LHW) power in the SOL has been accounted for, which together with fast electron losses can reproduce the LHCD experimental observations in different discharges of EAST. We have also analyzed EAST discharges where there is a significant ohmic contribution to the total current, and good agreement with experiment in terms of total current has been obtained. Also, the full-wave code TORLH has been used for the simulation of the LH physics in the EAST, including full-wave effects such as diffraction and focusing which may also play an important role in bridging the spectral gap. The comparisons between the GENRAY and the TORLH codes are done for both the Maxwellian and the quasi-linear electron Landau damping cases. In conclusion, these simulations represent an important addition to the validation studies of the GENRAY-CQL3D and TORLH models being used in weak absorption scenarios of tokamaks with large aspect ratio.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

SU-E-T-252: Developing a Pencil Beam Dose Calculation Algorithm for CyberKnife System

Purpose: Currently there are two dose calculation algorithms available in the Cyberknife planning system: ray-tracing and Monte Carlo, which is either not accurate or time-consuming for irregular field shaped by the MLC that was recently introduced. The purpose of this study is to develop a fast and accurate pencil beam dose calculation algorithm which can handle irregular field. Methods: A pencil beam dose calculation algorithm widely used in Linac system is modified. The algorithm models both primary (short range) and scatter (long range) components with a single input parameter: TPR{sub 20}/{sub 10}. The TPR{sub 20}/{sub 20}/{sub 10} value was first estimated to derive an initial set of pencil beam model parameters (PBMP). The agreement between predicted and measured TPRs for all cones were evaluated using the root mean square of the difference (RMSTPR), which was then minimized by adjusting PBMPs. PBMPs are further tuned to minimize OCR RMS (RMSocr) by focusing at the outfield region. Finally, an arbitrary intensity profile is optimized by minimizing RMSocr difference at infield region. To test model validity, the PBMPs were obtained by fitting to only a subset of cones (4) and applied to all cones (12) for evaluation. Results: With RMS values normalized to the dmax and all cones combined, the average RMSTPR at build-up and descending region is 2.3% and 0.4%, respectively. The RMSocr at infield, penumbra and outfield region is 1.5%, 7.8% and 0.6%, respectively. Average DTA in penumbra region is 0.5mm. There is no trend found in TPR or OCR agreement among cones or depths. Conclusion: We have developed a pencil beam algorithm for Cyberknife system. The prediction agrees well with commissioning data. Only a subset of measurements is needed to derive the model. Further improvements are needed for TPR buildup region and OCR penumbra. Experimental validations on MLC shaped irregular field needs to be performed. This work was partially supported by the National Natural Science Foundation of China (61171005) and the China Scholarship Council (CSC)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Paraxial ray-tracing approach for the simulation of ultrasonic inspection of welds

On-site inspection of bimetallic or austenitic welds can be very difficult to interpret owing to their internal structures. Skewing and splitting of the ultrasonic beam may occur due to the anisotropic and inhomogeneous properties of the welding material. In this paper, we present a ray-based method to simulate the propagation of ultrasonic waves in such structures. The formalism is based on dynamic ray tracing system in Cartesian coordinates along a reference ray. Standard ray tracing consists in the solution of a system of linear ordinary differential equations of the first order and is used to determine the trajectory of the ray. Likewise, dynamic ray tracing (DRT) also called paraxial ray tracing consists in the solution of an additional system of linear ordinary differential equations along the ray allowing paraxial quantities to be computed. It is used to evaluate the geometrical spreading and amplitude along the ray and in its vicinity. DRT is applied on a smooth representation of the elastic properties of the weld obtained thanks to an image processing technique applied on a macrograph of the weld. Simulation results are presented and compared to finite elements and experimental results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

A digitally reconstructed radiograph algorithm calculated from first principles

Purpose: To develop an algorithm for computing realistic digitally reconstructed radiographs (DRRs) that match real cone-beam CT (CBCT) projections with no artificial adjustments. Methods: The authors used measured attenuation data from cone-beam CT projection radiographs of different materials to obtain a function to convert CT number to linear attenuation coefficient (LAC). The effects of scatter, beam hardening, and veiling glare were first removed from the attenuation data. Using this conversion function the authors calculated the line integral of LAC through a CT along rays connecting the radiation source and detector pixels with a ray-tracing algorithm, producing raw DRRs. The effects of scatter, beam hardening, and veiling glare were then included in the DRRs through postprocessing. Results: The authors compared actual CBCT projections to DRRs produced with all corrections (scatter, beam hardening, and veiling glare) and to uncorrected DRRs. Algorithm accuracy was assessed through visual comparison of projections and DRRs, pixel intensity comparisons, intensity histogram comparisons, and correlation plots of DRR-to-projection pixel intensities. In general, the fully corrected algorithm provided a small but nontrivial improvement in accuracy over the uncorrected algorithm. The authors also investigated both measurement- and computation-based methods for determining the beam hardening correction, and found the computation-based method to be superior, as it accounted for nonuniform bowtie filter thickness. The authors benchmarked the algorithm for speed and found that it produced DRRs in about 0.35 s for full detector and CT resolution at a ray step-size of 0.5 mm. Conclusions: The authors have demonstrated a DRR algorithm calculated from first principles that accounts for scatter, beam hardening, and veiling glare in order to produce accurate DRRs. The algorithm is computationally efficient, making it a good candidate for iterative CT reconstruction techniques that require a data fidelity term based on the matching of DRRs and projections.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗