Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Performance Rating Method”

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 937 records · Page 52

Efficient method for the measurement of lifetime optical damage performance of thin film coatings from laser damage size analysis

Here, a laser damage test method based on damage size analysis (DSA) is described that simplifies the derivation of the lifetime optical damage threshold of film materials critical in the design of devices used in high-repetition-rate, high-power laser systems. The DSA method presented here is solely based on imaging to measure the damage site size produced from exposure to a known Gaussian-shaped beam with a fixed, systematically selected fluence well above the ablation threshold. The method locates the damage boundary produced from repeated exposures, using images with a high contrast, and maps it to the beam profile to extract a lifetime laser damage fluence threshold value. We validate the DSA approach using a few relevant transparent film material systems and by comparing it to the standard S /1 laser damage test method. The DSA method can be more efficient and accelerate materials development and validation necessary to support the design of high-power repetition-rated lasers and optoelectronic devices.

36 MATERIALS SCIENCE↗

Working fluid selection for organic Rankine cycle power generation using hot produced supercritical CO 2 from a geothermal reservoir

Geothermal heat mining simulations using supercritical CO 2 (sCO 2 ) were performed in this research. Working fluid selection criteria for power generation using sCO 2 from a geothermal reservoir are then presented for subcritical, superheated and supercritical organic Rankine cycles (ORCs). Meanwhile, method of working fluid classification for ORC is proposed. To get the most feasible ORC design, this study introduces the concept of “turning point” for isentropic and dry working fluids, as well as minimum turbine inlet temperature for wet working fluids. A thermodynamic model was developed with capabilities to obtain the optimal working fluid mass flow rate, evaporation temperature, superheated temperature, and supercritical pressure, to evaluate the thermal performance of the three ORC approaches using hot produced sCO 2 . With this model, thirty potential working fluids with critical temperatures in the range from 50 to 225 °C were screened for utilizing hot produced sCO 2 considering physical properties, environmental and safety impacts, and thermodynamic performances. Finally, the thermodynamic results were compared for all possible working fluids.

15 GEOTHERMAL ENERGY↗

SU-E-J-200: Operation of An Electron Accelerator On An Integrated MR-Linac System

Purpose: An integrated MRI guided radiotherapy system poses a challenge of operating a linear accelerator in the presence of a magnetic field as the magnetic force acting on the electrons could Result in radiation source displacement and subsequent reduction of dose output. It was the purpose of this work to test the performance of a linac in the presence of a 1.5T MRI system. Methods: The first experimental MRI guided radiotherapy system at UMC-Utrecht consisting of an Elekta linac rotating around a 1.5T Magnex magnet was examined. A passive magnetic shield was simulated, designed and installed to reduce the influence of the MRI magnet stray field on the electron beamline. The B field inside the shield was measured as a function of gantry angle and measurements of dose rate constancy upon gantry rotation were performed. Results: The magnitude of the magnetic field on the electron beam path without the shield was as high as 70G. It varied by up to 15G with gantry rotation due to the presence of metal beams in the bunker floor which resulted in dose output drop of up to 70% at certain gantry angles. With the prototype shield, field magnitude was reduced to well below 0.5G everywhere along the electron beam path. Field variation with gantry rotation was decreased to below 0.2G and enabled dose output of the linac to be recovered at all gantry angles. The homogeneity of the field inside the MRI magnet has not been compromised. Conclusion: It was demonstrated that the influence of the 1.5T magnet and the bunker design on the linac operation has been minimised. The performance will be further improved on the Elekta Atlantic system which incorporates a newly developed and optimised Philips magnet design and bunker construction. J Harasimowicz, D Roberts, I Shinton and S Sund are employed by Elekta Limited Crawley, H Wang and M Zhong are employed by Elekta Beijing Medical Systems Co. Ltd., J Overweg is employed by Philips Technologie GmbH Forschungslaboratorien.

60 APPLIED LIFE SCIENCES↗

