SEARCH · Engineering Papers

Results for “scattering”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 649 records · Page 36

All sources NASA NTRS DOE OSTI

DOE OSTI2023

A Refined Understanding of the Ice Cloud Longwave Scattering Effects in Climate Model

Abstract Because longwave (LW) absorption by greenhouse gases and clouds is more significant than the LW scattering effect by clouds, most climate models neglect cloud LW scattering to save computational costs. Ignoring cloud LW scattering directly overestimates outgoing longwave radiation (OLR). This study included ice‐cloud LW scattering treatment in the Exascale Energy Earth System Model (E3SM) version 2 and ran fully‐coupled simulations, prescribed sea surface temperature simulations, and offline radiative transfer calculations to comprehensively assess the impact of ice‐cloud LW scattering on global climate simulation. The instantaneous effect due to ice‐cloud LW scattering reduces the OLR by ∼1 W/m 2 on the global average and 2 W/m 2 on the tropical average. Tropospheric warming and high cloud amount reduction act to partially compensate for such instantaneous OLR reduction caused by the inclusion of LW scattering. When the simulation reaches the equilibrium, the surface warms by 0.66 K on average with respect to the simulation without LW scattering, with the Arctic surface temperature differences more than twice as large as that of the global mean. The impact of including LW scattering on the simulated climate change in response to 4 × CO 2 is also assessed. While including the cloud LW scattering does not significantly modify radiative forcing and total radiative feedback under such a scenario, it results in a 10% more positive cloud feedback.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2017

On-the-Fly Generation of Differential Resonance Scattering Probability Distribution Functions for Monte Carlo Codes

Current Monte Carlo codes use one of three models: (1) the asymptotic scattering model, (2) the free gas scattering model, or (3) the S(α,β) model, depending on the neutron energy and the specific Monte Carlo code. This thesis addresses the consequences of using the free gas scattering model, which assumes that the neutron interacts with atoms in thermal motion in a monatomic gas in thermal equilibrium at material temperature, T. Most importantly, the free gas model assumes the scattering cross section is constant over the neutron energy range, which is usually a good approximation for light nuclei, but not for heavy nuclei where the scattering cross section may have several resonances in the epithermal region. Several researchers in the field have shown that the exact resonance scattering model is temperaturedependent, and neglecting the resonances in the lower epithermal range can under-predict resonance absorption due to the upscattering phenomenon mentioned above, leading to an over-prediction of keff by several hundred pcm. Existing methods to address this issue involve changing the neutron weights or implementing an extra rejection scheme in the free gas sampling scheme, and these all involve performing the collision analysis in the center-of-mass frame, followed by a conversion back to the laboratory frame to continue the random walk of the neutron. The goal of this paper was to develop a sampling methodology that (1) accounted for the energydependent scattering cross sections in the collision analysis and (2) was performed in the laboratory frame,avoiding the conversion to the center-of-mass frame. The energy dependence of the scattering cross section was modeled with even-ordered polynomials (2nd and 4th order) to approximate the scattering cross section in Blackshaw’s equations for the moments of the differential scattering PDFs. These moments were used to sample the outgoing neutron speed and angle in the laboratory frame on-the-fly during the random walk of the neutron. Results for criticality studies on fuel pin and fuel assembly calculations using methods developed in this dissertation showed very close comparison to results using the reference Dopplerbroadened rejection correction (DBRC) scheme.

42 ENGINEERING↗

DOE OSTI2022

Anisotropic angular scattering models of elastic electron-neutral collisions for Monte Carlo plasma simulations

