Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Optimization”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Optimal Control Strategy With Efficiency and Reliability Improvement for Offshore DC Microgrids

Offshore microgrids, due to their remote location and lack of external energy support, face significant challenges in wide-range load operation and maintenance. Consequently, efficiency and reliability are critical concerns for converters in offshore dc microgrids. This article presents an optimal control strategy aimed at enhancing both efficiency and reliability. A normalized nonlinear relationship between power loss and thermal stress of a paralleled converter is first established. Based on this, a dual-objective optimization function with an active weight function as well as a system overall performance index is established. The active weight function dynamically adjusts the control priority based on converter efficiency and switching device thermal stress. Then, the optimal power-sharing strategy is derived by the Lagrange multiplier method with the proposed optimal function. Additionally, to accommodate a wide load range, an optimal selection strategy for operating converter combinations is proposed, requiring only low-bandwidth communication. Experiment verification is given to validate the effectiveness of the proposed control strategy. The experiment results demonstrate that the proposed control strategy can improve the overall performance of offshore microgrids by optimizing efficiency and reliability.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Uncertainty propagation and sensitivity analysis for constrained optimization of nuclear waste vitrification

Abstract The vitrification of high‐level waste (HLW) by heating a mixture of glass‐forming chemicals (GFCs) with the waste can be improved using a constrained optimization problem. This study explores how different uncertainty propagation (UP) methods implemented with the optimization process can affect the glass formulation of nuclear waste glasses. UP is the effort of propagating uncertain inputs through a system to understand and quantify output distributions. Uncertainty intervals are crafted from output distributions to inform the optimization algorithm. UP is often implemented with Monte Carlo (MC) sampling for large nonlinear systems, which can be difficult to implement within a constrained optimization algorithm that requires derivative information. Other UP methods often used for optimization under uncertainty (OUU) can be designed to work within an established constrained optimization framework. Methods of UP are evaluated in this study including iterative sampling approaches, first‐order approximations, and surrogate modeling with machine learning (ML). A method of dimensional reduction based on global sensitivity analysis is introduced to support the UP methods for the large dimensionality of the problem. Analytical UP methods able to achieve similar optimums 10 times faster than the baseline MC approach, and produce 93.9% similar output distributions are reported.

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

Optimization of a Cyclone Using Multiphase Flow Computational Fluid Dynamics

The U.S. Department of Energy National Energy Technology Laboratory's (NETL) 50 kW th chemical looping reactor (CLR) has an underperforming cyclone, which was designed using empirical correlations. To improve the performance of this cyclone using computational fluid dynamics (CFD)-based modeling simulations, four critical design parameters including the vortex tube radius and length, barrel radius, and the inlet width and height were optimized. Here, NETL's open source multiphase flow with interphase exchange (MFiX) CFD code has been used to model a series of cyclones by systematically varying the geometric design parameters. To perform the optimization process, the surrogate modeling and sensitivity analysis followed by the optimization capability in nodeworks was used. The basic methodology for the process is to employ a statistical design of experiments (DOE) method to generate sampling simulations that fill the design space. Corresponding CFD models are then created, executed, and postprocessed. A response surface is created to characterize the relationship between input parameters and the quantities of interest (QoI). Finally, the CFD-surrogate is used by an optimization method to find the optimal design condition based on the objective and constraints prescribed. The resulting optimal cyclone has a larger diameter and longer vortex tube, a larger diameter barrel, and a taller and narrower solids inlet. The improved design has a predicted pressure drop 11 times lower than the original design while reducing the mass loss by a factor of 2.3.

42 ENGINEERING↗
DOE OSTI2013

Proton beam scattering system optimization for clinical and research applications

