Search NASA⌕ Search

SEARCH · Search NASA

Results for “developable planes”

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 1,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Effects of Particle Mixing and Gravitational Settling on Charge Transport in Carbon Flow-Electrode Cells

Flowable carbon slurries are actively studied and under development for charge transport in various electrochemical systems including flow capacitors, capacitive deionization cells, semi-solid flow batteries, and lithium extraction. However, much less is known about in operando slurry flow dynamics and their corresponding effect on charge transport. We performed an experimental study of mixing and settling dynamics of slurry electrodes within an electrochemical flow cell during continuous operations. The electrochemical cell consisted of two horizontal co-flowing channels, separated by a cation-exchange membrane (CEM). Here we used high-speed optical imaging of planes parallel to gravity and simultaneous electrochemical measurements. At low flow rates, dense yet dynamic particle beds formed on the bottom electrode in each channel, which unexpectedly yielded the highest currents. This approach enables the operation of the flow cell at low system-average particle concentrations while leveraging gravity-driven particle settling to locally enhance carbon concentrations precisely at the current collector sites. Conversely, high flow rates were characterized by thin particle beds and well-mixed particle flows. In the latter case, the electrodes in closest proximity (located on either side of the CEM) achieved a current higher than the other electrode pairs. The observations have implications for slurry control and electrode designs in electrochemical systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

A Particle-Based Conditional Sampling Scheme for the Simulation of Transport in Fractured Rock With Diffusion Into Stagnant Water and Rock Matrix

In situ experiments and field-scale characterization studies have pointed out that, in fractured crystalline media, groundwater flow is highly channelized. This implies that, at the scale of a single fracture, only part of the fracture surface area is in contact with flowing water, while the rest of in-plane water is essentially stagnant and can be accessed by solutes via molecular diffusion. Despite their importance for contaminant retention, to date, there are no numerical or analytical approaches that could be used to assess the implication of stagnant water zones on solute transport in realistic large-scale Discrete Fracture Network-based models. In this study, we present an efficient and flexible algorithm for the simulation of transport in fractured rock with diffusion into stagnant water and rock matrix. The algorithm is a generalization of a previously developed numerical framework for time domain particle tracking in sparsely fractured rock. The key of the generalization is that total time in fracture (τ f ) is first evaluated using a Monte Carlo sampling and then a second sampling is performed conditioned on τ f . The algorithm has been successfully validated against existing independent solutions and the implication of diffusion into stagnant water and secondary diffusion into the matrix has been assessed for a realistic modeling scenario. The results show that, due to diffusion into stagnant water, contaminants are more strongly retarded. This increased retention is more significant for sorbing species, as a larger number of sorption sites is accessible. A high sensitivity to the flowing channel/stagnant water zone geometry has also been observed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Switching the spin cycloid in BiFeO3 with an electric field

Abstract Bismuth ferrite (BiFeO 3 ) is a multiferroic material that exhibits both ferroelectricity and canted antiferromagnetism at room temperature, making it a unique candidate in the development of electric-field controllable magnetic devices. The magnetic moments in BiFeO 3 are arranged into a spin cycloid, resulting in unique magnetic properties which are tied to the ferroelectric order. Previous understanding of this coupling has relied on average, mesoscale measurements. Using nitrogen vacancy-based diamond magnetometry, we observe the magnetic spin cycloid structure of BiFeO 3 in real space. This structure is magnetoelectrically coupled through symmetry to the ferroelectric polarization and this relationship is maintained through electric field switching. Through a combination of in-plane and out-of-plane electrical switching, coupled with ab initio studies, we have discovered that the epitaxy from the substrate imposes a magnetoelastic anisotropy on the spin cycloid, which establishes preferred cycloid propagation directions. The energy landscape of the cycloid is shaped by both the ferroelectric degree of freedom and strain-induced anisotropy, restricting the spin spiral propagation vector to changes to specific switching events.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The next-generation magnetic recoil spectrometer (MRSnext) on OMEGA and NIF for diagnosing ion temperature, yield, areal density, and alpha heating

