Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EXHALATION”

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 91 records · Page 5

Investigating Fracture Network Deformation Using Noble Gas Release

We investigate deformation mechanics of fracture networks in unsaturated fractured rocks from subsurface conventional detonation using dynamic noble gas measurements and changes in air permeability. We dynamically measured the noble gas isotopic composition and helium exhalation of downhole gas before and after a large subsurface conventional detonation. These noble gas measurements were combined with measurements of the subsurface permeability field from 64 discrete sampling intervals before and after the detonation and subsurface mapping of fractures in borehole walls before well completion. We saw no observable increase in radiogenic noble gas release from either an isotopic composition or a helium exhalation point of view. Large increases in permeability were observed in 13 of 64 discrete sampling intervals. Of the sampling intervals which saw large increases in flow, only two locations did not have preexisting fractures mapped at the site. Given the lack of noble gas release and a clear increase in permeability, we infer that most of the strain accommodation of the fractured media occurred along previously existing fractures, rather than the creation of new fractures, even for a high strain rate event. These results have significant implications for how we conceptualize the deformation of rocks with fracture networks above the percolation threshold, with application to a variety of geologic and geological engineering problems.

Gardner, W. Payton↗

Components of aircraft life support systems interact with each other and the user

The life support system in a tactical aircraft provides necessary supplemental oxygen to the aircrew. However, interactions among its various components may generate unexpected breathing loads. We focus here on the interactions between a regulator and breathing mask commonly used together in the U.S. Navy, the CRU-103 regulator and MBU 23/P mask, and some effects of the interactions on the user. The data reported were collected during a larger research effort examining potential physiological and cognitive effects of low regulator inlet pressures. Seventeen participants completed a series of tasks under mild exercise while breathing 40% O 2 (balance N 2 ) from an MBU-23/P mask supplied by a CRU-103 regulator with supply pressures 10, 6, 4, and 2 psig (CRU-103 specifications are for inlet pressures from 5 to 120 psig). Variables measured included flow to the mask and pressures at the regulator supply, in the hose to the mask, and in the mask. In addition to restricting inspiratory flow, low inlet pressure to the CRU-103 caused a counterintuitive overshoot in gas delivery pressure at end-inspiration, a mean increase of 1.5 cm H 2 O between the 10- and 2 psig conditions. The added pressure to the exhalation valve increased the expiratory threshold, the pressure to start expiratory flow, by approximately 2 cm H 2 O, increasing the effort needed to exhale.

60 APPLIED LIFE SCIENCES↗

Metabolic analyzer

An apparatus is described for the measurement of metabolic rate and breathing dynamics in which inhaled and exhaled breath are sensed by sealed, piston-displacement type spirometers. These spirometers electrically measure the volume of inhaled and exhaled breath. A mass spectrometer analyzes simultaneously for oxygen, carbon dioxide, nitrogen, and water vapor. Circuits responsive to the outputs of the spirometers, mass spectrometer, temperature, pressure, and timing signals compute oxygen consumption, carbon dioxide production, minute volume, and respiratory exchange ratio. A selective indicator provides for readout of these data at predetermined cyclic intervals.

Perry, C. L.↗

Metabolic analyzer

An apparatus is described for the measurement of metabolic rate and breathing dynamics in which inhaled and exhaled breath are sensed by sealed, piston-displacement type spirometers. These spirometers electrically measure the volume of inhaled and exhaled breath. A mass spectrometer analyzes simultaneously for oxygen, carbon dioxide, nitrogen and water vapor. Computation circuits are responsive to the outputs of the spirometers, mass spectrometer, temperature, pressure and timing signals and compute oxygen consumption, carbon dioxide production, minute volume and respiratory exchange ratio. A selective indicator provides for read-out of these data at predetermined cyclic intervals.

Rummel, J. A.↗

Metabolic analyzer

The metabolic analyzer was designed to support experiment M171. It operates on the so-called open circuit method to measure a subject's metabolic activity in terms of oxygen consumed, carbon dioxide produced, minute volume, respiratory exchange ratio, and tidal volume or vital capacity. The system operates in either of two modes. (1) In Mode I, inhaled respiratory volumes are actually measured by a piston spirometer. (2) In Mode II, inhaled volumes are calculated from the exhaled volume and the measured inhaled and exhaled nitrogen concentrations. This second mode was the prime mode for Skylab. Following is a brief description of the various subsystems and their operation.