Purpose: To develop and test a method for optimizing and constructing a dual scattering system in passively scattered proton therapy. Methods: A beam optics optimization algorithm was developed to optimize the thickness of the first scatterer (S1) and the profile (of both the high-Z material and Lexan) of the second scatterer (S2) to deliver a proton beam matching a given set of parameters, including field diameter, fluence, flatness, and symmetry. A new manufacturing process was also tested that allows the contoured second scattering foil to be created much more economically and quickly using Cerrobend casting. Two application-specific scattering systems were developed and tested using both experimental and Monte Carlo techniques to validate the optimization process described. Results: A scattering system was optimized and constructed to deliver large uniform irradiations of radiobiology samples at low dose rates. This system was successfully built and tested using film and ionization chambers. The system delivered a uniform radiation field of 50 cm diameter (to a dose of {+-}7% of the central axis) while the depth dose profile could be tuned to match the specifications of the particular investigator using modulator wheels and range shifters. A second scattering system for intermediate field size (4 cm < diameter < 10 cm) stereotactic radiosurgery and radiation therapy (SRS and SRT) treatments was also developed and tested using GEANT4. This system improved beam efficiency by over 70% compared with existing scattering systems while maintaining field flatness and depth dose profile. In both cases the proton range uniformity across the radiation field was maintained, further indicating the accuracy of the energy loss formalism in the optimization algorithm. Conclusions: The methods described allow for rapid prototyping of scattering foils to meet the demands of both research and clinical beam delivery applications in proton therapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Adaptation of the CVT algorithm for catheter optimization in high dose rate brachytherapy

Purpose: An innovative, simple, and fast method to optimize the number and position of catheters is presented for prostate and breast high dose rate (HDR) brachytherapy, both for arbitrary templates or template-free implants (such as robotic templates).Methods: Eight clinical cases were chosen randomly from a bank of patients, previously treated in our clinic to test our method. The 2D Centroidal Voronoi Tessellations (CVT) algorithm was adapted to distribute catheters uniformly in space, within the maximum external contour of the planning target volume. The catheters optimization procedure includes the inverse planning simulated annealing algorithm (IPSA). Complete treatment plans can then be generated from the algorithm for different number of catheters. The best plan is chosen from different dosimetry criteria and will automatically provide the number of catheters and their positions. After the CVT algorithm parameters were optimized for speed and dosimetric results, it was validated against prostate clinical cases, using clinically relevant dose parameters. The robustness to implantation error was also evaluated. Finally, the efficiency of the method was tested in breast interstitial HDR brachytherapy cases.Results: The effect of the number and locations of the catheters on prostate cancer patients was studied. Treatment plans with a better or equivalent dose distributions could be obtained with fewer catheters. A better or equal prostate V100 was obtained down to 12 catheters. Plans with nine or less catheters would not be clinically acceptable in terms of prostate V100 and D90. Implantation errors up to 3 mm were acceptable since no statistical difference was found when compared to 0 mm error (p > 0.05). No significant difference in dosimetric indices was observed for the different combination of parameters within the CVT algorithm. A linear relation was found between the number of random points and the optimization time of the CVT algorithm. Because the computation time decrease with the number of points and that no effects were observed on the dosimetric indices when varying the number of sampling points and the number of iterations, they were respectively fixed to 2500 and to 100. The computation time to obtain ten complete treatments plans ranging from 9 to 18 catheters, with the corresponding dosimetric indices, was 90 s. However, 93% of the computation time is used by a research version of IPSA. For the breast, on average, the Radiation Therapy Oncology Group recommendations would be satisfied down to 12 catheters. Plans with nine or less catheters would not be clinically acceptable in terms of V100, dose homogeneity index, and D90.Conclusions: The authors have devised a simple, fast and efficient method to optimize the number and position of catheters in interstitial HDR brachytherapy. The method was shown to be robust for both prostate and breast HDR brachytherapy. More importantly, the computation time of the algorithm is acceptable for clinical use. Ultimately, this catheter optimization algorithm could be coupled with a 3D ultrasound system to allow real-time guidance and planning in HDR brachytherapy.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

A DVH-guided IMRT optimization algorithm for automatic treatment planning and adaptive radiotherapy replanning

