SEARCH · Engineering Papers

Results for “Representatives”

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 1,333 records · Page 74

All sources NASA NTRS DOE OSTI

DOE OSTI2014

A 4H Silicon Carbide Gate Buffer for Integrated Power Systems

A gate buffer fabricated in a 2-mu m 4H silicon carbide (SiC) process is presented. The circuit is composed of an input buffer stage with a push-pull output stage, and is fabricated using enhancement mode N-channel FETs in a process optimized for SiC power switching devices. Simulation and measurement results of the fabricated gate buffer are presented and compared for operation at various voltage supply levels, with a capacitive load of 2 nF. Details of the design including layout specifics, simulation results, and directions for future improvement of this buffer are presented. In addition, plans for its incorporation into an isolated high-side/low-side gate-driver architecture, fully integrated with power switching devices in a SiC process, are briefly discussed. This letter represents the first reported MOSFET-based gate buffer fabricated in 4H SiC.

Ericson, N↗

DOE OSTI2021

Distribution System State Estimation With High Penetration of Demand Response Enabled Loads

Demand-side operations incentivize utility customers to take part in various grid services. A demand response enabled load (DREL) is a flexible grid asset that schedules electricity consumption in response to a time-of-use (TOU) energy price. Consequently, its energy profile differs from that of a conventional load that is insensitive to price. This difference may cause new challenges for distribution system state estimation (DSSE). It is well known that DSSE often needs to use pseudo-measurements based on historic load profiles to increase system observability. However, historic profiles of conventional loads are not representative of DREL behaviors. Additionally, the inaccuracy impacts DSSE results and other DSSE-dependent operations. In this paper, we propose an online pseudo-measurement generation approach for DSSE with DRELs. We formulate an optimization model to represent DRELs self-adjusting actions. Sampling-based stochastic optimization techniques are proposed to account for uncertainties in DRELs. A set of representative DREL behavior data corresponding to the samples are used to characterize DREL pseudo-measurements. Case studies with modified IEEE 123-bus test system verify the validity of the proposed work.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

A Hierarchical Local Electricity Market for a DER-rich Grid Edge

We report with increasing penetration of distributed energy resources (DER) in the distribution system, it is critical to design market structures that enable smooth integration of DERs. A hierarchical local electricity market (LEM) structure is proposed in this paper with a secondary market (SM) at the lower level representing secondary feeders and a primary market (PM) at the upper level, representing primary feeders, in order to effectively use DERs to increase grid efficiency and resilience. The lower level SM enforces budget, power balance, and flexibility constraints and accounts for costs related to consumers, such as their disutility, flexibility limits, and commitment reliability, while the upper level PM enforces grid physics constraints such as power balance and capacity limits, and also minimizes line losses. The hierarchical LEM is extensively evaluated using a modified IEEE-123 bus with high DER penetration, with each primary feeder consisting of at least three secondary feeders. Data from a GridLAB-D model is used to emulate realistic power injections and load profiles over the course of 24 hours. The performance of the LEM is illustrated by delineating the family of power-injection profiles across the primary and secondary feeders as well as corresponding local electricity tariffs that vary across the distribution grid. Through numerical simulations, the hierarchical LEM is shown to improve the efficiency of the market in terms of lowering overall costs, including both the distribution-level locational marginal prices (d-LMP) as well as retail tariffs paid by customers. Together, it represents an overall framework for a Distribution System Operator (DSO) who can provide the oversight for the entire LEM.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

Traveler: Navigating Task Parallel Traces for Performance Analysis

Understanding the behavior of software in execution is a key step in identifying and fixing performance issues. This is especially important in high performance computing contexts where even minor performance tweaks can translate into large savings in terms of computational resource use. To aid performance analysis, developers may collect an execution trace —a chronological log of program activity during execution. As traces represent the full history, developers can discover a wide array of possibly previously unknown performance issues, making them an important artifact for exploratory performance analysis. However, interactive trace visualization is difficult due to issues of data size and complexity of meaning. Traces represent nanosecond-level events across many parallel processes, meaning the collected data is often large and difficult to explore. The rise of asynchronous task parallel programming paradigms complicates the relation between events and their probable cause. Here, to address these challenges, we conduct a continuing design study in collaboration with high performance computing researchers. We develop diverse and hierarchical ways to navigate and represent execution trace data in support of their trace analysis tasks. Through an iterative design process, we developed Traveler , an integrated visualization platform for task parallel traces. Traveler provides multiple linked interfaces to help navigate trace data from multiple contexts. We evaluate the utility of Traveler through feedback from users and a case study, finding that integrating multiple modes of navigation in our design supported performance analysis tasks and led to the discovery of previously unknown behavior in a distributed array library.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2018

