Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Radiological Data”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

Noise spatial nonuniformity and the impact of statistical image reconstruction in CT myocardial perfusion imaging

Purpose: To achieve high temporal resolution in CT myocardial perfusion imaging (MPI), images are often reconstructed using filtered backprojection (FBP) algorithms from data acquired within a short-scan angular range. However, the variation in the central angle from one time frame to the next in gated short scans has been shown to create detrimental partial scan artifacts when performing quantitative MPI measurements. This study has two main purposes. (1) To demonstrate the existence of a distinct detrimental effect in short-scan FBP, i.e., the introduction of a nonuniform spatial image noise distribution; this nonuniformity can lead to unexpectedly high image noise and streaking artifacts, which may affect CT MPI quantification. (2) To demonstrate that statistical image reconstruction (SIR) algorithms can be a potential solution to address the nonuniform spatial noise distribution problem and can also lead to radiation dose reduction in the context of CT MPI. Methods: Projection datasets from a numerically simulated perfusion phantom and an in vivo animal myocardial perfusion CT scan were used in this study. In the numerical phantom, multiple realizations of Poisson noise were added to projection data at each time frame to investigate the spatial distribution of noise. Images from all datasets were reconstructed using both FBP and SIR reconstruction algorithms. To quantify the spatial distribution of noise, the mean and standard deviation were measured in several regions of interest (ROIs) and analyzed across time frames. In the in vivo study, two low-dose scans at tube currents of 25 and 50 mA were reconstructed using FBP and SIR. Quantitative perfusion metrics, namely, the normalized upslope (NUS), myocardial blood volume (MBV), and first moment transit time (FMT), were measured for two ROIs and compared to reference values obtained from a high-dose scan performed at 500 mA. Results: Images reconstructed using FBP showed a highly nonuniform spatial distribution of noise. This spatial nonuniformity led to large fluctuations in the temporal direction. In the numerical phantom study, the level of noise was shown to vary by as much as 87% within a given image, and as much as 110% between different time frames for a ROI far from isocenter. The spatially nonuniform noise pattern was shown to correlate with the source trajectory and the object structure. In contrast, images reconstructed using SIR showed a highly uniform spatial distribution of noise, leading to smaller unexpected noise fluctuations in the temporal direction when a short scan angular range was used. In the numerical phantom study, the noise varied by less than 37% within a given image, and by less than 20% between different time frames. Also, the noise standard deviation in SIR images was on average half of that of FBP images. In the in vivo studies, the deviation observed between quantitative perfusion metrics measured from low-dose scans and high-dose scans was mitigated when SIR was used instead of FBP to reconstruct images. Conclusions: (1) Images reconstructed using FBP suffered from nonuniform spatial noise levels. This nonuniformity is another manifestation of the detrimental effects caused by short-scan reconstruction in CT MPI. (2) Images reconstructed using SIR had a much lower and more uniform noise level and thus can be used as a potential solution to address the FBP nonuniformity. (3) Given the improvement in the accuracy of the perfusion metrics when using SIR, it may be desirable to use a statistical reconstruction framework to perform low-dose dynamic CT MPI.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

WE-G-18A-06: Sinogram Restoration in Helical Cone-Beam CT

Purpose: To extend CT sinogram restoration, which has been shown in 2D to reduce noise and to correct for geometric effects and other degradations at a low computational cost, from 2D to a 3D helical cone-beam geometry. Methods: A method for calculating sinogram degradation coefficients for a helical cone-beam geometry was proposed. These values were used to perform penalized-likelihood sinogram restoration on simulated data that were generated from the FORBILD thorax phantom. Sinogram restorations were performed using both a quadratic penalty and the edge-preserving Huber penalty. After sinogram restoration, Fourier-based analytical methods were used to obtain reconstructions. Resolution-variance trade-offs were investigated for several locations within the reconstructions for the purpose of comparing sinogram restoration to no restoration. In order to compare potential differences, reconstructions were performed using different groups of neighbors in the penalty, two analytical reconstruction methods (Katsevich and single-slice rebinning), and differing helical pitches. Results: The resolution-variance properties of reconstructions restored using sinogram restoration with a Huber penalty outperformed those of reconstructions with no restoration. However, the use of a quadratic sinogram restoration penalty did not lead to an improvement over performing no restoration at the outer regions of the phantom. Application of the Huber penalty to neighbors both within a view and across views did not perform as well as only applying the penalty to neighbors within a view. General improvements in resolution-variance properties using sinogram restoration with the Huber penalty were not dependent on the reconstruction method used or the magnitude of the helical pitch. Conclusion: Sinogram restoration for noise and degradation effects for helical cone-beam CT is feasible and should be able to be applied to clinical data. When applied with the edge-preserving Huber penalty, sinogram restoration leads to an improvement in resolution-variance tradeoffs.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Monitoring internal organ motion with continuous wave radar in CT

