Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Implementation”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Design Considerations for the Implementation of a High-Field-Side Transient CHI System on QUEST

Transient coaxial helicity injection (T-CHI), a method first developed on the small helicity injected torus (HIT-II) experiment and then validated on the much larger National Spherical Torus Experiment (NSTX) device, is a method to initiate an inductive-like tokamak plasma discharge without reliance on the central solenoid. A CHI discharge is initiated by driving current along magnetic flux that connects the inner and outer divertor plates on one end of the tokamak. To permit this, on both HIT-II and NSTX, toroidal ceramic insulators were used to electrically separate the inner and outer vessel components. The use of such large toroidal vacuum insulators may not be easy to implement in reactors. To address this issue, the Q-shu University experiments with steady-state spherical tokamak (ST) (QUEST) is developing a reactor-relevant CHI configuration in which one of the divertor plates is electrically insulated from the rest of the vessel. The first application of T-CHI on QUEST biased the CHI electrode to the outer vessel. While the CHI discharges could be easily generated, it was found that as the discharge filled the vessel, the separation distance between the injector magnetic flux footprints widened, a condition that is not favorable for the generation of closed flux surfaces. Biasing the electrode to the inner wall is a configuration similar to that used on NSTX and HIT-II, but initial testing in this configuration has proved to be challenging. The design described here overcomes the present limitation by locating the CHI electrode much closer to the CHI injector flux coil and using an NSTX-like gas injection manifold to enable high-field-side T-CHI startup on QUEST. Finally, the concepts described in this article should also benefit the future implementation of T-CHI systems in other tokamaks and spherical tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Comparative Implementation of High Performance Computing for Power System Dynamic Simulations

Dynamic simulation for transient stability assessment is one of the most important, but intensive, computations for power system planning and operation. Present commercial software is mainly designed for sequential computation to run a single simulation, which is very time consuming with a single processer. The application of High Performance Computing (HPC) to dynamic simulations is very promising in accelerating the computing process by parallelizing its kernel algorithms while maintaining the same level of computation accuracy. This paper describes the comparative implementation of four parallel dynamic simulation schemes in two state-of-the-art HPC environments: Message Passing Interface (MPI) and Open Multi-Processing (OpenMP). These implementations serve to match the application with dedicated multi-processor computing hardware and maximize the utilization and benefits of HPC during the development process.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

IMPLEMENTATION OF DETAILED CHEMISTRY MECHANISMS IN ENGINE SIMULATIONS

The stiffness of large chemistry mechanisms has been proved to be a major hurdle toward predictive engine simulations. As a result, detailed chemistry mechanisms with a few thousand species need to be reduced based on target conditions so that they can be accommodated within the available computational resources. The computational cost of simulations typically increases super-linearly with the number of species and reactions. This work aims to bring detailed chemistry mechanisms within the realm of engine simulations by coupling the framework of unsteady flamelets and fast chemistry solvers. A previously developed tabulated flamelet model (TFM) framework for nonpremixed combustion was used in this study. The flamelet solver consists of the traditional operator-splitting scheme with variable coefficient ordinary differential equation (ODE) solver (VODE) and a numerical Jacobian for solving the chemistry. In order to use detailed mechanisms with thousands of species, a new framework with the Livermore solver for ODEs in sparse form (LSODES) chemistry solver and an analytical Jacobian was implemented in this work. Results from 1D simulations show that with the new framework, the computational cost is linearly proportional to the number of species in a given chemistry mechanism. As a result, the new framework is 2-3 orders of magnitude faster than the conventional variable coefficient ODE (VODE) solver for large chemistry mechanisms. This new framework was used to generate unsteady flamelet libraries for ndodecane using a detailed chemistry mechanism with 2755 species and 11,173 reactions. The engine combustion network (ECN) spray A experiments, which consist of an igniting n-dodecane spray in turbulent, high-pressure engine conditions are simulated using large eddy simulations (LES) coupled with detailed mechanisms. A grid with 0.06 mm minimum cell size and 22 x 10(6) peak cell count was implemented. The framework is validated across a range of ambient temperatures against ignition delay and liftoff lengths (LOLs). Qualitative results from the simulations were compared against experimental OH and CH2O planar laser-induced fluorescence (PLIF) data. The models are able to capture the spatial and temporal trends in species compared to those observed in the experiments. Quantitative and qualitative comparisons between the predictions of the reduced and detailed mechanisms are presented in detail. The main goal of this study is to demonstrate that detailed reaction mechanisms (similar to 1000 species) can now be used in engine simulations with a linear increase in computation cost with number of species during the tabulation process and a small increase in the 3D simulation cost.