Model-Based Extracted Water Desalination System for Carbon Sequestration

Over the last 1.5 years, GE Global Research and Pennsylvania State University defined a model-based, scalable, and multi-stage extracted water desalination system that yields clean water, concentrated brine, and, optionally, salt. The team explored saline brines that ranged across the expected range for extracted water for carbon sequestration reservoirs (40,000 up to 220,000 ppm total dissolved solids, TDS). In addition, the validated the system performance at pilot scale with field-sourced water using GE’s pre-pilot and lab facilities. This project encompassed four principal tasks, in addition to Project Management and Planning: 1) identify a deep saline formation carbon sequestration site and a partner that are suitable for supplying extracted water; 2) conduct a techno-economic assessment and down-selection of pre-treatment and desalination technologies to identify a cost-effective system for extracted water recovery; 3) validate the downselected processes at the lab/pre-pilot scale; and 4) define the scope of the pilot desalination project. Highlights from each task are described below: Deep saline formation characterization The deep saline formations associated with the five DOE NETL 1260 Phase 1 projects were characterized with respect to their mineralogy and formation water composition. Sources of high TDS feed water other than extracted water were explored for high TDS desalination applications, including unconventional oil and gas and seawater reverse osmosis concentrate. Technoeconomic analysis of desalination technologies Techno-economic evaluations of alternate brine concentration technologies, including humidification-dehumidification (HDH), membrane distillation (MD), forward osmosis (FO), turboexpander-freeze, solvent extraction and high pressure reverse osmosis (HPRO), were conducted. These technologies were evaluated against conventional falling film-mechanical vapor recompression (FF-MVR) as a baseline desalination process. Furthermore, a quality function deployment (QFD) method was used to compare alternate high TDS desalination technologies to FF-MVR. High pressure reverse osmosis was found to a be a promising alternative desalination technology. A deep-dive technoeconomic analysis of HPRO was performed, including Capex and Opex estimates, for seawater RO (SWRO). Additionally, two additional cases were explored: 1) a comparison of a SWRO plus HPRO system to the option of doubling the size of a standard seawater RO system to achieve the same total pure water recovery rate; and 2) a flue gas desulfurization wastewater treatment zero-liquid discharge (ZLD) application, where preconcentration with RO (SWRO or SWRO + HPRO) before evaporation and crystallization was compared to FF-MVR and crystallization technologies without RO preconcentration. Pre-pilot process validation Pre-pilot-scale tests were conducted using field production water to validate key process steps for extracted water pretreatment. Approximately 5,000 gallons of field produced water was processed through, microfiltration, ultrafiltration, and steam regenerable sorbent operations. Smaller quantities were processed through microclarification. In addition, analytical methods (purge-and-trap gas chromatography and Hach TOC analytical methods) were validated. Lab-scale HPRO elements were constructed and tested at high pressures, to identify and mitigate technical risks of the technology. Lastly, improvements in RO membrane materials were identified as the necessary next step to achieve further improvement in element performance at high pressure. Scope of Field Pilot A field pilot for extracted water pretreatment was designed.

20 FOSSIL-FUELED POWER PLANTS↗

Accelerating transients with NekRS: GPU overlapping domain implementation and multi-rate timestepping

