Engineering topics

Ferreira, C

Publications and source records attributed to Ferreira, C.

10 records

DOE OSTI2016

SU-F-I-70: Investigation of Gafchromic EBT3 Film Energy Dependence Using Proton, Photon, and Electron Beams

Purpose: To investigate the energy dependence of Gafchromic EBT3 film over a range of clinically used proton, photon and electron energies. Methods: Proton beam energies of 117 and 204 MeV, corresponding respectively to ranges in water of 10 cm and 27 cm from a Mevion S250 double scatter system unit were used. Electron energies of 6 and 20 MeV and photon energies of 6 and 18 MV from a Varian Clinac 21EX Linac were used. Two pieces of film (5×5 cm{sup 2}) were irradiated sequentially for doses of 100, 500, and 1000 cGy for all energies and modalities. Films were placed on the central beam axis for a 10×10 cm{sup 2} field size in the middle of spread out Bragg peak (SOBP) for proton and in respective dmax for photon and electron energies. Films were scanned on a flatbed Epson Expression 10000 XL scanner on the central region of the scanning window using 48-bit, 300 dpi, and landscape orientation after 48 hours post-irradiation of film to account for optical density (OD) stabilization. Film analysis of the red channel was performed using ImageJ 1.48v (National Institutes of Health). Results: The energy dependency of EBT3 among all energies and modalities for all doses studied was small within measurement uncertainties (1σ = ± 4.1%). The mean net OD in red channel for films receiving the same dose in the same energy modality had standard deviations within 0.9% for photons, 4.9% for electrons and 1.8% for protons. It was observed that film pieces were activated during proton irradiation, e.g., 7 mR/hr at surface after 30 minutes of irradiation, lasting for 2 hours post irradiation. Conclusion: EBT3 energy dependency was evaluated for clinically used proton, photon, and electron energies. The film self-activation may have contributed to fog and negligible dose.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-F-T-48: Clinical Implementation of Brachytherapy Planning System for COMS Eye Plaques

Purpose: To commission the Brachytherapy Planning (BP) system (Varian, Palo Alto, CA) for the Collaborative Ocular Melanoma Study (COMS) eye plaques by evaluating dose differences against original plans from Nucletron Planning System (NPS). Methods: NPS system is the primary planning software for COMS-plaques at our facility; however, Brachytherapy Planning 11.0.47 (Varian Medical Systems) is used for secondary check and for seed placement configurations not originally commissioned. Dose comparisons of BP and NPS plans were performed for prescription of 8500 cGy at 5 mm depth and doses to normal structures: opposite retina, inner sclera, macula, optic disk and lens. Plans were calculated for Iodine-125 seeds (OncoSeeds, Model 6711) using COMS-plaques of 10, 12, 14, 16, 18 and 20 mm diameters. An in-house program based on inverse-square was utilized to calculate point doses for comparison as well. Results: The highest dose difference between BP and NPS was 3.7% for the prescription point for all plaques. Doses for BP were higher than doses reported by NPS for all points. The largest percent differences for apex, opposite retina, inner sclera, macula, optic disk, and lens were 3.2%, 0.9%, 13.5%, 20.5%, 15.7% and 2.2%, respectively. The dose calculated by the in-house program was 1.3% higher at the prescription point, and were as high as 42.1%, for points away from the plaque (i.e. opposite retina) when compared to NPS. Conclusion: Doses to the tumor, lens, retina, and optic nerve are paramount for a successful treatment and vision preservation. Both systems are based on TG-43 calculations and assume water medium tissue homogeneity (ρe=1, water medium). Variations seen may result from the different task group versions and/or mathematical algorithms of the software. BP was commissioned to serve as a backup system and it also enables dose calculation in cases where seeds don’t follow conventional placement configuration.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-F-T-143: Implementation of a Correction-Based Output Model for a Compact Passively Scattered Proton Therapy System