Kundu, Prithwish↗
DOE OSTI2021

Cyber-Physical System Implementation for Manufacturing With Analytics in the Cloud Layer

Effective and efficient modern manufacturing operations require the acceptance and incorporation of the fourth industrial revolution, also known as Industry 4.0. Traditional shop floors are evolving their production into smart factories. To continue this trend, a specific architecture for the cyber-physical system is required, as well as a systematic approach to automate the application of algorithms and transform the acquired data into useful information. This work makes use of an approach that distinguishes three layers that are part of the existing Industry 4.0 paradigm: edge, fog, and cloud. Each of the layers performs computational operations, transforming the data produced in the smart factory into useful information. Trained or untrained methods for data analytics can be incorporated into the architecture. A case study is presented in which a real-time statistical control process algorithm based on control charts was implemented. The algorithm automatically detects changes in the material being processed in a computerized numerical control (CNC) machine. The algorithm implemented in the proposed architecture yielded short response times. The performance was effective since it automatically adapted to the machining of aluminum and then detected when the material was switched to steel. The data were backed up in a database that would allow traceability to the line of g-code that performed the machining.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Potential scenarios for broadening stakeholder involvement in the implementing geological disposal technology platform

This paper analyses the potential for the involvement of different types of stakeholders in the Implementing Geological Disposal Technology Platform (IGD-TP). This analysis was conducted as part of the InSOTEC project, a three-year (2011- 2014) collaborative research project funded under the 7. Euratom Framework Programme (Grant Agreement nr. 269906). In our analysis, we consider the extent to which the IGDTP's practice as regards to stakeholder involvement matches its discourse, and what potential for improvement exists given its structural organisation as a European Technology Platform (ETPs). Technology Platforms (TPs) can be understood as knowledge networks, deliberately set up to influence (research) policy in a specific domain. We therefore use knowledge networks as a conceptual approach and look at the IGD-TP as a complex network which includes actors, knowledge and practices across different countries, focusing on a very specific topic (i.e. implementing geological disposal). We compare the way different stakeholders are involved in the IGD-TP to the practice of other ETPs, and explore how the IGD-TP is viewed by its members and by outsiders to the platform Applying Callon's framework of knowledge co-production (1999) we come to define different degrees of interaction between science, society and policy in view of defining research and development (R and D) priorities [1]. Subsequently we describe how these interactions could be conceptualised and interpreted for the IGD-TP. The current approach of the IGDTP can be mainly understood as classical model involving mainly expert stakeholders and scientists. Where there seems to be a good representation among IGD-TP members of industry, research institutes, and some members of the academic community this is not the case for other types of stakeholders, such as public authorities or civil society. At this stage, the overall approach of the IGD-TP would seem to restrict the scope of stakeholder involvement, as it narrows participation down to uniquely technology experts, hindering socio-technical manifestations. Our analysis nevertheless shows that there is room for engaging with a broader range of stakeholders in the field of radioactive waste management, if this is the intention of the IGD-TP. However, this would require a commitment to developing a common knowledge base including other stakeholders through a process of mutual adjustment and negotiation. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Implementation and Verification of Cable Bending Stiffness in MoorDyn