The simulation of nuclear transients using Computational Fluid Dynamics (CFD) presents significant computational challenges due to the inherent complexity and the wide separation in temporal scales between various flow physical phenomena. These disparities lead to high computational costs, often making the simulation of transients impractical without advanced techniques. Consequently, multiple research initiatives are being pursued by the NEAMS thermal-hydraulic area, some driven by academic institutions and some by national laboratories. Overall, they are exploring novel methods to make transient simulations more feasible and efficient. This report delves into recent advancements within the CFD code NekRS, specifically those achieved in Fiscal Year 2024 under the CONNECT effort, aimed at improving the performance and feasibility of transient simulations. The first major advancement involves the porting of NekRS to Aurora, one of the Department of Energy’s (DOE) most powerful supercomputers. Additionally, the report discusses the implementation of an overlapping domain capability within NekRS. This novel GPU-accelerated capability allows different spatial regions of the domain to be solved independently, enhancing the code’s efficiency, particularly when running large-scale simulations in complex domains. The scalability of this approach is demonstrated, highlighting its potential to transform how transients are approached in CFD simulations. Lastly, the report focuses on how this overlapping domain capability specifically accelerates transient simulations through multi-rate timestepping. By decoupling different regions and facilitating faster computations, this method offers a promising pathway to making nuclear transient simulations more computationally feasible, addressing one of the critical bottlenecks in the field. Together, these advancements represent a significant leap forward in transient simulation technology, bringing closer the possibility of handling highly complex nuclear scenarios with greater efficiency and accuracy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

SU-F-BRE-16: VMAT Commissioning and Quality Assurance (QA) of An Elekta Synergy-STM Linac Using ICOM Test HarnessTM

Purpose: To establish a set of tests based on the iCOM software that can be used to commission and perform periodic QA of VMAT delivery on the Elekta Synergy-S, commonly known as the Beam Modulator (BM). Methods: iCOM is used to create and deliver customized treatment fields to characterize the system in terms of 1) MLC positioning accuracy under static and dynamic delivery with full gantry rotation, 2) MLC positioning with known errors, 3) Maximum dose rate, 4) Maximum MLC speed, 5) Maximum gantry speed, 6) Synchronization: gantry speed versus dose rate, and 7) Synchronization: MLC speed versus dose rate. The resulting images were captured on the iView GT and exported in DICOM format to Dosimetry Check™ system for visual and quantitative analysis. For the initial commissioning phase, the system tests described should be supplemented with extensive patient QAs covering all clinically relevant treatment sites. Results: The system performance test suite showed that on our Synergy-S, MLC positioning was accurate under both static and dynamic deliveries. Intentional errors of 1 mm were also easily identified on both static and dynamic picket fence tests. Maximum dose rate was verified with stop watch to be consistently between 475-480 MU/min. Maximum gantry speed and MLC speed were 5.5 degree/s and 2.5 cm/s respectively. After accounting for beam flatness, both synchronization tests, gantry versus dose rate and MLC speed versus dose rate, were successful as the fields were uniform across the strips and there were no obvious cold/hot spots. Conclusion: VMAT commissioning and quality assurance should include machine characterization tests in addition to patient QAs. Elekta iCOM is a valuable tool for the design of customized VMAT field with specific MU, MLC leaf positions, dose rate, and indirect control of MLC and gantry speed at each of its control points.

07 ISOTOPE AND RADIATION SOURCES↗

Carbon-Binder Weight Loading Optimization for Improved Lithium-Ion Battery Rate Capability

Battery performance is strongly correlated with electrode microstructure and weight loading of the electrode components. Among them are the carbon-black and binder additives that enhance effective conductivity and provide mechanical integrity. However, these both reduce effective ionic transport in the electrolyte phase and reduce energy density. Therefore, an optimal additive loading is required to maximize performance, especially for fast charging where ionic transport is essential. Such optimization analysis is however challenging due to the nanoscale imaging limitations that prevent characterizing this additive phase and thus quantifying its impact on performance. Herein, an additive-phase generation algorithm has been developed to remedy this limitation and identify percolation threshold used to define a minimal additive loading. Improved ionic transport coefficients from reducing additive loading has been then quantified through homogenization calculation, macroscale model fitting, and experimental symmetric cell measurement, with good agreement between the methods. Rate capability test demonstrates capacity improvement at fast charge at the beginning of life, from 37% to 55%, respectively for high and low additive loading during 6C CC charging, in agreement with macroscale model, and attributed to a combination of lower cathode impedance, reduced electrode tortuosity and cathode thickness.