Purpose: To commission and investigate the accuracy of an output (cGy/MU) prediction model for a compact passively scattered proton therapy system. Methods: A previously published output prediction model (Sahoo et al, Med Phys, 35, 5088–5097, 2008) was commissioned for our Mevion S250 proton therapy system. This model is a correction-based model that multiplies correction factors (d/MUwnc=ROFxSOBPF xRSFxSOBPOCFxOCRxFSFxISF). These factors accounted for changes in output due to options (12 large, 5 deep, and 7 small), modulation width M, range R, off-center, off-axis, field-size, and off-isocenter. In this study, the model was modified to ROFxSOBPFxRSFxOCRxFSFxISF-OCFxGACF by merging SOBPOCF and ISF for simplicity and introducing a gantry angle correction factor (GACF). To commission the model, outputs over 1,000 data points were taken at the time of the system commissioning. The output was predicted by interpolation (1D for SOBPF, FSF, and GACF; 2D for RSF and OCR) with inverse-square calculation (ISF-OCR). The outputs of 273 combinations of R and M covering total 24 options were measured to test the model. To minimize fluence perturbation, scattered dose from range compensator and patient was not considered. The percent differences between the predicted (P) and measured (M) outputs were calculated to test the prediction accuracy ([P-M]/Mx100%). Results: GACF was required because of up to 3.5% output variation dependence on the gantry angle. A 2D interpolation was required for OCR because the dose distribution was not radially symmetric especially for the deep options. The average percent differences were −0.03±0.98% (mean±SD) and the differences of all the measurements fell within ±3%. Conclusion: It is concluded that the model can be clinically used for the compact passively scattered proton therapy system. However, great care should be taken when the field-size is less than 5×5 cm{sup 2} where a direct output measurement is required due to substantial output change by irregular block shape.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-F-T-160: Commissioning of a Single-Room Double-Scattering Proton Therapy System

Purpose: To report the detailed commissioning experience for a compact double-scattering Mevion S250 proton therapy system at a University Cancer Center site. Methods: The commissioning of the proton therapy system mainly consisted of ensuring integrity of mechanical and imaging system, beam data collection, and commissioning of a treatment planning system (TPS). First, mechanical alignment and imaging were tested including safety, interlocks, positional accuracy of couch and gantry, image quality, mechanical and imaging isocenter and so on. Second, extensive beam data (outputs, PDDs, and profiles) were collected and analyzed through effective sampling of range (R) and modulation width (M) from 24 beam options. Three different output (cGy/MU) prediction models were also commissioned as primary and secondary MU calculation tool. Third, the Varian Eclipse TPS was commissioned through five sets of data collections (in-water Bragg peak scans, in-air longitudinal fluence scans, in-air lateral profiles, in-air half-beam profiles, and an HU-to-stopping-power conversion curve) and accuracy of TPS calculation was tested using in-water scans and dose measurements with a 2D array detector with block and range compensator. Finally, an anthropomorphic phantom was scanned and heterogeneity effects were tested by inserting radiochromic films in the phantom and PET activation scans for range verification in conjunction with end-to-end test. Results: Beam characteristics agreed well with the vendor specifications; however, minor mismatches in R and M were found in some measurements during the beam data collection. These were reflected into the TPS commissioning such that the TPS could accurately predict the R and M within tolerance levels. The output models had a good agreement with measured outputs (<3% error). The end-to-end test using the film and PET showed reasonably the TPS predicted dose, R and M in heterogeneous medium. Conclusion: The proton therapy system was successfully commissioned and was released for clinical use.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-F-T-219: Verification of the Accuracy of a Relative Stopping Power (RSP) to Hounsfield Unit (HU) Calibration Curve Produced for Proton Treatment Planning Using EBT3-Gafchromic Film Dosimetry