Many laboratory and industrial plasma applications require accurate modeling techniques to understand the interplay between microscopic and macroscopic processes. A prime example of this interplay is how particle and Monte Carlo (MC) simulation codes describe angular scattering of electrons following elastic scattering events. The forward peaked nature of high energy electron elastic scattering is relatively trivial to accurately describe in plasma simulations. However, for lower energy collisions, which produce near isotropic or backward peaked differential cross sections, there is not a strong consensus among the plasma modeling community on how to best describe these angular scattering trends. Here, in this study, we propose a systematic method to approximate the aforementioned non-trivial angular scattering behavior with a formula that can be readily implemented in particle-in-cell (PIC) and/or MC plasma simulation codes. The present approach is specifically applied to fusion relevant atomic hydrogen and helium, as well as for molecular hydrogen, and results are also applicable to the atomic isotopes and homonuclear molecular isotopologues of these species. Comparisons between the present angular distribution function and benchmark scattering data were used to validate the proposed models. In addition, two-term Boltzmann calculations and PIC direct simulation MC simulations revealed that the proposed angular distribution function is accurate, agreeing very well with benchmark convergent close-coupling scattering calculations, and electron transport measurements. These studies confirmed that the present angular distribution function model can be utilized without the need of renormalization to the momentum transfer cross section (as opposed to using the elastic scattering integrated cross section), which has been suggested by several studies in order to correct for deficient angular scattering models, and to agree with transport measurements. Hence, the present anisotropic angular scattering model can be utilized to accurately model the momentum transfer as well as the electron trajectories of elastic collisions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Ab initio nucleon-nucleus elastic scattering with chiral effective field theory uncertainties

Effective interactions for nucleon-nucleus elastic scattering from first principles require the use of the same nucleon-nucleon interaction in the structure and reaction calculations and a consistent treatment of the relevant operators at each order. Systematic investigations of the effect of truncation uncertainties of chiral nucleon-nucleon ( N N ) forces have been carried out for scattering observables in the two- and three-nucleons system as well as for bound-state properties of light nuclei. Here we extend this type of study to proton and neutron elastic scattering for 16 O and 12 C. Using the frameworks of the spectator expansion of multiple scattering theory as well as the no-core shell model, we employ one specific chiral interaction from the LENPIC collaboration and consistently calculate the leading-order effective nucleon-nucleus interaction up to the third chiral order, from which we extract elastic scattering observables. Then we apply pointwise as well as correlated uncertainty quantification for the estimation of the chiral truncation error. We calculate and analyze proton elastic scattering observables for 16 O and 12 C as well as neutron elastic scattering observables for 12 C between 65 and 185 MeV projectile kinetic energy. We find in all cases qualitatively similar results for the chiral truncation uncertainties as in few-body systems and assess them using similar diagnostic tools. The order-by-order convergence of the scattering observables for 16 O and 12 C is very reasonable around 100 MeV, while for higher energies the chiral expansion parameter becomes too large for convergence. Comparing proton and neutron differential cross sections reveals that their description is comparable up to around 100 MeV. Here, we also find a nearly perfect correlation between the differential cross section for neutron scattering and the N N Wolfenstein amplitudes for small momentum transfers. The diagnostic tools for studying order-by-order convergence in observables in few-body systems can be employed for observables in nucleon-nucleus scattering with only minor modifications provided the momentum scale in the problems is not too large. We also find that the chiral N N interaction on which our study is based gives a very good description of differential cross sections for 16 O and 12 C as low as 65-MeV projectile energy, particularly in the forward direction. In addition, the very forward direction of the neutron differential cross section mirrors the behavior of the N N interaction amazingly well.

6 ≤ A ≤ 19↗

DOE OSTI2013

Proton beam scattering system optimization for clinical and research applications

Purpose: To develop and test a method for optimizing and constructing a dual scattering system in passively scattered proton therapy. Methods: A beam optics optimization algorithm was developed to optimize the thickness of the first scatterer (S1) and the profile (of both the high-Z material and Lexan) of the second scatterer (S2) to deliver a proton beam matching a given set of parameters, including field diameter, fluence, flatness, and symmetry. A new manufacturing process was also tested that allows the contoured second scattering foil to be created much more economically and quickly using Cerrobend casting. Two application-specific scattering systems were developed and tested using both experimental and Monte Carlo techniques to validate the optimization process described. Results: A scattering system was optimized and constructed to deliver large uniform irradiations of radiobiology samples at low dose rates. This system was successfully built and tested using film and ionization chambers. The system delivered a uniform radiation field of 50 cm diameter (to a dose of {+-}7% of the central axis) while the depth dose profile could be tuned to match the specifications of the particular investigator using modulator wheels and range shifters. A second scattering system for intermediate field size (4 cm < diameter < 10 cm) stereotactic radiosurgery and radiation therapy (SRS and SRT) treatments was also developed and tested using GEANT4. This system improved beam efficiency by over 70% compared with existing scattering systems while maintaining field flatness and depth dose profile. In both cases the proton range uniformity across the radiation field was maintained, further indicating the accuracy of the energy loss formalism in the optimization algorithm. Conclusions: The methods described allow for rapid prototyping of scattering foils to meet the demands of both research and clinical beam delivery applications in proton therapy.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2017