The relatively large motions experienced by floating wind turbines and wave energy converters pose a challenge for power cables, whose internal components provide significant bending resistance and are sensitive to deformation. The behavior and associated design considerations of power cables in these highly dynamic applications make coupled analysis relevant for design. Bending stiffness capabilities have recently been added to the lumped-mass mooring dynamics model MoorDyn to enable simulation of dynamic power cables. MoorDyn is a common modeling choice for floating wind energy simulation (often coupled with OpenFAST) and floating wave energy converter simulation (often coupled with WEC-Sim) but the model’s previous line elasticity formulation only considered axial stiffness. To properly capture the dynamics of power cables, a bending stiffness model has been added that approximates cable curvature based on the difference in tangent vectors of adjacent elements. The resulting bending moment is realized by applying forces on adjacent nodes, enabling cable modeling while leaving the underlying lumped-mass formulation unchanged. In this paper, the new bending stiffness implementation is verified in static conditions against analytical solutions and then in a dynamic power cable scenario in comparison with the commercial simulator OrcaFlex. The dynamic scenario uses prescribed motions and includes wave loadings on the cable. Results indicate correct implementation of bending stiffness and show close agreement with OrcaFlex.

50 EE - Wind and Water Power Program - Water (EE-4↗
DOE OSTI2014

SU-D-18C-06: Initial Experience with Implementing MRI Safety Guidelines for Patients with Pacemakers - Medical Physicist Perspective

Purpose: Several institutions have developed MRI guidelines for patients with MR-unsafe or MR-conditional pacemakers. Here we highlight the role of a medical physicist in implementing these guidelines for non-pacemaker dependent patients. Guidelines: Implementing these guidelines requires involvement from several medical specialties and a strong collaboration with the site MRI supervisor to develop a structured workflow. A medical physicist is required to be present during the scan to supervise the MR scanning and to maintain a safety checklist that ensures: 1) uninterrupted patient communication with the technologist, 2) continuous patient physiologic monitoring (e.g. blood pressure and electrocardiography) by a trained nurse, 3) redundant patient vitals monitoring (e.g. pulse oximetry) due to the possibility of in vivo electrocardiography reading fluctuations during image acquisition. A radiologist is strongly recommended to be available to review the images before patients are discharged from the scanner. Pacemaker MRI should be restricted to 1.5T field strength. The MRI sequences should be optimized by the physicist with regards to: a) SAR: limited to <1.5 W/Kg for MR-unsafe pacemakers in normal operating mode, b) RF exposure time: <30 min, c) Coils: use T/R coils but not restricted to such, d) Artifacts: further optimization of sequences whenever image quality is compromised due to the pacemaker. In particular, cardiac, breast and left-shoulder MRIs are most susceptible to these artifacts. Possible strategies to lower the SAR include: a) BW reduction, 2) echo-train-length reduction, 3) increase TR, 4) decrease number of averages, 5) decrease flip angle, 6) reduce slices and/or a combination of all the options. Conclusion: A medical physicist in collaboration with the MR supervisor plays an important role in the supervision/implementation of safe MR scanning of pacemaker patients. Developing and establishing a workflow has enabled our institution to scan over 30 patients with pacemakers without complications, including 3 cardiac MR exams.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-P-03: Implementing a Low Dose Lung Screening CT Program Meeting Regulatory Requirements

Purpose: To provide information pertaining to IROC Houston QA Center's (RPC) credentialing process for institutions participating in NCI-sponsored clinical trials. Purpose: Provide guidance to the Radiology Departments with the intent of implementing a Low Dose CT Screening Program using different CT Scanners with multiple techniques within the framework of the required state regulations. Method: State Requirements for the purpose of implementing a Low Dose CT Lung Protocol required working with the Radiology and Pulmonary Department in setting up a Low Dose Screening Protocol designed to reduce the radiation burden to the patients enrolled. Radiation dose measurements (CTDIvol) for various CT manufacturers (Siemens16, Siemens 64, Philips 64, and Neusoft128) for three different weight based protocols. All scans were reviewed by the Radiologist. Prior to starting a low dose lung screening protocol, information had to be submitted to the state for approval. Performing a Healing Arts protocol requires extensive information. This not only includes name and address of the applicant but a detailed description of the disease, the x-ray examination and the population to be examined. The unit had to be tested by a qualified expert using the technique charts. The credentials of all the operators, the supervisors and the Radiologists had to be submitted to the state. Results: All the appropriate documentation was sent to the state for review. The measured results between the Low Dose Protocol versus the default Adult Chest Protocol showed that there was a dose reduction of 65% for small (100-150 lb.) patient, 75% for the Medium patient (151-250 lbs.), and a 55% reduction for the Large patient ( over 250 lbs.). Conclusion: Measured results indicated that the Low Dose Protocol indeed lowered the screening patient's radiation dose and the institution was able to submit the protocol to the State's regulators.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