Purpose: To avoid motion artifacts in medical imaging or to minimize the exposure of healthy tissues in radiation therapy, medical devices are often synchronized with the patient's respiratory motion. Today's respiratory motion monitors require additional effort to prepare the patients, e.g., mounting a motion belt or placing an optical reflector on the patient's breast. Furthermore, they are not able to measure internal organ motion without implanting markers. An interesting alternative to assess the patient's organ motion is continuous wave radar. The aim of this work is to design, implement, and evaluate such a radar system focusing on application in CT.Methods: The authors designed a radar system operating in the 860 MHz band to monitor the patient motion. In the intended application of the radar system, the antennas are located close to the patient's body inside the table of a CT system. One receive and four transmitting antennas are used to avoid the requirement of exact patient positioning. The radar waves propagate into the patient's body and are reflected at tissue boundaries, for example at the borderline between muscle and adipose tissue, or at the boundaries of organs. At present, the authors focus on the detection of respiratory motion. The radar system consists of the hardware mentioned above as well as of dedicated signal processing software to extract the desired information from the radar signal. The system was evaluated using simulations and measurements. To simulate the radar system, a simulation model based on radar and wave field equations was designed and 4D respiratory-gated CT data sets were used as input. The simulated radar signals and the measured data were processed in the same way. The radar system hardware and the signal processing algorithms were tested with data from ten volunteers. As a reference, the respiratory motion signal was recorded using a breast belt simultaneously with the radar measurements.Results: Concerning the measurements of the test persons, there is a very good correlation (ρ= 0.917) between the respiratory motion phases received by the radar system and the external motion monitor. Our concept of using an array of transmitting antennas turned out to be widely insensitive to the positioning of the test persons. A time shift between the respiratory motion curves recorded with the radar system and the motion curves from the external respiratory monitor was observed which indicates a slight difference between internal organ motion and motion detected by the external respiratory monitor. The simulations were in good accordance with the measurements.Conclusions: A continuous wave radar operating in the near field of the antennas can be used to determine the respiratory motion of humans accurately. In contrast to trigger systems used today, the radar system is able to measure motion inside the body. If such a monitor was routinely available in clinical CT, it would be possible optimizing the scan start with respect to the respiratory state of the patient. Breathing commands would potentially widely be avoided, and as far as uncooperative patients or children are concerned, less sedation might be necessary. Further applications of the radar system could be in radiation therapy or interventional imaging for instance.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Improved UTE-based attenuation correction for cranial PET-MR using dynamic magnetic field monitoring

