Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “volume estimation”

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 595 records · Page 33

Feature detection and matching on atmospheric nuclear detonation video

Automated feature matching of nuclear detonations (NUDETs) enables three‐dimensional point cloud reconstruction, and establishment a volume‐based model to reduce uncertainty in estimating the yield of NUDETs solely from video. Establishing a volume‐based model requires feature correspondences between wide viewpoints of 58°–110° that are larger than scale‐invariant feature transform‐based techniques can reliably match. The presented technique detects relative bright features in the NUDET known as ‘hotspots,’ and matches them across wide viewpoints using a spherical based object model. Results show that hotspots can be detected with a 71.95% hit rate and 86.03% precision. Hotspots are matched to films from different viewpoints with 76.6% correctness and a standard deviation of 16.4%. Hotspot descriptors are also matched in time sequence with 99.6% correctness and a standard deviation of 1.07%. The results demonstrate that a spherical model can serve as a viable descriptor model for matching across wide viewpoints when the object is known to be spherical. It also demonstrates an automated feature detection and matching combination that enables features to be matched from unsynchronised video across wide viewpoints of 58°–110° on spherical objects where state‐of‐the‐art techniques are insufficient.

Schmitt, Daniel T.↗

Sound Speed as a Candidate for Internal Temperature Monitoring During Solid Phase Processing of Materials

Internal temperature is an important parameter for understanding the solid phase processing (SPP) of materials and improving feedback control of these processes. Friction stir processes, shear assisted processing and extrusion (ShAPE ™ ), and ultrahigh velocity (UHV) cold spray are examples of such processes. Simple thermal models based on temperatures measured at the boundary of the deforming body cannot capture the dynamic nature of temperature generation within the deforming solid. One possibly way of understanding internal temperature is by measuring sound speed, which varies as a function of temperature for materials and directing ultrasonic wave propagation through selected material volumes of interest. In this project, longitudinal wave speed is estimated for aluminum 6061. Piezoelectric materials bonded to a specimen provided data from which longitudinal wave speed was estimated. A discussion evaluated the linear relationship of wave speed versus temperature and postulated how wave speed measurements could be used to estimate the internal temperature of selected material volumes for feedback during SPP of materials. Instrumentation, data and analysis, and a discussion are reviewed.

36 MATERIALS SCIENCE↗

Physiological cross-sectional area of human leg muscles based on magnetic resonance imaging

Magnetic resonance imaging techniques were used to determine the physiological cross-sectional areas (PCSAs) of the major muscles or muscle groups of the lower leg. For 12 healthy subjects, the boundaries of each muscle or muscle group were digitized from images taken at 1-cm intervals along the length of the leg. Muscle volumes were calculated from the summation of each anatomical CSA (ACSA) and the distance between each section. Muscle length was determined as the distance between the most proximal and distal images in which the muscle was visible. The PCSA of each muscle was calculated as muscle volume times the cosine of the angle of fiber pinnation divided by fiber length, where published fiber length:muscle length ratios were used to estimate fiber lengths. The mean volumes of the major plantarflexors were 489, 245, and 140 cm3 for the soleus and medial (MG) and lateral (LG) heads of the gastrocnemius. The mean PCSA of the soleus was 230 cm2, about three and eight times larger than the MG (68 cm2) and LG (28 cm2), respectively. These PCSA values were eight (soleus), four (MG), and three (LG) times larger than their respective maximum ACSA. The major dorsiflexor, the tibialis anterior (TA), had a muscle volume of 143 cm2, a PCSA of 19 cm2, and an ACSA of 9 cm2. With the exception of the soleus, the mean fiber length of all subjects was closely related to muscle volume across muscles. The soleus fibers were unusually short relative to the muscle volume, thus potentiating its force potential.(ABSTRACT TRUNCATED AT 250 WORDS).

NASA Discipline Musculoskeletal↗

Thermodynamic estimation: Ionic materials

