Search NASA⌕ Search

SEARCH · Search NASA

Results for “matched field processing”

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

SU-F-BRF-09: A Non-Rigid Point Matching Method for Accurate Bladder Dose Summation in Cervical Cancer HDR Brachytherapy

Purpose: To propose and validate a deformable point matching scheme for surface deformation to facilitate accurate bladder dose summation for fractionated HDR cervical cancer treatment. Method: A deformable point matching scheme based on the thin plate spline robust point matching (TPSRPM) algorithm is proposed for bladder surface registration. The surface of bladders segmented from fractional CT images is extracted and discretized with triangular surface mesh. Deformation between the two bladder surfaces are obtained by matching the two meshes' vertices via the TPS-RPM algorithm, and the deformation vector fields (DVFs) characteristic of this deformation is estimated by B-spline approximation. Numerically, the algorithm is quantitatively compared with the Demons algorithm using five clinical cervical cancer cases by several metrics: vertex-to-vertex distance (VVD), Hausdorff distance (HD), percent error (PE), and conformity index (CI). Experimentally, the algorithm is validated on a balloon phantom with 12 surface fiducial markers. The balloon is inflated with different amount of water, and the displacement of fiducial markers is benchmarked as ground truth to study TPS-RPM calculated DVFs' accuracy. Results: In numerical evaluation, the mean VVD is 3.7(±2.0) mm after Demons, and 1.3(±0.9) mm after TPS-RPM. The mean HD is 14.4 mm after Demons, and 5.3mm after TPS-RPM. The mean PE is 101.7% after Demons and decreases to 18.7% after TPS-RPM. The mean CI is 0.63 after Demons, and increases to 0.90 after TPS-RPM. In the phantom study, the mean Euclidean distance of the fiducials is 7.4±3.0mm and 4.2±1.8mm after Demons and TPS-RPM, respectively. Conclusions: The bladder wall deformation is more accurate using the feature-based TPS-RPM algorithm than the intensity-based Demons algorithm, indicating that TPS-RPM has the potential for accurate bladder dose deformation and dose summation for multi-fractional cervical HDR brachytherapy. This work is supported in part by the National Natural ScienceFoundation of China (no 30970866 and no 81301940)

60 APPLIED LIFE SCIENCES↗

Revisiting single inclusive jet production: timelike factorization and reciprocity

Factorization theorems for single inclusive jet production play a crucial role in the study of jets and their substructure. In the case of small radius jets, the dynamics of the jet clustering can be factorized from both the hard production dynamics, and the dynamics of the low scale jet substructure measurement, and is described by a matching coefficient that can be computed in perturbative Quantum Chromodynamics (QCD). A proposed factorization formula describing this process has been previously presented in the literature, and is referred to as the semi-inclusive, or fragmenting jets formalism. By performing an explicit two-loop calculation, we show the inconsistency of this factorization formula, in agreement with another recent result in the literature. Building on recent progress in the factorization of single logarithmic observables, and the understanding of reciprocity, we then derive a new all-order factorization theorem for inclusive jet production. The use of a jet algorithm, being only a modification of the infrared structure of the measurement, modifies the structure of convolutions in the factorization theorem, as compared to inclusive fragmentation, but maintains the universality of the inclusive hard function and its associated Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution, which are ultraviolet properties. However, the non-trivial structure of convolutions in the factorization theorem implies that the jet functions exhibit a modified evolution. We perform an explicit two-loop calculation of the jet function in both N = 4 super Yang-Mills (SYM), and for all color channels in QCD, finding exact agreement with the structure derived from our renormalization group equations. In addition, we derive several new results, including an extension of our factorization formula to jet substructure observables, a jet algorithm definition of a generating function for the energy correlators, and new results for exclusive jet functions. Our results are a key ingredient for achieving precision jet substructure at colliders.

Effective Field Theories↗

Adaptation of existing facilities to isentropic compression experiments

We demonstrate that the established pulsed power infrastructure at the National High Magnetic Field Laboratory - Pulsed Field Facility (NHMFL-PFF) at the Los Alamos National Laboratory can be adapted to obtain high quality isentropic compression experiment (ICE) data on materials in extreme conditions of dynamic high pressure. Experiments utilized a single-turn magnet pulsed power system at the NHMFL-PFF that was originally designed to measure actinide samples in extremes of high magnetic field (to 300 Tesla). A simple modification to the single-turn magnet has converted it to a fast turnaround dynamic high pressure measurement system. This paper details the work done including important background details that indicate that much more can be accomplished with optimization of the load characteristics in terms of ultimate peak pressures. To match the rise time of the NHMFL capacitor bank ({approx}2 {mu}s versus {approx}0.5 {mu}s for the Sandia Z-machine) the sample dimensions can be relatively large, i.e., up to 5 mm thickness. The maximum stresses are {approx}50GPa (0.5 Mbar) at the maximum bank voltage (60 kV) and higher pressures may be possible if the sample is tamped. For the design and predictions of performance of the NHMFL-ICE experiment it is important to have good predictive models. A SPICE code simulation was chosen to model all aspects of the experiment, electrical and physical. To this end, accurate dynamic load models were developed to simulate the compression and expansion of the dynamic load at high pressures using shock physics principles. A series experiments have been performed which demonstrated the feasibility of the NHMFL-ICE technique. The results will be shown and discussed. The NHMFL-ICE technique is an excellent method for measuring equations of state (EOS) at megabar pressures. Because a complete EOS can be obtained in one experiment from zero to the peak pressure, and because many shots can be fired in one day, the technique promises to provide high quality EOS data at relatively low cost.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Enabling high-temperature nanophotonics for energy applications

The nascent field of high-temperature nanophotonics could potentially enable many important solid-state energy conversion applications, such as thermophotovoltaic energy generation, selective solar absorption, and selective emission of light. However, special challenges arise when trying to design nanophotonic materials with precisely tailored optical properties that can operate at high-temperatures (> 1,100 K). These include proper material selection and purity to prevent melting, evaporation, or chemical reactions; severe minimization of any material interfaces to prevent thermomechanical problems such as delamination; robust performance in the presence of surface diffusion; and long-range geometric precision over large areas with severe minimization of very small feature sizes to maintain structural stability. We report an approach for high-temperature nanophotonics that surmounts all of these difficulties. It consists of an analytical and computationally guided design involving high-purity tungsten in a precisely fabricated photonic crystal slab geometry (specifically chosen to eliminate interfaces arising from layer-by-layer fabrication) optimized for high performance and robustness in the presence of roughness, fabrication errors, and surface diffusion. It offers near-ultimate short-wavelength emittance and low, ultra-broadband long-wavelength emittance, along with a sharp cutoff offering 4:1 emittance contrast over 10% wavelength separation. This is achieved via Q-matching, whereby the absorptive and radiative rates of the photonic crystal’s cavity resonances are matched. Strong angular emission selectivity is also observed, with short-wavelength emission suppressed by 50% at 75° compared to normal incidence. Finally, a precise high-temperature measurement technique is developed to confirm that emission at 1,225 K can be primarily confined to wavelengths shorter than the cutoff wavelength.

