Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “motion”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Ground Motion Models (GMMs) Improvements Using Earthquake Simulations on High Performance Computers

A computationally efficient simulation platform was developed that can provide representative synthetic ground motions from crustal earthquakes in the Stable Continental Regions of Central and Eastern US (CEUS), using 3D modeling and high-performance computing. The main objective was to use synthetic ground motion to provide constrains to refinements of exiting ergodic Ground Motion Models (GMMs), for large magnitude earthquakes and near-fault distances, for which these models are less reliable. Physics-based broadband (0-5Hz) ground motion simulations were used to estimate the near-fault ground motion amplitudes and within event and between-event variabilities associated with fault rupture characteristics. As part of a strategy for selecting a reginal velocity model and validation of developed rupture modeling technique, ground motions from the moment magnitude Mw5.0 November 7, 2016, Cushing Oklahoma, and Mw5.8 September 3, 2016, Pawnee Oklahoma earthquakes were simulated. In our simulations we used a 3D regional velocity model that was based on Saikia’s 1D velocity model. Saikia’s model demonstrated better performance in modelling high frequency regional wave propagation for CEUS region. The proposed 3D model includes lateral variations added to the 1D background model using the stochastic scheme of Pitarka and Mellors. Comparisons of the simulations with recordings of both earthquakes demonstrated the reliability of our deterministic simulation approach while emphasizing the importance of including small-scale variability in the regional velocity model needed to reproduce the observed high-frequency wave scattering effects. As part of validation analysis, comparisons with different GMMs for a Mw6.5 earthquake in the CESUS region resulted in a very good match between the simulated and empirical ground motion models. Initial investigations of within-event and between-event ground motion variabilities for Mw6.5 scenario earthquakes on a strike-slip fault, suggest that they are strongly related to spatial slip and slip rate variations, average rupture velocity, rupture area and rupture initiation location. For certain scenarios we found that the ground motion variability observed at near-fault distances (< 5 km) also persists at longer distances. Regardless of the rupture scenario, the simulated ground motion tends to fully saturate at short distances and for all periods. The near-fault saturation has to do with the attenuation of waves propagating along the fault and local rupture radiation pattern that also contribute to stronger ground motion variation at such distances. Analysis of effects of rupture initiation location suggest that the peak ground motion (PGV) and spectral acceleration (SA) can be quite variable due to rupture directivity effects. Such effects are stronger at periods longer than 1s. The effect of the 1D velocity models and surface topography on simulated ground motion were investigated by comparing three component synthetic seismograms computed at selected sites. Effect of surface topography was considered using the ratio between spectral accelerations simulated for two 1D models with flat surface topography and realistic model with surface topography. Overall, the topography slightly amplifies (by ~30%) the ground motion amplitude in the frequency range 1-3Hz. The effect of topography is more visible in the surface and coda waves portion of the seismograms.

58 GEOSCIENCES↗
NASA NTRS1984

Oculometric indices of simulator and aircraft motion

The effects on eye scan behavior of both simulator and aircraft motion and sensitivity of an oculometric measure to motion effects was demonstrated. It was found that fixation time is sensitive to motion effects. Differences between simulator motion and no motion conditions during a series of simulated ILS approaches were studied. The mean fixation time for the no motion condition was found to be significantly longer than for the motion conditions. Eye scan parameters based on data collected in flight, and in fixed base simulation were investigated. Motion effects were evident when the subject was viewing a display supplying attitude and flight path information. The nature of the information provided by motion was examined. The mean fixation times for the no motion condition were significantly longer than for either motion condition, while the two motion conditions did not differ. It is shown that motion serves an alerting function, providing a cue or clue to the pilot that something happened. It is suggested that simulation without motion cues may represent an understatement of the true capacity of the pilot.

Comstock, J. R.↗
NASA NTRS1992

Role of orientation reference selection in motion sickness