25 ENERGY STORAGE↗

High-order shock-fitted detonation propagation in high explosives

In this paper, a highly accurate numerical shock and material interface fitting scheme composed of fifth-order spatial and third- or fifth-order temporal discretizations is applied to the two-dimensional reactive Euler equations in both slab and axisymmetric geometries. High rates of convergence are not typically possible with shock-capturing methods as the Taylor series analysis breaks down in the vicinity of discontinuities. Furthermore, for typical high explosive (HE) simulations, the effects of material interfaces at the charge boundary can also cause significant computational errors. Fitting a computational boundary to both the shock front and material interface (i.e. streamline) alleviates the computational errors associated with captured shocks and thus opens up the possibility of high rates of convergence for multi-dimensional shock and detonation flows. Several verification tests, including a Sedov blast wave, a Zel'dovich–von Neumann–Döring (ZND) detonation wave, and Taylor–Maccoll supersonic flow over a cone, are utilized to demonstrate high rates of convergence to nontrivial shock and reaction flows. Comparisons to previously published shock-capturing multi-dimensional detonations in a polytropic fluid with a constant adiabatic exponent (PF-CAE) are made, demonstrating significantly lower computational error for the present shock and material interface fitting method. For an error on the order of 10 m / s , which is similar to that observed in experiments, shock-fitting offers a computational savings on the order of 1000. In addition, the behavior of the detonation phase speed is examined for several slab widths to evaluate the detonation performance of PBX 9501 while utilizing the Wescott–Stewart–Davis (WSD) model, which is commonly used in HE modeling. It is found that the thickness effect curve resulting from this equation of state and reaction model using published values is dramatically more steep than observed in recent experiments. Utilizing the present fitting strategy, in conjunction with a nonlinear optimizer, a new set of reaction rate parameters improves the correlation of the model to experimental results. In conclusion, this new model is tested against two dimensional slabs as a validation test.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

The Variable Saturation Hydraulic Conductivity Method for Improving Soil Water Content Simulation in EPIC and APEX Models

Core Ideas We propose VSHC for calculating soil water percolation in continuous models. VSHC relates permeability to daily hydraulic properties to estimate hydraulic conductivity. VSHC performed well at field scale and improved estimates of peak streamflows. Soil water percolation is a key process in the life cycle of water in fields, watersheds, and river basins. The Environmental Policy Integrated Climate (EPIC) and the Agricultural Policy/Environmental eXtender (APEX) are continuous models developed for evaluating the environmental effects of agricultural management. Traditionally, these models have simulated soil water percolation processes using a tipping‐bucket approach, with the rate of flow limited by the saturated hydraulic conductivity. This simple approach often leads to inaccuracy in simulating elevated soil water conditions where soil water content (SWC) levels may remain above field capacity under prolonged wet weather periods or limited drainage. To overcome this deficiency, a new submodel, the variable saturation hydraulic conductivity (VSHC) method, was developed for simulating soil water percolation processes using a nonlinear equation to estimate the effective hydraulic conductivity as a function of the SWC and soil properties. The VSHC method was evaluated at three sites in the United States and two sites in Europe. In addition, a numerical solution of the Richards equation was used as a benchmark for SWC comparison. Results show that the VSHC method substantially improves the accuracy of the SWC simulation in long‐term simulations, particularly during wet periods. At the watershed scale, results on the Riesel Y2 watershed indicate that the VSHC method enhances model performance in the high‐flow regime of channel peak flows because of the improved estimation of SWC, which implies that the improved SWC simulation at the field scale is beneficial to hydrologic modeling at the watershed scale.

Doro, Luca↗

Rates of Ligand Exchange around the Bis-Oxalato Complex [NpO 2 (C 2 O 4 ) 2 ] 3- Measured by Using Multinuclear NMR Spectroscopy under Neutral to Semi-Alkaline Conditions [Rates of ligand exchange around the bis-oxalato complex NpO2(C2O4)23- measured using multinuclear NMR at neutral to semi-alkaline conditions]