Lem, J. D.↗

Metabolic rate meter and method

A method is described for measuring the dynamic metabolic rate of a human or animal. The ratio of the exhaled carbon dioxide to a known amount of C(13)02 introduced into the exhalation is determined by mass spectrometry. This provides an instantaneous measurement of the carbon dioxide generated.

Taylor, T. I.↗

A source of experimental underestimation of aerosol bolus deposition

We examined the measurement error in inhaled and exhaled aerosol concentration resulting from the bolus delivery system when small volumes of monodisperse aerosols are inspired to different lung depths. A laser photometer that illuminated approximately 75% of the breathing path cross section recorded low inhaled bolus half-widths (42 ml) and negative deposition values for shallow bolus inhalation when the inhalation path of a 60-ml aerosol was straight and unobstructed. We attributed these results to incomplete mixing of the inhaled aerosol bolus over the breathing path cross section, on the basis of simultaneous recordings of the photometer with a particle-counter sampling from either the center or the edge of the breathing path. Inserting a 90 degrees bend into the inhaled bolus path increased the photometer measurement of inhaled bolus half-width to 57 ml and yielded positive deposition values. Dispersion, which is predominantly affected by exhaled bolus half-width, was not significantly altered by the 90 degrees bend. We conclude that aerosol bolus-delivery systems should ensure adequate mixing of the inhaled bolus to avoid error in measurement of bolus deposition.

Non-NASA Center↗

Electronic Nose Development and Preliminary Human Breath Testing for Rapid, Non-Invasive COVID-19 Detection

We adapted an existing, spaceflight-proven, robust “electronic nose” (E-Nose) that uses an array of electrical resistivity-based nanosensors mimicking aspects of mammalian olfaction to conduct on-site, rapid screening for COVID-19 infection by measuring the pattern of sensor responses to volatile organic compounds (VOCs) in exhaled human breath. We built and tested multiple copies of a hand-held prototype E-Nose sensor system, composed of 64 chemically sensitive nanomaterial sensing elements tailored to COVID-19 VOC detection; data acquisition electronics; a smart tablet with software (App) for sensor control, data acquisition and display; and a sampling fixture to capture exhaled breath samples and deliver them to the sensor array inside the E-Nose. The sensing elements detect the combination of VOCs typical in breath at parts-per-billion (ppb) levels, with repeatability of 0.02% and reproducibility of 1.2%; the measurement electronics in the E-Nose provide measurement accuracy and signal-to-noise ratios comparable to benchtop instrumentation. Preliminary clinical testing at Stanford Medicine with 63 participants, their COVID-19-positive or COVID-19-negative status determined by concomitant RT-PCR, discriminated between these two categories of human breath with a 79% correct identification rate using “leave-one-out” training-and-analysis methods. Analyzing the E-Nose response in conjunction with body temperature and other non-invasive symptom screening using advanced machine learning methods, with a much larger database of responses from a wider swath of the population, is expected to provide more accurate on-the-spot answers. Additional clinical testing, design refinement, and a mass manufacturing approach are the main steps toward deploying this technology to rapidly screen for active infection in clinics and hospitals, public and commercial venues, or at home.

COVID-19↗

Interfractional Dose Variations in Intensity-Modulated Radiotherapy With Breath-Hold for Pancreatic Cancer