Purpose: Accuracy of a RSP-HU calibration curve produced for proton treatment planning is tested by comparing the treatment planning system dose grid to physical doses delivered on film by a Mevion S250 double-scattering proton unit. Methods: A single batch of EBT3 Gafchromic film was used for calibration and measurements. The film calibration curve was obtained using Mevion proton beam reference option 20 (15cm range, 10cm modulation). Paired films were positioned at the center of the spread out Bragg peak (SOBP) in solid water. The calibration doses were verified with an ion chamber, including background and doses from 20cGy to 350cGy. Films were scanned in a flatbed Epson-Expression 10000-XL scanner, and analyzed using the red channel. A Rando phantom was scanned with a GE LightSpeed CT Simulator. A single-field proton plan (Eclipse, Varian) was calculated to deliver 171cGy to the pelvis section (heterogeneous region), using a standard 4×4cm aperture without compensator, 7.89cm beam range, and 5.36cm SOBP. Varied depths of the calculated distal 90% isodose-line were recorded and compared. The dose distribution from film irradiated between Rando slices was compared with the calculated plans using RIT v.6.2. Results: Distal 90% isodose-line depth variation between CT scans was 2mm on average, and 4mm at maximum. Fine calculation of this variation was restricted by the dose calculation grid, as well as the slice thickness. Dose differences between calibrated film measurements and calculated doses were on average 5.93cGy (3.5%), with the large majority of differences forming a normal distribution around 3.5cGy (2%). Calculated doses were almost entirely greater than those measured. Conclusion: RSP to HU calibration curve is shown to produce distal depth variation within the margin of tolerance (±4.3mm) across all potential scan energies and protocols. Dose distribution calculation is accurate to 2–4% within the SOBP, including areas of high tissue heterogeneity.

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2016

SU-G-201-12: Investigation of Beta-Emitter 90Sr-90Y Dose Distribution Using Gafchromic EBT3 Film for Application On Conformal Skin Brachytherapy Device

Purpose: To investigate {sup 90}Sr-{sup 90}Y as a high dose rate (HDR) source for application in a conformal skin brachytherapy (CSBT) device. The CSBT device has been previously developed to provide patient specific treatment for small inoperable lesions and irregular surfaces. Methods: A popular beta emitter, {sup 90}Sr-{sup 90}Y was tested for feasibility in a CSBT device. A 1 cm diameter plaque was used to deliver dose to a solid water phantom containing EBT3 Gafchromic films arranged at the surface and perpendicular to it. Additionally, a 1 cm diameter 6 MeV electron beam was used to irradiate EBT3 film at 100 cm SSD with a 0.5 cm bolus. Films were digitized with an Epson Expression 10000 XL scanner and calibrated with a 6 MeV electron specific dose curve. Normalized percent depth doses (PDD) and dose profiles for both techniques were analyzed using ImageJ. Results: Dose distributions achieved with the {sup 90}Sr-{sup 90}Y sources were compared with those of external electron beam radiation therapy (EBRT). Penumbra (20%- 80%) for EBRT and 90Sr-90Y were 4.3 mm and 1.6 mm, respectively. PDD values of 50% (normalized to 2 mm) were 10.1 mm and 2.8 mm for electron and {sup 90}Sr-{sup 90}Y, respectively. Flatness (80% of the central beam profile) was 14.1% at a 5 mm depth for EBRT and 4.0% at surface for the {sup 90}Sr-{sup 90}Y. Conclusion: As expected, the PDDs of {sup 90}Sr-{sup 90}Y in water are shallower than that of external electron beams for the same field size. {sup 90}Sr-{sup 90}Y can be used in CSBT to provide patient specific treatment where shallower depth doses than that provided by electron external beams may be required: e.g. eyelids, nose, lips, ears, etc. The customizability of EBRT could be replicated by using multiple adjacent {sup 90}Sr-{sup 90}Y plaque placements.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2015

SU-E-T-424: Feasibility of 3D Printed Radiological Equivalent Customizable Tissue Like Materials