The overall objective of this proposal is to understand the relationship between human orientation control and motion sickness susceptibility. Three areas related to orientation control will be investigated. These three areas are (1) reflexes associated with the control of eye movements and posture, (2) the perception of body rotation and position with respect to gravity, and (3) the strategies used to resolve sensory conflict situations which arise when different sensory systems provide orientation cues which are not consistent with one another or with previous experience. Of particular interest is the possibility that a subject may be able to ignore an inaccurate sensory modality in favor of one or more other sensory modalities which do provide accurate orientation reference information. We refer to this process as sensory selection. This proposal will attempt to quantify subjects' sensory selection abilities and determine if this ability confers some immunity to the development of motion sickness symptoms. Measurements of reflexes, motion perception, sensory selection abilities, and motion sickness susceptibility will concentrate on pitch and roll motions since these seem most relevant to the space motion sickness problem. Vestibulo-ocular (VOR) and oculomotor reflexes will be measured using a unique two-axis rotation device developed in our laboratory over the last seven years. Posture control reflexes will be measured using a movable posture platform capable of independently altering proprioceptive and visual orientation cues. Motion perception will be quantified using closed loop feedback technique developed by Zacharias and Young (Exp Brain Res, 1981). This technique requires a subject to null out motions induced by the experimenter while being exposed to various confounding sensory orientation cues. A subject's sensory selection abilities will be measured by the magnitude and timing of his reactions to changes in sensory environments. Motion sickness susceptibility will be measured by the time required to induce characteristic changes in the pattern of electrogastrogram recordings while exposed to various sensory environments during posture and motion perception tests. The results of this work are relevant to NASA's interest in understanding the etiology of space motion sickness. If any of the reflex, perceptual, or sensory selection abilities of subjects are found to correlate with motion sickness susceptibility, this work may be an important step in suggesting a method of predicting motion sickness susceptibility. If sensory selection can provide a means to avoid sensory conflict, then further work may lead to training programs which could enhance a subject's sensory selection ability and therefore minimize motion sickness susceptibility.

Peterka, Robert J.↗
NASA NTRS1991

Relationship of area postrema to three putative measures of motion sickness

Although the rat has an incomplete emetic reflex, several species-specific responses to motion were proposed as measures of 'motion sickness' in rats. The purpose was to determine the dependence of these responses on one of several neural structures known to be essential to motion-induced vomiting in species with a complete emetic reflex. The Area Postrema (AP) was shown to play an important role in the production of motion sickness in vomiting species. The effects of thermo-cautery ablations of the AP on three different responses supposedly reflecting motion sickness in the rat were compared: conditioned taste aversion (CTA); drinking suppression; and fecal boli. Efficacy of the ablations was determined by subjecting ablated, sham-operated, and unoperated control animals to a CTA test which is known to require a functional AP. Animals with AP ablations failed to form CTA when 0.15 M LiCl was paired with a 10 percent sucrose solution, while sham-operated control subjects conditioned as well as the unoperated control subjects. The extent of the ablations was evaluated histologically at the end of the experiment. To determine the effects of the ablations on the measures of motion sickness, all animals were subjected to rotation for 30 min or 90 min on a platform displaced 20 deg from earth horizontal. Results indicate that ablation of AP in the rat has no effect on the formation of CTA to a 4 percent solution of cider paired with motion, on the suppression of drinking immediately after exposure to motion, or on the frequency of fecal boli during exposure to motion. This failure of AP ablations to eliminate the effects of motion on any of these responses discourages their use as equivalents of motion-induced vomiting. The appropriateness of other suggested measures, e.g., pica, remains untested but the dependence of such measures on stimulation more severe than commonly used in motion sickness research and the absence of a demonstration of their dependence on neural structures essential to motion sickness in vomiting species, suggest caution in the use of such responses. Further, until more is known about the neural structures underlying these putative measures, the rat will remain a questionable subject in which to study motion sickness.

R Sutton↗
NASA NTRS2016

Effects of Different Heave Motion Components on Pilot Pitch Control Behavior

The study described in this paper had two objectives. The first objective was to investigate if a different weighting of heave motion components decomposed at the center of gravity, allowing for a higher fidelity of individual components, would result in pilot manual pitch control behavior and performance closer to that observed with full aircraft motion. The second objective was to investigate if decomposing the heave components at the aircraft's instantaneous center of rotation rather than at the center of gravity could result in additional improvements in heave motion fidelity. Twenty-one general aviation pilots performed a pitch attitude control task in an experiment conducted on the Vertical Motion Simulator at NASA Ames under different hexapod motion conditions. The large motion capability of the Vertical Motion Simulator also allowed for a full aircraft motion condition, which served as a baseline. The controlled dynamics were of a transport category aircraft trimmed close to the stall point. When the ratio of center of gravity pitch heave to center of gravity heave increased in the hexapod motion conditions, pilot manual control behavior and performance became increasingly more similar to what is observed with full aircraft motion. Pilot visual and motion gains significantly increased, while the visual lead time constant decreased. The pilot visual and motion time delays remained approximately constant and decreased, respectively. The neuromuscular damping and frequency both decreased, with their values more similar to what is observed with real aircraft motion when there was an equal weighting of the heave of the center of gravity and heave due to rotations about the center of gravity. In terms of open- loop performance, the disturbance and target crossover frequency increased and decreased, respectively, and their corresponding phase margins remained constant and increased, respectively. The decomposition point of the heave components only had limited effects on pilot manual control behavior and performance.

manual control↗
DOE OSTI2015