Traditional plant functional groups explain variation in economic but not size-related traits across the tundra biome

Aim: Plant functional groups are widely used in community ecology and earth system modelling to describe trait variation within and across plant communities. However, this approach rests on the assumption that functional groups explain a large proportion of trait variation among species. We test whether four commonly used plant functional groups represent variation in six ecologically important plant traits. Location: Tundra biome. Time period: Data collected between 1964 and 2016. Major taxa studied: 295 tundra vascular plant species. Methods: We compiled a database of six plant traits (plant height, leaf area, specific leaf area, leaf dry matter content, leaf nitrogen, seed mass) for tundra species. We examined the variation in species-level trait expression explained by four traditional functional groups (evergreen shrubs, deciduous shrubs, graminoids, forbs), and whether variation explained was dependent upon the traits included in analysis. We further compared the explanatory power and species composition of functional groups to alternative classifications generated using post-hoc clustering of species-level traits. Results: Traditional functional groups explained significant differences in trait expression, particularly among traits associated with resource economics, which were consistent across sites and at the biome scale. However, functional groups explained 19% of overall trait variation and poorly represented differences in traits associated with plant size. Post-hoc classification of species did not correspond well with traditional functional groups, and explained twice as much variation in species-level trait expression. Main conclusion: Traditional functional groups only coarsely represent variation in well-measured traits within tundra plant communities, and better explain resource economic traits than size-related traits. We recommend caution when using functional group approaches to predict tundra ecosystem change, or ecosystem functions relating to plant size, such as albedo or carbon storage. We argue that alternative classifications or direct use of specific plant traits could provide new insight for ecological prediction and modelling.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Continental United States may lose 1.8 petagrams of soil organic carbon under climate change by 2100

Abstract Aims High‐resolution information on soils’ vulnerability to climate‐induced soil organic carbon (SOC) loss can enable environmental scientists, land managers, and policy makers to develop targeted mitigation strategies. This study aims to estimate baseline and decadal changes in continental US surface SOC stocks under future emission scenarios. Location Continental United States. Time period 2014–2100. Methods We used recent SOC field observations ( n = 6,213 sites), environmental factors ( n = 32), and an ensemble machine learning (ML) approach to estimate baseline SOC stocks in surface soils across the continental United States at 100‐m spatial resolution, and decadal changes under the projected climate scenarios of Coupled Model Intercomparison Project Phase Six (CMIP6) earth system models (ESMs). Results Baseline SOC projections from ML approaches captured more than 50% of variability in SOC observations, whereas ESMs represented only 6–16% of observed SOC variability. ML estimates showed a mean total loss of 1.8 Pg C from US surface soils under the high‐emission scenario by 2100, whereas ESMs showed no significant change in SOC stocks with wide variation among ESMs. Both ML and ESM predictions agree on the direction of SOC change (net emissions or sequestration) across 46–51% of continental US land area. These differences are attributable to the high‐resolution site‐specific data used in the ML models compared to the relatively coarse grid represented in CMIP6 ESMs. Main conclusions Our high‐resolution estimates of baseline SOC stocks, identification of key environmental controllers, and projection of SOC changes from US land cover types under future climate scenarios suggest the need for high‐resolution simulations of SOC in ESMs to represent the heterogeneity of SOC. We found that the SOC change is sensitive to key soil related factors (e.g. soil drainage and soil order) that have not been historically considered as input parameters in ESMs, because currently more than 95% variability in the SOC of CMIP6 ESMs is controlled by net primary productivity, temperature, and precipitation. Using additional environmental factors to estimate the baseline SOC stocks and predict the future trajectory of SOC change can provide more accurate results.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

In-Situ Species Concentration Measurements In Ammonia-Mix Flames Using Ftir Spectroscopy