14 SOLAR ENERGY↗

Gulf of Mexico miocene CO₂ site characterization mega transect

This project characterized the Miocene-age sub-seafloor stratigraphy in the near-offshore portion of the Gulf of Mexico adjacent to the Texas coast. The large number of industrial sources of carbon dioxide (CO₂) in coastal counties and the high density of onshore urbanization and environmentally sensitive areas make this offshore region extremely attractive for long-term storage of carbon dioxide emissions from industrial sources (CCS). The study leverages dense existing geologic data from decades of hydrocarbon exploration in and around the study area to characterize the regional geology for suitability and storage capacity. Primary products of the study include: regional static storage capacity estimates, sequestration “leads” and prospects with associated dynamic capacity estimates, experimental studies of CO₂-brine-rock interaction, best practices for site characterization, a large-format ‘Atlas’ of sequestration for the study area, and characterization of potential fluid migration pathways for reducing storage risks utilizing novel high-resolution 3D (HR3D) seismic surveys. In addition, three subcontracted studies address source-to-sink matching optimization, offshore well bore management and environmental aspects. The various geologic data and interpretations are integrated and summarized in a series of cross-sections and maps, which represent a primary resource for any near-term commercial deployment of CCS in the area. The regional study characterized and mapped important geologic features (e.g., Clemente-Tomas fault zone, the regionally extensive Marginulina A and Amphistegina B confining systems, etc.) that provided an important context for regional static capacity estimates and specific sequestration prospects of the study. A static capacity estimate of the majority of the Study area (14,467 mi 2 ) was estimated at 86 metric Gigatonnes. While local capacity estimates are likely to be lower due to reservoir-scale characteristics, the offshore Miocene interval is a storage resource of National interest for providing CO₂ storage as an atmospheric emissions abatement strategy. The natural petroleum system was used as an analog to infer seal quality and predict possible migration pathways of fluids in an engineered system of anthropogenic CO₂ injection and storage. The regional structural features (e.g., Clemente-Tomas fault zone) that exert primary control on the trapping and distribution of Miocene hydrocarbons are expected to perform similarly for CCS. Industrial-scale CCS will require storage capacity utilizing well-documented Miocene hydrocarbon (dominantly depleted gas) fields and their larger structural closures, as well as barren (unproductive, brine-filled) closures. No assessment was made of potential for CO₂ utilization for enhanced oil and gas recovery. The use of 3D numerical fluid flow simulations have been used in the study to greatly assist in characterizing the potential storage capacity of a specific reservoir. Due to the complexity of geologic systems (stratigraphic heterogeneity) and inherent limitations on producing a 3D geologic model, these simulations are typically simplified scenarios that explore the influence of model property variability (sensitivity study). A specific site offshore San Luis Pass (southern Galveston Island) was undertaken successfully, indicating stacked storage potential. Downscaling regional capacity estimates to the local scale (and the inverse) has proven challenging, and remains an outstanding gap in capacity assessments. In order to characterize regional seal performance and identify potential brine and CO₂ leakage pathways, results from three high-resolution 3D (HR3D) seismic datasets acquired by the study using novel HR3D (P-Cable) acquisition system showed steady and significant improvements in data quality because of improved acquisition and processing technique. Finely detailed faults and stratigraphy in the shallowest 1000 milliseconds (~800 m) of data allowed for the identification and mapping of unconformable surfaces including what is probably a surface associated with the last Pleistocene glacial lowstand. The identification of a previously unrecognized (in commercial seismic data) gas chimney that was clearly defined in the 2013 HR3D survey, indicates that HR3D surveys may be useful as both a characterization tool for the overburden of a potential carbon sequestration site and as an additional monitoring tool for future engineered injection sites. Geochemical modeling indicated that injection of CO₂ would result in minor dissolution of calcite, K-feldspar and albite. In addition, modeling of typical brines in Miocene age rocks indicate that approximately 5% of injection capacity would result from CO₂ dissolution into the brine. After extensive searches, no rock samples of the Marginulina A and Amphistegina B seals (“caprocks”) were obtained, but analyses of available core samples of other Miocene age mudrocks (seals or caprocks) indicate that they have sealing ability sufficient for potential CO 2 storage in underlying sandstone units.

01 COAL, LIGNITE, AND PEAT↗

Measuring linac photon beam energy through EPID image analysis of physically wedged fields

Purpose: Electronic portal imaging devices (EPIDs) have proven to be useful tools for measuring several parameters of interest in linac quality assurance (QA). However, a method for measuring linac photon beam energy using EPIDs has not previously been reported. In this report, such a method is devised and tested, based on fitting a second order polynomial to the profiles of physically wedged beams, where the metric of interest is the second order coefficientα. The relationship between α and the beam quality index [percentage depth dose at 10 cm depth (PDD{sub 10})] is examined to produce a suitable calibration curve between these two parameters. Methods: Measurements were taken in a water-tank for beams with a range of energies representative of the local QA tolerances about the nominal value 6 MV. In each case, the beam quality was found in terms of PDD{sub 10} for 100 × 100 mm{sup 2} square fields. EPID images of 200 × 200 mm{sup 2} wedged fields were then taken for each beam and the wedge profile was fitted in MATLAB 2010b (The MathWorks, Inc., Natick, MA). α was then plotted against PDD{sub 10} and fitted with a linear relation to produce the calibration curve. The uncertainty in α was evaluated by taking five repeat EPID images of the wedged field for a beam of 6 MV nominal energy. The consistency of measuring α was found by taking repeat measurements on a single linac over a three month period. The method was also tested at 10 MV by repeating the water-tank crosscalibration for a range of energies centered approximately about a 10 MV nominal value. Finally, the calibration curve from the test linac and that from a separate clinical machine were compared to test consistency of the method across machines in a matched fleet. Results: The relationship betweenα and PDD{sub 10} was found to be strongly linear (R{sup 2} = 0.979) while the uncertainty in α was found to be negligible compared to that associated with measuring PDD{sub 10} in the water-tank (±0.3%). The repeat measurements over a three month period showed the method to be reasonably consistent (i.e., well within the limits defined by local QA tolerances). The measurements were repeated on a matched machine and the same linear relationship between α and PDD{sub 10} was observed. The results for both machines were found to be indistinguishable across the energy range of interest (i.e., across and close to the thresholds defined by local QA tolerances), hence a single relation could be established across a matched fleet. Finally, the experiment was repeated on both linacs at 10 MV, where the linear relationship between α and PDD{sub 10} was again observed. Conclusions: The authors conclude that EPID image analysis of physically wedged beam profiles can be used to measure linac photon beam energy. The uncertainty in such a measurement is dominated by that associated with measuring PDD{sub 10} in the water-tank; hence, the accuracies of these two methods are directly comparable. This method provides a useful technique for quickly performing energy constancy measurements while saving significant clinical downtime for QA.