Development of an ejecta particle size measurement diagnostic based on Mie scattering

The goal of this work is to determine the feasibility of extracting the size of particles ejected from shocked metal surfaces (ejecta) from the angular distribution of light scattered by a cloud of such particles. The basis of the technique is the Mie theory of scattering, and implicit in this approach are the assumptions that the scattering particles are spherical and that single scattering conditions prevail. The meaning of this latter assumption, as far as experimental conditions are concerned, will become clear later. The solution to Maxwell’s equations for spherical particles illuminated by a plane electromagnetic wave was derived by Gustav Mie more than 100 years ago, but several modern treatises discuss this solution in great detail. The solution is a complicated series expansion of the scattered electric field, as well as the field within the particle, from which the total scattering and absorption cross sections as well as the angular distribution of scattered intensity can be calculated numerically. The detailed nature of the scattering is determined by the complex index of refraction of the particle material as well as the particle size parameter, x, which is the product of the wavenumber of the incident light and the particle radius, i.e. x = 2rπ= λ. Figure 1 shows the angular distribution of scattered light for different particle size parameters and two orthogonal incident light polarizations as calculated using the Mie solution. It is obvious that the scattering pattern is strongly dependent on the particle size parameter, becoming more forward-directed and less polarizationdependent as the particle size parameter increases. This trend forms the basis for the diagnostic design.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2017

Scattering matrix for the interaction between solar acoustic waves and sunspots. I. Measurements

Assessing the interaction between solar acoustic waves and sunspots is a scattering problem. The scattering matrix elements are the most commonly used measured quantities to describe scattering problems. We use the wavefunctions of scattered waves of NOAAs 11084 and 11092 measured in the previous study to compute the scattering matrix elements, with plane waves as the basis. The measured scattered wavefunction is from the incident wave of radial order n to the wave of another radial order n′, for n=0--5. For a time-independent sunspot, there is no mode mixing between different frequencies. An incident mode is scattered into various modes with different wavenumbers but the same frequency. Working in the frequency domain, we have the individual incident plane-wave mode, which is scattered into various plane-wave modes with the same frequency. This allows us to compute the scattering matrix element between two plane-wave modes for each frequency. Each scattering matrix element is a complex number, representing the transition from the incident mode to another mode. The amplitudes of diagonal elements are larger than those of the off-diagonal elements. The amplitude and phase of the off-diagonal elements are detectable only for n−1⩽n{sup ′}⩽n+1 and −3Δk⩽δk{sub x}⩽3Δk, where δk{sub x} is the change in the transverse component of the wavenumber and Δk = 0.035 rad Mm{sup −1}.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2023

Effect of Scattering Angle on DSCOVR/EPIC Observations