Hydrogen and ammonia represent two carbon-free fuel sources that could be used in place of current fossil energy sources in combustion systems. To develop optimized ammonia combustion systems, validated modeling tools are needed. In the open literature, it has been shown that the complex chemistry associated with fuel-bound nitrogen contained in ammonia differs greatly from natural gas or hydrogen combustion. As a result, several new chemical kinetic mechanisms have been developed. Many of these mechanisms have been validated experimentally, however this has primarily focused on bulk parameters such as laminar flame speed and ignition delay time. Critically, high quality measurements of species concentrations are needed under controlled conditions which are easily represented by simple models. In this paper, direct, in-situ measurements of species concentrations and gas temperature are performed in a laminar flat-flame burner. This arrangement enables comparison with 1D model predictions, better isolating chemical kinetics from the fluid dynamics. Quantitative species concentrations are determined by absorption spectroscopy using an FTIR spectrometer. Fuel compositions representative of cracked ammonia (NH 3 /H 2 ) and ammonia-natural gas (NH 3 /CH 4 ) are considered for rich and lean equivalence ratios. A major focus of the paper is on the selection of spectral features for nitric oxide and ammonia and correcting for large amounts of baseline H 2 O absorption.

36 MATERIALS SCIENCE↗

DOE OSTI2012

Special Section Guest Editorial: Laser Damage

Laser damage of optical materials, first reported in 1964, continues to limit the output energy and power of pulsed and continuous-wave laser systems. In spite of some 48 years of research in this area, interest from the international laser community to laser damage issues remains at a very high level and does not show any sign of decreasing. Moreover, it grows with the development of novel laser systems, for example, ultrafast and short-wavelength lasers that involve new damage effects and specific mechanisms not studied before. This interest is evident from the high level of attendance and presentations at the annual SPIE Laser Damage Symposium (aka, Boulder Damage Symposium) that has been held in Boulder, Colorado, since 1969. This special section of Optical Engineering is the first one devoted to the entire field of laser damage rather than to a specific part. It is prepared in response to growing interest from the international laser-damage community. Some papers in this special section were presented at the Laser Damage Symposium; others were submitted in response to the general call for papers for this special section. The 18 papers compiled into this special section represent many sides of the broad field of laser-damage research. They consider theoretical studies of the fundamental mechanisms of laser damage including laser-driven electron dynamics in solids (O. Brenk and B. Rethfeld; A. Nikiforov, A. Epifanov, and S. Garnov; T. Apostolova et al.), modeling of propagation effects for ultrashort high-intensity laser pulses (J. Gulley), an overview of mechanisms of inclusion-induced damage (M. Koldunov and A. Manenkov), the formation of specific periodic ripples on a metal surface by femtosecond laser pulses (M. Ahsan and M. Lee), and the laser-plasma effects on damage in glass (Y. Li et al). Material characterization is represented by the papers devoted to accurate and reliable measurements of absorption with special emphasis on thin films (C. Mühlig and S. Bublitz; B. Cho, E. Danielewicz, and J. Rudisill; W. Palm et al; and J. Lu et al.). Statistical treatment of measurements of the laser-damage threshold (J. Arenberg) and the relationship to damage mechanisms (F. Wagner et al.) represent the large subfield of laser-damage measurements. Various aspects of multilayer coating and thin-film characterization are considered in papers by B. Cho, J. Rudisill, and E. Danielewicz (spectral shift in multilayer mirrors) and R. Weber et al. (novel approach to damage studies based on third-harmonic generation microscopy). Of special interest for readers is the paper by C. Stolz that summarizes the results of four “thin-film damage competitions” organized as a part of the Laser Damage Symposium. Another paper is devoted to thermal annealing of damage precursors (N. Shen et al.). Finally, the influence of nano-size contamination on initiation of laser damage by ultrashort pulses is considered in paper of V. Komolov et al.

42 ENGINEERING↗

DOE OSTI2013

The need for application-based adaptation of deformable image registration