Magnetic Resonance Imaging Assessment of Spinal Cord and Cauda Equina Motion in Supine Patients With Spinal Metastases Planned for Spine Stereotactic Body Radiation Therapy

Purpose: To assess motion of the spinal cord and cauda equina, which are critical neural tissues (CNT), which is important when evaluating the planning organ-at-risk margin required for stereotactic body radiation therapy. Methods and Materials: We analyzed CNT motion in 65 patients with spinal metastases (11 cervical, 39 thoracic, and 24 lumbar spinal segments) in the supine position using dynamic axial and sagittal magnetic resonance imaging (dMRI, 3T Verio, Siemens) over a 137-second interval. Motion was segregated according to physiologic cardiorespiratory oscillatory motion (characterized by the average root mean square deviation) and random bulk shifts associated with gross patient motion (characterized by the range). Displacement was evaluated in the anteroposterior (AP), lateral (LR), and superior-inferior (SI) directions by use of a correlation coefficient template matching algorithm, with quantification of random motion measure error over 3 separate trials. Statistical significance was defined according to P<.05. Results: In the AP, LR, and SI directions, significant oscillatory motion was observed in 39.2%, 35.1%, and 10.8% of spinal segments, respectively, and significant bulk motions in all cases. The median oscillatory CNT motions in the AP, LR, and SI directions were 0.16 mm, 0.17 mm, and 0.44 mm, respectively, and the maximal statistically significant oscillatory motions were 0.39 mm, 0.41 mm, and 0.77 mm, respectively. The median bulk displacements in the AP, LR, and SI directions were 0.51 mm, 0.59 mm, and 0.66 mm, and the maximal statistically significant displacements were 2.21 mm, 2.87 mm, and 3.90 mm, respectively. In the AP, LR, and SI directions, bulk displacements were greater than 1.5 mm in 5.4%, 9.0%, and 14.9% of spinal segments, respectively. No significant differences in axial motion were observed according to cord level or cauda equina. Conclusions: Oscillatory CNT motion was observed to be relatively minor. Our results support the importance of controlling bulk patient motion and the practice of applying a planning organ-at-risk margin.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

MR-based motion correction for PET imaging using wired active MR microcoils in simultaneous PET-MR: Phantom study

Purpose: Artifacts caused by head motion present a major challenge in brain positron emission tomography (PET) imaging. The authors investigated the feasibility of using wired active MR microcoils to track head motion and incorporate the measured rigid motion fields into iterative PET reconstruction. Methods: Several wired active MR microcoils and a dedicated MR coil-tracking sequence were developed. The microcoils were attached to the outer surface of an anthropomorphic{sup 18}F-filled Hoffman phantom to mimic a brain PET scan. Complex rotation/translation motion of the phantom was induced by a balloon, which was connected to a ventilator. PET list-mode and MR tracking data were acquired simultaneously on a PET-MR scanner. The acquired dynamic PET data were reconstructed iteratively with and without motion correction. Additionally, static phantom data were acquired and used as the gold standard. Results: Motion artifacts in PET images were effectively removed by wired active MR microcoil based motion correction. Motion correction yielded an activity concentration bias ranging from −0.6% to 3.4% as compared to a bias ranging from −25.0% to 16.6% if no motion correction was applied. The contrast recovery values were improved by 37%–156% with motion correction as compared to no motion correction. The image correlation (mean ± standard deviation) between the motion corrected (uncorrected) images of 20 independent noise realizations and static reference was R{sup 2} = 0.978 ± 0.007 (0.588 ± 0.010, respectively). Conclusions: Wired active MR microcoil based motion correction significantly improves brain PET quantitative accuracy and image contrast.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Evaluation of respiratory and cardiac motion correction schemes in dual gated PET/CT cardiac imaging