The kinetics of ligand exchange between the free oxalate ion, C 2 O 4 2- , and the bis-oxalato Np V complex, [NpO 2 (C 2 O 4 ) 2 ] 3- , in aqueous solution are reported by using 13 C and 17 O NMR spectroscopy methods. Rates of exchange were measured in the pH regime of 6.5–9.0, at which speciation is shown to be suitably simple. Because the neptunium(V) complex is paramagnetic, the rates of ligand exchange were estimated by following the width of the 13 C and 17 O signals assigned to the free oxalate ion in solution and by applying the Swift–Connick method for measuring rates of exchange. A set of experiments were conducted in which pH and total oxalate concentration were varied, and the linear dependence of the rate on these parameters was demonstrated. Variable-temperature NMR spectroscopy was also performed to measure activation parameters of complexation. Finally, at pH<8.0, ΔH ≠ =16.9 ±4.9 kJ mol -1 and ΔS≠=-116.3 ±17.1kJ mol -1 K -1 , whereas at pH>8.0 there is almost no dependence on temperature, which is interpreted to indicate that hydrolysis is coupled to ligand exchange under these conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A data-centric weak supervised learning for highway traffic incident detection

Using the data from loop detector sensors for near-real-time detection of traffic incidents on highways is crucial to averting major traffic congestion. While recent supervised machine learning methods offer solutions to incident detection by leveraging human-labeled incident data, the false alarm rate is often too high to be used in practice. Specifically, the inconsistency in the human labeling of the incidents significantly affects the performance of supervised learning models. To that end, we focus on a data-centric approach to improve the accuracy and reduce the false alarm rate of traffic incident detection on highways. We develop a weak supervised learning workflow to generate high-quality training labels for the incident data without the ground truth labels, and we use those generated labels in the supervised learning setup for final detection. This approach comprises three stages. First, we introduce a data preprocessing and curation pipeline that processes traffic sensor data to generate high-quality training data through leveraging labeling functions, which can be domain knowledge-related or simple heuristic rules. Second, we evaluate the training data generated by weak supervision using three supervised learning models-random forest, k-nearest neighbors, and a support vector machine ensemble-and long short-term memory classifiers. The results show that the accuracy of all of the models improves significantly after using the training data generated by weak supervision. Third, we develop an online real-time incident detection approach that leverages the model ensemble and the uncertainty quantification while detecting incidents. Finally, we show that our proposed weak supervised learning workflow achieves a high incident detection rate (0.90) and low false alarm rate (0.08).

97 MATHEMATICS AND COMPUTING↗

A method for including external feed in depletion calculations with CRAM and implementation into ORIGEN

A method for including external feed with polynomial time dependence in depletion calculations with the Chebyshev Rational Approximation Method (CRAM) is presented and the implementation of CRAM to the ORIGEN module of the SCALE suite is described. In addition to being able to handle time-dependent feed rates, the new solver also adds the capability to perform adjoint calculations. Results obtained with the new CRAM solver and the original depletion solver of ORIGEN are compared to high precision reference calculations, which shows the new solver to be orders of magnitude more accurate. Lastly, in most cases, the new solver is up to several times faster due to not requiring similar substepping as the original one.

97 MATHEMATICS AND COMPUTING↗

Turbine scale and siting considerations in wind plant layout optimization and implications for capacity density