Purpose: To develop a novel algorithm that incorporates prior treatment knowledge into intensity modulated radiation therapy optimization to facilitate automatic treatment planning and adaptive radiotherapy (ART) replanning. Methods: The algorithm automatically creates a treatment plan guided by the DVH curves of a reference plan that contains information on the clinician-approved dose-volume trade-offs among different targets/organs and among different portions of a DVH curve for an organ. In ART, the reference plan is the initial plan for the same patient, while for automatic treatment planning the reference plan is selected from a library of clinically approved and delivered plans of previously treated patients with similar medical conditions and geometry. The proposed algorithm employs a voxel-based optimization model and navigates the large voxel-based Pareto surface. The voxel weights are iteratively adjusted to approach a plan that is similar to the reference plan in terms of the DVHs. If the reference plan is feasible but not Pareto optimal, the algorithm generates a Pareto optimal plan with the DVHs better than the reference ones. If the reference plan is too restricting for the new geometry, the algorithm generates a Pareto plan with DVHs close to the reference ones. In both cases, the new plans have similar DVH trade-offs as the reference plans. Results: The algorithm was tested using three patient cases and found to be able to automatically adjust the voxel-weighting factors in order to generate a Pareto plan with similar DVH trade-offs as the reference plan. The algorithm has also been implemented on a GPU for high efficiency. Conclusions: A novel prior-knowledge-based optimization algorithm has been developed that automatically adjust the voxel weights and generate a clinical optimal plan at high efficiency. It is found that the new algorithm can significantly improve the plan quality and planning efficiency in ART replanning and automatic treatment planning.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

SU-E-T-23: A Novel Two-Step Optimization Scheme for Tandem and Ovoid (T and O) HDR Brachytherapy Treatment for Locally Advanced Cervical Cancer

Purpose: To present a novel method allowing fast, true volumetric optimization of T and O HDR treatments and to quantify its benefits. Materials and Methods: 27 CT planning datasets and treatment plans from six consecutive cervical cancer patients treated with 4–5 intracavitary T and O insertions were used. Initial treatment plans were created with a goal of covering high risk (HR)-CTV with D90 > 90% and minimizing D2cc to rectum, bladder and sigmoid with manual optimization, approved and delivered. For the second step, each case was re-planned adding a new structure, created from the 100% prescription isodose line of the manually optimized plan to the existent physician delineated HR-CTV, rectum, bladder and sigmoid. New, more rigorous DVH constraints for the critical OARs were used for the optimization. D90 for the HR-CTV and D2cc for OARs were evaluated in both plans. Results: Two-step optimized plans had consistently smaller D2cc's for all three OARs while preserving good D90s for HR-CTV. On plans with “excellent” CTV coverage, average D90 of 96% (range 91–102), sigmoid D2cc was reduced on average by 37% (range 16–73), bladder by 28% (range 20–47) and rectum by 27% (range 15–45). Similar reductions were obtained on plans with “good” coverage, with an average D90 of 93% (range 90–99). For plans with inferior coverage, average D90 of 81%, an increase in coverage to 87% was achieved concurrently with D2cc reductions of 31%, 18% and 11% for sigmoid, bladder and rectum. Conclusions: A two-step DVH-based optimization can be added with minimal planning time increase, but with the potential of dramatic and systematic reductions of D2cc for OARs and in some cases with concurrent increases in target dose coverage. These single-fraction modifications would be magnified over the course of 4–5 intracavitary insertions and may have real clinical implications in terms of decreasing both acute and late toxicity.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-T-367: Optimization of DLG Using TG-119 Test Cases and a Weighted Mean Approach