The Earth Polychromatic Imaging Camera (EPIC) on the Deep Space Climate Observatory (DSCOVR) routinely captures reflected radiation from the whole sunlit side of the Earth in the near backward direction to monitor the changing planet. The instrument had routinely operated until 27 June 2019, when the spacecraft was placed in an extended safe hold due to degradation of an inertial navigation unit. DSCOVR returned to full operations on 2 March 2020. Since then, the range of scattering angles between the incident sunlight and sensor direction has been larger than before and the largest scattering angle reaches ∼178°, only 2° from perfect backscattering, proving a unique opportunity to study the top-of-atmosphere (TOA) reflectance under such extreme conditions. In the paper, we compare EPIC global spectral reflectances in 2021–2016. We found that there are four occasions when the scattering angle reaches about 178° and associated with them enhanced global daily average spectral reflectances in 2021. The scattering angle related reflectance enhancements are not found in 2016 data when the maximum scattering angle is about 174.5°. CERES data do not show such occasions in global daily reflected shortwave flux. As a result, those enhanced reflectance occasions are primarily due to the change in scattering angle. The enhancement due to changes in scattering angle depends strongly on wavelength, primarily because of wavelength dependence of cloud scattering phase function. Radiative transfer calculations show that the change in scattering angles has the largest impact on reflectance in the red and NIR channels at 680 nm and 780 nm and the smallest influence on reflectance in the UV channel at 388 nm, consistent with EPIC observations. The change of global average cloud amount also plays an important role in the reflectance enhancement. The influence of the cloud effect depends on whether the change is in phase or not with the change of scattering angle.

DSCOVR EPIC↗

NASA NTRS2018

"Independent" and "Dependent" Scattering by Particles in a Multi-Particle Group

The terms "independent" and "dependent" scattering are ubiquitous in the phenomenological discipline of light scattering by particulate media. Yet there is a wide range of ad hoc definitions of these terms, many of which are vague and conceptually inconsequential. In this paper we perform a first-principles analysis of these terms based on the rigorous volume-integral-equation formulation of electromagnetic scattering. We argue that scattering by a multi-particle group can be called independent if certain optical observables for the entire group can be expressed in appropriate single-particle observables. Otherwise one deals with the dependent scattering regime. The prime (and perhaps the only) examples of independent scattering are scattering scenarios described by the first-order scattering approximation and the first-principles radiative transfer theory.

Particle scattering↗

DOE OSTI2025

Molecular beam scattering of neon from flat jets of cold salty water

Molecular beam scattering experiments are carried out to study collisions between Ne atoms (E i = 24.3 kJ mol-1) and the surface of a cold salty water (8 m LiBr(aq), 230 K) flat jet. Translational energy distributions are collected as a function of scattering angle using a rotatable mass spectrometer. Impulsive scattering and thermal desorption contribute to the overall scattering distributions, but impulsive scattering dominates at all three incidence angles explored. Highly super-specular scattering is observed in the impulsive scattering channel that is attributed to anisotropic momentum transfer to the liquid surface. The thermal desorption channel exhibits a cos θ angular distribution. Compared to Ne scattering from dodecane, fractional energy loss in the impulsive scattering channel is much larger across a wide range of deflection angles. A soft-sphere model is applied to investigate the kinematics of energy transfer between the scatterer and liquid surface. Fitting to this model yields an effective surface mass of 250-60 +100 amu and internal excitation of 11.8 ± 1.6 kJ mol-1, both of which are considerably larger than for Ne/dodecane. It thus appears that energy transfer to cold salty water is more efficient than to a dodecane liquid surface, a result attributed to the extensive hydrogen-bonded network of liquid water and roughness of the liquid surface.

Yang, Walt↗

DOE OSTI2023

Evaporation and scattering of neon, methane, and water from a dodecane flat liquid jet

The evaporation and scattering of Ne, CD 4 , and D 2 O from a dodecane flat liquid jet are investigated in a molecular beam apparatus. The experiment yields translational energy distributions as a function of scattering angle by means of a rotatable mass spectrometer. In the evaporation experiments, one observes a Maxwell–Boltzmann distribution with a cos $θ$ angular distribution superimposed on a weak, isotropic background. The scattering experiments show contributions from impulsive scattering and thermal desorption. At select incident angles for the three systems, angular distributions show super-specular scattering for the impulsive scattering channel, an effect attributed to anisotropic momentum transfer to the liquid surface. The impulsive scattering channel is analyzed with a soft-sphere model to explore energy transfer between the scatterer and liquid as a function of deflection angle. In conclusion, compared to Ne scattering, the polyatomic gases exhibit more thermal desorption and, in the impulsive scattering channel, a higher degree of internal excitation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Measurements of attenuation and scattering properties of water-based liquid scintillator