TH-A-19A-11: Validation of GPU-Based Monte Carlo Code (gPMC) Versus Fully Implemented Monte Carlo Code (TOPAS) for Proton Radiation Therapy: Clinical Cases Study

Purpose: For proton radiation therapy, Monte Carlo simulation (MCS) methods are recognized as the gold-standard dose calculation approach. Although previously unrealistic due to limitations in available computing power, GPU-based applications allow MCS of proton treatment fields to be performed in routine clinical use, on time scales comparable to that of conventional pencil-beam algorithms. This study focuses on validating the results of our GPU-based code (gPMC) versus fully implemented proton therapy based MCS code (TOPAS) for clinical patient cases. Methods: Two treatment sites were selected to provide clinical cases for this study: head-and-neck cases due to anatomical geometrical complexity (air cavities and density heterogeneities), making dose calculation very challenging, and prostate cases due to higher proton energies used and close proximity of the treatment target to sensitive organs at risk. Both gPMC and TOPAS methods were used to calculate 3-dimensional dose distributions for all patients in this study. Comparisons were performed based on target coverage indices (mean dose, V90 and D90) and gamma index distributions for 2% of the prescription dose and 2mm. Results: For seven out of eight studied cases, mean target dose, V90 and D90 differed less than 2% between TOPAS and gPMC dose distributions. Gamma index analysis for all prostate patients resulted in passing rate of more than 99% of voxels in the target. Four out of five head-neck-cases showed passing rate of gamma index for the target of more than 99%, the fifth having a gamma index passing rate of 93%. Conclusion: Our current work showed excellent agreement between our GPU-based MCS code and fully implemented proton therapy based MC code for a group of dosimetrically challenging patient cases.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

TH-E-BRE-07: Development of Dose Calculation Error Predictors for a Widely Implemented Clinical Algorithm

Purpose: Several shortcomings of the current implementation of the analytic anisotropic algorithm (AAA) may lead to dose calculation errors in highly modulated treatments delivered to highly heterogeneous geometries. Here we introduce a set of dosimetric error predictors that can be applied to a clinical treatment plan and patient geometry in order to identify high risk plans. Once a problematic plan is identified, the treatment can be recalculated with more accurate algorithm in order to better assess its viability. Methods: Here we focus on three distinct sources dosimetric error in the AAA algorithm. First, due to a combination of discrepancies in smallfield beam modeling as well as volume averaging effects, dose calculated through small MLC apertures can be underestimated, while that behind small MLC blocks can overestimated. Second, due the rectilinear scaling of the Monte Carlo generated pencil beam kernel, energy is not properly transported through heterogeneities near, but not impeding, the central axis of the beamlet. And third, AAA overestimates dose in regions very low density (< 0.2 g/cm{sup 3}). We have developed an algorithm to detect the location and magnitude of each scenario within the patient geometry, namely the field-size index (FSI), the heterogeneous scatter index (HSI), and the lowdensity index (LDI) respectively. Results: Error indices successfully identify deviations between AAA and Monte Carlo dose distributions in simple phantom geometries. Algorithms are currently implemented in the MATLAB computing environment and are able to run on a typical RapidArc head and neck geometry in less than an hour. Conclusion: Because these error indices successfully identify each type of error in contrived cases, with sufficient benchmarking, this method can be developed into a clinical tool that may be able to help estimate AAA dose calculation errors and when it might be advisable to use Monte Carlo calculations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-C-BRD-02: A Team Focused Clinical Implementation and Failure Mode and Effects Analysis of HDR Skin Brachytherapy Using Valencia and Leipzig Surface Applicators