Thermodynamics establishes equilibrium relations among thermodynamic parameters (“properties”) and delineates the effects of variation of the thermodynamic functions (typically temperature and pressure) on those parameters. However, classical thermodynamics does not provide values for the necessary thermodynamic properties, which must be established by extra-thermodynamic means such as experiment, theoretical calculation, or empirical estimation. While many values may be found in the numerous collected tables in the literature, these are necessarily incomplete because either the experimental measurements have not been made or the materials may be hypothetical. The current paper presents a number of simple and relible estimation methods for thermodynamic properties, principally for ionic materials. The results may also be used as a check for obvious errors in published values. The estimation methods described are typically based on addition of properties of individual ions, or sums of properties of neutral ion groups (such as “double” salts, in the Simple Salt Approximation), or based upon correlations such as with formula unit volumes (Volume-Based Thermodynamics). - Graphical abstract: Thermodynamic properties of ionic materials may be readily estimated by summation of the properties of individual ions, by summation of the properties of ‘double salts’, and by correlation with formula volume. Such estimates may fill gaps in the literature, and may also be used as checks of published values. This simplicity arises from exploitation of the fact that repulsive energy terms are of short range and very similar across materials, while coulombic interactions provide a very large component of the attractive energy in ionic systems. Display Omitted - Highlights: • Estimation methods for thermodynamic properties of ionic materials are introduced. • Methods are based on summation of single ions, multiple salts, and correlations. • Heat capacity, entropy, lattice energy, enthalpy, Gibbs energy values are available.

36 MATERIALS SCIENCE↗

Reduced order models for prediction of groundwater quality impacts from CO₂ and brine leakage

A careful assessment of the risk associated with geologic CO₂ storage is critical to the deployment of large-scale storage projects. A potential risk is the deterioration of groundwater quality caused by the leakage of CO₂ and brine leakage from deep subsurface reservoirs. In probabilistic risk assessment studies, numerical modeling is the primary tool employed to assess risk. However, the application of traditional numerical models to fully evaluate the impact of CO₂ leakage on groundwater can be computationally complex, demanding large processing times and resources, and involving large uncertainties. As an alternative, reduced order models (ROMs) can be used as highly efficient surrogates for the complex process-based numerical models. In this study, we represent the complex hydrogeological and geochemical conditions in a heterogeneous aquifer and subsequent risk by developing and using two separate ROMs. The first ROM is derived from a model that accounts for the heterogeneous flow and transport conditions in the presence of complex leakage functions for CO₂ and brine. The second ROM is obtained from models that feature similar, but simplified flow and transport conditions, and allow for a more complex representation of all relevant geochemical reactions. To quantify possible impacts to groundwater aquifers, the basic risk metric is taken as the aquifer volume in which the water quality of the aquifer may be affected by an underlying CO₂ storage project. The integration of the two ROMs provides an estimate of the impacted aquifer volume taking into account uncertainties in flow, transport and chemical conditions. These two ROMs can be linked in a comprehensive system level model for quantitative risk assessment of the deep storage reservoir, wellbore leakage, and shallow aquifer impacts to assess the collective risk of CO₂ storage projects.

54 ENVIRONMENTAL SCIENCES↗

Performance Modeling of Materials-Based Hydrogen Storage Systems for Automotive Applications

One of the challenges of fuel cell vehicles is storage of the hydrogen fuel. In an effort to reduce the temperature and pressure limitations of storing neat hydrogen either as a high pressure gas or cryogenic liquid, material-based hydrogen storage is being developed. Hydrogen stored either as an adsorbent or as a chemical compound not only provides more reasonable pressures and temperatures but also improves the storage density of the hydrogen. Modeling tools have been developed as part of the U.S. Department of Energy's Hydrogen Storage Engineering Center of Excellence to assess the performance of these materials in light duty vehicle applications. The purpose of these modeling tools is to evaluate newly developed hydrogen storage materials and ensure that they meet customer expectations for driving range and performance before further investment. Recent advances in these modeling tools have been made to allow direct material property inputs. Such material properties as hydrogen capacity, isotherms and other thermodynamic and kinetic properties can reasonably be measured by materials researchers. Using these properties as inputs, the models can then provide estimates of the mass and volume of the storage system and the driving range and fuel efficiency of the vehicle during transient operation with four different drive cycles. The system sizing estimator and the framework model used to perform this analysis will be described and results presented for adsorbents and chemical based storage materials.

25 ENERGY STORAGE↗

Magnetohydrodynamics MHD Engineering Test Facility ETF 200 MWe power plant. Conceptual Design Engineering Report CDER. Volume 3: Costs and schedules