Purpose: Optimization of the dosimetric leaf gap (DLG) is an important step in commissioning the Eclipse treatment planning system for sliding window intensity-modulated radiation therapy (SW-IMRT) and RapidArc. Often the values needed for optimal dose delivery differ markedly from those measured at commissioning. We present a method to optimize this value using the AAPM TG-119 test cases. Methods: For SW-IMRT and RapidArc, TG-119 based test plans were created using a water-equivalent phantom. Dose distributions measured on film and ion chamber (IC) readings taken in low-gradient regions within the targets were analyzed separately. Since DLG is a single value per energy, SW-IMRT and RapidArc must be considered simultaneously. Plans were recalculated using a linear sweep from 0.02cm (the minimum DLG) to 0.3 cm. The calculated point doses were compared to the measured doses for each plan, and based on these comparisons an optimal DLG value was computed for each plan. TG-119 cases are designed to push the system in various ways, thus, a weighted mean of the DLG was computed where the relative importance of each type of plan was given a score from 0.0 to 1.0. Finally, SW-IMRT and RapidArc are assigned an overall weight based on clinical utilization. Our routine patient-QA (PQA) process was performed as independent validation. Results: For a Varian TrueBeam, the optimized DLG varied with σ = 0.044cm for SW-IMRT and σ = 0.035cm for RapidArc. The difference between the weighted mean SW-IMRT and RapidArc value was 0.038cm. We predicted utilization of 25% SW-IMRT and 75% RapidArc. The resulting DLG was ~1mm different than that found by commissioning and produced an average error of <1% for SW-IMRT and RapidArc PQA test cases separately. Conclusion: The weighted mean method presented is a useful tool for determining an optimal DLG value for commissioning Eclipse.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-T-593: Clinical Evaluation of Direct Aperture Optimization in Head/Neck and Prostate IMRT Treatment

Purpose: To investigate the planning and dosimetric advantages of direct aperture optimization (DAO) over beam-let optimization in IMRT treatment of head and neck (H/N) and prostate cancers. Methods: Five Head and Neck as well as five prostate patients were planned using the beamlet optimizer in Elekta-Xio ver 4.6 IMRT treatment planning system. Based on our experience in beamlet IMRT optimization, PTVs in H/N plans were prescribed to 70 Gy delivered by 7 fields. While prostate PTVs were prescribed to 76 Gy with 9 fields. In all plans, fields were set to be equally spaced. All cases were re-planed using Direct Aperture optimizer in Prowess Panther ver 5.01 IMRT planning system at same configurations and dose constraints. Plans were evaluated according to ICRU criteria, number of segments, number of monitor units and planning time. Results: For H/N plans, the near maximum dose (D2) and the dose that covers 95% D95 of PTV has improved by 4% in DAO. For organs at risk (OAR), DAO reduced the volume covered by 30% (V30) in spinal cord, right parotid, and left parotid by 60%, 54%, and 53% respectively. This considerable dosimetric quality improvement achieved using 25% less planning time and lower number of segments and monitor units by 46% and 51% respectively. In DAO prostate plans, Both D2 and D95 for the PTV were improved by only 2%. The V30 of the right femur, left femur and bladder were improved by 35%, 15% and 3% respectively. On the contrary, the rectum V30 got even worse by 9%. However, number of monitor units, and number of segments decreased by 20% and 25% respectively. Moreover the planning time reduced significantly too. Conclusion: DAO introduces considerable advantages over the beamlet optimization in regards to organs at risk sparing. However, no significant improvement occurred in most studied PTVs.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

SU-F-BRD-13: Quantum Annealing Applied to IMRT Beamlet Intensity Optimization