Purpose: To investigate the interfractional dose variations for intensity-modulated radiotherapy (RT) combined with breath-hold (BH) at end-exhalation (EE) for pancreatic cancer. Methods and Materials: A total of 10 consecutive patients with pancreatic cancer were enrolled. Each patient was fixed in the supine position on an individualized vacuum pillow with both arms raised. Computed tomography (CT) scans were performed before RT, and three additional scans were performed during the course of chemoradiotherapy using a conventional RT technique. The CT data were acquired under EE-BH conditions (BH-CT) using a visual feedback technique. The intensity-modulated RT plan, which used five 15-MV coplanar ports, was designed on the initial BH-CT set with a prescription dose of 39 Gy at 2.6 Gy/fraction. After rigid image registration between the initial and subsequent BH-CT scans, the dose distributions were recalculated on the subsequent BH-CT images under the same conditions as in planning. Changes in the dose-volume metrics of the gross tumor volume (GTV), clinical target volume (CTV = GTV + 5 mm), stomach, and duodenum were evaluated. Results: For the GTV and clinical target volume (CTV), the 95th percentile of the interfractional variations in the maximal dose, mean dose, dose covering 95% volume of the region of structure, and percentage of the volume covered by the 90% isodose line were within {+-}3%. Although the volume covered by the 39 Gy isodose line for the stomach and duodenum did not exceed 0.1 mL at planning, the volume covered by the 39 Gy isodose line for these structures was up to 11.4 cm{sup 3} and 1.8 cm{sup 3}, respectively. Conclusions: Despite variations in the gastrointestinal state and abdominal wall position at EE, the GTV and CTV were mostly ensured at the planned dose, with the exception of 1 patient. Compared with the duodenum, large variations in the stomach volume receiving high-dose radiation were observed, which might be beyond the negligible range in achieving dose escalation with intensity-modulated RT combined with BH at EE.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Respiratory Amplitude Guided 4-Dimensional Magnetic Resonance Imaging

Purpose: To evaluate the feasibility of prospectively guiding 4-dimensional (4D) magnetic resonance imaging (MRI) image acquisition using triggers at preselected respiratory amplitudes to achieve T{sub 2} weighting for abdominal motion tracking. Methods and Materials: A respiratory amplitude-based triggering system was developed and integrated into a commercial turbo spin echo MRI sequence. Initial feasibility tests were performed on healthy human study participants. Four respiratory states, the middle and the end of inhalation and exhalation, were used to trigger 4D MRI image acquisition of the liver. To achieve T{sub 2} weighting, the echo time and repetition time were set to 75 milliseconds and 4108 milliseconds, respectively. Single-shot acquisition, together with parallel imaging and partial k-space imaging techniques, was used to improve image acquisition efficiency. 4D MRI image sets composed of axial or sagittal slices were acquired. Results: Respiratory data measured and logged by the MRI scanner showed that the triggers occurred at the appropriate respiratory levels. Liver motion could be easily observed on both 4D MRI image datasets by sensing either the change of liver in size and shape (axial) or diaphragm motion (sagittal). Both 4D MRI image datasets were T{sub 2}-weighted as expected. Conclusions: This study demonstrated the feasibility of achieving T{sub 2}-weighted 4D MRI images using amplitude-based respiratory triggers. With the aid of the respiratory amplitude-based triggering system, the proposed method is compatible with most MRI sequences and therefore has the potential to improve tumor-tissue contrast in abdominal tumor motion imaging.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

SU-E-J-187: Individually Optimized Contrast-Enhancement 4D-CT for Pancreatic Adenocarcinoma in Radiotherapy Simulation

Purpose: To study the feasibility of individually optimized contrastenhancement (CE) 4D-CT for pancreatic adenocarcinoma (PDA) in radiotherapy simulation. To evaluate the image quality and contrast enhancement of tumor in the CE 4D-CT, compared to the clinical standard of CE 3D-CT and 4D-CT. Methods: In this IRB-approved study, each of the 7 PDA patients enrolled underwent 3 CT scans: a free-breathing 3D-CT with contrast (CE 3D-CT) followed by a 4D-CT without contrast (4D-CT) in the first study session, and a 4D-CT with individually synchronized contrast injection (CE 4D-CT) in the second study session. In CE 4D-CT, the time of full contrast injection was determined based on the time of peak enhancement for the test injection, injection rate, table speed, and longitudinal location and span of the pancreatic region. Physicians contoured both the tumor (T) and the normal pancreatic parenchyma (P) on the three CTs (end-of-exhalation for 4D-CT). The contrast between the tumor and normal pancreatic tissue was computed as the difference of the mean enhancement level of three 1 cm3 regions of interests in T and P, respectively. Wilcoxon rank sum test was used to statistically compare the scores and contrasts. Results: In qualitative evaluations, both CE 3D-CT and CE 4D-CT scored significantly better than 4D-CT (4.0 and 3.6 vs. 2.6). There was no significant difference between CE 3D-CT and CE 4D-CT. In quantitative evaluations, the contrasts between the tumor and the normal pancreatic parenchyma were 0.6±23.4, −2.1±8.0, and −19.6±28.8 HU, in CE 3D-CT, 4D-CT, and CE 4D-CT, respectively. Although not statistically significant, CE 4D-CT achieved better contrast enhancement between the tumor and the normal pancreatic parenchyma than both CE 3D-CT and 4DCT. Conclusion: CE 4D-CT achieved equivalent image quality and better contrast enhancement between tumor and normal pancreatic parenchyma than the clinical standard of CE 3D-CT and 4D-CT. This study was supported in part by Philips Healthcare.