Purpose: Ultrashort echo time (UTE) MRI has been proposed as a way to produce segmented attenuation maps for PET, as it provides contrast between bone, air, and soft tissue. However, UTE sequences require samples to be acquired during rapidly changing gradient fields, which makes the resulting images prone to eddy current artifacts. In this work it is demonstrated that this can lead to misclassification of tissues in segmented attenuation maps (AC maps) and that these effects can be corrected for by measuring the true k-space trajectories using a magnetic field camera. Methods: The k-space trajectories during a dual echo UTE sequence were measured using a dynamic magnetic field camera. UTE images were reconstructed using nominal trajectories and again using the measured trajectories. A numerical phantom was used to demonstrate the effect of reconstructing with incorrect trajectories. Images of an ovine leg phantom were reconstructed and segmented and the resulting attenuation maps were compared to a segmented map derived from a CT scan of the same phantom, using the Dice similarity measure. The feasibility of the proposed method was demonstrated inin vivo cranial imaging in five healthy volunteers. Simulated PET data were generated for one volunteer to show the impact of misclassifications on the PET reconstruction. Results: Images of the numerical phantom exhibited blurring and edge artifacts on the bone–tissue and air–tissue interfaces when nominal k-space trajectories were used, leading to misclassification of soft tissue as bone and misclassification of bone as air. Images of the tissue phantom and thein vivo cranial images exhibited the same artifacts. The artifacts were greatly reduced when the measured trajectories were used. For the tissue phantom, the Dice coefficient for bone in MR relative to CT was 0.616 using the nominal trajectories and 0.814 using the measured trajectories. The Dice coefficients for soft tissue were 0.933 and 0.934 for the nominal and measured cases, respectively. For air the corresponding figures were 0.991 and 0.993. Compared to an unattenuated reference image, the mean error in simulated PET uptake in the brain was 9.16% when AC maps derived from nominal trajectories was used, with errors in the SUV{sub max} for simulated lesions in the range of 7.17%–12.19%. Corresponding figures when AC maps derived from measured trajectories were used were 0.34% (mean error) and −0.21% to +1.81% (lesions). Conclusions: Eddy current artifacts in UTE imaging can be corrected for by measuring the true k-space trajectories during a calibration scan and using them in subsequent image reconstructions. This improves the accuracy of segmented PET attenuation maps derived from UTE sequences and subsequent PET reconstruction.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

3-Dimensional Magnetic Resonance Spectroscopic Imaging at 3 Tesla for Early Response Assessment of Glioblastoma Patients During External Beam Radiation Therapy

Purpose: To evaluate the utility of 3-dimensional magnetic resonance (3D-MR) proton spectroscopic imaging for treatment planning and its implications for early response assessment in glioblastoma multiforme. Methods and Materials: Eighteen patients with newly diagnosed, histologically confirmed glioblastoma had 3D-MR proton spectroscopic imaging (MRSI) along with T2 and T1 gadolinium-enhanced MR images at simulation and at boost treatment planning after 17 to 20 fractions of radiation therapy. All patients received standard radiation therapy (RT) with concurrent temozolomide followed by adjuvant temozolomide. Imaging for response assessment consisted of MR scans every 2 months. Progression-free survival was defined by the criteria of MacDonald et al. MRSI images obtained at initial simulation were analyzed for choline/N-acetylaspartate ratios (Cho/NAA) on a voxel-by-voxel basis with abnormal activity defined as Cho/NAA ≥2. These images were compared on anatomically matched MRSI data collected after 3 weeks of RT. Changes in Cho/NAA between pretherapy and third-week RT scans were tested using Wilcoxon matched-pairs signed rank tests and correlated with progression-free survival, radiation dose and location of recurrence using Cox proportional hazards regression. Results: After a median follow-up time of 8.6 months, 50% of patients had experienced progression based on imaging. Patients with a decreased or stable mean or median Cho/NAA values had less risk of progression (P<.01). Patients with an increase in mean or median Cho/NAA values at the third-week RT scan had a significantly greater chance of early progression (P<.01). An increased Cho/NAA at the third-week MRSI scan carried a hazard ratio of 2.72 (95% confidence interval, 1.10-6.71; P=.03). Most patients received the prescription dose of RT to the Cho/NAA ≥2 volume, where recurrence most often occurred. Conclusion: Change in mean and median Cho/NAA detected at 3 weeks was a significant predictor of early progression. The potential impact for risk-adaptive therapy based on early spectroscopic findings is suggested.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Relationship of Imaging Frequency and Planning Margin to Account for Intrafraction Prostate Motion: Analysis Based on Real-Time Monitoring Data