Purpose: We report on the first application of quantum annealing (QA) to the process of beamlet intensity optimization for IMRT. QA is a new technology, which employs novel hardware and software techniques to address various discrete optimization problems in many fields. Methods: We apply the D-Wave Inc. proprietary hardware, which natively exploits quantum mechanical effects for improved optimization. The new QA algorithm, running on this hardware, is most similar to simulated annealing, but relies on natural processes to directly minimize the free energy of a system. A simple quantum system is slowly evolved into a classical system, representing the objective function. To apply QA to IMRT-type optimization, two prostate cases were considered. A reduced number of beamlets were employed, due to the current QA hardware limitation of ∼500 binary variables. The beamlet dose matrices were computed using CERR, and an objective function was defined based on typical clinical constraints, including dose-volume objectives. The objective function was discretized, and the QA method was compared to two standard optimization Methods: simulated annealing and Tabu search, run on a conventional computing cluster. Results: Based on several runs, the average final objective function value achieved by the QA was 16.9 for the first patient, compared with 10.0 for Tabu and 6.7 for the SA. For the second patient, the values were 70.7 for the QA, 120.0 for Tabu, and 22.9 for the SA. The QA algorithm required 27–38% of the time required by the other two methods. Conclusion: In terms of objective function value, the QA performance was similar to Tabu but less effective than the SA. However, its speed was 3–4 times faster than the other two methods. This initial experiment suggests that QA-based heuristics may offer significant speedup over conventional clinical optimization methods, as quantum annealing hardware scales to larger sizes.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2014

Poster — Thur Eve — 61: A new framework for MPERT plan optimization using MC-DAO

This work combines the inverse planning technique known as Direct Aperture Optimization (DAO) with Intensity Modulated Radiation Therapy (IMRT) and combined electron and photon therapy plans. In particular, determining conditions under which Modulated Photon/Electron Radiation Therapy (MPERT) produces better dose conformality and sparing of organs at risk than traditional IMRT plans is central to the project. Presented here are the materials and methods used to generate and manipulate the DAO procedure. Included is the introduction of a powerful Java-based toolkit, the Aperture-based Monte Carlo (MC) MPERT Optimizer (AMMO), that serves as a framework for optimization and provides streamlined access to underlying particle transport packages. Comparison of the toolkit's dose calculations to those produced by the Eclipse TPS and the demonstration of a preliminary optimization are presented as first benchmarks. Excellent agreement is illustrated between the Eclipse TPS and AMMO for a 6MV photon field. The results of a simple optimization shows the functioning of the optimization framework, while significant research remains to characterize appropriate constraints.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Asymmetric dose–volume optimization with smoothness control for rotating-shield brachytherapy

Purpose: It is important to reduce fluence map complexity in rotating-shield brachytherapy (RSBT) inverse planning to improve delivery efficiency while maintaining plan quality. This study proposes an efficient and effective RSBT dose optimization method which enables to produce smooth fluence maps. Methods: Five cervical cancer patients each with a high-risk clinical-target-volume (HR-CTV) larger than 40 cm{sup 3} were considered as the test cases. The RSBT source was a partially shielded electronic brachytherapy source (Xoft Axxent™). The anchor RSBT plans generated by the asymmetric dose–volume optimization with smoothness control (ADOS) method were compared against those produced by the dose–surface optimization (DSO) method and inverse-planning with simulated annealing (IPSA). Either L{sub 1}-norm or L{sub 2}-norm was used to measure the smoothness of a fluence map in the proposed ADOS method as one weighted term of the objective function. Uniform dwell-time scaling was applied to all plans such that HR-CTV D{sub 90} was maximized without violating the D{sub 2cc} tolerances of the rectum, bladder, and sigmoid colon. The quality of the anchor plans was measured with HR-CTV D{sub 90} of the anchor plans. Single-shielded RSBT [(S-RSBT), RSBT with single, fix sized delivery window] and dynamic-sheilded RSBT [(D-RSBT), RSBT with dynamically varying sized delivery window] delivery plans generated based on the anchor plans were also measured, with delivery time constraints of 10, 20, and 30 min/fraction (fx). Results: The average HR-CTV D{sub 90} values of the anchor plans achieved by the ADOS, DSO, and IPSA methods were 111.5, 94.2, and 107.4 Gy, respectively, where the weighting parameter β used in ADOS with L{sub 2}-norm was set to be 100. By using S-RSBT sequencing and 20 min/fx delivery time, the corresponding D{sub 90} values were 88.8, 81.9, and 83.4 Gy; while using D-RSBT sequencing with 20 min/fx delivery time, the corresponding D{sub 90} values were 91.4, 88.3, and 78.9 Gy, respectively. The average optimization times for ADOS, DSO, and IPSA were, respectively, 77, 4, and 1800 s. By using L{sub 1}-norm instead of L{sub 2}-norm in the ADOS method, the optimization time was increased by 20 s, while the D{sub 90} was reduced by 6.8 Gy on average. ADOS-L1 was found to be more sensitive to the weighting parameter than ADOS-L2. If β was increased to 10 000, the D{sub 90} drop with ADOS-L1 was 38 Gy, while the drop with ADOS-L2 is 13 Gy. Conclusions: The ADOS method had a reasonable optimization time cost, while achieving comparable RSBT dose plans as the IPSA method, which is of much higher time complexity. Compared to the DSO and IPSA methods, ADOS is able to generate anchor plans which are more suitable for RSBT delivery while preserving the high quality of the original plans. Compared to ADOS-L1, ADOS-L2 is able to achieve better quality of anchor plans more efficiently.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-07: 4DCT Robust Optimization for Esophageal Cancer Using Intensity Modulated Proton Therapy

