SEARCH · Engineering Papers
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.
This paper discusses the efforts to create, using two NASA SP2 supercomputers, a 'Metacenter' which includes the capability to transparently and dynamically distribute the SP2 workload across the geographically separated systems. Functional components of the Phase 1 Metacenter are identified, outstanding issues are discussed, and the plan for the second phase of the project is outlined.
Bhutan's electricity system is already low-carbon through the use of hydropower. However, generation from hydroelectric plants is reduced during the dry winter months and will be impacted by climate change and other extreme weather events. Bhutan is mitigating these risks by diversifying their energy mix with clean energy resources. Bhutan is exploring photovoltaic (PV) solar energy development to enhance its energy system's overall resilience. To ensure efficient grid planning and solar integration, Bhutan's power generator, Druk Green Power Corporation, and the transmission and distribution utility, Bhutan Power Corporation, are partnering with the South Asia Group for Energy (SAGE) to develop their solar capabilities, integrate solar energy into the grid, and strengthen Bhutan's energy infrastructure and efficiency.
DISCO is a Python-based NREL software tool for conducting scalable, repeatable distribution analyses. Although DISCO was originally developed to support photovoltaic impact analyses, it can also be used to understand the impact of other distributed energy resources and load changes on distribution systems. This presentation took place on April 9, 2024, and included NREL power grid researchers Sherin Ann Abraham and Shibani Ghosh.
Purpose: Intensity-modulated radiotherapy (IMRT) treatment for lung cancer is difficult due to the heterogeneous dose distribution and excessive dose to the organs at risk (OARs). We introduce a simple method based on the base dose function (BDF) in Eclipse treatment planning system to overcome the difficulties. Methods: Thirteen patients suffered from stage III non-small cell lung cancer (NSCLC) were enrolled in the study. Three kinds of approaches were applied to obtain clinically acceptable treatment plans: 1) conventionally optimizing method with hot and cold spots re-optimization (CO); 2) target-divided optimizing method (TDO) in which the optimization objective in the lung density of planning target volume (PTV) was set to 2 to 4 Gy higher than in the soft tissue density; 3) base dose function (BDF) in which the treatment plan was produced based on the original plan for re-optimization. CO, TDO and BDF methods were then compared in terms of conformity index (CI), homogeneity index (HI), OARs sparing and monitor units (MUs). Additionally, delta4, portal dosimetry and IMSure were used to measure the dose delivering accuracy. Results: The BDF technique provided more superior CI and HI than the other two methods. Moreover, the new method also reduced the lung, esophagus, heart and spinal cord dose. However, the BDF plans needed extra 15% and 10% MUs than the CO and TDO methods. Dose verification results demonstrated good and comparable γ pass rates among the three methods. Conclusion: The proposed BDF method greatly improves the dose homogeneity and OARs sparing in the IMRT treatment for lung cancer.
Purpose: Real-time adaptive planning and treatment has been infeasible due in part to its high computational complexity. There have been many recent efforts to utilize graphics processing units (GPUs) to accelerate the computational performance and dose accuracy in radiation therapy. Data structure and memory access patterns are the key GPU factors that determine the computational performance and accuracy. In this paper, the authors present a nonvoxel-based (NVB) approach to maximize computational and memory access efficiency and throughput on the GPU. Methods: The proposed algorithm employs a ray-tracing mechanism to restructure the 3D data sets computed from the CT anatomy into a nonvoxel-based framework. In a process that takes only a few milliseconds of computing time, the algorithm restructured the data sets by ray-tracing through precalculated CT volumes to realign the coordinate system along the convolution direction, as defined by zenithal and azimuthal angles. During the ray-tracing step, the data were resampled according to radial sampling and parallel ray-spacing parameters making the algorithm independent of the original CT resolution. The nonvoxel-based algorithm presented in this paper also demonstrated a trade-off in computational performance and dose accuracy for different coordinate system configurations. In order to find the best balance between the computed speedup and the accuracy, the authors employed an exhaustive parameter search on all sampling parameters that defined the coordinate system configuration: zenithal, azimuthal, and radial sampling of the convolution algorithm, as well as the parallel ray spacing during ray tracing. The angular sampling parameters were varied between 4 and 48 discrete angles, while both radial sampling and parallel ray spacing were varied from 0.5 to 10 mm. The gamma distribution analysis method (γ) was used to compare the dose distributions using 2% and 2 mm dose difference and distance-to-agreement criteria, respectively. Accuracy was investigated using three distinct phantoms with varied geometries and heterogeneities and on a series of 14 segmented lung CT data sets. Performance gains were calculated using three 256 mm cube homogenous water phantoms, with isotropic voxel dimensions of 1, 2, and 4 mm. Results: The nonvoxel-based GPU algorithm was independent of the data size and provided significant computational gains over the CPU algorithm for large CT data sizes. The parameter search analysis also showed that the ray combination of 8 zenithal and 8 azimuthal angles along with 1 mm radial sampling and 2 mm parallel ray spacing maintained dose accuracy with greater than 99% of voxels passing the γ test. Combining the acceleration obtained from GPU parallelization with the sampling optimization, the authors achieved a total performance improvement factor of >175 000 when compared to our voxel-based ground truth CPU benchmark and a factor of 20 compared with a voxel-based GPU dose convolution method. Conclusions: The nonvoxel-based convolution method yielded substantial performance improvements over a generic GPU implementation, while maintaining accuracy as compared to a CPU computed ground truth dose distribution. Such an algorithm can be a key contribution toward developing tools for adaptive radiation therapy systems.
Purpose: To construct expected treatment dose for adaptive inverse planning optimization, and evaluate it on head and neck (h and n) cancer adaptive treatment modification. Methods: Adaptive inverse planning engine was developed and integrated in our in-house adaptive treatment control system. The adaptive inverse planning engine includes an expected treatment dose constructed using the daily cone beam (CB) CT images in its objective and constrains. Feasibility of the adaptive inverse planning optimization was evaluated retrospectively using daily CBCT images obtained from the image guided IMRT treatment of 19 h and n cancer patients. Adaptive treatment modification strategies with respect to the time and the number of adaptive inverse planning optimization during the treatment course were evaluated using the cumulative treatment dose in organs of interest constructed using all daily CBCT images. Results: Expected treatment dose was constructed to include both the delivered dose, to date, and the estimated dose for the remaining treatment during the adaptive treatment course. It was used in treatment evaluation, as well as in constructing the objective and constraints for adaptive inverse planning optimization. The optimization engine is feasible to perform planning optimization based on preassigned treatment modification schedule. Compared to the conventional IMRT, the adaptive treatment for h and n cancer illustrated clear dose-volume improvement for all critical normal organs. The dose-volume reductions of right and left parotid glands, spine cord, brain stem and mandible were (17 {+-} 6)%, (14 {+-} 6)%, (11 {+-} 6)%, (12 {+-} 8)%, and (5 {+-} 3)% respectively with the single adaptive modification performed after the second treatment week; (24 {+-} 6)%, (22 {+-} 8)%, (21 {+-} 5)%, (19 {+-} 8)%, and (10 {+-} 6)% with three weekly modifications; and (28 {+-} 5)%, (25 {+-} 9)%, (26 {+-} 5)%, (24 {+-} 8)%, and (15 {+-} 9)% with five weekly modifications. Conclusions: Adaptive treatment modification can be implemented including the expected treatment dose in the adaptive inverse planning optimization. The retrospective evaluation results demonstrate that utilizing the weekly adaptive inverse planning optimization, the dose distribution of h and n cancer treatment can be largely improved.
Purpose: To investigate the trade-off between vertebral column sparing and thecal-sac target coverage in craniospinal irradiation (CSI) of pediatric patients treated with passive-scattering (PS) and intensity modulated (IMPT) proton therapy. Methods: We selected 2 pediatric patients treated with PS CSI for medulloblastoma. Spinal irradiation was re-planned with IMPT. For all cases, we assumed prescription dose of 23.4 Gy(RBE), with the spinal canal receiving at least 95% of 23.4 Gy(RBE). PS planning was performed using the commercial system XiO. IMPT planning was done using the Astroid planning system. Beam arrangements consisted of (a) PS posterior-anterior (PA) field, PS-PA, (b) IMPT PA field, IMPT-PA, and (c) two posterior oblique IMPT fields, IMPT2 (-35°, 35°). Dose distributions were re-calculated using TOPAS Monte Carlo, along with LET distributions, to investigate LET variations within the target and vertebra anatomy. Variable RBE-weighed dose distributions were also calculated based on a dose and LET-dependent biophysical model. Dosimetric data were compared among the plans for the target volume, spinal cord and adjacent critical organs (thecal-sac and cauda equina). Results: IMPT2 resulted in better sparing of the posterior vertebral column (entrance region posterior to thecal-sac), where planned dose was approximately 6–8Gy(RBE). For IMPT-PA and PS-PA the MC-calculated dose to the posterior vertebral column was, on average, 20Gy and 18Gy respectively. For IMPT2 higher mean-LET (5keV/µm/(g/cm3)) values were observed in anterior vertebral column (beyond the thecal-sac) relative to IMPT-PA and PS-PA, where mean-LET was 3.5keV/µm/(g/cm3) and 2.5keV/µm/(g/cm3) respectively. The higher LET region observed for both IMPT plans was in the distal end of treatment fields, where dose delivered was less 5Gy(RBE). Conclusion: The two-oblique proton beams IMPT2 best spared the spinal column, while reducing the dose to the posterior spinal column from 18–20 to 6–8 Gy(RBE). The best LET distribution was obtained with the PS-PA fields.
Purpose: To compare dose distributions in stereotactic radiation surgery of brain lesions using gamma Knife, VMAT, conformal arcs, and IMRT in order to provide an optimal treatment. Methods: Dose distributions from single shot of 4C model of Gamma Knife at the helmet collimation sizes of 4, 8, 14, and 18 mm in diameter were compared with full arcs with the square shapes of 4×4 (or 5×5), 8×8 (or 10×10), and spherical shapes of 16 or 20 mm in diameter using EDR3 films in the same gamma knife QA phantom. Plans for ten SRS cases with single and multiple lesions were created in gamma knife plans and Pinnacle plans. The external beam plans had enlarged field size by 2-mm and used single conformal full circle arc for solitary lesion and none coplanar arcs/beams for multiple lesions. Coverage, conformity index, dose to critical organs, and integral dose to the brain and nearby critical structures were compared on all plans. Structures and dose matrices were registered in a Velocity deformable image registration system. Results: Single full circle arc from Elekta beam-modulate MLC (4-mm leaf thickness) and agility MLC (5-mm leaf thickness) have larger penumbra and less flatness than that of Gamma Knife single shot. None-coplanar arcs or beams were required to achieve similar dose distribution. In general, Gamma Knife plans provided significant less integral dose than that of linac-based plans. Benefits of IMRT and VMAT versus gamma Knife and conformal arcs were not significant. Conclusion: Our dose measurement and treatment planning evaluation clearly demonstrated dose distribution differences amount current popular SRS modalities for small solitary and multiple brain lesions. The trend of using MLC shape beams or arcs to replace conventional cones should be revisited in order to keep lower integral dose if the late correlates with some radiation-induced side effects. Pilot grant from Elekta LLC.
Purpose: A treatment planning system (TPS) was commissioned for clinical use with a fixed beam line proton delivery system. An outline of the data collection, modeling, and verification is provided. Methods: Beam data modeling for proton spot scanning in CMS Xio TPS requires the following measurements: (i) integral depth dose curves (IDDCs); (ii) absolute dose calibration; and (iii) beam spot characteristics. The IDDCs for 18 proton energies were measured using an integrating detector in a single spot field in a water phantom. Absolute scaling of the IDDCs were performed based on ion chamber measurements in mono-energetic 10×10 cm{sup 2} fields in water. Beam spot shapes were measured in air using a flat panel scintillator detector at multiple planes. For beam model verification, more than 45 uniform dose phantom and patient plans were generated. These plans were used to measure range, point dose, and longitudinal and lateral profiles. Tolerances employed for verification are: point dose and longitudinal profiles, ±2%; range, ±1 mm; FWHM for lateral profiles, ±2 mm; and patient plan dose distribution, gamma index of >90% at 3%/3 mm criteria. Results: More than 97% of the point dose measurements out of 115 were within +/-2% with maximum deviation of 3%. 98% of the ranges measured were within 1 mm with maximum deviation of 1.4mm. The normalized depth doses were within 2% at all depths. The maximum error in FWHM of lateral profiles was found to be less than 2mm. For 5 patient plans representing different anatomic sites, a total of 38 planes for 12 beams were analyzed for gamma index with average value of 99% and minimum of 94%. Conclusions: The planning system is successfully commissioned and can be safely deployed for clinical use. Measurements of IDDCs on user beam are highly recommended instead of using standard beam IDDCs.
Using technologies such as solar photovoltaics (PV) to provide power to a site during a large-scale grid outage can be an effective solution for enhancing resilience, if designed appropriately. This paper discusses storm-hardening strategies including wind-loading standards, secure racking systems, balance of electrical systems and pre-storm preparations.
Purpose: The use of respiratory gating for management of breathing motion during stereotactic body radiation therapy (SBRT) relies on a consistent relationship between the breathing signal and the actual position of the internal target. This relationship was investigated in patients treated for pancreatic cancer. Methods: Four patients with pancreatic cancer undergoing SBRT that had implanted fiducials in the tumor were included in this study. Treatment plans were generated based on the exhale phases (30–70%) from the pre-treatment 4DCT. The margin between the internal target volume (ITV) and the planning target volume was three mm. After patient setup using cone-beam CT, simultaneous fluoroscopic imaging and breathing motion monitoring were used during at least three breathing cycles to verify the fiducial position and to optimize the gating window. After treatment, fluoroscopic images were acquired for verification purposes and exported for retrospective analyses. Fiducial positions were determined using a template-matching algorithm. For each dataset, we established a linear relationship between the fiducial position and the anterior-posterior (AP) breathing signal. The relationships before and after treatment were compared and the dose distribution impact evaluated. Results: Seven pre- and post-treatment fluoroscopic pairs were available for fiducial position analyses in the superior-inferior (SI) and left-right (LR) directions, and five in the AP direction. Time between image acquisitions was typically six to eight minutes. An average absolute change of 1.2±0.7 mm (range: 0.1–1.7) of the SI fiducial position relative to the external signal was found. Corresponding numbers for the LR and AP fiducial positions were 0.9±1.0 mm (range: 0.2–3.0) and 0.5±0.4 mm (range: 0.2–1.2), respectively. The dose distribution impact was small in both the ITV and organs-at-risk. Conclusion: The relationship change between fiducial position and external breathing signal has been observed to be about 1 mm in four pancreas SBRT patients, leading to small dose distribution impact. Pettersson and Cervino are funded by a Varian Medical Systems grant.
Purpose: To evaluate a patient-specific QA program and system for constancy checking of a scanning delivery system developed at the National Institute of Radiological Sciences.Methods: For the patient-specific QA, all the planned beams are recalculated on a water phantom with treatment planning software (TPS). The recalculated dose distributions are compared with the measured distributions using a 2D ionization chamber array at several depths, and evaluated using gamma index analysis with criteria of 3% and 3 mm and a pass rate of 90%. For the constancy check, the authors developed the multiwire proportional chamber (MWPC), which can record the delivered 2D fluence images in a slice-by-slice manner. During irradiation for dosimetric QA with the 2D ionization chamber array and an accordion-type water phantom, the 2D fluence images are recorded using the MWPC in the delivery system. These recorded images are then compared to those taken in the treatment session to check the constancy check. This analysis also employs gamma index analysis using the same criteria as in the patient-specific QA. These patient-specific QA and constancy check evaluations were performed using the data of 122 patients.Results: In the patient-specific QA, the measured dose distributions agreed well with those calculated by the TPS, and the QA criteria were satisfied in all measurements. The additional check of the fluence comparison ensured the constancy of the delivered field during each treatment irradiation.Conclusions: The authors established a patient-specific QA program and additional check of delivery constancy in every treatment session. Fluence comparison is a strong tool for constancy checking of the delivery system.
The Glory program is an Earth and Solar science mission designed to broaden science community knowledge of the environment. The causes and effects of global warming have become a concern in recent years and Glory aims to contribute to the knowledge base of the science community. Glory is designed for two functions: one is solar viewing to monitor the total solar irradiance and the other is observing the Earth s atmosphere for aerosol composition. The former is done with an active cavity radiometer, while the latter is accomplished with an aerosol polarimeter sensor to discern atmospheric particles. The Glory program is managed by NASA Goddard Space Flight Center (GSFC) with Orbital Sciences in Dulles, VA as the prime contractor for the spacecraft bus, mission operations, and ground system. This paper will describe some of the more unique features of the Glory program including the integration and testing of the satellite and instruments as well as the science data processing. The spacecraft integration and test approach requires extensive analysis and additional planning to ensure existing components are successfully functioning with the new Glory components. The science mission data analysis requires development of mission unique processing systems and algorithms. Science data analysis and distribution will utilize our national assets at the Goddard Institute for Space Studies (GISS) and the University of Colorado's Laboratory for Atmospheric and Space Physics (LASP). The Satellite was originally designed and built for the Vegetation Canopy Lidar (VCL) mission, which was terminated in the middle of integration and testing due to payload development issues. The bus was then placed in secure storage in 2001 and removed from an environmentally controlled container in late 2003 to be refurbished to meet the Glory program requirements. Functional testing of all the components was done as a system at the start of the program, very different from a traditional program. The plan for Glory is to minimize any changes to the spacecraft in order to meet the Glory requirements. This means that the instrument designs must adhere to the existing interfaces and capabilities as much as possible. Given Glory's unique history and the potential science return, the program is one of significant value to both the science community and the world. The findings Glory promises will improve our understanding of the drivers for global climate change for a minimal investment. The program hopes to show that reuse of existing government assets can result in a lower cost, and fully successful mission.
Purpose: Treatment delivery monitoring is important for radiotherapy, which enables catching dosimetric error at the earliest possible opportunity. This project develops a virtual delivery system to monitor the dose delivery process of photon radiotherapy in real-time using GPU-based Monte Carlo (MC) method. Methods: The simulation process consists of 3 parallel CPU threads. A thread T1 is responsible for communication with a linac, which acquires a set of linac status parameters, e.g. gantry angles, MLC configurations, and beam MUs every 20 ms. Since linac vendors currently do not offer interface to acquire data in real time, we mimic this process by fetching information from a linac dynalog file at the set frequency. Instantaneous beam fluence map (FM) is calculated. A FM buffer is also created in T1 and the instantaneous FM is accumulated to it. This process continues, until a ready signal is received from thread T2 on which an inhouse developed MC dose engine executes on GPU. At that moment, the accumulated FM is transferred to T2 for dose calculations, and the FM buffer in T1 is cleared. Once the calculation finishes, the resulting 3D dose distribution is directed to thread T3, which displays it in three orthogonal planes overlaid on the CT image for treatment monitoring. This process continues to monitor the 3D dose distribution in real-time. Results: An IMRT and a VMAT cases used in our patient-specific QA are studied. Maximum dose differences between our system and treatment planning system are 0.98% and 1.58% for the two cases, respectively. The average time per MC calculation is 0.1sec with <2% relative uncertainty. The update frequency of ∼10Hz is considered as real time. Conclusion: By embedding a GPU-based MC code in a novel data/work flow, it is possible to achieve real-time MC dose calculations to monitor delivery process.
Electron spin resonance (ESR) and magnetic resonance imaging (MRI) can be used to measure radiation dose deposited in different milieu through its effects. Radiation can break chemical bonds and if they produce stable free radicals, ESR can measure their concentration through their spins and a dose can be inferred. Ionizing radiation can also promote polymerization and in this case proton relaxation times can be measured and an image weighed by T2 can be produced giving spatial information about dose. A review of the basics of these applications is presented concluding with an end-to-end test using a composite Gel-Alanine phantom to validate 3-dimensionally dose distribution delivered in a simulation of Volume Modulated Arch Therapy on the simultaneous treatment of multiple brain metastases. The results obtained with the gel and alanine dosimeters are consistent with the expected by the treatment planning system, showing the potential of this multidosimetric approach and validating dosimetrically the multiple brain metastases treatment using VMAT.
Purpose: To develop the detector for the four-dimensional dose distribution measurement. Methods: We made the prototype detector for four-dimensional dose distribution measurement using a cylindrical plastic scintillator (5 cm diameter) and a conical reflection grass. The plastic scintillator is used as a phantom. When the plastic scintillator is irradiated, the scintillation light was emitted according to absorbed dose distribution. The conical reflection grass was arranged to surround the plastic scintillator, which project to downstream the projection images of the scintillation light. Then, the projection image was reflected to 45 degree direction by flat reflection grass, and was recorded by camcorder. By reconstructing the three-dimensional dose distribution from the projection image recorded in each frame, we could obtain the four-dimensional dose distribution. First, we tested the characteristic according to the amount of emitted light. Then we compared of the light profile and the dose profile calculated with the radiotherapy treatment planning system. Results: The dose dependency of the amount of light showed linearity. The pixel detecting smaller amount of light had high sensitivity than the pixel detecting larger amount of light. However the difference of the sensitivity could be corrected from the amount of light detected in each pixel. Both of the depth light profile through the conical reflection grass and the depth dose profile showed the same attenuation in the region deeper than peak depth. In lateral direction, the difference of the both profiles was shown at outside field and penumbra region. We consider that the difference is occurred due to the scatter of the scintillation light in the plastic scintillator block. Conclusion: It was possible to obtain the amount of light corresponding to the absorbed dose distribution from the prototype detector. Four-dimensional dose distributions can be reconstructed with high accuracy by the correction of the scattered light.
This presentation provides a comprehensive overview of end-of-performance and decommissioning considerations for large-scale Battery Energy Storage Systems (BESS). It outlines expected system lifespans, midterm assessment needs, and pathways for extending operational life through augmentation or repowering. The presentation details regulatory requirements that govern decommissioning plans, cost estimates, financial assurance mechanisms, and performance obligations across multiple jurisdictions. It further examines end-of-life equipment management, including recycling, waste handling, transportation, and environmental compliance. Designed to support Malawi's electricity-sector institutions, the presentation highlights how planning for decommissioning and environmental stewardship can be integrated early in project development to ensure safe, financially accountable, and environmentally responsible BESS system retirement.
Approaches to automating the radiofrequency system used in the deep space network (DSN) are discussed. The goal is to operate DSN, which communicates with U.S. space probes that travel beyond the moon, in an unattended mode in order to reduce operating costs. The philosophy of a hierachical automation plan, which emphasizes distributed control and the utilization of DSN hardware development engineers to develop automation software, is discussed, and software development as well as common software and common microcontroller hardware are considered. It is planned that a station controller will serve as the operator interface and control and coordinate the operation of three system controllers (RF system, digital system, and antenna system). These systems are described.