The next-generation magnetic recoil spectrometer (MRSnext) is being designed to replace the current MRS at the National Ignition Facility and OMEGA for measurements of the neutron spectrum from an inertial confinement fusion implosion. The MRSnext will provide a far-superior performance and faster data turnaround than the current MRS systems, i.e., a 2× and 6× improvement in energy resolution at the NIF and OMEGA, respectively, and 20× improvement in data turnaround time. The substantially improved performance of the MRSnext is enabled by using electromagnets that provide a short focal plane (12–16 cm) and unprecedented flexibility for a wide range of applications. In addition to being able to measure neutron yield, apparent ion temperature, areal density, and plasma-flow velocity over a wide range of yields, the NIF MRSnext will be able to directly, uniquely assess the alpha heating of the fuel ions through measurements of the alpha knock-on tail in the neutron spectrum. The goal is to implement a radiation-hard electronic detection system capable of providing rapid data acquisition and analysis. The development of the MRSnext will also set the foundation for the more advanced, time-resolving MRSt and serve as a testbed for its implementation on the NIF.

47 OTHER INSTRUMENTATION↗
DOE OSTI2013

THE THOMSON SURFACE. II. POLARIZATION

The solar corona and heliosphere are visible via sunlight that is Thomson-scattered off free electrons, yielding a radiance against the celestial sphere. In this second part of a three-article series, we discuss linear polarization of this scattered light parallel and perpendicular to the plane of scatter in the context of heliospheric imaging far from the Sun. The difference between these two radiances (pB) varies quite differently with scattering angle, compared to the sum that would be detected in unpolarized light (B). The difference between these two quantities has long been used in a coronagraphic context for background subtraction and to extract some three-dimensional information about the corona; we explore how these effects differ in the wider-field heliospheric imaging case where small-angle approximations do not apply. We develop an appropriately simplified theory of polarized Thomson scattering in the heliosphere, discuss signal-to-noise considerations, invert the scattering equations analytically to solve the three-dimensional object location problem for small objects, discuss exploiting polarization for background subtraction, and generate simple forward models of several classes of heliospheric feature. We conclude that pB measurements of heliospheric material are much more localized to the Thomson surface than are B measurements, that the ratio pB/B can be used to track solar wind features in three dimensions for scientific and space weather applications better in the heliosphere than corona, and that,as an independent measurement of background signal, pB measurements may be used to reduce the effect of background radiances including the stably polarized zodiacal light.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Diverted Negative Triangularity plasmas on DIII–D: The benefit of high confinement without the liability of an edge pedestal

Diverted discharges at negative triangularity on the DIII-D tokamak sustain normalized confinement and pressure levels typical of standard H-mode scenarios (H 98y,2 ≅1, β N ≅3) without developing an edge pressure pedestal, despite the auxiliary power far exceeding the L → H power threshold expected from conventional scaling laws. The power degradation of confinement is substantially weaker than the ITER-89P scaling, resulting in a confinement factor that improves with increasing auxiliary power. Furthermore, the absence of the edge pedestal is beneficial in several aspects, such as eliminating the need for active mitigation or suppression of edge localized modes, low impurity retention and a reconstructed scrape-off layer heat flux width at the mid-plane that exceeds the ITPA multi-machine scaling law by up to 50%. Together with technological advantages granted by placing the divertor at larger radii, plasmas at Negative Triangularity without an edge pedestal feature both core confinement and power handling characteristics that are potentially suitable for operation in future fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Experimental and Monte Carlo studies of fluence corrections for graphite calorimetry in low- and high-energy clinical proton beams

Purpose: The aim of this study was to determine fluence corrections necessary to convert absorbed dose to graphite, measured by graphite calorimetry, to absorbed dose to water. Fluence corrections were obtained from experiments and Monte Carlo simulations in low- and high-energy proton beams. Methods: Fluence corrections were calculated to account for the difference in fluence between water and graphite at equivalent depths. Measurements were performed with narrow proton beams. Plane-parallel-plate ionization chambers with a large collecting area compared to the beam diameter were used to intercept the whole beam. High- and low-energy proton beams were provided by a scanning and double scattering delivery system, respectively. A mathematical formalism was established to relate fluence corrections derived from Monte Carlo simulations, using the FLUKA code [A. Ferrari et al., “FLUKA: A multi-particle transport code,” in CERN 2005-10, INFN/TC 05/11, SLAC-R-773 (2005) and T. T. Böhlen et al., “The FLUKA Code: Developments and challenges for high energy and medical applications,” Nucl. Data Sheets 120, 211–214 (2014)], to partial fluence corrections measured experimentally. Results: A good agreement was found between the partial fluence corrections derived by Monte Carlo simulations and those determined experimentally. For a high-energy beam of 180 MeV, the fluence corrections from Monte Carlo simulations were found to increase from 0.99 to 1.04 with depth. In the case of a low-energy beam of 60 MeV, the magnitude of fluence corrections was approximately 0.99 at all depths when calculated in the sensitive area of the chamber used in the experiments. Fluence correction calculations were also performed for a larger area and found to increase from 0.99 at the surface to 1.01 at greater depths. Conclusions: Fluence corrections obtained experimentally are partial fluence corrections because they account for differences in the primary and part of the secondary particle fluence. A correction factor, F(d), has been established to relate fluence corrections defined theoretically to partial fluence corrections derived experimentally. The findings presented here are also relevant to water and tissue-equivalent-plastic materials given their carbon content.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Improved Regional Moment Tensor Inversion for Moderately Large Earthquakes in the Western United States Using a 3D Earth Model Based on Full Waveform Tomography