Purpose: To utilize a deformable phantom to objectively evaluate the accuracy of 11 different deformable image registration (DIR) algorithms. Methods: The phantom represents an axial plane of the pelvic anatomy. Urethane plastic serves as the bony anatomy and urethane rubber with three levels of Hounsfield units (HU) is used to represent fat and organs, including the prostate. A plastic insert is placed into the phantom to simulate bladder filling. Nonradiopaque markers reside on the phantom surface. Optical camera images of these markers are used to measure the positions and determine the deformation from the bladder insert. Eleven different DIR algorithms are applied to the full and empty-bladder computed tomography images of the phantom (fixed and moving volumes, respectively) to calculate the deformation. The algorithms include those from MIM Software (MIM) and Velocity Medical Solutions (VEL) and nine different implementations from the deformable image registration and adaptive radiotherapy toolbox for Matlab. These algorithms warp one image to make it similar to another, but must utilize a method for regularization to avoid physically unrealistic deformation scenarios. The mean absolute difference (MAD) between the HUs at the marker locations on one image and the calculated location on the other serves as a metric to evaluate the balance between image similarity and regularization. To demonstrate the effect of regularization on registration accuracy, an additional beta version of MIM was created with a variable smoothness factor that controls the emphasis of the algorithm on regularization. The distance to agreement between the measured and calculated marker deformations is used to compare the overall spatial accuracy of the DIR algorithms. This overall spatial accuracy is also utilized to evaluate the phantom geometry and the ability of the phantom soft-tissue heterogeneity to represent patient data. To evaluate the ability of the DIR algorithms to accurately transfer anatomical contours, the rectum is delineated on both the fixed and moving images. A Dice similarity coefficient is then calculated between the contour on the fixed image and that transferred, via the calculated deformation, from the moving to the fixed image. Results: The phantom possesses sufficient soft-tissue heterogeneity to act as a proxy for patient data. Large discrepancies appear between the algorithms and the measured ground-truth deformation. VEL yields the smallest mean spatial error and a Dice coefficient of 0.90. MIM produces the lowest MAD value and the highest Dice coefficient of 0.96, but creates the largest spatial errors. Increasing the MIM smoothness factor above the default value improves the overall spatial accuracy, but the factor associated with the lowest mean error decreases the Dice coefficient to 0.85. Conclusions: Different applications of DIR require disparate balances between image similarity and regularization. A DIR algorithm that is optimized only for its ability to transfer anatomical contours will yield large deformation errors in homogeneous regions, which is problematic for dose mapping. For this reason, these algorithms must be tested for their overall spatial accuracy. The developed phantom is an objective tool for this purpose.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2013

The effect of irregular breathing patterns on internal target volumes in four-dimensional CT and cone-beam CT images in the context of stereotactic lung radiotherapy

Purpose: Stereotactic lung radiotherapy is complicated by tumor motion from patient respiration. Four-dimensional CT (4DCT) imaging is a motion compensation method used in treatment planning to generate a maximum intensity projection (MIP) internal target volume (ITV). Image guided radiotherapy during treatment may involve acquiring a volumetric cone-beam CT (CBCT) image and visually aligning the tumor to the planning 4DCT MIP ITV contour. Moving targets imaged with CBCT can appear blurred and currently there are no studies reporting on the effect that irregular breathing patterns have on CBCT volumes and their alignment to 4DCT MIP ITV contours. The objective of this work was therefore to image a phantom moving with irregular breathing patterns to determine whether any configurations resulted in errors in volume contouring or alignment. Methods: A Perspex thorax phantom was used to simulate a patient. Three wooden 'lung' inserts with embedded Perspex 'lesions' were moved up to 4 cm with computer-generated motion patterns, and up to 1 cm with patient-specific breathing patterns. The phantom was imaged on 4DCT and CBCT with the same acquisition settings used for stereotactic lung patients in the clinic and the volumes on all phantom images were contoured. This project assessed the volumes for qualitative and quantitative changes including volume, length of the volume, and errors in alignment between CBCT volumes and 4DCT MIP ITV contours. Results: When motion was introduced 4DCT and CBCT volumes were reduced by up to 20% and 30% and shortened by up to 7 and 11 mm, respectively, indicating that volume was being under-represented at the extremes of motion. Banding artifacts were present in 4DCT MIP images, while CBCT volumes were largely reduced in contrast. When variable amplitudes from patient traces were used and CBCT ITVs were compared to 4DCT MIP ITVs there was a distinct trend in reduced ITV with increasing amplitude that was not seen when compared to true ITVs. Breathing patterns with a rest period following expiration resulted in well-defined superior edges and were better aligned using an edge-to-edge alignment technique. In most cases, sinusoidal motion patterns resulted in the closest agreements to true values and the smallest misalignments. Conclusions: Strategies are needed to compensate for volume losses at the extremes of motion for both 4DCT MIP and CBCT images for larger and varied amplitudes, and for patterns with rest periods following expiration. Lesions moving greater than 2 cm would warrant larger treatment margins added to the 4DCT MIP ITV to account for the volume being under-represented at the extremes of motion. Lesions moving with a rest period following expiration would be better aligned using an edge-to-edge alignment technique. Sinusoidal patterns represented the ideal clinical scenario, reinforcing the importance of investigating clinically relevant motions and their effects on 4DCT MIP and CBCT volumes. Since most patients do not breathe sinusoidally this may lead to misinterpretation of previous studies using only sinusoidal motion.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2014