61 RADIATION PROTECTION AND DOSIMETRY↗

SU-E-J-115: Correlation of Displacement Vector Fields Calculated by Deformable Image Registration Algorithms with Motion Parameters of CT Images with Well-Defined Targets and Controlled-Motion

Purpose: To investigate correlation of displacement vector fields (DVF) calculated by deformable image registration algorithms with motion parameters in helical axial and cone-beam CT images with motion artifacts. Methods: A mobile thorax phantom with well-known targets with different sizes that were made from water-equivalent material and inserted in foam to simulate lung lesions. The thorax phantom was imaged with helical, axial and cone-beam CT. The phantom was moved with a cyclic motion with different motion amplitudes and frequencies along the superior-inferior direction. Different deformable image registration algorithms including demons, fast demons, Horn-Shunck and iterative-optical-flow from the DIRART software were used to deform CT images for the phantom with different motion patterns. The CT images of the mobile phantom were deformed to CT images of the stationary phantom. Results: The values of displacement vectors calculated by deformable image registration algorithm correlated strongly with motion amplitude where large displacement vectors were calculated for CT images with large motion amplitudes. For example, the maximal displacement vectors were nearly equal to the motion amplitudes (5mm, 10mm or 20mm) at interfaces between the mobile targets lung tissue, while the minimal displacement vectors were nearly equal to negative the motion amplitudes. The maximal and minimal displacement vectors matched with edges of the blurred targets along the Z-axis (motion-direction), while DVF’s were small in the other directions. This indicates that the blurred edges by phantom motion were shifted largely to match with the actual target edge. These shifts were nearly equal to the motion amplitude. Conclusions: The DVF from deformable-image registration algorithms correlated well with motion amplitude of well-defined mobile targets. This can be used to extract motion parameters such as amplitude. However, as motion amplitudes increased, image artifacts increased significantly and that limited image quality and poor correlation between the motion amplitude and DVF was obtained.

60 APPLIED LIFE SCIENCES↗

22” ADP Fan Rig Liners Design Report