60 APPLIED LIFE SCIENCES↗

MO-G-18C-03: Evaluation of Deformable Image Registration for Lung Motion Estimation Using Hyperpolarized Gas Tagging MRI

Purpose: Hyperpolarized gas (HP) tagging MRI is a novel imaging technique for direct measurement of lung motion during breathing. This study aims to quantitatively evaluate the accuracy of deformable image registration (DIR) in lung motion estimation using HP tagging MRI as references. Methods: Three healthy subjects were imaged using the HP MR tagging, as well as a high-resolution 3D proton MR sequence (TrueFISP) at the end-of-inhalation (EOI) and the end-of-exhalation (EOE). Ground truth of lung motion and corresponding displacement vector field (tDVF) was derived from HP tagging MRI by manually tracking the displacement of tagging grids between EOI and EOE. Seven different DIR methods were applied to the high-resolution TrueFISP MR images (EOI and EOE) to generate the DIR-based DVFs (dDVF). The DIR methods include Velocity (VEL), MIM, Mirada, multi-grid B-spline from Elastix (MGB) and 3 other algorithms from DIRART toolbox (Double Force Demons (DFD), Improved Lucas-Kanade (ILK), and Iterative Optical Flow (IOF)). All registrations were performed by independent experts. Target registration error (TRE) was calculated as tDVF – dDVF. Analysis was performed for the entire lungs, and separately for the upper and lower lungs. Results: Significant differences between tDVF and dDVF were observed. Besides the DFD and IOF algorithms, all other dDVFs showed similarity in deformation magnitude distribution but away from the ground truth. The average TRE for entire lung ranged 2.5−23.7mm (mean=8.8mm), depending on the DIR method and subject's breathing amplitude. Larger TRE (13.3–23.7mm) was found in subject with larger breathing amplitude of 45.6mm. TRE was greater in lower lung (2.5−33.9 mm, mean=12.4mm) than that in upper lung (2.5−11.9 mm, mean=5.8mm). Conclusion: Significant differences were observed in lung motion estimation between the HP gas tagging MRI method and the DIR methods, especially when lung motion is large. Large variation among different DIR methods was also observed.

60 APPLIED LIFE SCIENCES↗

Comparison of an alternative and existing binning methods to reduce the acquisition duration of 4D PET/CT