SU-E-T-51: Bayesian Network Models for Radiotherapy Error Detection

Purpose: To develop a probabilistic model of radiotherapy plans using Bayesian networks that will detect potential errors in radiation delivery. Methods: Semi-structured interviews with medical physicists and other domain experts were employed to generate a set of layered nodes and arcs forming a Bayesian Network (BN) which encapsulates relevant radiotherapy concepts and their associated interdependencies. Concepts in the final network were limited to those whose parameters are represented in the institutional database at a level significant enough to develop mathematical distributions. The concept-relation knowledge base was constructed using the Web Ontology Language (OWL) and translated into Hugin Expert Bayes Network files via the the RHugin package in the R statistical programming language. A subset of de-identified data derived from a Mosaiq relational database representing 1937 unique prescription cases was processed and pre-screened for errors and then used by the Hugin implementation of the Estimation-Maximization (EM) algorithm for machine learning all parameter distributions. Individual networks were generated for each of several commonly treated anatomic regions identified by ICD-9 neoplasm categories including lung, brain, lymphoma, and female breast. Results: The resulting Bayesian networks represent a large part of the probabilistic knowledge inherent in treatment planning. By populating the networks entirely with data captured from a clinical oncology information management system over the course of several years of normal practice, we were able to create accurate probability tables with no additional time spent by experts or clinicians. These probabilistic descriptions of the treatment planning allow one to check if a treatment plan is within the normal scope of practice, given some initial set of clinical evidence and thereby detect for potential outliers to be flagged for further investigation. Conclusion: The networks developed here support the use of probabilistic models into clinical chart checking for improved detection of potential errors in RT plans.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2014

SU-F-BRE-03: Consideration of a Track-Interaction Model for Radiochromic Film Response

Purpose: Conventional methods for characterizing the energy response of radiochromic film (RCF) typically involve assessing changes in response when exposed to various beam qualities and use Monte Carlo to determine absorbed dose. These methods represent RCF as a dose integrator of a homogeneous energy deposition volume. Apparent film saturation, nonlinearity, and intrinsic energy dependence are unpredicted with conventional methods. Recent work has shown significant RCF intrinsic energy dependence, which limits its use in absolute dosimetry. This work introduces a track-interaction model (TIM) for RCF and assesses its ability to predict total energy response. Methods: A TIM based on Katz single-hit theory was developed to accumulate energy flux along particle tracks within active crystals, represented as (1×1×20)um{sup 3} prisms about the Gafchromic™ EBT3 active volume using MCNP5 and Matlab. Energy flux contributed to film response only if near the threshold energy for polymerization in polydiacetylenes (2.5eV/monomer). Energy deposition in excess of maximum efficiency represented crystal saturation and did not contribute to film response. The TIM was applied to RCF exposed in air to various monoenergetic photon beams and Co-60. Geometric distribution of energy flux was found for each beam quality in a (1×1)mm{sup 2} RCF area. RCF response relative to Co-60 absorbed dose-to-water (S-TIM) was determined and compared to published values (S-PUB). Results: TIM successfully predicted that lower energy radiation is less effective at inducing polymerization, though the magnitude of the phenomenon was overpredicted. S-TIM was −29% and +20% for 20 and 40 keV, respectively. This agreed qualitatively with S-PUB of −27% and +16%. TIM-generated sensitometric curves contained the non-linearity and saturation apparent in RCF. Conclusion: This work indicates the possibility for TIMs to predict changes in RCF response to various energies. Future work will refine TIM by considering size distributions of active elements, other RCF formulations, and response in complex treatment fields. Acknowledgements: The authors wish to thank Dr. Christopher Soares at NIST (Gaithersburg, MD) for performing film measurements and Dr. David Lewis of Ashland, Inc. (Wayne, NJ) for making the EBT3 film available to us.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2014