Improvements in wind energy technology, reduced costs, and ambitious clean energy goals have led to projections of high wind contribution in coming years. Developing methodologies to design wind plants with a variety of siting constraints and turbine sizes helps enable high wind penetration, and gain a better understanding of how wind plants are sensitive to setback constraints and turbine design. In this paper, we present a two-step optimization method to simultaneously determine the optimal number of turbines and their locations in a wind plant domain divided into many small, discrete parcels. We present the optimized performance metrics of a wind plant optimized with different turbine sizes and ratings, and with different siting restrictions within the wind plant. Our results indicate that taller and larger turbines are more sensitive to increasing siting constraints. We also compare the optimal wind plant layouts and performance for wind plants optimized for minimum COE and maximum profit. Wind plants optimized for profit had 130%-190% of the capacity of plants optimized for COE, which demonstrates that the optimal results are greatly affected by the objective function, which should be carefully considered. Finally, in this paper we demonstrate the effect of increasing siting constraints on wind plant capacity density, and how the results change when different land areas are used to calculate capacity density. When using the entire wind plant boundary area to determine capacity density, increasing siting constraints decreases the capacity density. However, when we only use the available area (the area left after removing the siting constraints) to calculate the capacity density, increasing the siting constraints increases capacity density. This is a critical insight because of how capacity density is typically defined and used in research, and has important implications for assessment of technical potential and capacity expansion modeling, as well as future wind deployment potential.

17 WIND ENERGY↗

Electronic nose with a new feature reduction method and a multi-linear classifier for Chinese liquor classification

An electronic nose (e-nose) was designed to classify Chinese liquors of the same aroma style. A new method of feature reduction which combined feature selection with feature extraction was proposed. Feature selection method used 8 feature-selection algorithms based on information theory and reduced the dimension of the feature space to 41. Kernel entropy component analysis was introduced into the e-nose system as a feature extraction method and the dimension of feature space was reduced to 12. Classification of Chinese liquors was performed by using back propagation artificial neural network (BP-ANN), linear discrimination analysis (LDA), and a multi-linear classifier. The classification rate of the multi-linear classifier was 97.22%, which was higher than LDA and BP-ANN. Finally the classification of Chinese liquors according to their raw materials and geographical origins was performed using the proposed multi-linear classifier and classification rate was 98.75% and 100%, respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Realistic operation of two residential cordwood-fired outdoor hydronic heater appliances—Part 2: Particle number and size

The use of wood as a fuel for home heating is a concern from an environmental health and safety perspective as biomass combustion appliances emit high concentrations of particulate matter. Wood burning significantly contributes to wintertime particulate matter concentrations in many states in the northern United States. Of particular concern are outdoor wood-fired hydronic heaters. These devices are concerning as they tend to have very large combustion chambers and typical use patterns can result in long periods of low output, which result in an increased particulate matter emission rate relative to high heat output operating conditions. In this study, the performance of two hydronic heaters operating under different combustion conditions, including four different heat output categories approximately corresponding to categories I–IV denoted in Environmental Protection Agency Method 28 Outdoor Wood-fired Hydronic Heaters, and during start-up and reloading events were investigated. Measurements of flue gas particulate number concentration and size for particles with aerodynamic diameters between 0.006 and 10 µm were made using a dilution sampling system. The measured particle number concentration in the flue gas was between 0.71 and 420 million particles per cubic centimeter and was dependent on fuel loading and heat output. For each hydronic heater tested, the highest average particle concentration was found at the beginning of each test during the cold-start condition. Additionally, the majority of the particles had aerodynamic diameters less than 0.100 µm (particles of this size made up between 64% and 97% of all particles) and less than 1% of all particles had aerodynamic diameters greater than 1 µm for all phases. For particles in the accumulation mode, between 0.100 and 1 µm, the mean particle diameter was dependent on fuel loading and heat output.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Network Detection of Radiation Sources Using Localization-Based Approaches

Radiation source detection is a critical problem in homeland security-related applications. Deploying a network of detectors is expected to provide improved detection due to the combined, albeit dispersed, capture area of multiple detectors. Recently, localization-based detection algorithms provided performance gains beyond the simple “aggregated” area as a result of localization being enabled by the networked detectors. We introduce the following three localization-based detection approaches: 1) source-attractor radiation detection (SRD); 2) triangulation-based radiation source detection (TriRSD); and 3) the ratio of square distance-based radiation source detection (ROSD-RSD). We use canonical datasets from Domestic Nuclear Detection Office's intelligence radiation sensors systems tests to assess the performance of these methods. Extensive results demonstrate that SRD outperforms TriRSD and ROSD-RSD, and other existing detection algorithms based on the sequential probability ratio test and maximum likelihood estimation in terms of both false alarm and detection rates.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Deep Reinforcement Learning From Demonstrations to Assist Service Restoration in Islanded Microgrids