Purpose: Correction for intrafraction prostate motion becomes important for hypofraction treatment of prostate cancer. The purpose of this study was to estimate an ideal planning margin to account for intrafraction prostate motion as a function of imaging and repositioning frequency in the absence of continuous prostate motion monitoring. Methods and Materials: For 31 patients receiving intensity modulated radiation therapy treatment, prostate positions sampled at 10 Hz during treatment using the Calypso system were analyzed. Using these data, we simulated multiple, less frequent imaging protocols, including intervals of every 10, 15, 20, 30, 45, 60, 90, 120, 180, and 240 seconds. For each imaging protocol, the prostate displacement at the imaging time was corrected by subtracting prostate shifts from the subsequent displacements in that fraction. Furthermore, we conducted a principal component analysis to quantify the direction of prostate motion. Results: Averaging histograms of every 240 and 60 seconds for all patients, vector displacements of the prostate were, respectively, within 3 and 2 mm for 95% of the treatment time. A vector margin of 1 mm achieved 91.2% coverage of the prostate with 30 second imaging. The principal component analysis for all fractions showed the largest variance in prostate position in the midsagittal plane at 54° from the anterior direction, indicating that anterosuperior to inferoposterior is the direction of greatest motion. The smallest prostate motion is in the left-right direction. Conclusions: The magnitudes of intrafraction prostate motion along the superior-inferior and anterior-posterior directions are comparable, and the smallest motion is in the left-right direction. In the absence of continuous prostate motion monitoring, and under ideal circumstances, 1-, 2-, and 3-mm vector planning margins require a respective imaging frequency of every 15, 60, and 240 to account for intrafraction prostate motion while achieving adequate geometric target coverage for 95% of the time.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRS2001

Shielded Heavy-Ion Environment Linear Detector (SHIELD): an experiment for the Radiation and Technology Demonstration (RTD) Mission

Radiological assessment of the many cosmic ion species of widely distributed energies requires the use of theoretical transport models to accurately describe diverse physical processes related to nuclear reactions in spacecraft structures, planetary atmospheres and surfaces, and tissues. Heavy-ion transport models that were designed to characterize shielded radiation fields have been validated through comparison with data from thick-target irradiation experiments at particle accelerators. With the RTD Mission comes a unique opportunity to validate existing radiation transport models and guide the development of tools for shield design. For the first time, transport properties will be measured in free-space to characterize the shielding effectiveness of materials that are likely to be aboard interplanetary space missions. Target materials composed of aluminum, advanced composite spacecraft structure and other shielding materials, helium (a propellant) and tissue equivalent matrices will be evaluated. Large solid state detectors will provide kinetic energy and charge identification for incident heavy-ions and for secondary ions created in the target material. Transport calculations using the HZETRN model suggest that 8 g cm -2 thick targets would be adequate to evaluate the shielding effectiveness during solar minimum activity conditions for a period of 30 days or more.

NASA Center LaRC↗
DOE OSTI2013

Eliminating Inconsistencies in Simulation and Treatment Planning Orders in Radiation Therapy

Purpose: To identify deficiencies with simulation and treatment planning orders and to develop corrective measures to improve safety and quality. Methods and Materials: At Washington University, the DMAIIC formalism is used for process management, whereby the process is understood as comprising Define, Measure, Analyze, Improve, Implement, and Control activities. Two complementary tools were used to provide quantitative assessments: failure modes and effects analysis and reported event data. The events were classified by the user according to severity. The event rates (ie, number of events divided by the number of opportunities to generate an event) related to simulation and treatment plan orders were determined. Results: We analyzed event data from the period 2008-2009 to design an intelligent SIMulation and treatment PLanning Electronic (SIMPLE) order system. Before implementation of SIMPLE, event rates of 0.16 (420 of 2558) for a group of physicians that were subsequently used as a pilot group and 0.13 (787 of 6023) for all physicians were obtained. An interdisciplinary group evaluated and decided to replace the Microsoft Word-based form with a Web-based order system. This order system has mandatory fields and context-sensitive logic, an ability to create templates, and enables an automated process for communication of orders through an enterprise management system. After the implementation of the SIMPLE order, the event rate decreased to 0.09 (96 of 1001) for the pilot group and to 0.06 (145 of 2140) for all physicians (P<.0001). The average time to complete the SIMPLE form was 3 minutes, as compared with 7 minutes for the Word-based form. The number of severe events decreased from 10.7% (45 of 420) and 12.1% (96 of 787) to 6.2% (6 of 96) and 10.3% (15 of 145) for the pilot group and all physicians, respectively. Conclusions: There was a dramatic reduction in the total and the number of potentially severe events through use of the SIMPLE system. In addition, the order process has become more efficient and reliable.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