The nature of seismic sources for moderately large (moment magnitude, M w 5.0–6.5) events are commonly characterized by their moment tensor (MT) solutions and obtained by inversion of regional distance (200–1600 km) long‐period (20–50 s) waveforms. Regional MT estimates are often calculated from average plane‐layered, one‐dimensional (1D) velocity models. However, 1D model calculations can produce misfits in the arrival times and waveform shapes that introduce errors, particularly at longer distances or for shorter periods, which are necessary for analyzing lower magnitude events. Approximate Earth models (e.g., 1D) representing broad areas may be inadequate, particularly in the crust and uppermost mantle of tectonically complex regions. In this study, we show how a three‐dimensional (3D) Earth model obtained from full waveform inversion tomography can improve waveform fits and decrease phase errors. We developed a platform and workflow to perform routine 3D MT inversions and inverted MTs for 25 earthquakes in the western United States and seven nuclear explosions using an average 1D and a recent 3D Earth model, WUS256 (Rodgers et al., 2022). Using the 3D model improves waveform fits (variance reduction and phase time shifts) compared with the 1D model, and the 3D MT solutions are stable across large distances. This study shows that 3D models obtained from full waveform tomography can improve MTs and source characterization especially at far regional distances (>800 km).

Geosciences↗
DOE OSTI2019

The Measurement of Electrical and Thermal Conductivities in Warm Dense Matter

This document serves as the Final Technical Report for grant DE-SC0012333, The Measurement of Electrical and Thermal Conductivities in Warm Dense Matter . This work proposed to make the first systematic measurements of the reflectivity of warm dense matter (WDM) at multiple wavelengths by controlling and measuring the experimental factors necessary to use the reflectivity measurements to determine the electrical and thermal conductivity as a function of temperature and, to a lesser degree, density. A new technique was developed to measure the real and imaginary components of the electrical conductivity by making only scalar measurements of the relative reflectivity different angles. A series of preliminary experiments were designed to use the second, third, and fourth harmonics of the 800 nm driving laser. Reflectivity at each harmonic for two different angles would provide the scalar information needed to extract the real and imaginary parts of the index of refraction. Experiments were made using continuous wave lasers to determine the limits of accuracy on the measured reflectance for surrogate plane mirrors. It was determined, that while the technique was sound, the practical application of it was limited for the conditions on Scarlet. Due to significant challenges with the Scarlet Laser Facility at the time, progress was delayed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Blunt Impact Brain Injury using Cellular Injury Criterion