The estimated plant capital cost for a coal fired 200 MWE electric generating plant with open cycle magnetohydrodynamics is divided into principal accounts based on Federal Energy Regulatory Commision account structure. Each principal account is defined and its estimated cost subdivided into identifiable and major equipment systems. The cost data sources for compiling the estimates, cost parameters, allotments, assumptions, and contingencies, are discussed. Uncertainties associated with developing the costs are quantified to show the confidence level acquired. Guidelines established in preparing the estimated costs are included. Based on an overall milestone schedule related to conventional power plant scheduling experience and starting procurement of MHD components during the preliminary design phase there is a 6 1/2-year construction period. The duration of the project from start to commercial operation is 79 months. The engineering phase of the project is 4 1/2 years; the construction duration following the start of the man power block is 37 months.

Source record↗

DACC program cost and work breakdown structure-dictionary. General purpose aft cargo carrier study, volume 2

Results of detailed cost estimates and economic analysis performed on the updated 201 configuration of the dedicated Aft Cargo Carrier (DACC) are given. The objective of this economic analysis is to provide the National Aeronautics and Space Administration (NASA) with information on the economics of using the DACC on the Space Transportation System (STS). The detailed cost estimates for the DACC are presented by a work breakdown structure (WBS) to ensure that all elements of cost are considered in the economic analysis and related subsystem trades. Costs reported by WBS provide NASA with a basis for comparing competing designs and provide detailed cost information that can be used to forecast phase C/D planning for new projects or programs derived from preliminary conceptual design studies. The scope covers all STS and STS/DACC launch vehicle cost impacts for delivering an orbital transfer vehicle to a 120 NM low Earth orbit (LEO).

Source record↗

Design Tool for Estimating Chemical Hydrogen Storage System Characteristics for Light-Duty Fuel Cell Vehicles

The U.S. Department of Energy (DOE) has developed a vehicle framework model to simulate fuel cell-based light-duty vehicle operation for various hydrogen storage systems. This transient model simulates the performance of the storage system, fuel cell, and vehicle for comparison to DOE’s Technical Targets using four drive cycles/profiles. Chemical hydrogen storage models have been developed for the Framework model for both exothermic and endothermic materials. Despite the utility of such models, they require that material researchers input system design specifications that cannot be easily estimated. To address this challenge, a design tool has been developed that allows researchers to directly enter kinetic and thermodynamic chemical hydrogen storage material properties into a simple sizing module that then estimates the systems parameters required to run the storage system model. Additionally, this design tool can be used as a standalone executable file to estimate the storage system mass and volume outside of the framework model and compare it to the DOE Technical Targets. These models will be explained and exercised with existing hydrogen storage materials.

hydrogen storage materials, Fuel Cell Vehicles, de↗

TH-CD-206-12: Image-Based Motion Estimation for Plaque Visualization in Coronary Computed Tomography Angiography

Purpose: Visualization and quantification of coronary artery calcification and atherosclerotic plaque benefits from coronary artery motion (CAM) artifact elimination. This work applies a rigid linear motion model to a Volume of Interest (VoI) for estimating motion estimation and compensation of image degradation in Coronary Computed Tomography Angiography (CCTA). Methods: In both simulation and testbench experiments, translational CAM was generated by displacement of the imaging object (i.e. simulated coronary artery and explanted human heart) by ∼8 mm, approximating the motion of a main coronary branch. Rotation was assumed to be negligible. A motion degraded region containing a calcification was selected as the VoI. Local residual motion was assumed to be rigid and linear over the acquisition window, simulating motion observed during diastasis. The (negative) magnitude of the image gradient of the reconstructed VoI was chosen as the motion estimation objective and was minimized with Covariance Matrix Adaptation Evolution Strategy (CMAES). Results: Reconstruction incorporated the estimated CAM yielded signification recovery of fine calcification structures as well as reduced motion artifacts within the selected local region. The compensated reconstruction was further evaluated using two image similarity metrics, the structural similarity index (SSIM) and Root Mean Square Error (RMSE). At the calcification site, the compensated data achieved a 3% increase in SSIM and a 91.2% decrease in RMSE in comparison with the uncompensated reconstruction. Conclusion: Results demonstrate the feasibility of our image-based motion estimation method exploiting a local rigid linear model for CAM compensation. The method shows promising preliminary results for the application of such estimation in CCTA. Further work will involve motion estimation of complex motion corrupted patient data acquired from clinical CT scanner.