Water-based liquid scintillator (WbLS) is a hybrid detector medium which has been proposed as a fill material for future large-volume rare event searches, including ν − detectors. This family of scintillating suspensions promises waterlike attenuation and scattering behavior, while offering scintillation light yield for sub-Cherenkov-threshold events. These features may allow for improved vertex and energy resolution, and concomitant improvements in background rejection and particle identification. While subscale measurements of light yield, timing and pulse-shape have been performed in several WbLS formulations, detailed measurements of the attenuation and scattering properties of WbLS remain a critical, unresolved step along the pathway to deployment. In pursuit of a better understanding of these parameters, a “long-arm” attenuation and scattering instrument has been developed at Lawrence Livermore National Laboratory, dubbed LASE (Livermore Attenuation and Scattering Experiment). Optical property measurements have been performed using LASE in DI water, Gd-water, WbLS and Gd-WbLS. Measurements of the optical properties of WbLS (1% LAB-PPO scintillator) provided by Brookhaven National Laboratory demonstrate an attenuation minimum of 4.87 × 10 − 4 ± 1.42 × 10 − 4 cm − 1 at 450 nm, while Gd-loaded WbLS (1% LAB-PPO scintillator, 0.1% Gd) demonstrated a minimum attenuation coefficient of 4.77 × 10 − 4 ± 1.27 × 10 − 4 cm − 1 at the same wavelength. Measurements of scattering in WbLS show a scattering coefficient of 3.39 × 10 − 4 ± 1.98 × 10 − 5 cm − 1 at the 450 nm attenuation minimum, while Gd-WbLS has a scattering coefficient of 2.8 × 10 − 4 ± 1.63 × 10 − 5 cm − 1 at that wavelength. These scattering and attenuation coefficients are significantly larger than measured for DI water at similar wavelengths ( 8.8 × 10 − 5 cm − 1 ± 3.8 × 10 − 5 attenuation at 430 nm, 2.2 × 10 − 5 ± 1.28 × 10 − 6 cm − 1 scattering). Though the WbLS attenuation and scattering coefficients measured are lower than the corresponding values published for ultrapure oils and LAB-PPO liquid scintillator, the material tested attenuates significantly more than water across the blue and green portions of the visible spectrum. This attenuation is markedly stronger at wavelengths under 430 nm. Published by the American Physical Society 2025

47 OTHER INSTRUMENTATION↗

DOE OSTI2011

A software-based x-ray scatter correction method for breast tomosynthesis

Purpose: To develop a software-based scatter correction method for digital breast tomosynthesis (DBT) imaging and investigate its impact on the image quality of tomosynthesis reconstructions of both phantoms and patients. Methods: A Monte Carlo (MC) simulation of x-ray scatter, with geometry matching that of the cranio-caudal (CC) view of a DBT clinical prototype, was developed using the Geant4 toolkit and used to generate maps of the scatter-to-primary ratio (SPR) of a number of homogeneous standard-shaped breasts of varying sizes. Dimension-matched SPR maps were then deformed and registered to DBT acquisition projections, allowing for the estimation of the primary x-ray signal acquired by the imaging system. Noise filtering of the estimated projections was then performed to reduce the impact of the quantum noise of the x-ray scatter. Three dimensional (3D) reconstruction was then performed using the maximum likelihood-expectation maximization (MLEM) method. This process was tested on acquisitions of a heterogeneous 50/50 adipose/glandular tomosynthesis phantom with embedded masses, fibers, and microcalcifications and on acquisitions of patients. The image quality of the reconstructions of the scatter-corrected and uncorrected projections was analyzed by studying the signal-difference-to-noise ratio (SDNR), the integral of the signal in each mass lesion (integrated mass signal, IMS), and the modulation transfer function (MTF). Results: The reconstructions of the scatter-corrected projections demonstrated superior image quality. The SDNR of masses embedded in a 5 cm thick tomosynthesis phantom improved 60%-66%, while the SDNR of the smallest mass in an 8 cm thick phantom improved by 59% (p < 0.01). The IMS of the masses in the 5 cm thick phantom also improved by 15%-29%, while the IMS of the masses in the 8 cm thick phantom improved by 26%-62% (p < 0.01). Some embedded microcalcifications in the tomosynthesis phantoms were visible only in the scatter-corrected reconstructions. The visibility of the findings in two patient images was also improved by the application of the scatter correction algorithm. The MTF of the images did not change after application of the scatter correction algorithm, indicating that spatial resolution was not adversely affected. Conclusions: Our software-based scatter correction algorithm exhibits great potential in improving the image quality of DBT acquisitions of both phantoms and patients. The proposed algorithm does not require a time-consuming MC simulation for each specific case to be corrected, making it applicable in the clinical realm.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2014

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↗