Purpose: Cardiac imaging suffers from both respiratory and cardiac motion. One of the proposed solutions involves double gated acquisitions. Although such an approach may lead to both respiratory and cardiac motion compensation there are issues associated with (a) the combination of data from cardiac and respiratory motion bins, and (b) poor statistical quality images as a result of using only part of the acquired data. The main objective of this work was to evaluate different schemes of combining binned data in order to identify the best strategy to reconstruct motion free cardiac images from dual gated positron emission tomography (PET) acquisitions. Methods: A digital phantom study as well as seven human studies were used in this evaluation. PET data were acquired in list mode (LM). A real-time position management system and an electrocardiogram device were used to provide the respiratory and cardiac motion triggers registered within the LM file. Acquired data were subsequently binned considering four and six cardiac gates, or the diastole only in combination with eight respiratory amplitude gates. PET images were corrected for attenuation, but no randoms nor scatter corrections were included. Reconstructed images from each of the bins considered above were subsequently used in combination with an affine or an elastic registration algorithm to derive transformation parameters allowing the combination of all acquired data in a particular position in the cardiac and respiratory cycles. Images were assessed in terms of signal-to-noise ratio (SNR), contrast, image profile, coefficient-of-variation (COV), and relative difference of the recovered activity concentration. Results: Regardless of the considered motion compensation strategy, the nonrigid motion model performed better than the affine model, leading to higher SNR and contrast combined with a lower COV. Nevertheless, when compensating for respiration only, no statistically significant differences were observed in the performance of the two motion models considered. Superior image SNR and contrast were seen using the affine respiratory motion model in combination with the diastole cardiac bin in comparison to the use of the whole cardiac cycle. In contrast, when simultaneously correcting for cardiac beating and respiration, the elastic respiratory motion model outperformed the affine model. In this context, four cardiac bins associated with eight respiratory amplitude bins seemed to be adequate. Conclusions: Considering the compensation of respiratory motion effects only, both affine and elastic based approaches led to an accurate resizing and positioning of the myocardium. The use of the diastolic phase combined with an affine model based respiratory motion correction may therefore be a simple approach leading to significant quality improvements in cardiac PET imaging. However, the best performance was obtained with the combined correction for both cardiac and respiratory movements considering all the dual-gated bins independently through the use of an elastic model based motion compensation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Adaptation of the modified Bouc–Wen model to compensate for hysteresis in respiratory motion for the list-mode binning of cardiac SPECT and PET acquisitions: Testing using MRI

Purpose: Binning list-mode acquisitions as a function of a surrogate signal related to respiration has been employed to reduce the impact of respiratory motion on image quality in cardiac emission tomography (SPECT and PET). Inherent in amplitude binning is the assumption that there is a monotonic relationship between the amplitude of the surrogate signal and respiratory motion of the heart. This assumption is not valid in the presence of hysteresis when heart motion exhibits a different relationship with the surrogate during inspiration and expiration. The purpose of this study was to investigate the novel approach of using the Bouc–Wen (BW) model to provide a signal accounting for hysteresis when binning list-mode data with the goal of thereby improving motion correction. The study is based on the authors’ previous observations that hysteresis between chest and abdomen markers was indicative of hysteresis between abdomen markers and the internal motion of the heart. Methods: In 19 healthy volunteers, they determined the internal motion of the heart and diaphragm in the superior–inferior direction during free breathing using MRI navigators. A visual tracking system (VTS) synchronized with MRI acquisition tracked the anterior–posterior motions of external markers placed on the chest and abdomen. These data were employed to develop and test the Bouc–Wen model by inputting the VTS derived chest and abdomen motions into it and using the resulting output signals as surrogates for cardiac motion. The data of the volunteers were divided into training and testing sets. The training set was used to obtain initial values for the model parameters for all of the volunteers in the set, and for set members based on whether they were or were not classified as exhibiting hysteresis using a metric derived from the markers. These initial parameters were then employed with the testing set to estimate output signals. Pearson’s linear correlation coefficient between the abdomen, chest, average of chest and abdomen markers, and Bouc–Wen derived signals versus the true internal motion of the heart from MRI was used to judge the signals match to the heart motion. Results: The results show that the Bouc–Wen model generated signals demonstrated strong correlation with the heart motion. This correlation was slightly larger on average than that of the external surrogate signals derived from the abdomen marker, and average of the abdomen and chest markers, but was not statistically significantly different from them. Conclusions: The results suggest that the proposed model has the potential to be a unified framework for modeling hysteresis in respiratory motion in cardiac perfusion studies and beyond.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Accelerated acquisition of tagged MRI for cardiac motion correction in simultaneous PET-MR: Phantom and patient studies

Purpose: Degradation of image quality caused by cardiac and respiratory motions hampers the diagnostic quality of cardiac PET. It has been shown that improved diagnostic accuracy of myocardial defect can be achieved by tagged MR (tMR) based PET motion correction using simultaneous PET-MR. However, one major hurdle for the adoption of tMR-based PET motion correction in the PET-MR routine is the long acquisition time needed for the collection of fully sampled tMR data. In this work, the authors propose an accelerated tMR acquisition strategy using parallel imaging and/or compressed sensing and assess the impact on the tMR-based motion corrected PET using phantom and patient data. Methods: Fully sampled tMR data were acquired simultaneously with PET list-mode data on two simultaneous PET-MR scanners for a cardiac phantom and a patient. Parallel imaging and compressed sensing were retrospectively performed by GRAPPA and kt-FOCUSS algorithms with various acceleration factors. Motion fields were estimated using nonrigid B-spline image registration from both the accelerated and fully sampled tMR images. The motion fields were incorporated into a motion corrected ordered subset expectation maximization reconstruction algorithm with motion-dependent attenuation correction. Results: Although tMR acceleration introduced image artifacts into the tMR images for both phantom and patient data, motion corrected PET images yielded similar image quality as those obtained using the fully sampled tMR images for low to moderate acceleration factors (<4). Quantitative analysis of myocardial defect contrast over ten independent noise realizations showed similar results. It was further observed that although the image quality of the motion corrected PET images deteriorates for high acceleration factors, the images were still superior to the images reconstructed without motion correction. Conclusions: Accelerated tMR images obtained with more than 4 times acceleration can still provide relatively accurate motion fields and yield tMR-based motion corrected PET images with similar image quality as those reconstructed using fully sampled tMR data. The reduction of tMR acquisition time makes it more compatible with routine clinical cardiac PET-MR studies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