A liner design study was conducted as part of a cooperative effort between five (5) government/industry teams that together seek to demonstrate the technology for designing and manufacturing acoustic liners that are twenty five percent (25%) more efficient than 1992 technology liners. The study emphasized teaming collaboration. The improved liners were designed by Pratt & Whitney and Boeing Airplane Company and built by Rohr Inc, and will be tested in the NASA/P& W 22-inch ADP fan rig at NASA Lewis Research Center's 9' x 15' wind tunnel. Design guidelines and decisions were made collectively during monthly design review telecons and at the formal final design review. The tools that were used to design the new improved liners were not new, but the process that was developed as part of this study that led into the evaluation and selection of the final and best liner designs is new. Until now, the liner design process as practiced by different industry teams varied greatly. Some procedures were based on empirical liner attenuation databases which could not adequately account for engine-to-engine hardwall fan noise spectral differences, while others were entirely theoretical. And regardless of which method one chose, the major difficulty was still the lack of understanding and knowledge of the actual hardwall fan source noise modal structure and farfield SPL spectra. The NASA-led effort to attempt actual measurements of the fan source noise modal structure by means of the rotating microphone array is expected to contribute significantly to the understanding of the nature of the fan tone noise modes, but the application to broadband noise is still a long way off. The new design process included consideration of the measured fan tone modes, but for the majority of the spectra a separate, systematic process was used. It is a common practice in the engine/nacelle industry that acoustic liners for new products are designed before actual measured hardwall engine far-field noise spectra are available. Target noise spectra are normally derived from existing engine noise databases with some adjustments to absolute SPLs and frequencies. This practice has been acceptable as long as the new engine was a derivative of the base engine. However, in the case of the ADP, the transition from a current engine base is too - great, and the adjustments could not adequately account for the quantum changes in SPL spectral differences. Examples were the 1992 single degree of freedom (SDOF) inlet and aft liners (designed by P&W) that would be used as the baseline liners against which the new improved liners' efficien¬cies would -be measured. As will be shown, these baseline liners have been found to be deeper than the desire optimum depths. The new liner design process begins with the selection of the target hard wall fan noise spectra. These spectra were obtained by scaling up (5.91 scale factor) the measured hardwall fan spectral data from the 22” ADP fan rig. Next, the modal energy contents for each 1/3-octave band center frequency of these hardwall fan noise spectra are estimated. (Within each 1/3-octave band center frequency, the model energy distribution approximation is for both tone and broadband noise). For the inlet noise, P&W uses a derivative of the NASA Lewis (Ed Rice's) method of classifying; and grouping propagating modes by their cutoff ratios. The next step is to assign energy level to each group of modes having the same cutoff ratio. In Ed Rice's model, the modes were grouped according to ten (10) cutoff ratio intervals with center values located at 1.026, 1.085, 1.155, 1.24, 1.35, 1.49, 1.69, 2.0, and 4.47 (with equal number of modes in each interval). The modes were assumed to have equal energy. P&W’s model expands the cutoff ratio-mode grouping into two hundred (200) smaller cutoff ratio intervals with center values located f1".'m 1.003 to 11.5 in 199 increasing incremental intervals. Next, P&W's model uses a "2-parameter" normal distribution as a template to assign energy levels to these 200 pre-determined cutoff ratio values (for each 1/3-octave band center frequency). In the past, P&W had conducted an extensive study to determine what "2-parameter" values are appropriate for fullscale inlet liners, and had developed a set of twenty-four (24) "2-pararneter" values (i.e. one for each 1/3-octave frequency band) that when used in Ed Rice's Inlet Attenuation Prediction Method produced predicted liner attenuations that closely matched measured liner attenuations from several P&W's engines. In the absence of actual measured tone and broadband modal data from the 22" ADP rig, the process will use P&W's proprietary set of "2-parameter" values for this liner design study. For the aft noise, Boeing uses a modal energy approximation that the "transport energy" of each propagating mode is equal. This approximation is almost the same as the equal energy per mode approximation, except for modes that are near cutoff. Boeing's model forces these modes to have lower energy levels. Both P&W and Boeing agreed that the "transport energy" approximation should work well for the ADP aft fan noise which appears to be dominated by broadband noise. The design process then proceeds to calculate the optimum liner impedances for each frequency in both the inlet and the aft. These optimum impedances represent the target impedances that the designed liners should have. Liners with impedances matching the optimum impedances at all frequencies are "ideal" liners. These ideal liners are theoretically the best liners. Unfortunately, it has been showed that it is impossible to design and build such ideal liners. The next best liners are ones that have impedances matching the optimum impedances over some frequency range (not all frequencies as for the ideal). This is accomplished by the use of "frequency weightings". Several of P&W's and Boeing's existing computer decks were used for optimizing and matching the designed liner impedances to the target optimum values (with the various frequency weightings specified). The optimization produces liner candidates with predicted liner impedances and descriptions of their physical liner characteristics. These candidate liner impedances are then used to predict their spectral attenuation characteristics which are then used together with the target hardwall fan noise spectra to determine the resulting treated noise spectra and PNLT values. Further optimization around the selected candidate designs yield final designs that are best in PNLT attenuations. Use of the optimization decks allowed a large number of liner candidates to be screened in a relatively short period of time. This design process is systematic and is efficient. The new inlet SDOF liner design obtained from the improved process was predicted to be 34% more effective (per unit area) than the 1992 baseline inlet liner. This inlet liner design is a 112 rayl wovenwiremesh facesheet over a 0.312-inch deep honeycomb core. The new aft SDOF liner design was predicted to be 52% more efficient than the aft baseline liner. The new aft SDOF liner is "segmented" with a shallower liner on the core cowl (inner duct wall) and a deeper liner on the fan cowl (outer duct wall). The shallower liner is a 70.6 rayl woven-wiremesh facesheet over a 0.141-inch deep honeycomb core. The deeper liner is a 68 rayl woven-wiremesh facesheet over a 0.309-inch deep honeycomb core. All liner dimensions are for the 22-inch model-scale ADP fan rig liners. The selected advanced liners are "segmented" double-layer (DDOF) liners for the inlet and aft locations. Also, for the inlet, a bulk liner with ceramic foam for wider broadband noise absorption was also selected. Triple-layer liner designs were not considered since the model scaled liners ( 1/5. 91) were dimensionally too small to be built correctly, and irrin earlier concept study, Boeing found a triple-layer to have only very small benefits over a double-layer. The inlet DDOF liner was predicted to be 83% more effective than the baseline. The inlet DDOF design is a 78 rayl facesheet over a 0.080 top cavity depth, a 68 rayl septum and a 0.227-inch bottom cavity depth. The inlet bulk liner is a 60 rayl facesheet over a 0.33-inch deep honeycomb filled with high temperature (HTP) ceramic foam with a density of 4.8 lb/cu.ft and a flow resistivity of 167 rayl/cm. The inlet bulk liner was predicted to be 83% more effective than the baseline liner. The aft DDOF liners are segmented with a shallower liner on the core cowl and a deeper liner on the fan cowl. The shallow DDOF liner is a 49.8 rayl facesheet over a 0.093-inch top cavity depth, a septum of 88.2 rayls over a 0.181-inch bottom cavity depth. The deep DDOF liner is a 12.9 rayl facesheet over a 0.140-inch top cavity depth, a septum of 53.1 ray ls over a 0.258-inch botom cavity depth. The segmented aft DDOF liners were predicted to be 86% more efficient than the baseline liner.

turbofan acoustic treatment↗

Kansas City, Missouri, Streetlight Electric Vehicle Charging: Strategies and challenges for site selection of streetlight electric vehicle infrastructure in Kansas City, Missouri (Final Report)

Public streetlight charging, whether on streets in central business districts or residential areas, provides easy charging access for apartment residents and homeowners alike. While most electric vehicle (EV) drivers charge at home, they do so in garages or on driveways they own. For renters and residents of multifamily housing (MFH), however, this may not be an option. EVs have a lower cost of ownership compared to conventional vehicles, and a used EV may be an affordable option for a lower-income household. But without easy access to charging, even a low-cost used EV may not be an option for a prospective buyer. An affordable curbside charging network has the potential to expand EV adoption into neighborhoods that have to date seen minimal interest and uptake of the technology and associated charging infrastructure. Streetlight charging networks can provide an economical, scalable, and effective approach to providing equitable and convenient charging. Metropolitan Energy Center (MEC) is dedicated to the mission of creating resource efficiency, environmental health, and economic vitality in the Kansas City region and beyond. Since 1983, MEC has provided resources, outreach, and training to make alternative fuels and energy efficiency commonplace. MEC led a streetlight charging pilot project that installed limited EV charging infrastructure on the streetlight system in Kansas City, Missouri, to demonstrate and test the benefits of curbside charging for EVs at existing on-street parking locations. The project aimed to cost-effectively expand the charging network in Kansas City to support residential charging and provide infrastructure in one or more charging deserts throughout the city. This pilot evaluates the impact and overall success of streetlight charging based on community feedback, utilization of charging infrastructure, technical feasibility, and cost. The project has pursued a data- and community-driven site selection process designed to identify sites with high demand and high opportunity for EV charging. This project was funded by the U.S. Department of Energy (DOE) and awarded to MEC through a competitive proposal process. The novelty and complexity of this project required an organization that could facilitate collaboration across levels of government, community members, and industry partners. For the past 25 years, through Kansas City Regional Clean Cities, MEC has worked with numerous public and private fleets on a variety of projects to improve the environmental performance and efficiency of the regional vehicle fleet. To advance affordable, efficient, and clean transportation efforts, DOE Clean Cities and Communities coalitions create local networks of public and private sector stakeholders and engage communities. Rooted within their local communities, the coalitions serve as experts and ambassadors, bringing to bear the collective knowledge, experience, and practical know-how of the entire network from within DOE, its national laboratories, and diverse stakeholders in the field. MEC and its project partners made in-kind contributions to leverage federal dollars for the benefit of the Kansas City community. Findings from this project will help determine the best applications for streetlight charging technologies to maximize funding impact and serve community needs. The team evaluated locations based on expected charging demand, technical feasibility, safety considerations, and enhanced charging network siting needs. Throughout the project, the team gathered feedback and evaluated ways to make public charging for EVs available to all community members. The insights will help Kansas City and other communities streamline future efforts to support EV drivers through public charging in the city right-of-way. Furthermore, this project will inform citywide guidance for future installations. MEC is committed to a transparent and publicly accessible approach that encourages the collaborative evaluation of streetlight charging. The project has engaged the community to proactively identify and evaluate the benefits and impacts of streetlight charging. It was a priority for the project to ensure the benefits of this pilot are distributed equitably to all members of the Kansas City community and that new charging opportunities and associated resources are available in diverse neighborhoods across the city. The charging infrastructure supports an affordable curbside charging network that will enable more drivers to choose EVs and provide easy charging access for all community members interested in driving an EV. The community feedback received through this project informed future resources and opportunities to make EVs more accessible to all members of the Kansas City community. MEC worked with several community partners on this project, including Missouri University of Science and Technology (MST), Pennsylvania State University (Penn State), the National Renewable Energy Laboratory (NREL); the city of Kansas City, Missouri; Evergy; Black and McDonald (B&M); LilyPad EV; EVNoire; and Westside Housing Organization (WHO). Project partners contributed to the cost match required for DOE grants through capital expenditures, personnel, and other in-kind contributions. Detailed descriptions of project team organizations can be found in Appendix A. Project Partners. Analysts at NREL and MST/PennState developed site maps based on demand and equity considerations. MEC conducted outreach to community members to garner input on project design and site selection, and received approval from the Missouri Public Service Commission (PSC) for Evergy’s EV charging station ownership. MEC worked with all partners to gather additional siting criteria and developed a site selection evaluation checklist, and partners conducted site visits to proposed installation sites. Next, B&M, Evergy, and the city executed all site agreements, conducted site-specific engineering design, acquired associated permits, and issued notices to proceed site by site or in small batches. Finally, from January to April 2023, the project team installed 23 EV charging stations built on Kansas City’s streetlight system in six council districts. Evergy will own, operate, and monitor the stations for 10 years, sharing charging data with MEC for at least 1 year.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

WE-D-9A-01: A Novel Mesh-Based Deformable Surface-Contour Registration

Purpose: Initial guess is vital for 3D-2D deformable image registration (DIR) while dealing with large deformations for adaptive radiation therapy. A fast procedure has been developed to deform body surface to match 2D body contour on projections. This surface-contour DIR will provide an initial deformation for further complete 3D DIR or image reconstruction. Methods: Both planning CT images and come-beam CT (CBCT) projections are preprocessed to create 0–1 binary mask. Then the body surface and CBCT projection body contours are extracted by Canny edge detector. A finite element modeling system was developed to automatically generate adaptive meshes based on the image surface. After that, the projections of the CT surface voxels are computed and compared with corresponding 2D projection contours from CBCT scans. As a result, the displacement vector field (DVF) on mesh vertices around the surface was optimized iteratively until the shortest Euclidean distance between the pixels on the projections of the deformed CT surface and the corresponding CBCT projection contour is minimized. With the help of the tetrahedral meshes, we can smoothly diffuse the deformation from the surface into the interior of the volume. Finally, the deformed CT images are obtained by the optimal DVF applied on the original planning CT images. Results: The accuracy of the surface-contour registration is evaluated by 3D normalized cross correlation increased from 0.9176 to 0.9957 (sphere-ellipsoid phantom) and from 0.7627 to 0.7919 (H and N cancer patient data). Under the GPU-based implementation, our surface-contour-guided method on H and N cancer patient data takes 8 seconds/iteration, about 7.5 times faster than direct 3D method (60 seconds/iteration), and it needs fewer optimization iterations (30 iterations vs 50 iterations). Conclusion: The proposed surface-contour DIR method can substantially improve both the accuracy and the speed of reconstructing volumetric images, which is helpful for applying in adaptive radiotherapy. This research is supported by CPRIT individual investigator award RP110329.

60 APPLIED LIFE SCIENCES↗

Joint EM-NE-International Study of Glass Behavior over Geologic Time Scales - 12303

Vitrification has been chosen as the best demonstrated available technology for waste immobilization worldwide. To date, the contributions of physical and chemical processes controlling the long-term glass dissolution rate in geologic disposal remain uncertain, leading to a lack of international consensus on a rate law for glass corrosion. Existing rate laws have overcome uncertainty through conservatism, but a thorough mechanistic understanding of waste form durability in geologic environments would improve public and regulator confidence. If it is possible to take credit for the true durability of the waste form in repository system evaluations, then it is possible to design the repository with less conservatism with concomitant cost savings. To gain a fundamental understanding of the dissolution rate law, six nations have joined together to formulate a joint plan for collaborative research into the mechanisms controlling the long-term corrosion of glass. This report highlights the technical program plan behind the US portion of this effort, with an emphasis on the current understanding (and limitations) of several mechanistic theories for glass corrosion. Some recent results are presented to provide an example of the ongoing research. Atom probe tomography has been used to provide a high-resolution analysis of elemental concentration gradients present at the hydrated glass / pristine glass interface in SON68 after 25.75 years of corrosion in a simulated granitic groundwater at 90 deg. C. The most valuable result of these initial studies is the success of the technique. Characterization by APT had never been previously demonstrated for glass corrosion layers. The resolution of APT is a powerful addition to the tools with which we can investigate the mechanisms dominating glass corrosion. Some other key results of this study include the observation that the elemental interfacial width between the hydrated glass and pristine glass appears to be much sharper (∼2 nm for B, Na and Al) than had been previously measured using nanoSIMS (∼240 nm). It is not clear whether the APT analysis and nanoSIMS characterizations were possibly performed on topographically unique regions, or whether nanoSIMS overestimated the elemental width. However, the APT data seems very convincing that the elemental width can be much sharper than was previously thought. This result calls into question some of the assumptions made for the diffusion-control models of glass dissolution, since such a sharp profile would not match the diffusion coefficients used to date. Other results, such as the observation of apparently layered concentration profiles, show that gel evolution is not as simple as is currently assumed in nearly every model. This task is a good example of the collaborative nature of the I-TEAM effort. Based on experimental needs and differences in expertise, scientists from DOE and CEA worked together to change the level of understanding in the field. These types of interactions are nearly ubiquitous among the tasks in the technical program plan. With the excellence of the team in place and the willingness of the participants to work together for a common understanding, the stated goal of consensus on the mechanistic basis for radionuclide release from glass is well within reach. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Distortions in charged-particle images of laser direct-drive inertial confinement fusion implosions

Energetic charged particles generated by inertial confinement fusion (ICF) implosions encode information about the spatial morphology of the hotspot and dense fuel during the time of peak fusion reactions. The knock-on deuteron imager (KoDI) was developed at the Omega Laser Facility to image these particles in order to diagnose low-mode asymmetries in the hotspot and dense fuel layer of cryogenic deuterium–tritium ICF implosions. However, the images collected are distorted in several ways that prevent reconstruction of the deuteron source. In this paper, we describe these distortions and a series of attempts to mitigate or compensate for them. We present several potential mechanisms for the distortions, including a new model for scattering of charged particles in filamentary electric or magnetic fields surrounding the implosion. Particle-tracing is used to create synthetic KoDI data based on the filamentary field model that reproduces the main experimentally observed image distortions. We conclude that the filamentary scattering model best matches the observed image distortions. Finally, we discuss potential impacts of filamentary fields on other charged-particle diagnostics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigating CT to CBCT image registration for head and neck proton therapy as a tool for daily dose recalculation

Purpose: Intensity modulated proton therapy (IMPT) of head and neck (H and N) cancer patients may be improved by plan adaptation. The decision to adapt the treatment plan based on a dose recalculation on the current anatomy requires a diagnostic quality computed tomography (CT) scan of the patient. As gantry-mounted cone beam CT (CBCT) scanners are currently being offered by vendors, they may offer daily or weekly updates of patient anatomy. CBCT image quality may not be sufficient for accurate proton dose calculation and it is likely necessary to perform CBCT CT number correction. In this work, the authors investigated deformable image registration (DIR) of the planning CT (pCT) to the CBCT to generate a virtual CT (vCT) to be used for proton dose recalculation. Methods: Datasets of six H and N cancer patients undergoing photon intensity modulated radiation therapy were used in this study to validate the vCT approach. Each dataset contained a CBCT acquired within 3 days of a replanning CT (rpCT), in addition to a pCT. The pCT and rpCT were delineated by a physician. A Morphons algorithm was employed in this work to perform DIR of the pCT to CBCT following a rigid registration of the two images. The contours from the pCT were deformed using the vector field resulting from DIR to yield a contoured vCT. The DIR accuracy was evaluated with a scale invariant feature transform (SIFT) algorithm comparing automatically identified matching features between vCT and CBCT. The rpCT was used as reference for evaluation of the vCT. The vCT and rpCT CT numbers were converted to stopping power ratio and the water equivalent thickness (WET) was calculated. IMPT dose distributions from treatment plans optimized on the pCT were recalculated with a Monte Carlo algorithm on the rpCT and vCT for comparison in terms of gamma index, dose volume histogram (DVH) statistics as well as proton range. The DIR generated contours on the vCT were compared to physician-drawn contours on the rpCT. Results: The DIR accuracy was better than 1.4 mm according to the SIFT evaluation. The mean WET differences between vCT (pCT) and rpCT were below 1 mm (2.6 mm). The amount of voxels passing 3%/3 mm gamma criteria were above 95% for the vCT vs rpCT. When using the rpCT contour set to derive DVH statistics from dose distributions calculated on the rpCT and vCT the differences, expressed in terms of 30 fractions of 2 Gy, were within [−4, 2 Gy] for parotid glands (D{sub mean}), spinal cord (D{sub 2%}), brainstem (D{sub 2%}), and CTV (D{sub 95%}). When using DIR generated contours for the vCT, those differences ranged within [−8, 11 Gy]. Conclusions: In this work, the authors generated CBCT based stopping power distributions using DIR of the pCT to a CBCT scan. DIR accuracy was below 1.4 mm as evaluated by the SIFT algorithm. Dose distributions calculated on the vCT agreed well to those calculated on the rpCT when using gamma index evaluation as well as DVH statistics based on the same contours. The use of DIR generated contours introduced variability in DVH statistics.

60 APPLIED LIFE SCIENCES↗

Development and clinical evaluation of automatic fiducial detection for tumor tracking in cine megavoltage images during volumetric modulated arc therapy

Purpose: Real-time tracking of implanted fiducials in cine megavoltage (MV) imaging during volumetric modulated arc therapy (VMAT) delivery is complicated due to the inherent low contrast of MV images and potential blockage of dynamic leaves configurations. The purpose of this work is to develop a clinically practical autodetection algorithm for motion management during VMAT. Methods: The expected field-specific segments and the planned fiducial position from the Eclipse (Varian Medical Systems, Palo Alto, CA) treatment planning system were projected onto the MV images. The fiducials were enhanced by applying a Laplacian of Gaussian filter in the spatial domain for each image, with a blob-shaped object as the impulse response. The search of implanted fiducials was then performed on a region of interest centered on the projection of the fiducial when it was within an open field including the case when it was close to the field edge or partially occluded by the leaves. A universal template formula was proposed for template matching and normalized cross correlation was employed for its simplicity and computational efficiency. The search region for every image was adaptively updated through a prediction model that employed the 3D position of the fiducial estimated from the localized positions in previous images. This prediction model allowed the actual fiducial position to be tracked dynamically and was used to initialize the search region. The artifacts caused by electronic interference during the acquisition were effectively removed. A score map was computed by combining both morphological information and image intensity. The pixel location with the highest score was selected as the detected fiducial position. The sets of cine MV images taken during treatment were analyzed with in-house developed software written in MATLAB (The Mathworks, Inc., Natick, MA). Five prostate patients were analyzed to assess the algorithm performance by measuring their positioning accuracy during treatment. Results: The algorithm was able to accurately localize the fiducial position on MV images with success rates of more than 90% per case. The percentage of images in which each fiducial was localized in the studied cases varied between 23% and 65%, with at least one fiducial having been localized between 40% and 95% of the images. This depended mainly on the modulation of the plan and fiducial blockage. The prostate movement in the presented cases varied between 0.8 and 3.5 mm (mean values). The maximum displacement detected among all patients was of 5.7 mm. Conclusions: An algorithm for automatic detection of fiducial markers in cine MV images has been developed and tested with five clinical cases. Despite the challenges posed by complex beam aperture shapes, fiducial localization close to the field edge, partial occlusion of fiducials, fast leaf and gantry movement, and inherently low MV image quality, good localization results were achieved in patient images. This work provides a technique for enabling real-time accurate fiducial detection and tumor tracking during VMAT treatments without the use of extra imaging dose.

60 APPLIED LIFE SCIENCES↗

Hardware Implementation of a Bilateral Subtraction Filter

A bilateral subtraction filter has been implemented as a hardware module in the form of a field-programmable gate array (FPGA). In general, a bilateral subtraction filter is a key subsystem of a high-quality stereoscopic machine vision system that utilizes images that are large and/or dense. Bilateral subtraction filters have been implemented in software on general-purpose computers, but the processing speeds attainable in this way even on computers containing the fastest processors are insufficient for real-time applications. The present FPGA bilateral subtraction filter is intended to accelerate processing to real-time speed and to be a prototype of a link in a stereoscopic-machine- vision processing chain, now under development, that would process large and/or dense images in real time and would be implemented in an FPGA. In terms that are necessarily oversimplified for the sake of brevity, a bilateral subtraction filter is a smoothing, edge-preserving filter for suppressing low-frequency noise. The filter operation amounts to replacing the value for each pixel with a weighted average of the values of that pixel and the neighboring pixels in a predefined neighborhood or window (e.g., a 9 9 window). The filter weights depend partly on pixel values and partly on the window size. The present FPGA implementation of a bilateral subtraction filter utilizes a 9 9 window. This implementation was designed to take advantage of the ability to do many of the component computations in parallel pipelines to enable processing of image data at the rate at which they are generated. The filter can be considered to be divided into the following parts (see figure): a) An image pixel pipeline with a 9 9- pixel window generator, b) An array of processing elements; c) An adder tree; d) A smoothing-and-delaying unit; and e) A subtraction unit. After each 9 9 window is created, the affected pixel data are fed to the processing elements. Each processing element is fed the pixel value for its position in the window as well as the pixel value for the central pixel of the window. The absolute difference between these two pixel values is calculated and used as an address in a lookup table. Each processing element has a lookup table, unique for its position in the window, containing the weight coefficients for the Gaussian function for that position. The pixel value is multiplied by the weight, and the outputs of the processing element are the weight and pixel-value weight product. The products and weights are fed to the adder tree. The sum of the products and the sum of the weights are fed to the divider, which computes the sum of products the sum of weights. The output of the divider is denoted the bilateral smoothed image. The smoothing function is a simple weighted average computed over a 3 3 subwindow centered in the 9 9 window. After smoothing, the image is delayed by an additional amount of time needed to match the processing time for computing the bilateral smoothed image. The bilateral smoothed image is then subtracted from the 3 3 smoothed image to produce the final output. The prototype filter as implemented in a commercially available FPGA processes one pixel per clock cycle. Operation at a clock speed of 66 MHz has been demonstrated, and results of a static timing analysis have been interpreted as suggesting that the clock speed could be increased to as much as 100 MHz.