Purpose: To investigate the feasibility of 3D printing CT# specific radiological equivalent tissue like materials. Methods: A desktop 3D printer was utilized to create a series of 3 cm x 3 cm x 2 cm PLA plastic blocks of varying fill densities. The fill pattern was selected to be hexagonal (Figure 1). A series of blocks was filled with paraffin and compared to a series filled with air. The blocks were evaluated with a “GE Lightspeed” 16 slice CT scanner and average CT# of the centers of the materials was determined. The attenuation properties of the subsequent blocks were also evaluated through their isocentric irradiation via “TrueBeam” accelerator under six beam energies. Blocks were placed upon plastic-water slabs of 4 cm in thickness assuring electronic equilibrium and data was collected via Sun Nuclear “Edge” diode detector. Relative changes in dose were compared with those predicted by Varian “Eclipse” TPS. Results: The CT# of 3D printed blocks was found to be a controllable variable. The fill material was able to narrow the range of variability in each sample. The attenuation of the block tracked with the density of the total fill structure. Assigned CT values in the TPS were seen to fall within an expected range predicted by the CT scans of the 3D printed blocks. Conclusion: We have demonstrated that it is possible to 3D print materials of varying tissue equivalencies, and that these materials have radiological properties that are customizable and predictable.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

SU-E-T-652: Quantification of Dosimetric Uncertainty of I-125 COMS Eye Plaque

Purpose: To compare dosimetrically three plan calculation systems (Plato, Varian Brachytherapy, and in-house-made Excel) available for I-125 COMS eye plaque treatment with measurement. Methods: All systems assume homogeneous media and calculations are based on a three-dimensional Cartesian coordinates, Plato and Brachytherapy Planning are based on AAPM TG-43 and the in-house Excel program only on inverse square corrections. Doses at specific depths were measured with EBT3 Gafchromic film from a fully loaded and a partially loaded 16 mm plaque (13 and 8 seeds respectively, I-125, model 6711 GE, Oncura). Measurements took place in a water tank, utilizing solid water blocks and a 3D-printed plaque holder. Taking advantage that gafchromic film has low energy dependence, a dose step wedge was delivered with 6 MV photon beam from a Varian 2100 EX linac for calibration. The gray-scale to dose in cGy was obtained with an Epson Expression 10000 XL scanner in the green channel. Treatment plans were generated for doses of 2200 cGy to a depth of 7 mm, and measurements were taken on a sagittal plane. Results: The calculated dose at the prescription point was 2242, 2344, and 2211 cGy with Excel, Brachyvision and Plato respectively for a fully loaded plaque, for the partially loaded plaque the doses were 2266, 2477, and 2193 cGy respectively. At 5 mm depth the doses for Brachyvision and Plato were comparable (3399 and 3267 cGy respectively), however, the measured dose in film was 3180 cGy which was lower by as much as 6.4% in the fully loaded plaque and 7.6% in the partially loaded plaque. Conclusion: Careful methodology and calibration are essential when measuring doses at specific depth due to the sensitivity and rapid dose fall off of I-125.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2015

TU-AB-201-10: Novel Conformal Skin Brachytherapy Device