MagicPlate-512: A 2D silicon detector array for quality assurance of stereotactic motion adaptive radiotherapy

Purpose: Spatial and temporal resolutions are two of the most important features for quality assurance instrumentation of motion adaptive radiotherapy modalities. The goal of this work is to characterize the performance of the 2D high spatial resolution monolithic silicon diode array named “MagicPlate-512” for quality assurance of stereotactic body radiation therapy (SBRT) and stereotactic radiosurgery (SRS) combined with a dynamic multileaf collimator (MLC) tracking technique for motion compensation. Methods: MagicPlate-512 is used in combination with the movable platform HexaMotion and a research version of radiofrequency tracking system Calypso driving MLC tracking software. The authors reconstruct 2D dose distributions of small field square beams in three modalities: in static conditions, mimicking the temporal movement pattern of a lung tumor and tracking the moving target while the MLC compensates almost instantaneously for the tumor displacement. Use of Calypso in combination with MagicPlate-512 requires a proper radiofrequency interference shielding. Impact of the shielding on dosimetry has been simulated by GEANT4 and verified experimentally. Temporal and spatial resolutions of the dosimetry system allow also for accurate verification of segments of complex stereotactic radiotherapy plans with identification of the instant and location where a certain dose is delivered. This feature allows for retrospective temporal reconstruction of the delivery process and easy identification of error in the tracking or the multileaf collimator driving systems. A sliding MLC wedge combined with the lung motion pattern has been measured. The ability of the MagicPlate-512 (MP512) in 2D dose mapping in all three modes of operation was benchmarked by EBT3 film. Results: Full width at half maximum and penumbra of the moving and stationary dose profiles measured by EBT3 film and MagicPlate-512 confirm that motion has a significant impact on the dose distribution. Motion, no motion, and motion with MLC tracking profiles agreed within 1 and 0.4 mm, respectively, for all field sizes tested. Use of electromagnetic tracking system generates a fluctuation of the detector baseline up to 10% of the full scale signal requiring a proper shielding strategy. MagicPlate-512 is also able to reconstruct the dose variation pulse-by-pulse in each pixel of the detector. An analysis of the dose transients with motion and motion with tracking shows that the tracking feedback algorithm used for this experiment can compensate effectively only the effect of the slower transient components. The fast changing components of the organ motion can contribute only to discrepancy of the order of 15% in penumbral region while the slower components can change the dose profile up to 75% of the expected dose. Conclusions: MagicPlate-512 is shown to be, potentially, a valid alternative to film or 2D ionizing chambers for quality assurance dosimetry in SRS or SBRT. Its high spatial and temporal resolutions allow for accurate reconstruction of the profile in any conditions with motion and with tracking of the motion. It shows excellent performance to reconstruct the dose deposition in real time or retrospectively as a function of time for detailed analysis of the effect of motion in a specific pixel or area of interest.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-J-60: Evaluation of Temporal Lag in Radiotherapy Gating for Tumor Motion Trajectories