60 APPLIED LIFE SCIENCES↗

Current state and future projections of drying processes in the US food and pulp and paper sectors: Energy, economic, and environmental assessment

The pulp and paper (P/P) and food sectors are the third- and fifth-largest industrial energy consumers in the United States, with total on-site energy consumption of 2,039 TBtu and 1,144 TBtu, respectively. Thermal drying processes for moisture removal, which are energy-intensive, play a critical role in both industries. This study is the first to evaluate state- and national-level US drying energy demand for these sectors from 2020 to 2050. To complete this evaluation, we developed a thermodynamic modeling framework integrated with economic and environmental models to compute product-specific drying energy intensity and estimate the sector-specific costs and emissions profiles associated with drying operations. The model-predicted energy intensity was validated against the literature. Using current and projected annual production volumes in these sectors, we estimated total drying energy use. Results indicate that drying accounts for 22 % of total energy consumption in the P/P sector and 10 % in the food sector. The estimated annual energy cost (2020) to operate thermal dryers is $\$$919 M in the P/P sector and $\$$417 M in the food sector. Additionally, drying contributes to 25 % of total CO 2 e emissions in the P/P sector (including biogenic) and 15 % of emissions in the food sector. Regional performance shows that the Southern US is the leading energy consumer for P/P drying, whereas the Midwest leads in food drying. This study presents both potential solutions to enhance drying efficiency and barriers to implementation. Energy efficiency improvements, low-carbon fuels, and electrification are discussed as key pathways for reducing costs and optimizing industrial drying processes.

3E analysis↗

Subtask 1.1 – Advanced Characterization of Unconventional Oil and Gas Reservoirs to Enhance CO2 Storage Resource Estimates