DOE OSTI2014

SU-E-I-64: X-Ray Coherent Scatter Mammography Simulation

Purpose: Conventional mammography has poor contrast between healthy and cancerous tissues due to the small difference in attenuation properties. Coherent scatter potentially provides more information because interference of coherently scattered radiation depends on the average intermolecular spacing, and can be used to characterize tissue types. However, typical coherent scatter analysis techniques are not compatible with rapid low dose screening modalities. Coherent scatter slot scan imaging is a novel imaging technique which provides new information with higher contrast. In this work a simulation was performed of coherent scatter slot scan imaging to assess its performance and provide system optimization. Methods: In coherent scatter imaging, the coherent scatter is exploited using a conventional slot scan mammography system with anti-scatter grids tilted at the characteristic angle of cancerous tissues. A Monte Carlo simulation was used to simulate the coherent scatter imaging. System optimization was performed across several parameters, including source voltage and filtration, tilt angle, source and grid distances, grid ratio, and shielding geometry. Results: The contrast was high and increased as the grid tilt angle was increased beyond the characteristic angle for the modeled carcinoma, and as the source voltage increased. Source filtration did not have a significant effect on contrast. Increasing grid ratio improved contrast at the expense of decreasing the signal to noise ratio (SNR). As the tumor size is decreased, additional shields to block the fat signal are necessary in order to detect the cancer signal. Conclusion: Coherent scatter analysis using a wide slot setup is promising as an enhancement for screening mammography. This work was supported by National Institutes of Health, # 7 R01EB009715.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2014

SU-E-T-439: An Improved Formula of Scatter-To-Primary Ratio for Photon Dose Calculation