Purpose: To develop a 4DCT robust optimization method to reduce the dosimetric impact from respiratory motion in intensity modulated proton therapy (IMPT) for esophageal cancer. Methods: Four esophageal cancer patients were selected for this study. The different phases of CT from a set of 4DCT were incorporated into the worst-case dose distribution robust optimization algorithm. 4DCT robust treatment plans were designed and compared with the conventional non-robust plans. Result doses were calculated on the average and maximum inhale/exhale phases of 4DCT. Dose volume histogram (DVH) band graphic and ΔD95%, ΔD98%, ΔD5%, ΔD2% of CTV between different phases were used to evaluate the robustness of the plans. Results: Compare to the IMPT plans optimized using conventional methods, the 4DCT robust IMPT plans can achieve the same quality in nominal cases, while yield a better robustness to breathing motion. The mean ΔD95%, ΔD98%, ΔD5% and ΔD2% of CTV are 6%, 3.2%, 0.9% and 1% for the robustly optimized plans vs. 16.2%, 11.8%, 1.6% and 3.3% from the conventional non-robust plans. Conclusion: A 4DCT robust optimization method was proposed for esophageal cancer using IMPT. We demonstrate that the 4DCT robust optimization can mitigate the dose deviation caused by the diaphragm motion.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-574: Novel Chance-Constrained Optimization in Intensity-Modulated Proton Therapy Planning to Account for Range and Patient Setup Uncertainties

Purpose: We propose to apply a probabilistic framework, namely chanceconstrained optimization, in the intensity-modulated proton therapy (IMPT) planning subject to range and patient setup uncertainties. The purpose is to hedge against the influence of uncertainties and improve robustness of treatment plans. Methods: IMPT plans were generated for a typical prostate patient. Nine dose distributions are computed — the nominal one and one each for ±5mm setup uncertainties along three cardinal axes and for ±3.5% range uncertainty. These nine dose distributions are supplied to the solver CPLEX as chance constraints to explicitly control plan robustness under these representative uncertainty scenarios with certain probability. This probability is determined by the tolerance level. We make the chance-constrained model tractable by converting it to a mixed integer optimization problem. The quality of plans derived from this method is evaluated using dose-volume histogram (DVH) indices such as tumor dose homogeneity (D5% – D95%) and coverage (D95%) and normal tissue sparing like V70 of rectum, V65, and V40 of bladder. We also compare the results from this novel method with the conventional PTV-based method to further demonstrate its effectiveness Results: Our model can yield clinically acceptable plans within 50 seconds. The chance-constrained optimization produces IMPT plans with comparable target coverage, better target dose homogeneity, and better normal tissue sparing compared to the PTV-based optimization [D95% CTV: 67.9 vs 68.7 (Gy), D5% – D95% CTV: 11.9 vs 18 (Gy), V70 rectum: 0.0 % vs 0.33%, V65 bladder: 2.17% vs 9.33%, V40 bladder: 8.83% vs 21.83%]. It also simultaneously makes the plan more robust [Width of DVH band at D50%: 2.0 vs 10.0 (Gy)]. The tolerance level may be varied to control the tradeoff between plan robustness and quality. Conclusion: The chance-constrained optimization generates superior IMPT plan compared to the PTV-based optimization with explicit control of plan robustness. NIH/NCI K25CA168984, Eagles Cancer Research Career Development, The Lawrence W. and Marilyn W. Matteson Fund for Cancer Research, Mayo ASU Seed Grant, and The Kemper Marley Foundation.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