A stochastic model for tumor geometry evolution during radiation therapy in cervical cancer

Purpose: To develop mathematical models to predict the evolution of tumor geometry in cervical cancer undergoing radiation therapy. Methods: The authors develop two mathematical models to estimate tumor geometry change: a Markov model and an isomorphic shrinkage model. The Markov model describes tumor evolution by investigating the change in state (either tumor or nontumor) of voxels on the tumor surface. It assumes that the evolution follows a Markov process. Transition probabilities are obtained using maximum likelihood estimation and depend on the states of neighboring voxels. The isomorphic shrinkage model describes tumor shrinkage or growth in terms of layers of voxels on the tumor surface, instead of modeling individual voxels. The two proposed models were applied to data from 29 cervical cancer patients treated at Princess Margaret Cancer Centre and then compared to a constant volume approach. Model performance was measured using sensitivity and specificity. Results: The Markov model outperformed both the isomorphic shrinkage and constant volume models in terms of the trade-off between sensitivity (target coverage) and specificity (normal tissue sparing). Generally, the Markov model achieved a few percentage points in improvement in either sensitivity or specificity compared to the other models. The isomorphic shrinkage model was comparable to the Markov approach under certain parameter settings. Convex tumor shapes were easier to predict. Conclusions: By modeling tumor geometry change at the voxel level using a probabilistic model, improvements in target coverage and normal tissue sparing are possible. Our Markov model is flexible and has tunable parameters to adjust model performance to meet a range of criteria. Such a model may support the development of an adaptive paradigm for radiation therapy of cervical cancer.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

Imaging dose in breast radiotherapy: does breast size affect the dose to the organs at risk and the risk of secondary cancer to the contralateral breast?

Correct target positioning is crucial for accurate dose delivery in breast radiotherapy resulting in utilisation of daily imaging. However, the radiation dose from daily imaging is associated with increased probability of secondary induced cancer. The aim of this study was to quantify doses associated with three imaging modalities and investigate the correlation of dose and varying breast size in breast radiotherapy. Planning computed tomography (CT) data sets of 30 breast cancer patients were utilised to simulate the dose received by various organs from a megavoltage computed tomography (MV-CT), megavoltage electronic portal image (MV-EPI) and megavoltage cone-beam computed tomography (MV-CBCT). The mean dose to organs adjacent to the target volume (contralateral breast, lungs, spinal cord and heart) were analysed. Pearson correlation analysis was performed to determine the relationship between imaging dose and primary breast volume and the lifetime attributable risk (LAR) of induced secondary cancer was calculated for the contralateral breast. The highest contralateral breast mean dose was from the MV-CBCT (1.79 Gy), followed by MV-EPI (0.22 Gy) and MV-CT (0.11 Gy). A similar trend was found for all organs at risk (OAR) analysed. The primary breast volume inversely correlated with the contralateral breast dose for all three imaging modalities. As the primary breast volume increases, the likelihood of a patient developing a radiation-induced secondary cancer to the contralateral breast decreases. MV-CBCT showed a stronger relationship between breast size and LAR of developing a radiation-induced contralateral breast cancer in comparison with the MV-CT and MV-EPI. For breast patients, imaging dose to OAR depends on imaging modality and treated breast size. When considering the use of imaging during breast radiotherapy, the patient's breast size and contralateral breast dose should be taken into account.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

Incorporating Cancer Stem Cells in Radiation Therapy Treatment Response Modeling and the Implication in Glioblastoma Multiforme Treatment Resistance