Purpose: Evaluating timing differences between LINAC beam ON/OFF and the estimation of tumor positioning using gating systems is essential for establishing confidence when treating with gating during radiotherapy, and is an annual requirement of TG-142. Temporal discrepancies between the trajectories of external marker surrogates and beam delivery may vary depending upon the type of external marker motion, which is quantified in this work for several trajectories. Methods: A precise robotic 3D motion stage performed several trajectories typically used for gating phantoms, including sinusoidal and Lujan-type motion; a commercial respiratory motion simulator was also employed. The true motions were monitored using variable resistors. The beam ON/OFF was controlled separately by two RPM (Varian) systems, an integrated version delivered by a Varian Truebeam LINAC and version 1.6 delivered by a Varian Trilogy, and measured using a diode. The resistor and diode signals were read by a multichannel digital oscilloscope, and timing differences between beam ON/OFF as measured by the diode and the phantom motion were determined using a peak detection algorithm. Results: Timing differences between beam ON/OFF and 3D stage motion peaks (diode—true motion timing) were computed to be 79.4 & 57.7ms for sinusoidal motion and 109.1 & 63.6ms for Lujan-type motion on the Truebeam LINAC, for beam ON and OFF, respectively. Timing differences for the Trilogy LINAC were 34.4 & 55.2ms for the sinusoidal motion and 29.0 & 26.3ms for the Lujan-type motion, for beam ON and OFF, respectively. With the commercial motion simulator, the timing differences were found to be −9ms and −78ms for beam ON/OFF, respectively, with the Truebeam, and −97.6ms and −60.9ms for beam ON/OFF, respectively, with the Trilogy. Conclusion: Setup-dependent temporal lags were found using this methodology. These discrepancies have the potential to influence quality assurance on gating systems and ultimately the treatment dose during radiotherapy. NIH National Institute of Biomedical Imaging and Bioengineering grant T32 EB002103-21.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Characterizing spatiotemporal information loss in sparse-sampling-based dynamic MRI for monitoring respiration-induced tumor motion in radiotherapy

Purpose: Sparse-sampling and reconstruction techniques represent an attractive strategy to achieve faster image acquisition speeds, while maintaining adequate spatial resolution and signal-to-noise ratio in rapid magnetic resonance imaging (MRI). The authors investigate the use of one such sequence, broad-use linear acquisition speed-up technique (k-t BLAST) in monitoring tumor motion for thoracic and abdominal radiotherapy and examine the potential trade-off between increased sparsification (to increase imaging speed) and the potential loss of “true” information due to greater reliance on a priori information. Methods: Lung tumor motion trajectories in the superior–inferior direction, previously recorded from ten lung cancer patients, were replayed using a motion phantom module driven by an MRI-compatible motion platform. Eppendorf test tubes filled with water which serve as fiducial markers were placed in the phantom. The modeled rigid and deformable motions were collected in a coronal image slice using balanced fast field echo in conjunction with k-t BLAST. Root mean square (RMS) error was used as a metric of spatial accuracy as measured trajectories were compared to input data. The loss of spatial information was characterized for progressively increasing acceleration factor from 1 to 16; the resultant sampling frequency was increased approximately from 2.5 to 19 Hz when the principal direction of the motion was set along frequency encoding direction. In addition to the phantom study, respiration-induced tumor motions were captured from two patients (kidney tumor and lung tumor) at 13 Hz over 49 s to demonstrate the impact of high speed motion monitoring over multiple breathing cycles. For each subject, the authors compared the tumor centroid trajectory as well as the deformable motion during free breathing. Results: In the rigid and deformable phantom studies, the RMS error of target tracking at the acquisition speed of 19 Hz was approximately 0.3–0.4 mm, which was smaller than the reconstructed pixel resolution of 0.67 mm. In the patient study, the dynamic 2D MRI enabled the monitoring of cycle-to-cycle respiratory variability present in the tumor position. It was seen that the range of centroid motion as well as the area covered due to target motion during each individual respiratory cycle was underestimated compared to the entire motion range observed over multiple breathing cycles. Conclusions: The authors’ initial results demonstrate that sparse-sampling- and reconstruction-based dynamic MRI can be used to achieve adequate image acquisition speeds without significant information loss for the task of radiotherapy guidance. Such monitoring can yield spatial and temporal information superior to conventional offline and online motion capture methods used in thoracic and abdominal radiotherapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

An integrated model-driven method for in-treatment upper airway motion tracking using cine MRI in head and neck radiation therapy