The Advanced Combat Helmet (\ACH") military specification (\mil-spec") requires a helmeted magnesium (\Mg") Department of Transportation (\DOT") headform be dropped vertically, with an impact speed of 3.1 m/s (10 ft/s), onto a steel hemispherical target. The pass/fail criteria are based on translational acceleration (150 G) alone, absent of any rotational component. Without a rotational component, the specification's injury risk application is limited to skull fracture and peripheral hematomas (subdural, subarachnoid), since this translational acceleration injury risk assessment is based on the Wayne State Tolerance Curve (\WSTC"). To provide a more comprehensive view of injury for the entire brain, an alternative approach is needed. To meet this need, we worked with a larger group called PANTHER, a collaboration between national laboratories, industry, and academia. Collaborations specific to research and results presented here come from efforts led by Mr. Ron Szalkowski and Mr. Sushant Malave, Ms. Alice Fawzi, and Dr. Christian Franck. We have developed a prototypical injury risk criterion based on the neuronal response to abrupt changes in general motion (translation, rotation, or both). The cellular-based mild traumatic brain injury (\cbmTBI") criterion utilizes both the strain and strain rate of brain tissue to account for the stretch and rate of stretch that occurs throughout the brain as a result of blunt impact to the head. We conducted physical experiments of an ACH-fitted magnesium headform, which produced repeatable headform peak accelerations. Then, we developed a simulation of the experiment, and validated the simulation output with the experimental data. We then substituted the magnesium headform with a human headform, consisting of skin, muscle, bone, gray matter, white matter, cerebral-spinal fluid, membranes, vasculature, intravertebral discs, airway and sinus. We quantified brain injury risk using the cbmTBI criterion, using the current mil-spec test and a modified test. The modified mil-spec test used an inclined anvil target that was located posterior to the crown of the helmet in the axial plane. While the current mil-spec test produced brain deformation from head translation alone, the modified test produced brain deformation from head translation and rotation, which is closer to most real world and combat theater impacts (e.g., such as occur in tertiary blast exposure). Compared to the current mil-spec test, the modified test produced elevated strains in the human digital twin. These data, mapped to the cbmTBI criterion, suggest increased injury risk for blunt impacts that cause rotation and translation, rather than just translation alone. Moreover, these data may lead to a rotational performance metric, which is rooted in the biology and pathology of the brain's response to impact and blast, and which should be used to improve next-generation helmet designs.

42 ENGINEERING↗
DOE OSTI2017

Growth Oscillatory Zoning in Erythrite, Ideally Co3(AsO4)2·8H2O: Structural Variations in Vivianite-Group Minerals

The crystal structure of an oscillatory zoned erythrite sample from Aghbar mine, Bou Azzer, Morocco, was refined using synchrotron high-resolution powder X-ray diffraction (HRPXRD) data, Rietveld refinement, space group C2/m, and Z = 2. The crystal contains two sets of oscillatory zones that appear to have developed during epitaxial growth. The unit-cell parameters obtained are a = 10.24799(3) Å, b = 13.42490(7) Å, c = 4.755885(8) Å, β = 105.1116(3)°, and V = 631.680(4) Å3. The empirical formula for erythrite, obtained with electron-probe micro-analysis (EPMA), is [Co2.78Zn0.11Ni0.07Fe0.04]Σ3.00(AsO4)2·8H2O. Erythrite belongs to the vivianite-type structure that contains M1O2(H2O)4 octahedra and M22O6(H2O)4 octahedral dimers that are linked by TO4 (T5+ = As or P) tetrahedra to form complex layers parallel to the (010) plane. These layers are connected by hydrogen bonds. The average [6] = 2.122(1) Å and average [6] = 2.088(1) Å. With space group C2/m, there are two solid solutions: M3(AsO4)2·8H2O and M3(PO4)2·8H2O where M2+ = Mg, Fe, Co, Ni, or Zn. In these As- and P-series, using data from this study and from the literature, we find that their structural parameters evolve linearly with V and in a nearly parallel manner despite of the large difference in size between P5+ (0.170 Å) and As5+ (0.355 Å) cations. Average [4], [6], and [6] distances increase linearly with V. The average distance is affected by M atoms, whereas the average distance is unaffected because it contains shorter and stronger P–O bonds. Although As- and P-series occur naturally, there is no structural reason why similar V-series vivianite-group minerals do not occur naturally or cannot be synthesized.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Enhanced Spontaneous Polarization in Ultrathin SnTe Films with Layered Antipolar Structure

2D SnTe films with a thickness of as little as 2 atomic layers (ALs) have recently been shown to be ferroelectric with in–plane polarization. Remarkably, they exhibit transition temperatures (T c ) much higher than that of bulk SnTe. Here, combining molecular beam epitaxy, variable temperature scanning tunneling microscopy, and ab initio calculations, the underlying mechanism of the Tc enhancement is unveiled, which relies on the formation of γ–SnTe, a van der Waals orthorhombic phase with antipolar inter–layer coupling in few–AL thick SnTe films. In this phase, 4n – 2 AL (n = 1, 2, 3…) thick films are found to possess finite in–plane polarization (space group Pmn2 1 ), while 4n AL thick films have zero total polarization (space group Pnma). Above 8 AL, the γ–SnTe phase becomes metastable, and can convert irreversibly to the bulk rock salt phase as the temperature is increased. This finding unambiguously bridges experiments on ultrathin SnTe films with predictions of robust ferroelectricity in GeS–type monochalcogenide monolayers. Here, the observed high transition temperature, together with the strong spin–orbit coupling and van der Waals structure, underlines the potential of atomically thin γ–SnTe films for the development of novel spontaneous polarization–based devices.