Purpose: To perform a preliminary exploration with a simplistic mathematical cancer stem cell (CSC) interaction model to determine whether the tumor-intrinsic heterogeneity and dynamic equilibrium between CSCs and differentiated cancer cells (DCCs) can better explain radiation therapy treatment response with a dual-compartment linear-quadratic (DLQ) model. Methods and Materials: The radiosensitivity parameters of CSCs and DCCs for cancer cell lines including glioblastoma multiforme (GBM), non–small cell lung cancer, melanoma, osteosarcoma, and prostate, cervical, and breast cancer were determined by performing robust least-square fitting using the DLQ model on published clonogenic survival data. Fitting performance was compared with the single-compartment LQ (SLQ) and universal survival curve models. The fitting results were then used in an ordinary differential equation describing the kinetics of DCCs and CSCs in response to 2- to 14.3-Gy fractionated treatments. The total dose to achieve tumor control and the fraction size that achieved the least normal biological equivalent dose were calculated. Results: Smaller cell survival fitting errors were observed using DLQ, with the exception of melanoma, which had a low α/β = 0.16 in SLQ. Ordinary differential equation simulation indicated lower normal tissue biological equivalent dose to achieve the same tumor control with a hypofractionated approach for 4 cell lines for the DLQ model, in contrast to SLQ, which favored 2 Gy per fraction for all cells except melanoma. The DLQ model indicated greater tumor radioresistance than SLQ, but the radioresistance was overcome by hypofractionation, other than the GBM cells, which responded poorly to all fractionations. Conclusion: The distinct radiosensitivity and dynamics between CSCs and DCCs in radiation therapy response could perhaps be one possible explanation for the heterogeneous intertumor response to hypofractionation and in some cases superior outcome from stereotactic ablative radiation therapy. The DLQ model also predicted the remarkable GBM radioresistance, a result that is highly consistent with clinical observations. The radioresistance putatively stemmed from accelerated DCC regrowth that rapidly restored compartmental equilibrium between CSCs and DCCs.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

SU-E-T-546: Modeling and Validation of the a New Proton Therapy System Using a Monte-Carlo Environment Optimized for Protons

8 particle histories were run. Results: Range measurements of the Monte-Carlo simulations matched the measured data within 1mm. Distal fall-off of the simulated fields matched within <1mm. Lateral penumbra and field size measurements of the standard-sized square and half-beam blocked fields matched within 1mm at all three planes compared. A small difference was seen in the in-air profiles at doses <0%. The suspected cause of the difference was the aperture shape. The measured data utilized a divergent aperture. The Monte-Carlo calculation used a non-divergent aperture. Conclusion: The validation measurements indicate that we were able to accurately model the MEVION s250 Proton therapy system using Monte-Carlo Calculations. This may reduce the commissioning time for future users. Purpose: Monte-Carlo modeling is an important tool for understanding the behavior of therapeutic proton beams in a heterogeneous media such as the patient. To gain confidence that a Monte-Carlo model is accurate in complex geometries and media, it must first be compared with measurement in simple situations. This study documents the validation of our Monte-Carlo Model. Methods: A model of the MEVION s250 Proton therapy system was created in the TOPAS Monte-Carlo environment using machine geometry and field shaping system information provided by the vendor. For each of 24 options, validation of the TOPAS model was performed by comparing the dose scored by TOPAS to the dose measurements obtained during the commissioning of the treatment planning system. The measurements compared consisted of: pristine peak depth-dose profiles, in-air profiles for a standard-sized square field (20cm×20cm or 10cm×10cm depending on the maximum field size for each option) at isocenter and at 20cm upstream and downstream of isocenter, and in-air profiles with a half-beam blocked aperture at isocenter and at 20cm upstream and downstream of isocenter. For all Monte-Carlo simulations,.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2012

Can Drugs Enhance Hypofractionated Radiotherapy? A Novel Method of Modeling Radiosensitization Using In Vitro Data