SU-E-T-590: Optimizing Magnetic Field Strengths with Matlab for An Ion-Optic System in Particle Therapy Consisting of Two Quadrupole Magnets for Subsequent Simulations with the Monte-Carlo Code FLUKA

Purpose: Aim of this study was to optimize the magnetic field strengths of two quadrupole magnets in a particle therapy facility in order to obtain a beam quality suitable for spot beam scanning. Methods: The particle transport through an ion-optic system of a particle therapy facility consisting of the beam tube, two quadrupole magnets and a beam monitor system was calculated with the help of Matlab by using matrices that solve the equation of motion of a charged particle in a magnetic field and field-free region, respectively. The magnetic field strengths were optimized in order to obtain a circular and thin beam spot at the iso-center of the therapy facility. These optimized field strengths were subsequently transferred to the Monte-Carlo code FLUKA and the transport of 80 MeV/u C12-ions through this ion-optic system was calculated by using a user-routine to implement magnetic fields. The fluence along the beam-axis and at the iso-center was evaluated. Results: The magnetic field strengths could be optimized by using Matlab and transferred to the Monte-Carlo code FLUKA. The implementation via a user-routine was successful. Analyzing the fluence-pattern along the beam-axis the characteristic focusing and de-focusing effects of the quadrupole magnets could be reproduced. Furthermore the beam spot at the iso-center was circular and significantly thinner compared to an unfocused beam. Conclusion: In this study a Matlab tool was developed to optimize magnetic field strengths for an ion-optic system consisting of two quadrupole magnets as part of a particle therapy facility. These magnetic field strengths could subsequently be transferred to and implemented in the Monte-Carlo code FLUKA to simulate the particle transport through this optimized ion-optic system.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

SU-E-T-780: Use Robustness Optimization (RO) Method to Improve the Planning Efficiency for Pencil Beam Scanning Cranial Spinal Irradiation

Purpose: We evaluate the feasibility of using robustness optimization (RO) function to improve the planning efficiency of pencil beam scanning (PBS) craniospinal irradiation (CSI) with gradient matching technique. Methods: A CSI patient was planned with 2 lateral brain fields and 4 posterior fields to cover the entire spine to maximal field of 24 cm × 20 cm on a compact PBS gantry, ProteusONE. CSI plans were generated using traditional volumetric gradient dose optimization (VGDO) and robustness optimization (RO) method respectively. In traditional VGDO, besides the sectioned spine target volumes, gradient volume (GV) were generated as 4 equally spaced structures e.g. 80%, 60%, 40%, and 20% of prescription dose. In RO method, only sectioned spine target volumes with an overlap of 4cm were created. In the robustness optimization settings, 5mm uncertainty in superior and inferior direction was defined for auto gradient optimization. Dosimetric metrics of conformity number (CN), homogeneity index (HI), and maximal cord doses were compared in Raystation version 4.6.100.6. Results: In VGDO method, total 16 GV structures and five 100% dose level target structures were contoured compared to total 5 target structures in RO method which saves 30 min in contour. With the same PTV coverage (95% volume receive 30.6Gy prescription dose), maximum cord dose is 32.64Gy in VGDO and 31.94Gy in RO. HI is 1.03 and 1.04 for VGDO and RO respectively. CN is 0.93 and 0.94 for VGDO and RO respectively. Conclusions: The dosimetric comparison demonstrated both methods are equivalent in terms of plan quality. With robust optimization for CSI gradient matching, it efficiently reduces the amount of planning target contour structure by factor of 4 and thus improves the planning efficiency especially for 4 or more gradient junction area.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