Purpose: For the first time, MRI-guided radiation therapy systems can acquire cine images to dynamically monitor in-treatment internal organ motion. However, the complex head and neck (H&N) structures and low-contrast/resolution of on-board cine MRI images make automatic motion tracking a very challenging task. In this study, the authors proposed an integrated model-driven method to automatically track the in-treatment motion of the H&N upper airway, a complex and highly deformable region wherein internal motion often occurs in an either voluntary or involuntary manner, from cine MRI images for the analysis of H&N motion patterns. Methods: Considering the complex H&N structures and ensuring automatic and robust upper airway motion tracking, the authors firstly built a set of linked statistical shapes (including face, face-jaw, and face-jaw-palate) using principal component analysis from clinically approved contours delineated on a set of training data. The linked statistical shapes integrate explicit landmarks and implicit shape representation. Then, a hierarchical model-fitting algorithm was developed to align the linked shapes on the first image frame of a to-be-tracked cine sequence and to localize the upper airway region. Finally, a multifeature level set contour propagation scheme was performed to identify the upper airway shape change, frame-by-frame, on the entire image sequence. The multifeature fitting energy, including the information of intensity variations, edge saliency, curve geometry, and temporal shape continuity, was minimized to capture the details of moving airway boundaries. Sagittal cine MR image sequences acquired from three H&N cancer patients were utilized to demonstrate the performance of the proposed motion tracking method. Results: The tracking accuracy was validated by comparing the results to the average of two manual delineations in 50 randomly selected cine image frames from each patient. The resulting average dice similarity coefficient (93.28% ± 1.46%) and margin error (0.49 ± 0.12 mm) showed good agreement between the automatic and manual results. The comparison with three other deformable model-based segmentation methods illustrated the superior shape tracking performance of the proposed method. Large interpatient variations of swallowing frequency, swallowing duration, and upper airway cross-sectional area were observed from the testing cine image sequences. Conclusions: The proposed motion tracking method can provide accurate upper airway motion tracking results, and enable automatic and quantitative identification and analysis of in-treatment H&N upper airway motion. By integrating explicit and implicit linked-shape representations within a hierarchical model-fitting process, the proposed tracking method can process complex H&N structures and low-contrast/resolution cine MRI images. Future research will focus on the improvement of method reliability, patient motion pattern analysis for providing more information on patient-specific prediction of structure displacements, and motion effects on dosimetry for better H&N motion management in radiation therapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-119: Pilot Study On the Location-Based Lung Motion Assessment

Purpose: In most of lung treatment cases with various radiotherapy beam modalities, 4DCT images are obtained in order to define ITV. ITV is defined with the signal from motion monitoring system, e.g. RPM. However, the signal is not consistent with tumor motion because it varies with location, its size, age, gender, etc. In the present study, the location-based motion assessment is presented. Methods: 4DCT images of 70 patients were reviewed: 28-left-lung and 42-right-lung patients; 36-female and 34-male patients; the age range of 51.2–89.9; tumor-size range of 0.75–9.50cm with 25% of these adherent to bony-anatomy. Philips Big-Bore Simulation CT and RPM systems were used. The study was performed as follows. First, RPM signal and tumor motion in superior-inferior direction was compared. Second, the tumor size and its motion amplitude in all directions were measured at multiple locations. Third, the average tumor motion was calculated to assess general motion amplitudes at various locations. Results: RPM amplitude is not consistent with lung tumor motion amplitude. The tumors of similar sizes at similar location present various motion amplitude up to 1.1cm difference, but in average, the standard deviation was <0.5cm. Almost regardless of tumor sizes, the tumor motion was greatest at lower lobe location (>=1.0cm), and the smallest at upper lobe location and when adherent to bony-anatomy (<=0.5cm). Conclusion: The tumor size affects the motion amplitude less than does the tumor location. However, as the study results indicate that tumor motion has noticeable variation and so further study with more patient cases is needed. Also, for the same patient, the RPM signal presents instability of breathing, and clinically the patient with the instability of RPM breathing of <=10% is selected for respiratory-gated radiotherapy and ∼25% of patients under current study was treated. Patient-specific motion-uncertainty margins are considered to be added following further study.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-128: Dosimetric Impact of Esophagus Motion in Spine Stereotactic Body Radiotherapy

Purpose: Acute esophageal toxicity is a common side effect in spine stereotactic body radiotherapy (SBRT). The respiratory motion may alter esophageal position from the planning scan resulting in excessive esophageal dose. Here we assessed the dosimetric impact resulting from the esophageal motion using 4DCT. Methods: Nine patients treated to their thoracic spines in one fraction of 24 Gy were identified for this study. The original plan on a free breathing CT was copied to each phase image of a 4DCT scan, recalculated, scaled, and accumulated to the free breathing CT using deformable image registration. A segment of esophagus was contoured in the vicinity of treatment target. Esophagus dose volume histogram (DVH) was generated for both the original planned dose and the accumulated 4D dose for comparison. In parallel, we performed a chained deformable registration of 4DCT phase images to estimate the motion magnitude of the esophagus in a breathing cycle. We examined the correlation between the motion magnitude and the dosimetric deviation. Results: The esophageal motion mostly exhibited in the superior-inferior direction. The cross-sectional motion was small. Esophagus motion at T1 vertebra level (0.7 mm) is much smaller than that at T11 vertebra level (6.5 mm). The difference of Dmax between the original and 4D dose distributions ranged from 9.1 cGy (esophagus motion: 5.6 mm) to 231.1 cGy (esophagus motion: 3.1 mm). The difference of D(5cc) ranged from 5 cGy (esophagus motion: 3.1 mm) to 85 cGy (esophagus motion: 3.3 mm). There was no correlation between the dosimetric deviation and the motion magnitude. The V(11.9Gy)<5cc constraint was met for each patient when examining the DVH calculated from the 4D dose. Conclusion: Respiratory motion did not result in substantial dose increase to esophagus in spine SBRT. 4DCT simulation may not be necessary with regards to esophageal dose assessment.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-G-BRA-13: An Advanced Deformable Lung Phantom for Analyzing the Dosimetric Impact of Respiratory Motion