Purpose: Hypofractionated radiotherapy (hRT) is being explored for a number of malignancies. The potential benefit of giving concurrent chemotherapy with hRT is not known. We sought to predict the effects of combined modality treatments by using mathematical models derived from laboratory data. Methods and Materials: Data from 26 published clonogenic survival assays for cancer cell lines with and without the use of radiosensitizing chemotherapy were collected. The first three data points of the RT arm of each assay were used to derive parameters for the linear quadratic (LQ) model, the multitarget (MT) model, and the generalized linear quadratic (gLQ) model. For each assay and model, the difference between the predicted and observed surviving fractions at the highest tested RT dose was calculated. The gLQ model was fitted to all the data from each RT cell survival assay, and the biologically equivalent doses in 2-Gy fractions (EQD2s) of clinically relevant hRT regimens were calculated. The increase in cell kill conferred by the addition of chemotherapy was used to estimate the EQD2 of hRT along with a radiosensitizing agent. For comparison, this was repeated using conventionally fractionated RT regimens. Results: At a mean RT dose of 8.0 Gy, the average errors for the LQ, MT, and gLQ models were 1.63, 0.83, and 0.56 log units, respectively, favoring the gLQ model (p < 0.05). Radiosensitizing chemotherapy increased the EQD2 of hRT schedules by an average of 28% to 82%, depending on disease site. This increase was similar to the gains predicted for the addition of chemotherapy to conventionally fractionated RT. Conclusions: Based on published in vitro assays, the gLQ equation is superior to the LQ and MT models in predicting cell kill at high doses of RT. Modeling exercises demonstrate that significant increases in biologically equivalent dose may be achieved with the addition of radiosensitizing agents to hRT. Clinical study of this approach is warranted.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

Multileaf collimator tongue-and-groove effect on depth and off-axis doses: A comparison of treatment planning data with measurements and Monte Carlo calculations

To investigate how accurately treatment planning systems (TPSs) account for the tongue-and-groove (TG) effect, Monte Carlo (MC) simulations and radiochromic film (RCF) measurements were performed for comparison with TPS results. Two commercial TPSs computed the TG effect for Varian Millennium 120 multileaf collimator (MLC). The TG effect on off-axis dose profile at 3 depths of solid water was estimated as the maximum depth and the full width at half maximum (FWHM) of the dose dip at an interleaf position. When compared with the off-axis dose of open field, the maximum depth of the dose dip for MC and RCF ranged from 10.1% to 20.6%; the maximum depth of the dose dip gradually decreased by up to 8.7% with increasing depths of 1.5 to 10 cm and also by up to 4.1% with increasing off-axis distances of 0 to 13 cm. However, TPS results showed at most a 2.7% decrease for the same depth range and a negligible variation for the same off-axis distances. The FWHM of the dose dip was approximately 0.19 cm for MC and 0.17 cm for RCF, but 0.30 cm for Eclipse TPS and 0.45 cm for Pinnacle TPS. Accordingly, the integrated value of TG dose dip for TPS was larger than that for MC and RCF and almost invariant along the depths and off-axis distances. We concluded that the TG dependence on depth and off-axis doses shown in the MC and RCF results could not be appropriately modeled by the TPS versions in this study.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2021

MACCS Theory Manual

This report describes the models of the MACCS computer code as presented in MACCS Version 3.10.0. The purpose of MACCS is to simulate the impact of severe accidents at nuclear power plants on the surrounding environment. MACCS has been developed by Sandia National Laboratories for the U.S. Nuclear Regulatory Commission. From a given release of radioactive material into the atmosphere, MACCS estimates the extent and magnitude of radiological contamination, offsite doses, protective actions, socioeconomic impacts and costs, and health effects. Since the weather at the time of an accident is not predictable, MACCS supports various sampling options to run a representative set of simulations to evaluate weather variability. MACCS simulates atmospheric transport with a straight-line Gaussian plume segment model. From the estimated air and ground concentrations, MACCS models dose projections through several dose exposure pathways. These exposures can be offset by protective actions during the emergency response and long-term recovery of the accident. MACCS users directly specify the evacuation and sheltering area, while other protective actions (e.g., relocation, farmland restrictions, decontamination) are based on user-specified dose or concentration limits. While protective actions help reduce dose accumulation, they also cause social and economic impacts. MACCS models the extent of displaced individuals and land contamination, and the cost of offsite property damage, economic disruptions, and various accident expenditures caused by protective actions. Finally, from the dose accumulation, MACCS estimates early and stochastic health effects according to dose-response models. The purpose of consequence analyses is to be able to understand and estimate the impact of nuclear accidents. Consequence analysis is an essential tool to inform determinations of adequate protection of the public, to understand nuclear power hazards, to measure the value of regulations, and to help us appreciate the importance of nuclear safety. As such, MACCS has a variety of regulatory uses including environmental analyses (10 CFR 51.53, 52.47), regulatory cost-benefit analyses, backfit analyses (10 CFR 50.109), consequence analysis studies such as SOARCA (NUREG-1935), Level 3 PRA studies, and risk-informing of emergency planning (10 CFR 50 App. E and 50.47). This report updates the previous MACCS theory manual (NUREG/CR-4691 Vol. 2; Chanin, Sprung, Ritchie, & Jow, 1990) and accompanies the MACCS User's Guide (SAND-2021-1588) that describes the use and input requirements of the graphical user interface of MACCS known as WinMACCS. The MACCS User's Guide is also a reference guide that describes data input file formats, describes various software components in the MACCS code suite, and provides a set of example tutorials for running WinMACCS. Also, soon to be published is a MACCS input parameter guidance report (NUREG/CR-7270) that provides technical bases for commonly used MACCS input values. This page left blank

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2011