Purpose: and Leipzig applicators (VLAs) are single-channel brachytherapy surface applicators used to treat skin lesions up to 2cm diameter. Source dwell times can be calculated and entered manually after clinical set-up or ultrasound. This procedure differs dramatically from CT-based planning; the novelty and unfamiliarity could lead to severe errors. To build layers of safety and ensure quality, a multidisciplinary team created a protocol and applied Failure Modes and Effects Analysis (FMEA) to the clinical procedure for HDR VLA skin treatments. Methods: team including physicists, physicians, nurses, therapists, residents, and administration developed a clinical procedure for VLA treatment. The procedure was evaluated using FMEA. Failure modes were identified and scored by severity, occurrence, and detection. The clinical procedure was revised to address high-scoring process nodes. Results: Several key components were added to the clinical procedure to minimize risk probability numbers (RPN): -Treatments are reviewed at weekly QA rounds, where physicians discuss diagnosis, prescription, applicator selection, and set-up. Peer review reduces the likelihood of an inappropriate treatment regime. -A template for HDR skin treatments was established in the clinical EMR system to standardize treatment instructions. This reduces the chances of miscommunication between the physician and planning physicist, and increases the detectability of an error during the physics second check. -A screen check was implemented during the second check to increase detectability of an error. -To reduce error probability, the treatment plan worksheet was designed to display plan parameters in a format visually similar to the treatment console display. This facilitates data entry and verification. -VLAs are color-coded and labeled to match the EMR prescriptions, which simplifies in-room selection and verification. Conclusion: Multidisciplinary planning and FMEA increased delectability and reduced error probability during VLA HDR Brachytherapy. This clinical model may be useful to institutions implementing similar procedures.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Initial implementation of the conversion from the energy-subtracted CT number to electron density in tissue inhomogeneity corrections: An anthropomorphic phantom study of radiotherapy treatment planning

Purpose: To achieve accurate tissue inhomogeneity corrections in radiotherapy treatment planning, the authors had previously proposed a novel conversion of the energy-subtracted computed tomography (CT) number to an electron density (ΔHU–ρ{sub e} conversion), which provides a single linear relationship between ΔHU and ρ{sub e} over a wide range of ρ{sub e}. The purpose of this study is to present an initial implementation of the ΔHU–ρ{sub e} conversion method for a treatment planning system (TPS). In this paper, two example radiotherapy plans are used to evaluate the reliability of dose calculations in the ΔHU–ρ{sub e} conversion method. Methods: CT images were acquired using a clinical dual-source CT (DSCT) scanner operated in the dual-energy mode with two tube potential pairs and an additional tin (Sn) filter for the high-kV tube (80–140 kV/Sn and 100–140 kV/Sn). Single-energy CT using the same DSCT scanner was also performed at 120 kV to compare the ΔHU–ρ{sub e} conversion method with a conventional conversion from a CT number to ρ{sub e} (Hounsfield units, HU–ρ{sub e} conversion). Lookup tables for ρ{sub e} calibration were obtained from the CT image acquisitions for tissue substitutes in an electron density phantom (EDP). To investigate the beam-hardening effect on dosimetric uncertainties, two EDPs with different sizes (a body EDP and a head EDP) were used for the ρ{sub e} calibration. Each acquired lookup table was applied to two radiotherapy plans designed using the XiO TPS with the superposition algorithm for an anthropomorphic phantom. The first radiotherapy plan was for an oral cavity tumor and the second was for a lung tumor. Results: In both treatment plans, the performance of the ΔHU–ρ{sub e} conversion was superior to that of the conventional HU–ρ{sub e} conversion in terms of the reliability of dose calculations. Especially, for the oral tumor plan, which dealt with dentition and bony structures, treatment planning with the HU–ρ{sub e} conversion exhibited apparent discrepancies between the dose distributions and dose–volume histograms (DVHs) of the body-EDP and head-EDP calibrations. In contrast, the dose distributions and DVHs of the body-EDP and head-EDP calibrations coincided with each other almost perfectly in the ΔHU–ρ{sub e} conversion for 100–140 kV/Sn. The difference between the V{sub 100}’s (the mean planning target volume receiving 100% of the prescribed dose; a DVH parameter) of the body-EDP and head-EDP calibrations could be reduced to less than 1% using the ΔHU–ρ{sub e} conversion, but exceeded 11% for the HU–ρ{sub e} conversion. Conclusions: The ΔHU–ρ{sub e} conversion can be implemented for currently available TPS’s without any modifications or extensions. The ΔHU–ρ{sub e} conversion appears to be a promising method for providing an accurate and reliable inhomogeneity correction in treatment planning for any ill-conditioned scans that include (i) the use of a calibration EDP that is nonequivalent to the patient’s body tissues, (ii) a mismatch between the size of the patient and the calibration EDP, or (iii) a large quantity of high-density and high-atomic-number tissue structures.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-10: A Clinical Implementation and the Dosimetric Evidence in High Dose Rate Vaginal Multichannel Applicator Brachytherapy