MO-FG-204-08: Optimization-Based Image Reconstruction From Unevenly Distributed Sparse Projection Views

Purpose: Optimization-based reconstruction has been proposed and investigated for reconstructing CT images from sparse views, as such the radiation dose can be substantially reduced while maintaining acceptable image quality. The investigation has so far focused on reconstruction from evenly distributed sparse views. Recognizing the clinical situations wherein only unevenly sparse views are available, e.g., image guided radiation therapy, CT perfusion and multi-cycle cardiovascular imaging, we investigate the performance of optimization-based image reconstruction from unevenly sparse projection views in this work. Methods: The investigation is carried out using the FORBILD and an anthropomorphic head phantoms. In the study, 82 views, which are evenly sorted out from a full (360°) axial CT scan consisting of 984 views, form sub-scan I. Another 82 views are sorted out in a similar manner to form sub-scan II. As such, a CT scan with sparse (164) views at 1:6 ratio are formed. By shifting the two sub-scans relatively in view angulation, a CT scan with unevenly distributed sparse (164) views at 1:6 ratio are formed. An optimization-based method is implemented to reconstruct images from the unevenly distributed views. By taking the FBP reconstruction from the full scan (984 views) as the reference, the root mean square (RMS) between the reference and the optimization-based reconstruction is used to evaluate the performance quantitatively. Results: In visual inspection, the optimization-based method outperforms the FBP substantially in the reconstruction from unevenly distributed, which are quantitatively verified by the RMS gauged globally and in ROIs in both the FORBILD and anthropomorphic head phantoms. The RMS increases with increasing severity in the uneven angular distribution, especially in the case of anthropomorphic head phantom. Conclusion: The optimization-based image reconstruction can save radiation dose up to 12-fold while providing acceptable image quality for advanced clinical applications wherein only unevenly distributed sparse views are available. Research Grants: W81XWH-12-1-0138 (DoD), Sinovision Technologies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

TU-CD-304-10: Development and Optimization of “Compton Lens” Collimator Design for Increased Dose Rate in SRS

Purpose: To provide a proof of concept for a new collimator design to increase the dose rate at isocenter for stereotactic radiosurgery (SRS) treatments by taking advantage of off axis Compton scattered photons which are attenuated in current collimators. Methods: A fundamentally new collimator design was developed and optimized by introducing a series of slits to a standard block collimator. The introduction of slits allowed for initially off axis radiation that was scattered in the direction of the target to contribute to the target dose. For optimization, the design was broken into two parts: an upper interaction plate where primary scattering occurs, and a lower “Compton slit” region which allows for scattered photons traveling toward the target to reach isocenter. To optimize the design, a series of simulations were performed using MCNP6 in which several key parameters were adjusted and the output was compared to a standard collimator. Key parameters modified included the collimator material, cone size, and interaction plate thickness. The effects of using energies different than the traditional 6 MV beam were also explored. Results: An optimized collimator design utilizing a solid interaction plate with a Cesium-137 beam and a 4 mm cone size resulted in a dose rate increase on the order of 5% relative to standard collimators in use. Currently, designs incorporating a Cesium-137 source are the most feasible due to necessary size and weight concerns for 6 MV beams. Conclusion: Preliminary designs provide a proof of concept and indicate a potential to improve upon the dose rate of current collimators while not largely compromising the sharp dose falloff inherent to SRS. Further optimization into the geometry and positioning of the interaction plate, as well as slit optimization, will likely lead to further dose rate increases than were observed in this study. Funding for this work was provided by the Wisconsin Alumni Research Foundation (WARF). Authors have a conflict of interest as the technology discussed is patent pending, and thus there is a financial interest in the work.

60 APPLIED LIFE SCIENCES↗