Purpose: Respiratory motion is a source of artifacts that reduce image quality in PET. Four dimensional (4D) PET/CT is one approach to overcome this problem. Existing techniques to limiting the effects of respiratory motions are based on prospective phase binning which requires a long acquisition duration (15–25 min). This time is uncomfortable for the patients and limits the clinical exploitation of 4D PET/CT. In this work, the authors evaluated an existing method and an alternative retrospective binning method to reduce the acquisition duration of 4D PET/CT. Methods: The authors studied an existing mixed-amplitude binning (MAB) method and an alternative binning method by mixed-phases (MPhB). Before implementing MPhB, they analyzed the regularity of the breathing patterns in patients. They studied the breathing signal drift and missing CT slices that could be challenging for implementing MAB. They compared the performance of MAB and MPhB with current binning methods to measure the maximum uptake, internal volume, and maximal range of tumor motion. Results: MPhB can be implemented depending on an optimal phase (in average, the exhalation peak phase −4.1% of the entire breathing cycle duration). Signal drift of patients was in average 35% relative to the breathing amplitude. Even after correcting this drift, MAB was feasible in 4D CT for only 64% of patients. No significant differences appeared between the different binning methods to measure the maximum uptake, internal volume, and maximal range of tumor motion. The authors also determined the inaccuracies of MAB and MPhB to measure the maximum amplitude of tumor motion with three bins (less than 3 mm for movement inferior to 12 mm, up to 6.4 mm for a 21 mm movement). Conclusions: The authors proposed an alternative binning method by mixed-phase binning that halves the acquisition duration of 4D PET/CT. Mixed-amplitude binning was challenging because of signal drift and missing CT slices. They showed that more than three bins were necessary for a more accurate measurement of the maximum amplitude of the tumor motion. However, the current 4D-CT technology limits the increase of the number of bins in 4D PET/CT because of missing CT slices. One can reconstruct 4D PET images with more bins but without attenuation/scatter correction.

60 APPLIED LIFE SCIENCES↗

Rapid estimation of 4DCT motion-artifact severity based on 1D breathing-surrogate periodicity

Purpose: Motion artifacts are common in patient four-dimensional computed tomography (4DCT) images, leading to an ill-defined tumor volume with large variations for radiotherapy treatment and a poor foundation with low imaging fidelity for studying respiratory motion. The authors developed a method to estimate 4DCT image quality by establishing a correlation between the severity of motion artifacts in 4DCT images and the periodicity of the corresponding 1D respiratory waveform (1DRW) used for phase binning in 4DCT reconstruction. Methods: Discrete Fourier transformation (DFT) was applied to analyze 1DRW periodicity. The breathing periodicity index (BPI) was defined as the sum of the largest five Fourier coefficients, ranging from 0 to 1. Distortional motion artifacts (excluding blurring) of cine-scan 4DCT at the junctions of adjacent couch positions around the diaphragm were classified in three categories: incomplete, overlapping, and duplicate anatomies. To quantify these artifacts, discontinuity of the diaphragm at the junctions was measured in distance and averaged along six directions in three orthogonal views. Artifacts per junction (APJ) across the entire diaphragm were calculated in each breathing phase and phase-averaged APJ{sup ¯}, defined as motion-artifact severity (MAS), was obtained for each patient. To make MAS independent of patient-specific motion amplitude, two new MAS quantities were defined: MAS{sup D} is normalized to the maximum diaphragmatic displacement and MAS{sup V} is normalized to the mean diaphragmatic velocity (the breathing period was obtained from DFT analysis of 1DRW). Twenty-six patients’ free-breathing 4DCT images and corresponding 1DRW data were studied. Results: Higher APJ values were found around midventilation and full inhalation while the lowest APJ values were around full exhalation. The distribution of MAS is close to Poisson distribution with a mean of 2.2 mm. The BPI among the 26 patients was calculated with a value ranging from 0.25 to 0.93. The DFT calculation was within 3 s per 1DRW. Correlations were found between 1DRW periodicity and 4DCT artifact severity: −0.71 for MAS{sup D} and −0.73 for MAS{sup V}. A BPI greater than 0.85 in a 1DRW suggests minimal motion artifacts in the corresponding 4DCT images. Conclusions: The breathing periodicity index and motion-artifact severity index are introduced to assess the relationship between 1DRW and 4DCT. A correlation between 1DRW periodicity and 4DCT artifact severity has been established. The 1DRW periodicity provides a rapid means to estimate 4DCT image quality. The rapid 1DRW analysis and the correlative relationship can be applied prospectively to identify irregular breathers as candidates for breath coaching prior to 4DCT scan and retrospectively to select high-quality 4DCT images for clinical motion-management research.

60 APPLIED LIFE SCIENCES↗

SU-E-T-281: Dose Measurements of Modulated Spot-Scanning Particle Beams with Beam-Gating of Respiratory-Phase