Purpose: The multichannel cylindrical applicator has a distinctive modification of the traditional single channel cylindrical applicator. The novel multichannel applicator has additional peripheral channels that provide more flexibility both in treatment planning process and outcomes. To protect by reducing doses to adjacent organ at risk (OAR) while maintaining target coverage with inverse plan optimization are the goals for such novel Brachytherapy device. Through a series of comparison and analysis of reults in more than forty patients who received HDR Brachytherapy using multichannel vaginal applicator, this procedure has been implemented in our institution. Methods: Multichannel planning was CT image based. The CTV of 5mm vaginal cuff rind with prescribed length was well reconstructed as well as bladder and rectum. At least D95 of CTV coverage is 95% of prescribed dose. Multichannel inverse plan optimization algorithm not only shapes target dose cloud but set dose avoids to OAR’s exclusively. The doses of D2cc, D5cc and D5; volume of V2Gy in OAR’s were selected to compare with single channel results when sole central channel is only possibility. Results: Study demonstrates plan superiorly in OAR’s doe reduction in multi-channel plan. The D2cc of the rectum and bladder were showing a little lower for multichannel vs. single channel. The V2Gy of the rectum was 93.72% vs. 83.79% (p=0.007) for single channel vs. multichannel respectively. Absolute reduced mean dose of D5 by multichannel was 17 cGy (s.d.=6.4) and 44 cGy (s.d.=15.2) in bladder and rectum respectively. Conclusion: The optimization solution in multichannel was to maintain D95 CTV coverage while reducing the dose to OAR’s. Dosimetric advantage in sparing critical organs by using a multichannel applicator in HDR Brachytherapy treatment of the vaginal cuff is so promising and has been implemented clinically.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-121: Analyzing the Broadening Effect On the Bragg Peak Due to Heterogeneous Geometries and Implementing User-Routines in the Monte-Carlo Code FLUKA in Order to Reduce Computation Time

Purpose: Aim of this study was to analyze the modulating, broadening effect on the Bragg Peak due to heterogeneous geometries like multi-wire chambers in the beam path of a particle therapy beam line. The effect was described by a mathematical model which was implemented in the Monte-Carlo code FLUKA via user-routines, in order to reduce the computation time for the simulations. Methods: The depth dose curve of 80 MeV/u C12-ions in a water phantom was calculated using the Monte-Carlo code FLUKA (reference curve). The modulating effect on this dose distribution behind eleven mesh-like foils (periodicity ∼80 microns) occurring in a typical set of multi-wire and dose chambers was mathematically described by optimizing a normal distribution so that the reverence curve convoluted with this distribution equals the modulated dose curve. This distribution describes a displacement in water and was transferred in a probability distribution of the thickness of the eleven foils using the water equivalent thickness of the foil’s material. From this distribution the distribution of the thickness of one foil was determined inversely. In FLUKA the heterogeneous foils were replaced by homogeneous foils and a user-routine was programmed that varies the thickness of the homogeneous foils for each simulated particle using this distribution. Results: Using the mathematical model and user-routine in FLUKA the broadening effect could be reproduced exactly when replacing the heterogeneous foils by homogeneous ones. The computation time was reduced by 90 percent. Conclusion: In this study the broadening effect on the Bragg Peak due to heterogeneous structures was analyzed, described by a mathematical model and implemented in FLUKA via user-routines. Applying these routines the computing time was reduced by 90 percent. The developed tool can be used for any heterogeneous structure in the dimensions of microns to millimeters, in principle even for organic materials like lung tissue.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-329: Dosimetric Impact of Implementing Metal Artifact Reduction Methods and Metal Energy Deposition Kernels for Photon Dose Calculations