The Energy & Environmental Research Center (EERC) developed methods and approaches for advanced characterization of unconventional oil and gas reservoirs to enhance CO2 storage resource estimates using samples collected from organic-rich shales of the Bakken Formation. The project was funded under U.S. Department of Energy (DOE) Cooperative Agreement No. DE-FE0024233, with nonfederal cost share provided by Hitachi High Technologies America. The project participants included the EERC, Hitachi High Technologies America, the team at the National Energy Technology Laboratory’s (NETL’s) computed tomography (CT) scanning laboratory in Morgantown, West Virginia, and several NETL staff who have worked on CO2 storage resource estimation. Investigations on the use of subsurface geologic formations for CO2 storage have been ongoing for many years. In the last decade, improvements in directional drilling technology and hydraulic fracturing of organic-rich shale and other tight rock formations have opened new opportunities for CO2 storage, with the potential cobenefit of enhancing gas and oil recovery. Gas-bearing, organic-rich shales have been of particular interest as a carbon storage target because their ability to store large volumes of CO2 is implied by their ability to store methane; however, estimating the CO2 storage resource of organic-rich shales and other unconventional reservoirs is challenging because of the lack of detailed geologic and petrophysical data that are needed to develop improved volumetric equations. According to NETL’s method for estimating regional-scale CO2 storage resources, the key factors that control CO2 storage in geologic formations include void space within induced fractures, natural fractures and matrix pores, as well as the occurrence and distribution of organic matter (OM) and clay minerals. A key goal of this effort was to use field emission scanning electron microscopy (FESEM) analysis of fine-grain unconventional reservoir rocks coupled with advanced image analysis and machine learning techniques to quantify the occurrence and distribution of porosity, OM, and clay minerals. By better understanding the potential migration pathways of CO2 within tight rocks and quantifying the occurrence of sorptive components (OM and clays) along potential flow pathways, methods to estimate the CO2 storage resource of such rocks can be refined. The Bakken petroleum system (Bakken) was used as the source of rock samples for this effort because it includes organic-rich shales with contrasting degrees of thermal maturity as well as tight, nonshale clastic rocks with varying clay mineral content. To better understand the key components of interest in tight rock samples, the EERC worked closely with Hitachi to optimize the FESEM settings necessary to resolve the sample mineralogy, porosity and OM content at micro- to nanoscale resolution. The EERC also worked closely with Hitachi to apply and improve its advanced mineral identification and characterization system (AMICS) software so that it is better suited for mineral identification at the high resolution necessary to resolve features of interest in fine-grained rocks such as organic-rich shales. The EERC also applied advanced image analysis software packages, including Ilastik, to classify and quantify the key components of interest within the samples using FESEM imagery. These advanced techniques allowed for quantification and differentiation of the void space within fractures, pores and OM, as well as a technique to determine and quantify the key mineral and OM components surrounding pore spaces. Being able to better quantify the void space and the potential accessibility of CO2 to sorptive minerals or OM within tight rocks allows for improved estimation of the CO2 storage resource potential of organic-rich shales and other unconventional formations. The information derived from the image analysis techniques was used to refine NETL’s method for estimating regional-scale CO2 storage resources in collaboration with staff from NETL. An additional goal of this work was to better understand rates of CO2 permeation, flow pathways, and sorption in organic-rich shales and nonshale reservoir rocks and to evaluate potential changes in CO2 flow over time. The EERC conducted long-duration (several weeks) CO2 flow-through testing at reservoir conditions using Bakken shale and nonshale samples to assess CO2 permeation rates and retention/sorption within core plugs to represent bulk matrix behavior. Conventional isothermal adsorption testing over a range of pressure conditions was also performed on pulverized Bakken samples to better assess potential maximum sorptive capacity of the samples. The results of this testing illustrated the higher sorptive capacity of the Bakken shales when compared to the sorptive capacity of the nonshale reservoir rock, which contains very little OM when compared to the Bakken shales. The results also illustrate the differences in permeation behavior as a function of fracture flow versus matrix flow. In an attempt to visualize CO2 permeation and flow behavior in Bakken shale samples over time, the EERC collaborated with NETL’s CT Lab in Morgantown, West Virginia. A variety of tests were performed at the NETL lab on unfractured and fractured samples of the Bakken shale. The unfractured Bakken shale samples did not result in CO2 migration within the samples over the testing duration; thus imagery of CO2 migration pathways within the matrix of Bakken samples was not obtained. The results of the flow-through testing within fractured samples illustrate changes in CO2 permeation rates over time, with a trend of decreasing CO2 flow over time. The results of this effort suggest that FESEM analysis coupled with advanced image analysis and machine learning techniques are effective in characterizing and quantifying the mineralogy, porosity, fractures, and OM content of 2-D imagery at the high resolution needed to identify features of interest in tight rocks. In addition, the image analysis techniques allow for identification and quantification of different porosity types as well as the key mineral components and/or OM that CO2 would contact in the reservoir as it migrates through pores and fracture networks. This is key information needed to better estimate the CO2 storage resource potential of organic-rich shales and other tight reservoirs. Additional work is needed to determine the sorptive capacity of the individual sample components (OM and clay minerals) at reservoir conditions so that this approach and the refined CO2 storage resource equation can be applied to other shale and unconventional reservoirs. In addition, subsequent work should be performed to better understand the number of images needed to use stereological approaches to estimate the volume of sample components using 2-D imagery, to explore machine learning techniques to differentiate between induced vs. natural fractures, and to estimate the connectivity of pores and fracture networks using 2-D imagery. Being able to derive this information at the microscale would facilitate scaling up of the sample characterization data used for reservoir modeling and simulation to estimate the CO2 storage resource potential of tight formations.

03 NATURAL GAS↗

Application of the Redlich-Kister expansion for estimating the density of molten fluoride psuedo-ternary salt systems of nuclear industry interest