Huertas, Andres↗

A Transmission Electron Microscopy Study of a Refractory Metal Grain from a Calcium-Aluminum-Rich Inclusion in the Leoville CV3 Chondrite

Introduction: Calcium-aluminum-rich inclusions (CAIs) are an important component of chondritic meteorites. They can contain materials that are thermodynamically predicted and isotopically age dated to be among the first-formed solids in our solar system [1-5]. Observed in some CAIs are micron to sub-micron sized inclusions rich in Fe, Ni, and high-Z elements such as Pt, Os, Ir and W, in the form of refractory metal nuggets (RMNs), fremdlinges, and ‘nugget like objects’ (NLOs) [1,6]. Refractory siderophile elements such as Os, Ir and Ru are thermodynamically predicted to condense at temperatures well in excess of the major CAI phases such as melilite, perovskite, spinel and hibonite [2,7-9]. These refractory metal inclusions in CAIs can therefore serve as probes into the thermodynamic landscape of the early solar protoplanetary disk. Here we report on a refractory grain identified in a CAI of the Leoville CV3 chondrite. This work is part of an ongoing effort to gain insight into the thermochemistry of the early solar system through systematic analyses of the structure and chemistry of various components in CAIs [10-13]. Sample and Analytical Techniques: A fluffy type A CAI (Fig. 1A) was identified in a section of the Leoville, CV3 chondrite (Center for Meteorite Studies, Arizona State University collection, #821_C_3) using a JOEL-JXA 8530F electron microprobe at Arizona State University. Backscattered electron (BSE) imaging and energy-dispersive X-ray spectroscopy (EDS) were used to identify refractory metal grains in the CAI using a Thermo Fisher (formerly FEI) Helios NanoLab 660 G3 focused-ion-beam scanning-electron microscope (FIBSEM) located at the Kuiper Materials Imaging and Characterization Facility (KMICF) at the Lunar and Planetary Laboratory, University of Arizona. The FIB is equipped with an EDAX EDS system. We selected one of the larger (micron-sized) refractory metal grains, designated as ‘Spud’ (Fig. 1B) for further analysis. ‘Spud’ was extracted and thinned to electron transparency (<100 nm) using the FIB-SEM located in KMICF, following methods described by [14- 15]. The FIB section was analyzed using a 200 keV Hitachi HF5000 scanning transmission electron microscope (S/TEM) located at KMICF. The HF5000 is equipped with cold-field emission gun, 3rd-order spherical aberration corrector for STEM imaging, and an Oxford Instruments X-Max N 100 TLE energydispersive spectroscopy (EDS) system with dual 100 mm2 windowless silicon-drift detectors (Ω = 2.0 sr). Selected-area electron-diffraction (SAED) patterns were acquired to aid in determination of crystallinity and phase. Results: The mineralogy, texture, and morphology of the CAI are consistent with that of a fluffy type A (FTA) CAI [16]. BSE imaging at high magnifications revealed grains with high contrast, indicative of compositions rich in elements of higher atomic number relative to surrounding material. These high-Z grains have sizes that range from ∼250 nm to 4 µm. EDS analyses confirm that the bright grains are metal-rich inclusions. A minor fraction of the grains are composed of only Fe and Ni, but the majority (∼60%) of the identified inclusions also contained various refractory siderophiles including Os, Ru, Zr, Ir and Mo. EDS analysis on the FIB-SEM of Spud shows that it contains Fe, Ni, Mo and Ru. High-angle annular dark-field (HAAFD) imaging and EDS mapping in the TEM (Fig. 2) show that Spud occurs in melilite (Ca1.9Al1.99Si1.06O7). Spud contains a subhedral to anhedral morphology and is compositionally heterogenous (polyphasic, Fig. 2). Local spatial correlation occurs among Fe, Ni, and Pt, and also among Os, Ru, and Mo. SAED patterns show that the Fe-Ni-Pt, Fe-Os-Mo-Ru and Fe-Pt regions are crystalline. Discussion: CAIs can contain various types of inclusions rich in Fe, Ni and refractory siderophiles such as Os, Ru, W and Pt [1]. RMNs are micron-sized, single phase alloy grains and can contain Os, Ir, Ru and Rh [1,7,17]. NLOs are also micron-sized inclusions, but contain two phases, a refractory metal, and an oxide [6]. Fremdlinge are the largest of such inclusions (tens of microns in size) and are complex aggregates of Fe-Ni alloy, silicates, oxides, and sulfides [1,17]. While the size of Spud matches previous descriptions of RMNs and NLOs, Spud is neither a single-phase alloy like RMNs, nor does it contain one metal phase and one oxide like NLOs. Spud does not match the above described categories of refractory metal inclusions. The presence of refractory siderophiles such as Mo, Os, Ru, and Pt suggests a high-temperature origin. Thermodynamic modelling by [7] indicates condensation temperatures of 1917 K, 1693 K, 1613 K and 1415 K for Os, Mo, Ru and Pt respectively. These models also show that following the initial condensation of a refractory metal, alloying of solutes such as Fe, Ni and W, occurs in levels proportional to their partial pressures in the surrounding gas. Such alloying occurs at temperatures above the condensation temperatures of common CAI phases such as melilite (1529 K), perovskite (1441 K), spinel (197 K) and forsterite (1354 K) [2]. The polyphasic nature of Spud could be the result of such high-temperature alloying, possibly shortly after the condensation of Mo and Ru at 1693 K and 1613 K respectively. That Spud occurs as an inclusion is consistent with it having formed prior to and at temperature above that of its host melilite in this FTA CAI, which is qualitatively consistent with such prior thermodynamic modeling. Further, the polyphasic nature of Spud is similar to refractory grains from a FTA CAI in the Northwest Africa (NWA) 8323, CV3 chondrite [11-13]. These data suggest that such refractory metal grains could have been widespread in the inner and early solar protoplanetary disk and represent some of the earliest formed solids to have condensed. Acknowledgments: Research and instrumentation supported by NASA grants #NNX12AL47G, #NNX15AJ22G and #80NSSC19K0509, and NSF grants #1531243 and #0619599. Fig 2. STEM data on ‘Spud’. HAADF Image (Top) False-color EDS Maps (Bottom) References: [1] MacPherson G. J. (2014) T. of Geochem. Vol I: Met. And Cosmochem. Processes, 139-179. [2] Lodders K. (2003) ApJ, 591, 1220-1247. [3] Ebel D. S. (2006) Met. & the Early S. Sys. II., 253- 277. [4] Amelin Y. (2002) Science, 297, 1678-1683. [5] Connelly J.N. et.al. (2012) Science, 338, 651-655. [6] Schwander D. et al. (2015) GCA, 18, 70-87. [7] Palme H. and Wlotzka F. (1976) EPSL, 33, 45-60. [8] Berg T. et al. (2009) ApJ, 702, 172-176. [9] Liffman K. et al. (2021) Icarus, 221, 89-105. [10] Zega T.J. et al. (2021) PSJ, 2, 115. [11] Ramprasad T. et al. (2020) LPSC LI, Abstract #2472. [12] Ramprasad T. et al. (2021) Microscopy & Microanalysis, S1, 2792-2794. [13] Ramprasad T. et al. (2021) 84th MetSoc, Abstract #6123. [14] Zega T.J. et al. (2007) MAPS, 42, 1373-1386. [15] Ramprasad T. et al. (2022) MAPS, in revision. [16] Grossman L. (1975), GCA, 39, 433-454. [17] El Goresy A. et al. (1978) LPSC IX, Abstract#1100

T. Ramprasad↗

Optimization of a single module of CUPID

CUPID is the next generation experiment which will use scintillating cryogenic calorimeters to search for the neutrinoless double β decay. This unobserved process would shed light on the nature of the neutrino, which up to our knowledge could be a Majorana or a Dirac particle, and would give us an important hint to explain the lack of antimatter in the universe. This ambitious search needs a detector with unique characteristics such as an extremely low background level and an excellent energy resolution. CUPID is now in advanced R&D state to optimize the detector design in order to completely exploit the potentialities of scintillating cryogenic calorimeters. In the following I will describe the test performed at the LNGS (Laboratori Nazionali del Gran Sasso) of a single module of the future CUPID detector. In this contribution we present the performance obtained with a novel assembly concept, proving that it matches the requirements for CUPID.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

TH-CD-207A-03: A Surface Deformation Driven Respiratory Model for Organ Motion Tracking in Lung Cancer Radiotherapy

Purpose: To propose and validate a novel real-time surface-mesh-based internal organ-external surface motion and deformation tracking method for lung cancer radiotherapy. Methods: Deformation vector fields (DVFs) which characterizes the internal and external motion are obtained by registering the internal organ and tumor contours and external surface meshes to a reference phase in the 4D CT images using a recent developed local topology preserved non-rigid point matching algorithm (TOP). A composite matrix is constructed by combing the estimated internal and external DVFs. Principle component analysis (PCA) is then applied on the composite matrix to extract principal motion characteristics and finally yield the respiratory motion model parameters which correlates the internal and external motion and deformation. The accuracy of the respiratory motion model is evaluated using a 4D NURBS-based cardiac-torso (NCAT) synthetic phantom and three lung cancer cases. The center of mass (COM) difference is used to measure the tumor motion tracking accuracy, and the Dice’s coefficient (DC), percent error (PE) and Housdourf’s distance (HD) are used to measure the agreement between the predicted and ground truth tumor shape. Results: The mean COM is 0.84±0.49mm and 0.50±0.47mm for the phantom and patient data respectively. The mean DC, PE and HD are 0.93±0.01, 0.13±0.03 and 1.24±0.34 voxels for the phantom, and 0.91±0.04, 0.17±0.07 and 3.93±2.12 voxels for the three lung cancer patients, respectively. Conclusions: We have proposed and validate a real-time surface-mesh-based organ motion and deformation tracking method with an internal-external motion modeling. The preliminary results conducted on a synthetic 4D NCAT phantom and 4D CT images from three lung cancer cases show that the proposed method is reliable and accurate in tracking both the tumor motion trajectory and deformation, which can serve as a potential tool for real-time organ motion and deformation monitoring in lung cancer radiotherapy. This work is supported in part by grant from VARIAN MEDICAL SYSTEMS INC, the National Natural Science Foundation of China (no 81428019 and no 81301940), the Guangdong Natural Science Foundation (2015A030313302)and the 2015 Pearl River S&T Nova Program of Guangzhou (201506010096).

60 APPLIED LIFE SCIENCES↗