Purpose: To investigate two strategies for reducing dose calculation errors near metal implants: use of CT metal artifact reduction methods and implementation of metal-based energy deposition kernels in the convolution/superposition (C/S) method. Methods: Radiochromic film was used to measure the dose upstream and downstream of titanium and Cerrobend implants. To assess the dosimetric impact of metal artifact reduction methods, dose calculations were performed using baseline, uncorrected images and metal artifact reduction Methods: Philips O-MAR, GE’s monochromatic gemstone spectral imaging (GSI) using dual-energy CT, and GSI imaging with metal artifact reduction software applied (MARs).To assess the impact of metal kernels, titanium and silver kernels were implemented into a commercial collapsed cone C/S algorithm. Results: The CT artifact reduction methods were more successful for titanium than Cerrobend. Interestingly, for beams traversing the metal implant, we found that errors in the dimensions of the metal in the CT images were more important for dose calculation accuracy than reduction of imaging artifacts. The MARs algorithm caused a distortion in the shape of the titanium implant that substantially worsened the calculation accuracy. In comparison to water kernel dose calculations, metal kernels resulted in better modeling of the increased backscatter dose at the upstream interface but decreased accuracy directly downstream of the metal. We also found that the success of metal kernels was dependent on dose grid size, with smaller calculation voxels giving better accuracy. Conclusion: Our study yielded mixed results, with neither the metal artifact reduction methods nor the metal kernels being globally effective at improving dose calculation accuracy. However, some successes were observed. The MARs algorithm decreased errors downstream of Cerrobend by a factor of two, and metal kernels resulted in more accurate backscatter dose upstream of metals. Thus, these two strategies do have the potential to improve accuracy for patients with metal implants in certain scenarios. This work was supported by Public Health Service grants CA 180803 and CA 10953 awarded by the National Cancer Institute, United States of Health and Human Services, and in part by Mobius Medical Systems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-465: Implementation of An Automated Collision Detection Program Using Open Source Software for the Pinnacle Treatment Planning System

Purpose: Potential collisions between the gantry head and the patient or table assembly are difficult to detect in most treatment planning systems. We have developed and implemented a novel software package for the representation of potential gantry collisions with the couch assembly at the time of treatment planning. Methods: Physical dimensions of the Varian Edge linear accelerator treatment head were measured and reproduced using the Visual Python display package. A script was developed for the Pinnacle treatment planning system to generate a file with the relevant couch, gantry, and isocenter positions for each beam in a planning trial. A python program was developed to parse the information from the TPS and produce a representative model of the couch/gantry system. Using the model and the Visual Python libraries, a rendering window is generated for each beam that allows the planner to evaluate the possibility of a collision. Results: Comparison against heuristic methods and direct verification on the machine validated the collision model generated by the software. Encounters of <1 cm between the gantry treatment head and table were visualized as collisions in our virtual model. Visual windows were created depicting the angle of collision for each beam, including the anticipated table coordinates. Visual rendering of a 6 arc trial with multiple couch positions was completed in under 1 minute, with network bandwidth being the primary bottleneck. Conclusion: The developed software allows for quick examination of possible collisions during the treatment planning process and helps to prevent major collisions prior to plan approval. The software can easily be implemented on future planning systems due to the versatility and platform independence of the Python programming language. Further integration of the software with the treatment planning system will allow the possibility of patient-gantry collision detection for a range of treatment machines.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-467: Implementation of Monte Carlo Dose Calculation for a Multileaf Collimator Equipped Robotic Radiotherapy System