The development of next-generation molten salt reactors relies on accurate knowledge of the thermophysical properties of the candidate coolant and fueled molten salts. These thermophysical properties include density, viscosity, thermal conductivity, and heat capacity. Because of difficulties in measuring thermophysical properties of molten salts, there are many gaps in the current state of thermophysical property knowledge of these salts, particularly those that contain actinides or beryllium. Therefore, leveraging modeling techniques to estimate unknown molten salt thermophysical properties and guide future experimental measurements has high value for the nuclear industry. Here, the densities of molten fluoride pseudo-ternary salt systems, which are of interest to the nuclear industry, were estimated using Redlich-Kister expansion and Muggianu interpolation techniques. The pseudo-ternary systems considered for estimation in this study were NaF-LiF-ZrF 4 , LiF-BeF 2 -ZrF 4 , LiF-BeF 2 -ThF 4 , NaF-LiF-BeF 2 , NaF-KF-BeF 2 , NaF-ZrF 4 -UF 4 , and NaF-BeF 2 -UF 4 . This Redlich-Kister estimation approach accounts for nonideal mixing behavior based on pseudo-binary subsystem interaction parameters determined from experimentally measured pseudo-binary system density data sets. The Redlich-Kister estimation was compared with the method of additive molar volumes, which assumes ideal mixing. Additionally, the Redlich-Kister approach was used to determine previously unknown binary and ternary interaction parameters based on experimentally measured density data sets for select pseudo-ternary salt systems. The results of this study show improvement in density estimation using the Redlich-Kister approach for all systems considered compared with estimation by additive molar volumes. Furthermore, this analysis allowed for the estimation of nonideal density behavior in experimentally unstudied ZrF 4 -UF 4 and BeF 2 -UF 4 , as well as the quantification of ternary interaction in NaF-LiF-ZrF 4 , NaF-BeF 2 -UF 4 , and NaF-ZrF 4 -UF 4 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Estimates of Lava Eruption Rates at Alba Patera, Mars

The Martian volcanic complex Alba Patera exhibits a suite of well-defined, long and relatively narrow lava flows qualitatively resembling those found in Hawaii. Even without any information on the duration of the Martian flows, eruption rates (total volume discharge/duration of the extrusion) estimates are implied by the physical dimensions of the flows and the likely conjecture that Stephan-Boltzmann radiation is the dominating thermal loss mechanism. The ten flows in this analysis emanate radially from the central vent and were recently measured in length, plan areas, and average thicknesses by shadow measurement techniques. The dimensions of interest are shown. Although perhaps morphologically congruent to certain Hawaiian flows, the dramatically expanded physical dimensions of the Martian flows argues for some markedly distinct differences in lava flow composition for eruption characteristics.

Baloga, S. M.↗

Cytoplasm-to-myonucleus ratios following microgravity

The cytoplasmic volume-to-myonucleus ratio in the tibialis anterior and gastrocnemius muscles of juvenile rats after 5.4 days of microgravity was studied. Three groups of rats (n = 8 each) were used. The experimental group (space rats) was flown aboard the space shuttle Discovery (NASA, STS-48), while two ground-based groups, one hindlimb suspended (suspended rats), one non-suspended (control), served as controls. Single fibre analysis revealed a significant decrease in cross-sectional area (microns2) in the gastrocnemius for both the space and the suspended rats; in the tibialis anterior only the suspended rats showed a significant decrease. Myonuclei counts (myonuclei per mm) in both the tibialis anterior and gastrocnemius were significantly increased in the space rats but not in the suspended rats. The mean myonuclear volume (individual nuclei: microns3) in tibialis anterior fibres from the space rats, and in gastrocnemius fibres from both the space and the suspended rats, was significantly lower than that in the respective control group. Estimation of the total myonuclear volume (microns3 per.mm), however, revealed no significant differences between the three groups in either the tibialis anterior or gastrocnemius. The described changes in the cross-sectional area and myonuclei numbers resulted in significant decreases in the cytoplasmic volume-to-myonucleus ratio (microns3 x 10(3)) in both muscles and for both space and suspended rats (tibialis anterior; 15.6 +/- 0.6 (space), 17.2 +/- 1.0 (suspended), 20.8 +/- 0.9 (control): gastrocnemius; 13.4 +/- 0.4 (space) and 14.9 +/- 1.1 (suspended) versus 18.1 +/- 1.1 (control)). These results indicate that even short periods of unweighting due to microgravity or limb suspension result in changes in skeletal muscle fibres which lead to significant decreases in the cytoplasmic volume-to-myonucleus ratio.

NASA Discipline Cell Biology↗

TH-EF-BRA-08: A Novel Technique for Estimating Volumetric Cine MRI (VC-MRI) From Multi-Slice Sparsely Sampled Cine Images Using Motion Modeling and Free Form Deformation