MO-C-17A-05: A Three-Dimensional Head-And-Neck Phantom for Validation of Kilovoltage- and Megavoltage-Based Deformable Image Registration

Purpose: To develop a three-dimensional (3D) deformable head-and-neck (H and N) phantom with realistic tissue contrast for both kilovoltage and megavoltage computed tomography and use it to objectively evaluate deformable image registration (DIR) algorithms. Methods: The phantom represents H and N patient anatomy. It is constructed from thermoplastic, which becomes pliable in boiling water, and hardened epoxy resin. Using a system of additives, the Hounsfield unit (HU) values of these materials were tuned to mimic anatomy for both kilovoltage (kV) and megavoltage (MV) imaging. The phantom opened along a sagittal midsection to reveal nonradiopaque markers, which were used to characterize the phantom deformation. The deformed and undeformed phantom was scanned with kV and MV computed tomography. Additionally, a calibration curve was created to change the HUs of the MV scans to be similar to kV HUs, (MC). The extracted ground-truth deformation was then compared to the results of two commercially available DIR algorithms, from Velocity Medical Solutions and MIM Software. Results: The phantom produced a 3D deformation, representing neck flexion, with a magnitude of up to 8 mm and was able represent tissue HUs for both kV and MV imaging modalities. The two tested deformation algorithms yielded vastly different results. For kV-kV registration, MIM made the lowest mean error, and Velocity made the lowest maximum error. For MV-MV, kV-MV, and kV-MC Velocity produced both the lowest mean and lowest maximum errors. Conclusion: The application of DIR across different imaging modalities is particularly difficult, due to differences in tissue HUs and the presence of imaging artifacts. For this reason, DIR algorithms must be validated specifically for this purpose. The developed H and N phantom is an effective tool for this purpose.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2014

An initial study on the estimation of time-varying volumetric treatment images and 3D tumor localization from single MV cine EPID images

Purpose: In this work the authors develop and investigate the feasibility of a method to estimate time-varying volumetric images from individual MV cine electronic portal image device (EPID) images. Methods: The authors adopt a two-step approach to time-varying volumetric image estimation from a single cine EPID image. In the first step, a patient-specific motion model is constructed from 4DCT. In the second step, parameters in the motion model are tuned according to the information in the EPID image. The patient-specific motion model is based on a compact representation of lung motion represented in displacement vector fields (DVFs). DVFs are calculated through deformable image registration (DIR) of a reference 4DCT phase image (typically peak-exhale) to a set of 4DCT images corresponding to different phases of a breathing cycle. The salient characteristics in the DVFs are captured in a compact representation through principal component analysis (PCA). PCA decouples the spatial and temporal components of the DVFs. Spatial information is represented in eigenvectors and the temporal information is represented by eigen-coefficients. To generate a new volumetric image, the eigen-coefficients are updated via cost function optimization based on digitally reconstructed radiographs and projection images. The updated eigen-coefficients are then multiplied with the eigenvectors to obtain updated DVFs that, in turn, give the volumetric image corresponding to the cine EPID image. Results: The algorithm was tested on (1) Eight digital eXtended CArdiac-Torso phantom datasets based on different irregular patient breathing patterns and (2) patient cine EPID images acquired during SBRT treatments. The root-mean-squared tumor localization error is (0.73 ± 0.63 mm) for the XCAT data and (0.90 ± 0.65 mm) for the patient data. Conclusions: The authors introduced a novel method of estimating volumetric time-varying images from single cine EPID images and a PCA-based lung motion model. This is the first method to estimate volumetric time-varying images from single MV cine EPID images, and has the potential to provide volumetric information with no additional imaging dose to the patient.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2014

Poster — Thur Eve — 03: Application of the non-negative matrix factorization technique to [{sup 11}C]-DTBZ dynamic PET data for the early detection of Parkinson's disease