Purpose: To improve the treatment efficiency and capabilities for full-body treatment, a robotic radiosurgery system has equipped with a multileaf collimator (MLC) to extend its accuracy and precision to radiation therapy. To model the MLC and include it in the Monte Carlo patient dose calculation is the goal of this work. Methods: The radiation source and the MLC were carefully modeled to consider the effects of the source size, collimator scattering, leaf transmission and leaf end shape. A source model was built based on the output factors, percentage depth dose curves and lateral dose profiles measured in a water phantom. MLC leaf shape, leaf end design and leaf tilt for minimizing the interleaf leakage and their effects on beam fluence and energy spectrum were all considered in the calculation. Transmission/leakage was added to the fluence based on the transmission factors of the leaf and the leaf end. The transmitted photon energy was tuned to consider the beam hardening effects. The calculated results with the Monte Carlo implementation was compared with measurements in homogeneous water phantom and inhomogeneous phantoms with slab lung or bone material for 4 square fields and 9 irregularly shaped fields. Results: The calculated output factors are compared with the measured ones and the difference is within 1% for different field sizes. The calculated dose distributions in the phantoms show good agreement with measurements using diode detector and films. The dose difference is within 2% inside the field and the distance to agreement is within 2mm in the penumbra region. The gamma passing rate is more than 95% with 2%/2mm criteria for all the test cases. Conclusion: Implementation of Monte Carlo dose calculation for a MLC equipped robotic radiosurgery system is completed successfully. The accuracy of Monte Carlo dose calculation with MLC is clinically acceptable. This work was supported by Accuray Inc.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-500: Initial Implementation of GPU-Based Particle Swarm Optimization for 4D IMRT Planning in Lung SBRT

Purpose 4D-IMRT planning, combined with dynamic MLC tracking delivery, utilizes the temporal dimension as an additional degree of freedom to achieve improved OAR-sparing. The computational complexity for such optimization increases exponentially with increase in dimensionality. In order to accomplish this task in a clinically-feasible time frame, we present an initial implementation of GPU-based 4D-IMRT planning based on particle swarm optimization (PSO). Methods The target and normal structures were manually contoured on ten phases of a 4DCT scan of a NSCLC patient with a 54cm3 right-lower-lobe tumor (1.5cm motion). Corresponding ten 3D-IMRT plans were created in the Eclipse treatment planning system (Ver-13.6). A vendor-provided scripting interface was used to export 3D-dose matrices corresponding to each control point (10 phases × 9 beams × 166 control points = 14,940), which served as input to PSO. The optimization task was to iteratively adjust the weights of each control point and scale the corresponding dose matrices. In order to handle the large amount of data in GPU memory, dose matrices were sparsified and placed in contiguous memory blocks with the 14,940 weight-variables. PSO was implemented on CPU (dual-Xeon, 3.1GHz) and GPU (dual-K20 Tesla, 2496 cores, 3.52Tflops, each) platforms. NiftyReg, an open-source deformable image registration package, was used to calculate the summed dose. Results The 4D-PSO plan yielded PTV coverage comparable to the clinical ITV-based plan and significantly higher OAR-sparing, as follows: lung Dmean=33%; lung V20=27%; spinal cord Dmax=26%; esophagus Dmax=42%; heart Dmax=0%; heart Dmean=47%. The GPU-PSO processing time for 14940 variables and 7 PSO-particles was 41% that of CPU-PSO (199 vs. 488 minutes). Conclusion Truly 4D-IMRT planning can yield significant OAR dose-sparing while preserving PTV coverage. The corresponding optimization problem is large-scale, non-convex and computationally rigorous. Our initial results indicate that GPU-based PSO with further software optimization can make such planning clinically feasible. This work was supported through funding from the National Institutes of Health and Varian Medical Systems.

60 APPLIED LIFE SCIENCES↗