Purpose: Scatter-to-primary ratio (SPR) is an important dosimetric quantity that describes the contribution from the scatter photons in an external photon beam. The purpose of this study is to develop an improved analytical formula to describe SPR as a function of circular field size (r) and depth (d) using Monte Carlo (MC) simulation. Methods: MC simulation was performed for Mohan photon spectra (Co-60, 4, 6, 10, 15, 23 MV) using EGSNRC code. Point-spread scatter dose kernels in water are generated. The scatter-to-primary ratio (SPR) is also calculated using MC simulation as a function of field size for circular field size with radius r and depth d. The doses from forward scatter and backscatter photons are calculated using a convolution of the point-spread scatter dose kernel and by accounting for scatter photons contributing to dose before (z' d) reaching the depth of interest, d, where z' is the location of scatter photons, respectively. The depth dependence of the ratio of the forward scatter and backscatter doses is determined as a function of depth and field size. Results: We are able to improve the existing 3-parameter (a, w, d0) empirical formula for SPR by introducing depth dependence for one of the parameter d0, which becomes 0 for deeper depths. The depth dependence of d0 can be directly calculated as a ratio of backscatter-to-forward scatter doses for otherwise the same field and depth. With the improved empirical formula, we can fit SPR for all megavoltage photon beams to within 2%. Existing 3-parameter formula cannot fit SPR data for Co-60 to better than 3.1%. Conclusion: An improved empirical formula is developed to fit SPR for all megavoltage photon energies to within 2%.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

Fully iterative scatter corrected digital breast tomosynthesis using GPU-based fast Monte Carlo simulation and composition ratio update

Purpose: In digital breast tomosynthesis (DBT), scatter correction is highly desirable, as it improves image quality at low doses. Because the DBT detector panel is typically stationary during the source rotation, antiscatter grids are not generally compatible with DBT; thus, a software-based scatter correction is required. This work proposes a fully iterative scatter correction method that uses a novel fast Monte Carlo simulation (MCS) with a tissue-composition ratio estimation technique for DBT imaging. Methods: To apply MCS to scatter estimation, the material composition in each voxel should be known. To overcome the lack of prior accurate knowledge of tissue composition for DBT, a tissue-composition ratio is estimated based on the observation that the breast tissues are principally composed of adipose and glandular tissues. Using this approximation, the composition ratio can be estimated from the reconstructed attenuation coefficients, and the scatter distribution can then be estimated by MCS using the composition ratio. The scatter estimation and image reconstruction procedures can be performed iteratively until an acceptable accuracy is achieved. For practical use, (i) the authors have implemented a fast MCS using a graphics processing unit (GPU), (ii) the MCS is simplified to transport only x-rays in the energy range of 10–50 keV, modeling Rayleigh and Compton scattering and the photoelectric effect using the tissue-composition ratio of adipose and glandular tissues, and (iii) downsampling is used because the scatter distribution varies rather smoothly. Results: The authors have demonstrated that the proposed method can accurately estimate the scatter distribution, and that the contrast-to-noise ratio of the final reconstructed image is significantly improved. The authors validated the performance of the MCS by changing the tissue thickness, composition ratio, and x-ray energy. The authors confirmed that the tissue-composition ratio estimation was quite accurate under a variety of conditions. Our GPU-based fast MCS implementation took approximately 3 s to generate each angular projection for a 6 cm thick breast, which is believed to make this process acceptable for clinical applications. In addition, the clinical preferences of three radiologists were evaluated; the preference for the proposed method compared to the preference for the convolution-based method was statistically meaningful (p < 0.05, McNemar test). Conclusions: The proposed fully iterative scatter correction method and the GPU-based fast MCS using tissue-composition ratio estimation successfully improved the image quality within a reasonable computational time, which may potentially increase the clinical utility of DBT.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

Library based x-ray scatter correction for dedicated cone beam breast CT

Purpose: The image quality of dedicated cone beam breast CT (CBBCT) is limited by substantial scatter contamination, resulting in cupping artifacts and contrast-loss in reconstructed images. Such effects obscure the visibility of soft-tissue lesions and calcifications, which hinders breast cancer detection and diagnosis. In this work, we propose a library-based software approach to suppress scatter on CBBCT images with high efficiency, accuracy, and reliability. Methods: The authors precompute a scatter library on simplified breast models with different sizes using the GEANT4-based Monte Carlo (MC) toolkit. The breast is approximated as a semiellipsoid with homogeneous glandular/adipose tissue mixture. For scatter correction on real clinical data, the authors estimate the breast size from a first-pass breast CT reconstruction and then select the corresponding scatter distribution from the library. The selected scatter distribution from simplified breast models is spatially translated to match the projection data from the clinical scan and is subtracted from the measured projection for effective scatter correction. The method performance was evaluated using 15 sets of patient data, with a wide range of breast sizes representing about 95% of general population. Spatial nonuniformity (SNU) and contrast to signal deviation ratio (CDR) were used as metrics for evaluation. Results: Since the time-consuming MC simulation for library generation is precomputed, the authors’ method efficiently corrects for scatter with minimal processing time. Furthermore, the authors find that a scatter library on a simple breast model with only one input parameter, i.e., the breast diameter, sufficiently guarantees improvements in SNU and CDR. For the 15 clinical datasets, the authors’ method reduces the average SNU from 7.14% to 2.47% in coronal views and from 10.14% to 3.02% in sagittal views. On average, the CDR is improved by a factor of 1.49 in coronal views and 2.12 in sagittal views. Conclusions: The library-based scatter correction does not require increase in radiation dose or hardware modifications, and it improves over the existing methods on implementation simplicity and computational efficiency. As demonstrated through patient studies, the authors’ approach is effective and stable, and is therefore clinically attractive for CBBCT imaging.

60 APPLIED LIFE SCIENCES↗