Purpose: A novel conformal skin brachytherapy (CSBT) device was developed to provide patient specific treatment for small inoperable lesions and irregular surfaces that are not good candidates for electron external beam, e.g. eyelids, nose, lips, ears, etc. Methods: A prototype was built and tested using radioactive seeds spaced by 1.5 mm attached to a 3D printed template. Seeds were allowed to move independently at treatment surface to conform to target geometry. Beta-emitter Praseodymium-142 (half-life: 19.12 h, average beta energy: 809 keV) and Iodine-125 Model 6711 (half-life: 59.4 days, average gamma energy: 28 keV) seeds were used to test the feasibility of the device in solid water phantom for different prescription depths and surface areas. Percent depth doses and dose profiles were analyzed. MCNPX2.6 Monte Carlo Simulation and Gafchromic EBT3 film measurements were performed for a single Pr-142 seed. Treatment with nineteen I-125 seeds was planned in BrachyVision and measured experimentally for a total dose of 500 cGy at 5.0 mm depth. Results: For a single Pr-142 seed, dose was 500.0 cGy at surface and 17.4 cGy at 5.00 mm. BrachyVision planned and film doses for I-125 seed matrix were 500 cGy and 473.5 cGy at 5 mm depth and 171.0 cGy and 201.0 cGy at 10 mm depth. Total covered surface area for Pr-142 and I-125 were 2.35 mm{sup 2} and 88.24 mm{sup 2} respectively. Profile dose fall-off showed to be very conformal to contact area with 50% of the total surface dose from a single Pr-142 and array of I-125 seeds being deposited at 0.72 mm and 0.65 mm respectively from the central axis. Conclusion: CSBT device provides a high conformal dose to small surfaces. BrachyVision can be used to predict clinical dose distributions for multiple seeds matrix. Different radioactive seeds can be used to suit prescription depth and treatment area.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2011

Development and testing of improved statistical wind power forecasting methods.

Wind power forecasting (WPF) provides important inputs to power system operators and electricity market participants. It is therefore not surprising that WPF has attracted increasing interest within the electric power industry. In this report, we document our research on improving statistical WPF algorithms for point, uncertainty, and ramp forecasting. Below, we provide a brief introduction to the research presented in the following chapters. For a detailed overview of the state-of-the-art in wind power forecasting, we refer to [1]. Our related work on the application of WPF in operational decisions is documented in [2]. Point forecasts of wind power are highly dependent on the training criteria used in the statistical algorithms that are used to convert weather forecasts and observational data to a power forecast. In Chapter 2, we explore the application of information theoretic learning (ITL) as opposed to the classical minimum square error (MSE) criterion for point forecasting. In contrast to the MSE criterion, ITL criteria do not assume a Gaussian distribution of the forecasting errors. We investigate to what extent ITL criteria yield better results. In addition, we analyze time-adaptive training algorithms and how they enable WPF algorithms to cope with non-stationary data and, thus, to adapt to new situations without requiring additional offline training of the model. We test the new point forecasting algorithms on two wind farms located in the U.S. Midwest. Although there have been advancements in deterministic WPF, a single-valued forecast cannot provide information on the dispersion of observations around the predicted value. We argue that it is essential to generate, together with (or as an alternative to) point forecasts, a representation of the wind power uncertainty. Wind power uncertainty representation can take the form of probabilistic forecasts (e.g., probability density function, quantiles), risk indices (e.g., prediction risk index) or scenarios (with spatial and/or temporal dependence). Statistical approaches to uncertainty forecasting basically consist of estimating the uncertainty based on observed forecasting errors. Quantile regression (QR) is currently a commonly used approach in uncertainty forecasting. In Chapter 3, we propose new statistical approaches to the uncertainty estimation problem by employing kernel density forecast (KDF) methods. We use two estimators in both offline and time-adaptive modes, namely, the Nadaraya-Watson (NW) and Quantilecopula (QC) estimators. We conduct detailed tests of the new approaches using QR as a benchmark. One of the major issues in wind power generation are sudden and large changes of wind power output over a short period of time, namely ramping events. In Chapter 4, we perform a comparative study of existing definitions and methodologies for ramp forecasting. We also introduce a new probabilistic method for ramp event detection. The method starts with a stochastic algorithm that generates wind power scenarios, which are passed through a high-pass filter for ramp detection and estimation of the likelihood of ramp events to happen. The report is organized as follows: Chapter 2 presents the results of the application of ITL training criteria to deterministic WPF; Chapter 3 reports the study on probabilistic WPF, including new contributions to wind power uncertainty forecasting; Chapter 4 presents a new method to predict and visualize ramp events, comparing it with state-of-the-art methodologies; Chapter 5 briefly summarizes the main findings and contributions of this report.

17 WIND ENERGY↗