2d ferroelectrics↗
DOE OSTI2022

A reaction–diffusion model for grayscale digital light processing 3D printing

We report that digital light processing (DLP) 3D printing is an additive manufacturing process that utilizes light patterns to photopolymerize a liquid resin into a solid. Due to the accuracy of modern digital micromirror devices (DMD) and recent advances in resin chemistry, it is now possible to create functionally graded structures using different light intensity values, also known as grayscale DLP (g-DLP). Different intensities of light lead to differences in the polymer crosslinking density after curing, which ultimately produces a part with gradients of material properties. However, g-DLP is a complicated process. First, the DLP printing is a highly coupled chemical and physical process that involves light propagation, chemical reactions, species diffusion, heat transfer, volume shrinkage, and changes in mechanical behaviors of the curing resin. Second, in g-DLP, light gradients create strong in plane gradients of chemical species concentrations in the curing liquid resin due to the strong dependence of light intensity on the rate of monomer crosslinking. Furthermore, light gradients through the depth create concentration gradients due to the degree of cure dependent light absorption and the use of photoabsorbers. These complex physical features of the printing process must be understood in order to properly control printing parameters such as light exposure time, printing speed, and grayscale variations to achieve accurate mechanical properties. In this paper, a photopolymerization reaction–diffusion model is developed and used in conjunction with experiments to investigate the coupled effects of light propagation, chemical reaction rates, and species diffusion during g-DLP 3D printing. The model is implemented numerically utilizing the finite difference method and simulation results are compared to experimental findings of simple printed structures. The agreement between experimental and model predictions of simple quantities of interest, such as geometric feature sizes, shows that the model can capture the overcure due to free-radical and other species diffusion during printing when grayscale patterns are employed. This model lays the groundwork for future extensions that can incorporate more complex coupled physics such as heat transfer, volume shrinkage, and material property evolution, which are critically important in utilizing g-DLP 3D printing for the fabrication of high-performance parts which excellent geometric and material property tolerances.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Binary dislocation junction formation and strength in hexagonal close-packed crystals

This work examines binary dislocation interactions, junction formation and junction strengths in hexagonal close-packed ( hcp ) crystals. Through a line-tension model and dislocation dynamics (DD) simulations, the interaction and dissociation of different sets of binary junctions are investigated involving one dislocation on the (011¯0) prismatic plane and a second dislocation on one of the following planes: (0001) basal, (11¯00) prismatic, (11¯01) primary pyramidal, or (2¯112) secondary pyramidal. Varying pairs of Burgers vectors are chosen from among the common types the basal type < a > 1/3 < 112¯0 >, prismatic type < c > <0001>, and pyramidal type < a+c > 1/3 < 112¯3¯ >. For binary interaction due to dislocation intersection, both the analytical results and DD-simulations indicate a relationship between symmetry of interaction maps and the relative magnitude of the Burgers vectors that constitute the junction. Using analytical formulae, a simple regressive model is also developed to represent the junction yield surface. The equation is treated as a degenerated super elliptical equation to quantify the aspect ratio and tilting angle. Lastly, the results provide analytical insights on binary dislocation interactions that may occur in general hcp metals.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Beyond Epitaxy: Ion Implantation as a Tool for Orbital Engineering

Manipulating electronic orbital states in quantum materials provides a powerful means of controlling their physical properties and technological functionality. Here, we demonstrate that orbital populations in strongly correlated oxide thin films can be continuously and reversibly tuned through postsynthesis He ion implantation. Using LaNiO 3 as a model system, we show that the orbital preference can be systematically adjusted from favoring in-plane d x 2 –y 2 occupation toward out-of-plane d z 2 states through precise control of ion fluence. Unlike conventional heteroepitaxial approaches that lock in orbital configurations during growth, this strain-doping technique enables continuous orbital tuning and the selective modification of specific film regions after device fabrication. We demonstrate the practical impact of this control by achieving a 7-fold enhancement in oxygen reduction reaction catalysis. This work establishes ion implantation as a powerful approach for orbital engineering that complements existing synthesis-based strategies while offering unique advantages for both basic research and device development.

Doping↗
DOE OSTI2012