Purpose: To present dosimetric variations due to target movements with beam-gating windows of various respiratory-phase in homogeneous phantom irradiated by modulated spot-scanning carbon-ions and protons in Siemens IONTRIS system. Methods: For the safety and efficacy of proton therapy to treat three patients with lung cancer during our clinical trial, residual motion was required within 5mm. To study dosimetric variations due to respiratory movement, a target of 4.0cmx4.0cm with 2cm thick was moved with distances of 3.5mm, 4.4mm, 6.0mm, 8.3mm, and 11.0mm with beam-gating windows between various phases of inhalation or exhalation at a sin-wave breathing pattern. The target was irritated at 110mm depth by modulated carbon-ions with focus sizes of 4.1mm to 4.6mm, and a grid size of 1mm between spots. And, by modulated protons with focus sizes of 11.4mm to 13.6mm, and a grid size of 1.9mm between spots. A 4.0cmx4.0cm field size was used for both beams. EBT3 films was placed at the center of target for measurements. Delivery dose was 5.0 Gy with >1.0% uniformity over target. The uniformity and field size of each measured 2D lateral profiles were extracted. Results: By irradiating films to a doses at linear region of dose response, the uniformity and field size were extracted by measured optical density. The measured deviations from calculated field width and dose increase with increased motion amplitude. Larger non-uniformity was observed for carbonion with smaller focus size in comparing with protons. For movements of 4.4mm and 6.0mm, Optical Density uniformity of 3.80% and 5.66%were observed for carbon beam. But, is 3.46% for protons with 11.4mm movement. Conclusion: Our investigations showed that 5% optical density uniformity for carbon-ions and protons might be acceptable for treatments with 8.3mm movements in homogeneous phantoms. Dose variation introduced in complex anatomy of real patients need further investigation.

60 APPLIED LIFE SCIENCES↗

Lung motion estimation using dynamic point shifting: An innovative model based on a robust point matching algorithm

Purpose: Image-guided radiotherapy is an advanced 4D radiotherapy technique that has been developed in recent years. However, respiratory motion causes significant uncertainties in image-guided radiotherapy procedures. To address these issues, an innovative lung motion estimation model based on a robust point matching is proposed in this paper. Methods: An innovative robust point matching algorithm using dynamic point shifting is proposed to estimate patient-specific lung motion during free breathing from 4D computed tomography data. The correspondence of the landmark points is determined from the Euclidean distance between the landmark points and the similarity between the local images that are centered at points at the same time. To ensure that the points in the source image correspond to the points in the target image during other phases, the virtual target points are first created and shifted based on the similarity between the local image centered at the source point and the local image centered at the virtual target point. Second, the target points are shifted by the constrained inverse function mapping the target points to the virtual target points. The source point set and shifted target point set are used to estimate the transformation function between the source image and target image. Results: The performances of the authors’ method are evaluated on two publicly available DIR-lab and POPI-model lung datasets. For computing target registration errors on 750 landmark points in six phases of the DIR-lab dataset and 37 landmark points in ten phases of the POPI-model dataset, the mean and standard deviation by the authors’ method are 1.11 and 1.11 mm, but they are 2.33 and 2.32 mm without considering image intensity, and 1.17 and 1.19 mm with sliding conditions. For the two phases of maximum inhalation and maximum exhalation in the DIR-lab dataset with 300 landmark points of each case, the mean and standard deviation of target registration errors on the 3000 landmark points of ten cases by the authors’ method are 1.21 and 1.04 mm. In the EMPIRE10 lung registration challenge, the authors’ method ranks 24 of 39. According to the index of the maximum shear stretch, the authors’ method is also efficient to describe the discontinuous motion at the lung boundaries. Conclusions: By establishing the correspondence of the landmark points in the source phase and the other target phases combining shape matching and image intensity matching together, the mismatching issue in the robust point matching algorithm is adequately addressed. The target registration errors are statistically reduced by shifting the virtual target points and target points. The authors’ method with consideration of sliding conditions can effectively estimate the discontinuous motion, and the estimated motion is natural. The primary limitation of the proposed method is that the temporal constraints of the trajectories of voxels are not introduced into the motion model. However, the proposed method provides satisfactory motion information, which results in precise tumor coverage by the radiation dose during radiotherapy.