In this work, a class of non-negative matrix factorization (NMF) technique known as alternating non-negative least squares, combined with the projected gradient method, is used to analyze twenty-five [{sup 11}C]-DTBZ dynamic PET/CT brain data. For each subject, a two-factor model is assumed and two factors representing the striatum (factor 1) and the non-striatum (factor 2) tissues are extracted using the proposed NMF technique and commercially available factor analysis software “Pixies”. The extracted factor 1 and 2 curves represent the binding site of the radiotracer and describe the uptake and clearance of the radiotracer by soft tissues in the brain, respectively. The proposed NMF technique uses prior information about the dynamic data to obtain sample time-activity curves representing the striatum and the non-striatum tissues. These curves are then used for “warm” starting the optimization. Factor solutions from the two methods are compared graphically and quantitatively. In healthy subjects, radiotracer uptake by factors 1 and 2 are approximately 35–40% and 60–65%, respectively. The solutions are also used to develop a factor-based metric for the detection of early, untreated Parkinson's disease. The metric stratifies healthy subjects from suspected Parkinson's patients (based on the graphical method). The analysis shows that both techniques produce comparable results with similar computational time. The “semi-automatic” approach used by the NMF technique allows clinicians to manually set a starting condition for “warm” starting the optimization in order to facilitate control and efficient interaction with the data.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2014

A three-dimensional head-and-neck phantom for validation of multimodality deformable image registration for adaptive radiotherapy

Purpose: To develop a three-dimensional (3D) deformable head-and-neck (H and N) phantom with realistic tissue contrast for both kilovoltage (kV) and megavoltage (MV) imaging modalities and use it to objectively evaluate deformable image registration (DIR) algorithms. Methods: The phantom represents H and N patient anatomy. It is constructed from thermoplastic, which becomes pliable in boiling water, and hardened epoxy resin. Using a system of additives, the Hounsfield unit (HU) values of these materials were tuned to mimic anatomy for both kV and MV imaging. The phantom opens along a sagittal midsection to reveal radiotransparent markers, which were used to characterize the phantom deformation. The deformed and undeformed phantoms were scanned with kV and MV imaging modalities. Additionally, a calibration curve was created to change the HUs of the MV scans to be similar to kV HUs, (MC). The extracted ground-truth deformation was then compared to the results of two commercially available DIR algorithms, from Velocity Medical Solutions and MIM software. Results: The phantom produced a 3D deformation, representing neck flexion, with a magnitude of up to 8 mm and was able to represent tissue HUs for both kV and MV imaging modalities. The two tested deformation algorithms yielded vastly different results. For kV–kV registration, MIM produced mean and maximum errors of 1.8 and 11.5 mm, respectively. These same numbers for Velocity were 2.4 and 7.1 mm, respectively. For MV–MV, kV–MV, and kV–MC Velocity produced similar mean and maximum error values. MIM, however, produced gross errors for all three of these scenarios, with maximum errors ranging from 33.4 to 41.6 mm. Conclusions: The application of DIR across different imaging modalities is particularly difficult, due to differences in tissue HUs and the presence of imaging artifacts. For this reason, DIR algorithms must be validated specifically for this purpose. The developed H and N phantom is an effective tool for this purpose.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2015

Investigating the accuracy of microstereotactic-body-radiotherapy utilizing anatomically accurate 3D printed rodent-morphic dosimeters