Measurement of turbulent flow upstream and downstream of a circular pipe bend

We measured velocity distribution in cross sections of a fully developed turbulent pipe flow upstream and downstream of a 90 degree sign bend by synchronizing two sets of a particle image velocimetry (PIV) system. Unsteady undulation of Dean vortices formed downstream from the bend was characterized by the azimuthal position of the stagnation point found on the inner and outer sides of the bend. Linear stochastic estimation was applied to capture the upstream flow field conditioned by the azimuthal location of the stagnation point downstream from the bend. When the inner-side stagnation point stayed below (above) the symmetry plane, the conditional streamwise velocity upstream from the bend exhibited high-speed streaks extended in a quasi-streamwise direction on the outer side of the curvature above (below) the symmetry plane.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

Watching a signaling protein function in real time via 100-ps time-resolved Laue crystallography

To understand how signaling proteins function, it is necessary to know the time-ordered sequence of events that lead to the signaling state. We recently developed on the BioCARS 14-IDB beamline at the Advanced Photon Source the infrastructure required to characterize structural changes in protein crystals with near-atomic spatial resolution and 150-ps time resolution, and have used this capability to track the reversible photocycle of photoactive yellow protein (PYP) following trans-to-cis photoisomerization of its p-coumaric acid (pCA) chromophore over 10 decades of time. The first of four major intermediates characterized in this study is highly contorted, with the pCA carbonyl rotated nearly 90° out of the plane of the phenolate. A hydrogen bond between the pCA carbonyl and the Cys69 backbone constrains the chromophore in this unusual twisted conformation. Density functional theory calculations confirm that this structure is chemically plausible and corresponds to a strained cis intermediate. This unique structure is short-lived (~600 ps), has not been observed in prior cryocrystallography experiments, and is the progenitor of intermediates characterized in previous nanosecond time-resolved Laue crystallography studies. The structural transitions unveiled during the PYP photocycle include trans/cis isomerization, the breaking and making of hydrogen bonds, formation/relaxation of strain, and gated water penetration into the interior of the protein. This mechanistically detailed, near-atomic resolution description of the complete PYP photocycle provides a framework for understanding signal transduction in proteins, and for assessing and validating theoretical/computational approaches in protein biophysics.

Schotte, Friedrich↗
DOE OSTI2014

TH-A-18C-09: Ultra-Fast Monte Carlo Simulation for Cone Beam CT Imaging of Brain Trauma

Purpose: Application of cone-beam CT (CBCT) to low-contrast soft tissue imaging, such as in detection of traumatic brain injury, is challenged by high levels of scatter. A fast, accurate scatter correction method based on Monte Carlo (MC) estimation is developed for application in high-quality CBCT imaging of acute brain injury. Methods: The correction involves MC scatter estimation executed on an NVIDIA GTX 780 GPU (MC-GPU), with baseline simulation speed of ~1e7 photons/sec. MC-GPU is accelerated by a novel, GPU-optimized implementation of variance reduction (VR) techniques (forced detection and photon splitting). The number of simulated tracks and projections is reduced for additional speed-up. Residual noise is removed and the missing scatter projections are estimated via kernel smoothing (KS) in projection plane and across gantry angles. The method is assessed using CBCT images of a head phantom presenting a realistic simulation of fresh intracranial hemorrhage (100 kVp, 180 mAs, 720 projections, source-detector distance 700 mm, source-axis distance 480 mm). Results: For a fixed run-time of ~1 sec/projection, GPU-optimized VR reduces the noise in MC-GPU scatter estimates by a factor of 4. For scatter correction, MC-GPU with VR is executed with 4-fold angular downsampling and 1e5 photons/projection, yielding 3.5 minute run-time per scan, and de-noised with optimized KS. Corrected CBCT images demonstrate uniformity improvement of 18 HU and contrast improvement of 26 HU compared to no correction, and a 52% increase in contrast-tonoise ratio in simulated hemorrhage compared to “oracle” constant fraction correction. Conclusion: Acceleration of MC-GPU achieved through GPU-optimized variance reduction and kernel smoothing yields an efficient (<5 min/scan) and accurate scatter correction that does not rely on additional hardware or simplifying assumptions about the scatter distribution. The method is undergoing implementation in a novel CBCT dedicated to brain trauma imaging at the point of care in sports and military applications. Research grant from Carestream Health. JY is an employee of Carestream Health.

60 APPLIED LIFE SCIENCES↗