Microgrids can be operated in island mode during utility grid outages to support service restoration and improve system resilience. To schedule and dispatch distributed energy resources (DERs) in an islanded microgrid, conventional model-based methods rely on accurate distribution network models and lack generalization and adaptability. Data-driven methods are promising for DER coordination but face practical challenges such as potential hazards to microgrids during online training and insufficient online training opportunities due to low outage rates. This paper presents a novel two-stage learning framework that builds on the deep deterministic policy gradient from demonstrations to identify an optimal restoration strategy. At the pre-training stage, imitation learning is applied to equip the control agent with expert experiences to guarantee acceptable initial performance. At the online training stage, action clipping, reward shaping, and expert demonstrations are leveraged to ensure safe exploration while accelerating the training process. In conclusion, the proposed method is illustrated using the IEEE 123-node system and compared with a representative model-based method and the standard deep deterministic policy gradient method to prove solution accuracy and demonstrate increased computational efficiency.

Du, Yan↗

4D tumor centroid tracking using orthogonal 2D dynamic MRI: Implications for radiotherapy planning

Purpose: Current pretreatment, 4D imaging techniques are suboptimal in that they sample breathing motion over a very limited “snapshot” in time. Heretofore, long-duration, 4D motion characterization for radiotherapy planning, margin optimization, and validation have been impractical for safety reasons, requiring invasive markers imaged under x-ray fluoroscopy. To characterize 3D tumor motion and associated variability over durations more consistent with treatments, the authors have developed a practical dynamic MRI (dMRI) technique employing two orthogonal planes acquired in a continuous, interleaved fashion.Methods: 2D balanced steady-state free precession MRI was acquired continuously over 9–14 min at approximately 4 Hz in three healthy volunteers using a commercial 1.5 T system; alternating orthogonal imaging planes (sagittal, coronal, sagittal, etc.) were employed. The 2D in-plane pixel resolution was 2 × 2 mm{sup 2} with a 5 mm slice profile. Simultaneous with image acquisition, the authors monitored a 1D surrogate respiratory signal using a device available with the MRI system. 2D template matching-based anatomic feature registration, or tracking, was performed independently in each orientation. 4D feature tracking at the raw frame rate was derived using spline interpolation.Results: Tracking vascular features in the lung for two volunteers and pancreatic features in one volunteer, the authors have successfully demonstrated this method. Registration error, defined here as the difference between the sagittal and coronal tracking result in the SI direction, ranged from 0.7 to 1.6 mm (1σ) which was less than the acquired image resolution. Although the healthy volunteers were instructed to relax and breathe normally, significantly variable respiration was observed. To demonstrate potential applications of this technique, the authors subsequently explored the intrafraction stability of hypothetical tumoral internal target volumes and 3D spatial probability distribution functions. The surrogate respiratory information allowed the authors to show how this technique can be used to study correlations between internal and external (surrogate) information over these prolonged durations. However, compared against the gold standard of the time stamps in the dMRI frames, the temporal synchronization of the surrogate 1D respiratory information was shown to be likely unreliable.Conclusions: The authors have established viability of a novel and practical pretreatment, 4D tumor centroid tracking method employing a commercially available dynamic MRI sequence. Further developments from the vendor are likely needed to provide a reliably synchronized surrogate 1D respiratory signal, which will likely broaden the utility of this method in the pretreatment radiotherapy planning context.

62 RADIOLOGY AND NUCLEAR MEDICINE↗