Characterization of an x-ray phase contrast imaging system based on the miniature synchrotron MIRRORCLE-6X

Purpose: The implementation of in-line x-ray phase contrast imaging (PCI) for soft-tissue patient imaging is hampered by the lack of a bright and spatially coherent x-ray source that fits into the hospital environment. This article provides a quantitative characterization of the phase-contrast enhancement of a PCI system based on the miniature synchrotron technology MIRRORCLE-6X. Methods: The phase-contrast effect was measured using an edge response of a plexiglass plate as a function of the incident angle of radiation. We have developed a comprehensive x-ray propagation model based on the system's components, properties, and geometry in order to interpret the measurement data. Monte-Carlo simulations are used to estimate the system's spectral properties and resolution. Results: The measured ratio of the detected phase-contrast to the absorption contrast is currently in the range 100% to 200%. Experiments show that with the current implementation of the MIRRORCLE-6X, a target smaller than 30-40 {mu}m does not lead to a larger phase-contrast. The reason for this is that the fraction of x-rays produced by the material (carbon filament and glue) that is used for mounting the target in the electron beam is more than 25% of the total amount of x-rays produced. This increases the apparent source size. The measured phase-contrast is at maximum two times larger than the absorption contrast with the current set-up. Conclusions: Calculations based on our model of the present imaging system predict that the phase-contrast can be up to an order of magnitude larger than the absorption contrast in case the materials used for mounting the target in the electron beam do not (or hardly) produce x-rays. The methods described in this paper provide vital feedback for guiding future modifications to the design of the x-ray target of MIRRORCLE-type system and configuration of the in-line PCI systems in general.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2011

The Association of Inbreeding With Lung Fibrosis Incidence in Beagle Dogs That Inhaled 238PuO2 or 239PuO2.

Studies of health effects in animals after exposure to internally deposited radionuclides were intended to supplement observational studies in humans. Both nuclear workers and Beagle dogs have exhibited plutonium associated lung fibrosis; however, the dogs smaller gene pool may limit the applicability of findings to humans. Data on Beagles that inhaled either plutonium-238 dioxide (238PuO2) or plutonium-239 dioxide (239PuO2) were analyzed. Wright's Coefficient of Inbreeding was used to measure genetic or familial susceptibility and was assessed as an explanatory variable when modeling the association between lung fibrosis incidence and plutonium exposure. Lung fibrosis was diagnosed in approximately 80% of the exposed dogs compared with 23.7% of the control dogs. The maximum degree of inbreeding was 9.4%. Regardless of isotope, the addition of inbreeding significantly improved the model in female dogs but not in males. In female dogs an increased inbreeding coefficient predicted decreased hazard of a lung fibrosis diagnosis. Lung fibrosis was common in these dogs with inbreeding affecting models of lung fibrosis incidence in females but not in males. The apparent protective effect in females predicted by these models of lung fibrosis incidence is likely to be minimal given the small degree of inbreeding in these groups.

62 RADIOLOGY AND NUCLEAR MEDICINE↗