Purpose: To develop a technique to estimate on-board VC-MRI using multi-slice sparsely-sampled cine images, patient prior 4D-MRI, motion-modeling and free-form deformation for real-time 3D target verification of lung radiotherapy. Methods: A previous method has been developed to generate on-board VC-MRI by deforming prior MRI images based on a motion model(MM) extracted from prior 4D-MRI and a single-slice on-board 2D-cine image. In this study, free-form deformation(FD) was introduced to correct for errors in the MM when large anatomical changes exist. Multiple-slice sparsely-sampled on-board 2D-cine images located within the target are used to improve both the estimation accuracy and temporal resolution of VC-MRI. The on-board 2D-cine MRIs are acquired at 20–30frames/s by sampling only 10% of the k-space on Cartesian grid, with 85% of that taken at the central k-space. The method was evaluated using XCAT(computerized patient model) simulation of lung cancer patients with various anatomical and respirational changes from prior 4D-MRI to onboard volume. The accuracy was evaluated using Volume-Percent-Difference(VPD) and Center-of-Mass-Shift(COMS) of the estimated tumor volume. Effects of region-of-interest(ROI) selection, 2D-cine slice orientation, slice number and slice location on the estimation accuracy were evaluated. Results: VCMRI estimated using 10 sparsely-sampled sagittal 2D-cine MRIs achieved VPD/COMS of 9.07±3.54%/0.45±0.53mm among all scenarios based on estimation with ROI-MM-ROI-FD. The FD optimization improved estimation significantly for scenarios with anatomical changes. Using ROI-FD achieved better estimation than global-FD. Changing the multi-slice orientation to axial, coronal, and axial/sagittal orthogonal reduced the accuracy of VCMRI to VPD/COMS of 19.47±15.74%/1.57±2.54mm, 20.70±9.97%/2.34±0.92mm, and 16.02±13.79%/0.60±0.82mm, respectively. Reducing the number of cines to 8 enhanced temporal resolution of VC-MRI by 25% while maintaining the estimation accuracy. Estimation using slices sampled uniformly through the tumor achieved better accuracy than slices sampled non-uniformly. Conclusions: Preliminary studies showed that it is feasible to generate VC-MRI from multi-slice sparsely-sampled 2D-cine images for real-time 3D-target verification. This work was supported by the National Institutes of Health under Grant No. R01-CA184173 and a research grant from Varian Medical Systems.

60 APPLIED LIFE SCIENCES↗

TU-AB-BRA-09: A Novel Method of Generating Ultrafast Volumetric Cine MRI (VC-MRI) Using Prior 4D-MRI and On-Board Phase-Skipped Encoding Acquisition for Radiotherapy Target Localization

Purpose: To develop a technique generating ultrafast on-board VC-MRI using prior 4D-MRI and on-board phase-skipped encoding k-space acquisition for real-time 3D target tracking of liver and lung radiotherapy. Methods: The end-of-expiration (EOE) volume in 4D-MRI acquired during the simulation was selected as the prior volume. 3 major respiratory deformation patterns were extracted through the principal component analysis of the deformation field maps (DFMs) generated between EOE and all other phases. The on-board VC-MRI at each instant was considered as a deformation of the prior volume, and the deformation was modeled as a linear combination of the extracted 3 major deformation patterns. To solve the weighting coefficients of the 3 major patterns, a 2D slice was extracted from VC-MRI volume to match with the 2D on-board sampling data, which was generated by 8-fold phase skipped-encoding k-space acquisition (i.e., sample 1 phase-encoding line out of every 8 lines) to achieve an ultrafast 16–24 volumes/s frame rate. The method was evaluated using XCAT digital phantom to simulate lung cancer patients. The 3D volume of end-ofinhalation (EOI) phase at the treatment day was used as ground-truth onboard VC-MRI with simulated changes in 1) breathing amplitude and 2) breathing amplitude/phase change from the simulation day. A liver cancer patient case was evaluated for in-vivo feasibility demonstration. Results: The comparison between ground truth and estimated on-board VC-MRI shows good agreements. In XCAT study with changed breathing amplitude, the volume-percent-difference(VPD) between ground-truth and estimated tumor volumes at EOI was 6.28% and the Center-of-Mass-Shift(COMS) was 0.82mm; with changed breathing amplitude and phase, the VPD was 8.50% and the COMS was 0.54mm. The study of liver patient case also demonstrated a promising in vivo feasibility of the proposed method Conclusion: Preliminary results suggest the feasibility to estimate ultrafast VC-MRI for on-board target localization with phase skipped-encoding k-space acquisition. Research grant from NIH R01-184173.

60 APPLIED LIFE SCIENCES↗