60 APPLIED LIFE SCIENCES↗

Multivariate analysis for the estimation of target localization errors in fiducial marker-based radiotherapy

Purpose: To assess the target localization error (TLE) in terms of the distance between the target and the localization point estimated from the surrogates (|TMD|), the average of respiratory motion for the surrogates and the target (|aRM|), and the number of fiducial markers used for estimating the target (n). Methods: This study enrolled 17 lung cancer patients who subsequently underwent four fractions of real-time tumor tracking irradiation. Four or five fiducial markers were implanted around the lung tumor. The three-dimensional (3D) distance between the tumor and markers was at maximum 58.7 mm. One of the markers was used as the target (P{sub t}), and those markers with a 3D |TMD{sub n}| ≤ 58.7 mm at end-exhalation were then selected. The estimated target position (P{sub e}) was calculated from a localization point consisting of one to three markers except P{sub t}. Respiratory motion for P{sub t} and P{sub e} was defined as the root mean square of each displacement, and |aRM| was calculated from the mean value. TLE was defined as the root mean square of each difference between P{sub t} and P{sub e} during the monitoring of each fraction. These procedures were performed repeatedly using the remaining markers. To provide the best guidance on the answer with n and |TMD|, fiducial markers with a 3D |aRM ≥ 10 mm were selected. Finally, a total of 205, 282, and 76 TLEs that fulfilled the 3D |TMD| and 3D |aRM| criteria were obtained for n = 1, 2, and 3, respectively. Multiple regression analysis (MRA) was used to evaluate TLE as a function of |TMD| and |aRM| in each n. Results: |TMD| for n = 1 was larger than that for n = 3. Moreover, |aRM| was almost constant for all n, indicating a similar scale for the marker’s motion near the lung tumor. MRA showed that |aRM| in the left–right direction was the major cause of TLE; however, the contribution made little difference to the 3D TLE because of the small amount of motion in the left–right direction. The TLE calculated from the MRA ({sup MRA}TLE) increased as |TMD| and |aRM| increased and adversely decreased with each increment of n. The median 3D {sup MRA}TLE was 2.0 mm (range, 0.6–4.3 mm) for n = 1, 1.8 mm (range, 0.4–4.0 mm) for n = 2, and 1.6 mm (range, 0.3–3.7 mm) for n = 3. Although statistical significance between n = 1 and n = 3 was observed in all directions, the absolute average difference and the standard deviation of the {sup MRA}TLE between n = 1 and n = 3 were 0.5 and 0.2 mm, respectively. Conclusions: A large |TMD| and |aRM| increased the differences in TLE between each n; however, the difference in 3D {sup MRA}TLEs was, at most, 0.6 mm. Thus, the authors conclude that it is acceptable to continue fiducial marker-based radiotherapy as long as |TMD| is maintained at ≤58.7 mm for a 3D |aRM| ≥ 10 mm.

60 APPLIED LIFE SCIENCES↗

A Multiscale Bidirectional Coupling Framework

The lung is geometrically articulated across multiple scales from the trachea to the alveoli. A major computational challenge is to tightly link ODEs that describe lower scales to 3D finite element or finite volume models of airway mechanics using iterative communication between scales. In this study, we developed a novel multiscale computational framework for bidirectionally coupling 3D CFD models and systems of lower order ODEs. To validate the coupling framework, a four and eight generation Weibel lung model was constructed. For the coupled CFD-ODE simulations, the lung models were truncated at different generations and a RL circuit represented the truncated portion. The flow characteristics from the coupled models were compared to untruncated full 3D CFD models at peak inhalation and peak exhalation. Results showed that at no time or simulation was the difference in mass flux and/or pressure at a given location between uncoupled and coupled models was greater than 2.43%. The flow characteristics at prime locations for the coupled models showed good agreement to uncoupled models. Remarkably, due to reuse of the Krylov subspace, the cost of the ODE coupling is not much greater than uncoupled full 3D-CFD computations with simple prescribed pressure values at the outlets.

59 BASIC BIOLOGICAL SCIENCES↗