Purpose: Sophisticated small animal irradiators, incorporating cone-beam-CT image-guidance, have recently been developed which enable exploration of the efficacy of advanced radiation treatments in the preclinical setting. Microstereotactic-body-radiation-therapy (microSBRT) is one technique of interest, utilizing field sizes in the range of 1–15 mm. Verification of the accuracy of microSBRT treatment delivery is challenging due to the lack of available methods to comprehensively measure dose distributions in representative phantoms with sufficiently high spatial resolution and in 3 dimensions (3D). This work introduces a potential solution in the form of anatomically accurate rodent-morphic 3D dosimeters compatible with ultrahigh resolution (0.3 mm{sup 3}) optical computed tomography (optical-CT) dose read-out. Methods: Rodent-morphic dosimeters were produced by 3D-printing molds of rodent anatomy directly from contours defined on x-ray CT data sets of rats and mice, and using these molds to create tissue-equivalent radiochromic 3D dosimeters from Presage. Anatomically accurate spines were incorporated into some dosimeters, by first 3D printing the spine mold, then forming a high-Z bone equivalent spine insert. This spine insert was then set inside the tissue equivalent body mold. The high-Z spinal insert enabled representative cone-beam CT IGRT targeting. On irradiation, a linear radiochromic change in optical-density occurs in the dosimeter, which is proportional to absorbed dose, and was read out using optical-CT in high-resolution (0.5 mm isotropic voxels). Optical-CT data were converted to absolute dose in two ways: (i) using a calibration curve derived from other Presage dosimeters from the same batch, and (ii) by independent measurement of calibrated dose at a point using a novel detector comprised of a yttrium oxide based nanocrystalline scintillator, with a submillimeter active length. A microSBRT spinal treatment was delivered consisting of a 180° continuous arc at 225 kVp with a 20 × 10 mm field size. Dose response was evaluated using both the Presage/optical-CT 3D dosimetry system described above, and independent verification in select planes using EBT2 radiochromic film placed inside rodent-morphic dosimeters that had been sectioned in half. Results: Rodent-morphic 3D dosimeters were successfully produced from Presage radiochromic material by utilizing 3D printed molds of rat CT contours. The dosimeters were found to be compatible with optical-CT dose readout in high-resolution 3D (0.5 mm isotropic voxels) with minimal artifacts or noise. Cone-beam CT image guidance was possible with these dosimeters due to sufficient contrast between high-Z spinal inserts and tissue equivalent Presage material (CNR ∼10 on CBCT images). Dose at isocenter measured with optical-CT was found to agree with nanoscintillator measurement to within 2.8%. Maximum dose in line profiles taken through Presage and film dose slices agreed within 3%, with FWHM measurements through each profile found to agree within 2%. Conclusions: This work demonstrates the feasibility of using 3D printing technology to make anatomically accurate Presage rodent-morphic dosimeters incorporating spinal-mimicking inserts. High quality optical-CT 3D dosimetry is feasible on these dosimeters, despite the irregular surfaces and implanted inserts. The ability to measure dose distributions in anatomically accurate phantoms represents a powerful useful additional verification tool for preclinical microSBRT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2015

Development of realistic physical breast phantoms matched to virtual breast phantoms based on human subject data

Purpose: Physical phantoms are essential for the development, optimization, and evaluation of x-ray breast imaging systems. Recognizing the major effect of anatomy on image quality and clinical performance, such phantoms should ideally reflect the three-dimensional structure of the human breast. Currently, there is no commercially available three-dimensional physical breast phantom that is anthropomorphic. The authors present the development of a new suite of physical breast phantoms based on human data. Methods: The phantoms were designed to match the extended cardiac-torso virtual breast phantoms that were based on dedicated breast computed tomography images of human subjects. The phantoms were fabricated by high-resolution multimaterial additive manufacturing (3D printing) technology. The glandular equivalency of the photopolymer materials was measured relative to breast tissue-equivalent plastic materials. Based on the current state-of-the-art in the technology and available materials, two variations were fabricated. The first was a dual-material phantom, the Doublet. Fibroglandular tissue and skin were represented by the most radiographically dense material available; adipose tissue was represented by the least radiographically dense material. The second variation, the Singlet, was fabricated with a single material to represent fibroglandular tissue and skin. It was subsequently filled with adipose-equivalent materials including oil, beeswax, and permanent urethane-based polymer. Simulated microcalcification clusters were further included in the phantoms via crushed eggshells. The phantoms were imaged and characterized visually and quantitatively. Results: The mammographic projections and tomosynthesis reconstructed images of the fabricated phantoms yielded realistic breast background. The mammograms of the phantoms demonstrated close correlation with simulated mammographic projection images of the corresponding virtual phantoms. Furthermore, power-law descriptions of the phantom images were in general agreement with real human images. The Singlet approach offered more realistic contrast as compared to the Doublet approach, but at the expense of air bubbles and air pockets that formed during the filling process. Conclusions: The presented physical breast phantoms and their matching virtual breast phantoms offer realistic breast anatomy, patient variability, and ease of use, making them a potential candidate for performing both system quality control testing and virtual clinical trials.

60 APPLIED LIFE SCIENCES↗