Purpose: The difference between three-dimensional (3D) and four-dimensional (4D) dose is affected by factors such as tumor size and motion. To quantitatively analyze the effects of these factors, a phantom that can independently control for each factor is required. The purpose of this study is to develop a deformable lung phantom with the above attributes and evaluate characteristics. Methods: A phantom was designed to simulate diaphragm motion with amplitude in the range 1 to 7 cm and various periods of regular breathing. To simulate different size tumors, tumors were produced by pouring liquid silicone into custom molds created by a 3D printer. The accuracy of phantom diaphragm motion was assessed using calipers and protractor. To control tumor motion, tumor trajectories were evaluated using 4D computed tomography (CT), and diaphragm-tumor correlation curve was calculated by curve fitting method. Three-dimensional dose and 4D dose were calculated and compared according to tumor motion. Results: The accuracy of phantom diaphragm motion was less than 1 mm. Maximum tumor motion amplitudes in the left-right and anterior-posterior directions were 0.08 and 0.12 cm, respectively, in a 10 cm{sup 3} tumor, and 0.06 and 0.27 cm, respectively, in a 90 cm{sup 3} tumor. The diaphragm-tumor correlation curve showed that tumor motion in the superior-inferior direction was increased with increasing diaphragm motion. In the 10 cm{sup 3} tumor, the tumor motion was larger than the 90 cm{sup 3} tumor. According to tumor motion, variation of dose difference between 3D and 4D was identified. Conclusion: The developed phantom can independently control factors such as tumor size and motion. In potentially, this phantom can be used to quantitatively analyze the dosimetric impact of respiratory motion according to the factors that influence the difference between 3D and 4D dose. This research was supported by the Mid-career Researcher Program through NRF funded by the Ministry of Science, ICT & Future Planning of Korea (NRF-2014R1A2A1A10050270) and by the Radiation Technology R&D program through the National Research Foundation of Korea funded by the Ministry of Science, ICT & Future Planning (No. 2013M2A2A7038291)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

MO-FG-BRA-09: Towards an Optimal Breath-Holding Procedure for Radiotherapy: Differences in Organ Motion During Inhalation and Exhalation Breath-Holds

Purpose: Breath-holding (BH) is often used to reduce organ motion during radiotherapy. The aim of this study was to determine the differences in pancreatic and diaphragmatic motion during BH between inhalation and exhalation BHs with variable lung volumes and to investigate whether motion increases/decreases during BH. Methods: Sixteen healthy volunteers were asked to perform four different 60-second BHs, from fully inflated to fully deflated lungs (i.e. lung volumes of: 100%, ∼70%, ∼30% and 0% of inspiratory capacity) three times (total of 192 BHs). During each BH, we obtained single-slice (coronal) magnetic-resonance scans with spatial resolution 0.93×0.93×8.0 mm3 and temporal resolution 0.6 s. We used 2-dimensional image correlation to obtain the motion of pancreatic head and diaphragm during BH. Motion magnitude in inferior-superior direction was obtained by determining the maximum displacement during BH. Results: Pancreatic and diaphragmatic drifts occurred during BH and were mostly in the superior direction. We observed significantly smaller pancreatic and diaphragmatic motion magnitudes in inferior-superior direction during exhalation BHs (BH{sub 30%} and BH{sub 0%}) compared to inhalation BHs (BH{sub 100%} and BH{sub 70%}). The mean motion magnitudes of the pancreatic head were 7.0, 6.5, 4.4 and 4.2 mm during BH{sub 100%}, BH{sub 70%}, BH{sub 30%} and BH{sub 0%}, respectively, and mean BH durations were 59.9, 59.1, 59.0 and 52.7 s. For the diaphragm, mean motion magnitudes were 9.8, 9.0, 5.6 and 4.3 mm, respectively. When considering 30-second BHs, as often used in the clinic, the motion was most pronounced during the first 10 s and excluding these from the analysis (yielding an effective BH period of 20 s) significantly reduced (P≤0.002) organ motion. Conclusion: Organ motion was significantly smaller during exhalation BHs compared to inhalation BHs. Also, motion was largest at the start of BH. Hence, waiting for 10 s may significantly decrease motion of the pancreas and diaphragm during treatment.

60 APPLIED LIFE SCIENCES↗