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 1,657 records · Page 92

Uncertainty quantification for deep learning in particle accelerator applications

With the advent of increased computational resources and improved algorithms, machine learning-based models are being increasingly applied to complex problems in particle accelerators. However, such data-driven models may provide overly confident predictions with unknown errors and uncertainties. For reliable deployment of machine learning models in high-regret and safety-critical systems such as particle accelerators, estimates of prediction uncertainty are needed along with accurate point predictions. In this investigation, we evaluate Bayesian neural networks (BNN) as an approach that can provide accurate predictions along with reliably quantified uncertainties for particle accelerator problems, and compare their performance with bootstrapped ensembles of neural networks. We select three accelerator setups for this evaluation: a storage ring, a photoinjector, and a linac. The problems span different data volumes and dimensionalities (e.g., scalar predictions as well as image outputs). It is found that BNN provide accurate predictions of the mean along with reliable estimates of predictive uncertainty across the test cases. In this vein, BNN may offer an attractive alternative to deterministic deep learning tools to generate accurate predictions with quantified uncertainties in particle accelerator applications.

43 PARTICLE ACCELERATORS↗

2016 Billion-Ton Report: Advancing Domestic Resources for a Thriving Bioeconomy, Volume 2: Environmental Sustainability Effects of Select Scenarios from Volume 1

With the goal of understanding environmental effects of a growing bioeconomy, the U.S. Department of Energy (DOE), national laboratories, and U.S. Forest Service research laboratories, together with academic and industry collaborators, undertook a study to estimate environmental effects of potential biomass production scenarios in the United States, with an emphasis on agricultural and forest biomass. Potential effects investigated include changes in soil organic carbon (SOC), greenhouse gas (GHG) emissions, water quality and quantity, air emissions, and biodiversity. Effects of altered land-management regimes were analyzed based on select county-level biomass-production scenarios for 2017 and 2040 taken from the 2016 Billion-Ton Report: Advancing Domestic Resources for a Thriving Bioeconomy (BT16), volume 1, which assumes that the land bases for agricultural and forestry would not change over time. The scenarios reflect constraints on biomass supply (e.g., excluded areas; implementation of management practices; and consideration of food, feed, forage, and fiber demands and exports) that intend to address sustainability concerns. Nonetheless, both beneficial and adverse environmental effects might be expected. To characterize these potential effects, this research sought to estimate where and under what modeled scenarios or conditions positive and negative environmental effects could occur nationwide. The report also includes a discussion of land-use change (LUC) (i.e., land management change) assumptions associated with the scenario transitions (but not including analysis of indirect LUC [ILUC]), analyses of climate sensitivity of feedstock productivity under a set of potential scenarios, and a qualitative environmental effects analysis of algae production under carbon dioxide (CO2) co-location scenarios. Because BT16 biomass supplies are simulated independent of a defined end use, most analyses do not include benefits from displacing fossil fuels or other products, with the exception of including a few illustrative cases on potential reductions in GHG emissions and fossil energy consumption associated with using biomass supplies for fuel, power, heat, and chemicals.

09 BIOMASS FUELS↗

SU-E-T-399: Determination of the Radiobiological Parameters That Describe the Dose-Response Relations of Xerostomia and Disgeusia From Head and Neck Radiotherapy

Purpose: To estimate the radiobiological parameters that describe the doseresponse relations of xerostomia and disgeusia from head and neck cancer radiotherapy. To identify the organs that are best correlated with the manifestation of those clinical endpoints. Finally, to evaluate the goodnessof- fit by comparing the model predictions against the actual clinical results. Methods: In this study, 349 head and neck cancer patients were included. For each patient the dose volume histograms (DVH) of parotids (separate and combined), mandible, submandibular glands (separate and combined) and salivary glands were calculated. The follow-up of those patients was recorded at different times after the completion of the treatment (7 weeks, 3, 7, 12, 18 and 24 months). Acute and late xerostomia and acute disgeusia were the clinical endpoints examined. A maximum likelihood fitting was performed to calculate the best estimates of the parameters used by the relative seriality model. The statistical methods of the error distribution, the receiver operating characteristic (ROC) curve, the Pearson's test and the Akaike's information criterion were utilized to assess the goodness-of-fit and the agreement between the pattern of the radiobiological predictions with that of the clinical records. Results: The estimated values of the radiobiological parameters of salivary glands are D50 = 25.2 Gy, γ = 0.52, s = 0.001. The statistical analysis confirmed the clinical validity of those parameters (area under the ROC curve = 0.65 and AIC = 38.3). Conclusion: The analysis proved that the treatment outcome pattern of the patient material can be reproduced by the relative seriality model and the estimated radiobiological parameters. Salivary glands were found to have strong volume dependence (low relative seriality). Diminishing the biologically effective uniform dose to salivary glands below 30 Gy may significantly reduce the risk of complications to the patients irradiated for prostate cancer.

60 APPLIED LIFE SCIENCES↗

Sintering and collapse of synthetic coal ash and slag cones as observed through constant heating rate optical dilatometry

Ash agglomeration and slagging are responsible for costly maintenance problems in coal-fired boilers. Boiler design and operation rely on critical ash fusion temperatures defined by geometric criteria that explicitly ignore shrinkage or warping. Reproducibility of ash fusion test (AFT) results is known to be a major problem, given the subjective nature of some criteria used in critical temperature reporting. Decades of faulty predictions based on AFT data underscore the need for an objective measure of the likelihood for agglomeration and slagging as a function of coal ash composition and temperature. Heating microscopes and ash fusibility determinators are AFT instruments that are widely available and readily adapted for optical dilatometry of ash cones through digital image analysis. Optical dilatometry enables the continuous measurement of cone height and face area, and the calculation of changes in cone volume and density as a function of temperature during constant heating. A strong correlation was found between the sintering onset temperatures of two synthetic ash compositions measured through optical dilatometry and their corresponding solidus temperatures predicted through thermophysical modeling. The sintering onset temperatures of synthetic slag cones appears consistent with glass transition temperatures estimated using the Kauzmann-Beaman rule. Cone collapse is highly correlated with softening, hemispherical, and fluid temperatures for ash cones. Rapid changes in solid, liquid, and gas volume fractions near the collapse temperature suggest that a transition from capillary to slurry state may play an important role in triggering collapse. Slag cone collapse is better explained by a decrease in viscosity as a function of temperature in the glass working range. Furthermore, the sintering onset temperature and cone collapse temperature are readily understood through thermophysical modeling based on ash composition and provide two additional critical temperatures in the modeling of coal ash agglomeration and slagging.

01 COAL, LIGNITE, AND PEAT↗

In situ study of microstructure evolution and α → ω phase transition in annealed and pre-deformed Zr under hydrostatic loading

The detailed study of the effect of the initial microstructure on its evolution under hydrostatic compression before, during, and after the irreversible α → ω phase transformation and during pressure release in Zr using in situ x-ray diffraction is presented. Two samples were studied: one is plastically pre-deformed Zr with saturated hardness and the other is annealed. Phase transformation α → ω initiates at lower pressure for a pre-deformed sample but for a volume fraction of ω Zr, c > 0.7, a larger volume fraction is observed for the annealed sample. This implies that the proportionality between the athermal resistance to the transformation and the yield strength in the continuum phase transformation theory is invalid; an advanced version of the theory is outlined. Phenomenological plasticity theory under hydrostatic loading is outlined in terms of microstructural parameters, and plastic strain is estimated. During transformation, the first rule is suggested, i.e., the average domain size, microstrain, and dislocation density in ω Zr for c < 0.8 are functions of the volume fraction, c of ω Zr only, which are independent of the plastic strain tensor prior to transformation and pressure. The microstructure is not inherited during phase transformation. Surprisingly, for the annealed sample, the final dislocation density and the average microstrain after pressure release in the ω phase are larger than for the severely pre-deformed sample. The results suggest that an extended experimental basis is required for the predictive models for the combined pressure-induced phase transformations and microstructure evolutions.

36 MATERIALS SCIENCE↗

Medical Risk Estimates and Clinical Capability Needs for A Long Duration Artemis Mission

BACKGROUND Human exploration spaceflight missions to the Moon and Mars present unprecedented challenges for in-mission medical care. Compared with the ISS, the greater distance from Earth will mean increased mission durations, communication delays, limited to no resupply opportunities, and significant limitations on the evacuation of ill or injured crew. Spacecraft mass, volume, and power will be curtailed while higher demands will be placed on the crew’s knowledge, skills, and abilities. In this higher risk environment, it is important to: a) quantitatively estimate human system risk attributable to medical conditions, a process known as Probabilistic Risk Analysis, and b) use these estimates to inform medical system design. IMPACT (Informing Mission Planning via Analysis of Complex Tradespaces) is a PRA and medical trade space analysis tool developed by NASA to advance exploration mission medical system design. IMPACT v1.0 improves upon and will soon replace NASA’s existing tool, the Integrated Medical Model, with: a novel evidence base baselined to exploration environments; an expanded list of 119 medical conditions; a significant increase in the number of medical resources that can be utilized and in the flexibility of their use; and the modelling of time lost performing mission-specific tasks due to medical conditions. DISCUSSION: This abstract will present IMPACT estimates of medical system risk and clinical capability needs for an extended duration Artemis mission of approximately 9 months, with phases including outbound transit, 3 months on the Lunar Gateway space station, 3 months on the Lunar surface with EVAs, 3 months on Gateway (simulating the return phase of a Mars mission), and transit back to Earth. Medical system risk estimates include loss of crew life (LOCL), consideration of medical evacuation (known as return to definitive care – RTDC), and an estimate of crew time lost due to medical conditions (Task Time Lost – TTL). The presentation will also describe the medical conditions that are the greatest drivers of risk as well as the clinical capabilities and resources that have the largest effect on risk.

W Thompson↗

High-Pressure Mechanical Response of Two Vitreous Silicates

The improved understanding of the high-pressure response of vitreous silicates is relevant to sought-after improvements of transparent protective systems. In this study, the quasi-static high-pressure (> 1 GPa) mechanical responses of BOROFLOAT® and Starphire® vitreous silicates were examined using three methods.The first combined x-ray radiography with compression testing in a Paris-Edinburgh anvil at the Advanced Photon Source at Argonne National Laboratory, and is a method recently developed at the University of Alabama at Birmingham. An outcome of that testing is the portrayal of volume change as a function of (quasi) hydrostatic pressure or its equation of state. The second employed method was spherical indentation and is one that can be used to identify the onset of permanent deformation in brittle materials. The outcome of that testing is the estimation of an apparent yield stress. The last method involved the use of a double-toroid diamond anvil cell to compress a specimen. This test method subjects a relatively large volume of material to a quasi hydrostatic stress and whose volume enables material postmortem analysis. An outcome of that testing, when coupled with precision density measurements, is the identification of the apparent hydrostatic pressure when permanent densification initiates.The three test methods and results from the testing of BOROFLOAT® and Starphire® vitreous silicates are presented.

36 MATERIALS SCIENCE↗

Comparison of internal dose estimates obtained using organ-level, voxel S value, and Monte Carlo techniques

Purpose: The authors’ objective was to compare internal dose estimates obtained using the Organ Level Dose Assessment with Exponential Modeling (OLINDA/EXM) software, the voxel S value technique, and Monte Carlo simulation. Monte Carlo dose estimates were used as the reference standard to assess the impact of patient-specific anatomy on the final dose estimate. Methods: Six patients injected with{sup 99m}Tc-hydrazinonicotinamide-Tyr{sup 3}-octreotide were included in this study. A hybrid planar/SPECT imaging protocol was used to estimate {sup 99m}Tc time-integrated activity coefficients (TIACs) for kidneys, liver, spleen, and tumors. Additionally, TIACs were predicted for {sup 131}I, {sup 177}Lu, and {sup 90}Y assuming the same biological half-lives as the {sup 99m}Tc labeled tracer. The TIACs were used as input for OLINDA/EXM for organ-level dose calculation and voxel level dosimetry was performed using the voxel S value method and Monte Carlo simulation. Dose estimates for {sup 99m}Tc, {sup 131}I, {sup 177}Lu, and {sup 90}Y distributions were evaluated by comparing (i) organ-level S values corresponding to each method, (ii) total tumor and organ doses, (iii) differences in right and left kidney doses, and (iv) voxelized dose distributions calculated by Monte Carlo and the voxel S value technique. Results: The S values for all investigated radionuclides used by OLINDA/EXM and the corresponding patient-specific S values calculated by Monte Carlo agreed within 2.3% on average for self-irradiation, and differed by as much as 105% for cross-organ irradiation. Total organ doses calculated by OLINDA/EXM and the voxel S value technique agreed with Monte Carlo results within approximately ±7%. Differences between right and left kidney doses determined by Monte Carlo were as high as 73%. Comparison of the Monte Carlo and voxel S value dose distributions showed that each method produced similar dose volume histograms with a minimum dose covering 90% of the volume (D90) agreeing within ±3%, on average. Conclusions: Several aspects of OLINDA/EXM dose calculation were compared with patient-specific dose estimates obtained using Monte Carlo. Differences in patient anatomy led to large differences in cross-organ doses. However, total organ doses were still in good agreement since most of the deposited dose is due to self-irradiation. Comparison of voxelized doses calculated by Monte Carlo and the voxel S value technique showed that the 3D dose distributions produced by the respective methods are nearly identical.

60 APPLIED LIFE SCIENCES↗

Estimating Lossy Compressibility of Scientific Data Using Deep Neural Networks

Simulation based scientific applications generate increasingly large amounts of data on high-performance computing (HPC) systems. To allow data to be stored and analyzed efficiently, data compression is often utilized to reduce the volume and velocity of data. However, a question often raised by domain scientists is the level of compression that can be expected so that they can make more informed decisions, balancing between accuracy and performance. In this letter, we propose a deep neural network based approach for estimating the compressibility of scientific data. To train the neural network, we build both general features as well as compressor-specific features so that the characteristics of both data and lossy compressors are captured in training. Our approach is demonstrated to outperform a prior analytical model as well as a sampling based approach in the case of a biased estimation, i.e., for SZ. However, for the unbiased estimation (i.e., ZFP), the sampling based approach yields the best accuracy, despite the high overhead involved in sampling the target dataset.

97 MATHEMATICS AND COMPUTING↗

Analysis of nodal coverage utilizing image guided radiation therapy for primary gynecologic tumor volumes

Purpose: To evaluate radiation dose delivered to pelvic lymph nodes, if daily Image Guided Radiation Therapy (IGRT) was implemented with treatment shifts based on the primary site (primary clinical target volume [CTV]). Our secondary goal was to compare dosimetric coverage with patient outcomes. Materials and methods: A total of 10 female patients with gynecologic malignancies were evaluated retrospectively after completion of definitive intensity-modulated radiation therapy (IMRT) to their pelvic lymph nodes and primary tumor site. IGRT consisted of daily kilovoltage computed tomography (CT)-on-rails imaging fused with initial planning scans for position verification. The initial plan was created using Varian's Eclipse treatment planning software. Patients were treated with a median radiation dose of 45 Gy (range: 37.5 to 50 Gy) to the primary volume and 45 Gy (range: 45 to 64.8 Gy) to nodal structures. One IGRT scan per week was randomly selected from each patient's treatment course and re-planned on the Eclipse treatment planning station. CTVs were recreated by fusion on the IGRT image series, and the patient's treatment plan was applied to the new image set to calculate delivered dose. We evaluated the minimum, maximum, and 95% dose coverage for primary and nodal structures. Reconstructed primary tumor volumes were recreated within 4.7% of initial planning volume (0.9% to 8.6%), and reconstructed nodal volumes were recreated to within 2.9% of initial planning volume (0.01% to 5.5%). Results: Dosimetric parameters averaged less than 10% (range: 1% to 9%) of the original planned dose (45 Gy) for primary and nodal volumes on all patients (n = 10). For all patients, ≥99.3% of the primary tumor volume received ≥ 95% the prescribed dose (V95%) and the average minimum dose was 96.1% of the prescribed dose. In evaluating nodal CTV coverage, ≥ 99.8% of the volume received ≥ 95% the prescribed dose and the average minimum dose was 93%. In evaluating individual IGRT sessions, we found that 6 patients had an estimated minimal nodal CTV dose less than 90% (range: 78 to 99%) of that planned. With a median follow-up of 42.5 months, 2 patients experienced systemic disease progression at an average of 19.6 months. One patient was found to have a local or regional failure with an average follow-up of 42 months. Conclusion: Using only 3 dimensional IGRT corrections in gynecological radiation allows excellent coverage of the primary target volume and good average nodal CTV coverage. If IGRT corrections are based on alignment to the primary tumor volume, and is only able to be corrected in 3 degrees, this can create situations in which nodal volumes may be under dosed. Utilizing multiple IGRT sessions appears to average out dose discrepancies over the course of treatment. The implication of underdosing in a single IGRT session needs further evaluation in future studies. Based on the concern of minimum dose to a nodal target volume, these findings may signal caution when using IGRT and IMRT in gynecological radiation patients. Possible techniques to overcome this situation may include averaging shifts between tumor and nodal volume, use of a treatment couch with 6° of freedom, deformable registration, or adaptive planning.

61 RADIATION PROTECTION AND DOSIMETRY↗

SU-F-T-323: A Post-Mastectomy Radiation Therapy Dose Distribution Study Using Nanodots and Films

Purpose: In post-mastectomy radiation therapy (RT), skin dose must be accurately estimated to assess skin reactions such as erythema, desquamation and necrosis. Planning systems cannot always provide accurate dosimetry for target volumes distal to skin. Therefore, in-vivo dosimetry is necessary. A female anthropomorphic phantom was used with optically stimulated luminescence dosimeters (nanoDots) to measure dose to chest wall skin. In addition, EBT2 films was employed to measure dose to left lung and heart in post-mastectomy RT. Methods: Films and nanoDots were calibrated under full buildup conditions at 100cm SAD for 6MV photons. Five pieces of films were placed between slabs of Rando phantom to assess dose to left lung and heart. Two layers of 0.5cm thick bolus were used to cover the whole left chest. Six pairs of nanoDots were placed at medical and lateral aspects on the bolus surface, between the 0.5cm bolus layers, and under the bolus. Three control nanoDots were placed on chest wall to quantify imaging dose. The phantom was CT scanned with all dosimeters in place, and treatment planning was performed with tangential fields (200cGy). All dosimeters were contoured on CT and dose was extracted. NanoDots were read using nanoDot reader and films were scanned using film scanner. The measured and calculated doses were tabulated. Results: Dose to 12 nanoDots were evaluated. Dose variance for surface nanoDots were +3.8%, +2.7%, −5% and −9.8%. Those at lateral positions, with greater beam obliquity had larger variance than the medial positions. A similar trend was observed for other nanoDots (Table1). Point doses from films for heart and the left lung were 112.7cGy and 108.7cGy, with +10.2% and +9.04% deviation from calculated values, respectively. Conclusion: Dosimetry provided by the advanced planning system was verified using NanoDots and films. Both nanoDots and films provided good estimation of dose distribution in post-mastectomy RT.

60 APPLIED LIFE SCIENCES↗

TPSAS-NF1676L-28149-DND

An Optimal Estimation (OE) inversion method has been developed to retrieve vertical profiles of aerosol size distributions, volume/surface concentrations, complex refractive indices, and single scattering albedos using the High Spectral Resolution Lidar (HSRL) and polarimeter measurements. By combining the active and passive remote sensing data, we can take advantages of the high vertical resolution of the lidar measurements and the high information content of the combined observations. There are three components to the new retrieval system: 1) simultaneous aerosol vertical profile retrieval using HSRL data, 2) standalone aerosol retrieval using Research Scanning Polarimeter (RSP), and 3) combined HSRL and RSP optimal estimation retrieval. For the LIDAR retrieval, we solve for number concentrations, mode fractions, aerosol size distributions, and real and imaginary part of refractive indices simultaneously. Based on this information, other parameters such as effective particle radius, surface and volume concentration, single scattering albedo can be derived. The combined HSRL-polarimeter algorithm will provide improved single scattering albedo and refractive index retrievals. The newly developed retrieval algorithm is very fast, and it provides an optimal estimation solution with associated error estimates. By solving for the aerosol vertical profiles using HSRL extinctions and backscattering coefficients from different layers, the new OE method can effectively reduce the impact of measurement noise on the retrieved aerosol microphysical properties. It also enables us to effectively combing HSRL and RSP data in the single OE retrieval system since the passive polarimeter instruments measure path-integrated intensity and polarizations. We have performed sensitivity studies using simulated data and applied the OE algorithms to NASA airborne field data measured by the HSRL-2 and the RSP instruments.

Xu Liu↗

A complex approach to assessing porous structure of structured ceramics obtained by SPS technique

Highlights: • A complex approach to microstructural analysis of nanostructured SPS ceramics was presented. • A bunch of gas, intrusion and microscopic methods were applied for open/closed porosity determination in nanoceramics. • Kinetic TPR analysis for structural changes in nanoceramic was carried out. • 3D pore visualization for nanoceramics was provided via SEM/FIB tomography • Porosity calculations were performed via BET, BJH, t-plot, Washburn equation, FEI Avizo 9.1 software results. - Abstract: Spark plasma sintering (SPS) is a highly prospective technology providing rapid consolidation of nano-disperse compositions into structured ceramics with complex structural organization. Characterizing pore structure of such systems remains a challenge due to different types, size and shape of the pores as well as structural defects occurring in the material, thus requiring special methods to quantify certain porosity type. The paper presents a complex approach to studying structure of SPS hematite-based structured ceramics. The method involves a number of gas, intrusion and electron microscopy probing techniques to characterize open and closed pores in the ceramics bulk. Low temperature nitrogen sorption has been implemented to estimate open micro- and mesoporosity in the size range 0.7–50 nm using BET, BJH and t-plot methods to process experimental data. Additionally, structural peculiarities of consolidated structured ceramics have been studied using temperature-programmed reduction (TPR) with the kinetic analysis of the results being carried out. Volume of macropores sized above 50 nm has been measured with mercury intrusion. Closed pores and structural defects have been quantified with SEM and FIB/SEM tomography. By comparing experimental and estimated results we have shown that addressing all porosity types provides a highly precise structural characterization of the material.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Rainfall estimation over oceans from SMMR and SSM/I microwave data

The brightness temperature (T sub b) measured at 37 GHz shows fairly strong emission from rain, and only slight effects caused by scattering by ice above the rain clouds. At frequencies below 37 GHz, where the FOV is larger and the volume extinction coefficient is weaker, it is found that the observations do not yield appreciable additional information about rain. At 85 GHz (FOV = 15 km), where the volume extinction is considerably larger, direct information about rain below the clouds is usually masked. Based on the above idea, 37 GHz observations with a 30 km FOV from SMMR and SSM/I are selected to develop an empirical method for the estimation of rain rate. In this method, the statistics of the observed T sub b's at 37 GHz in a rain storm are related to the rain rate statistics in that storm. The underestimation of rain rate, arising from the inability of the radiometer to respond sensitively to rain rate above a given threshold is rectified in this technique with the aid of two parameters that depend on the total water vapor content in the atmosphere. The retrieved rain rates compare favorably with radar observations and monthly mean global maps of rain derived from this technique over the oceans.

Prabhakara, C.↗

SU-E-QI-11: Measurement of Renal Pyruvate-To-Lactate Exchange with Hyperpolarized 13C MRI

Purpose: Previous work [1] modeling the metabolic flux between hyperpolarized [1-13C]pyruvate and [1-13C]lactate in magnetic resonance spectroscopic imaging (MRSI) experiments failed to account for vascular signal artifacts. Here, we investigate a method to minimize the vascular signal and its impact on the fidelity of metabolic modeling. Methods: MRSI was simulated for renal metabolism in MATLAB both with and without bipolar gradients. The resulting data were fit to a two-site exchange model [1], and the effects of vascular partial volume artifacts on kinetic modeling were assessed. Bipolar gradients were then incorporated into a gradient echo sequence to validate the simulations experimentally. The degree of diffusion weighting (b = 32 s/mm{sup 2}) was determined empirically from 1H imaging of murine renal vascular signal. The method was then tested in vivo using MRSI with bipolar gradients following injection of hyperpolarized [1-{sup 13}C]pyruvate (∼80 mM at 20% polarization). Results: In simulations, vascular signal contaminated the renal metabolic signal at resolutions as high as 2 × 2 mm{sup 2} due to partial volume effects. The apparent exchange rate from pyruvate to lactate (k{sub p}) was underestimated in the presence of these artifacts due to contaminating pyruvate signal. Incorporation of bipolar gradients suppressed vascular signal and improved the accuracy of kp estimation. Experimentally, the in vivo results supported the ability of bipolar gradients to suppress vascular signal. The in vivo exchange rate increased, as predicted in simulations, from k{sub p} = 0.012 s-{sup 1} to k{sub p} = 0.020-{sup 1} after vascular signal suppression. Conclusion: We have demonstrated the limited accuracy of the two-site exchange model in the presence of vascular partial volume artifacts. The addition of bipolar gradients suppressed vascular signal and improved model accuracy in simulations. Bipolar gradients largely affected kp estimation in vivo. Currently, slow-flowing spins in small vessels and capillaries are only partially suppressed, so further improvement is possible. Funding support: Seed Grant from the Radiological Society of North America, GE Healthcare, University of Wisconsin Graduate School.

60 APPLIED LIFE SCIENCES↗

The Volume of Gas Hydrate-Bound Gas in the Northern Gulf of Mexico

The northern Gulf of Mexico is known to host gas hydrate in submarine sediments. Estimating the amount of gas hydrate for both carbon cycle and resource interest has been ongoing for more than three decades. A large range of estimates (from 0.2 to 680 trillion cubic meters (TCM)) for hydrate-bound natural gas at standard temperature and pressure exists. We bring a new perspective to resource estimates by using ~800 publicly available petroleum industry well logs assessed for natural gas hydrate. Our resulting probabilistic range of the gas hydrate-bound natural gas in the northern Gulf of Mexico ranges from 37 TCM (10th percentile estimate) to 78 TCM (90th percentile estimate) of gas hydrate-bound natural gas, with a mean estimate of 56 TCM. This suggests that the gas hydrate resource density of the northern Gulf of Mexico is 4 times lower than that previously estimated by the Bureau of Ocean Energy Management. Finally, our results also indicate that over half of the gas hydrate in the area is strategically hosted in sand reservoirs, which is significantly higher than that previously estimated.

58 GEOSCIENCES↗

The use of strain tensor to estimate thoracic tumors deformation

Purpose: Respiration-induced kinematics of thoracic tumors suggests a simple analogy with elasticity, where a strain tensor is used to characterize the volume of interests. The application of the biomechanical framework allows for the objective determination of tumor characteristics. Methods: Four-dimensional computed tomography provides the snapshots of the patient's anatomy at the end of inspiration and expiration. Image registration was used to obtain the displacement vector fields and deformation fields, which allows one for the determination of the strain tensor. Its departure from the identity matrix gauges the departure of the medium from rigidity. The tensorial characteristic of each GTV voxel was determined and averaged. To this end, the standard Euclidean matrix norm as well as the Log-Euclidean norm were employed. Tensorial anisotropy was gauged with the fractional anisotropy measure which is based on the normalized variance of the tensors eigenvalues. Anisotropy was also evaluated with the geodesic distance in the Log-Euclidean framework of a given strain tensor to its closest isotropic counterpart. Results: The averaged strain tensor was determined for each of the 15 retrospectively analyzed thoracic GTVs. The amplitude of GTV motion varied from 0.64 to 4.21 with the average of 1.20 cm. The GTV size ranged from 5.16 to 149.99 cc with the average of 43.19 cc. The tensorial analysis shows that deformation is inconsiderable and that the tensorial anisotropy is small. The Log-Euclidean distance of averaged strain tensors from the identity matrix ranged from 0.06 to 0.31 with the average of 0.19. The Frobenius distance from the identity matrix is similar and ranged from 0.06 to 0.35 with the average of 0.21. Their fractional anisotropy ranged from 0.02 to 0.12 with the average of 0.07. Their geodesic anisotropy ranged from 0.03 to 0.16 with the average of 0.09. These values also indicate insignificant deformation. Conclusions: The tensorial framework allows for direct measurements of tissue deformation. It goes beyond the evaluation of deformation via comparison of shapes. It is an independent and objective determination of tissue properties. This methodology can be used to determine possible changes in lung properties due to radiation therapy and possible toxicities.

60 APPLIED LIFE SCIENCES↗

Sci—Thur AM: YIS - 11: Estimation of Bladder-Wall Cumulative Dose in Multi-Fraction Image-Based Gynaecological Brachytherapy Using Deformable Point Set Registration

Dosimetric parameters based on dose-volume histograms (DVH) of contoured structures are routinely used to evaluate dose delivered to target structures and organs at risk. However, the DVH provides no information on the spatial distribution of the dose in situations of repeated fractions with changes in organ shape or size. The aim of this research was to develop methods to more accurately determine geometrically localized, cumulative dose to the bladder wall in intracavitary brachytherapy for cervical cancer. The CT scans and treatment plans of 20 cervical cancer patients were used. Each patient was treated with five high-dose-rate (HDR) brachytherapy fractions of 600cGy prescribed dose. The bladder inner and outer surfaces were delineated using MIM Maestro software (MIM Software Inc.) and were imported into MATLAB (MathWorks) as 3-dimensional point clouds constituting the “bladder wall”. A point-set registration toolbox for MATLAB, Coherent Point Drift (CPD), was used to non-rigidly transform the bladder-wall points from four of the fractions to the coordinate system of the remaining (reference) fraction, which was chosen to be the emptiest bladder for each patient. The doses were accumulated on the reference fraction and new cumulative dosimetric parameters were calculated. The LENT-SOMA toxicity scores of these patients were studied against the cumulative dose parameters. Based on this study, there was no significant correlation between the toxicity scores and the determined cumulative dose parameters.

60 APPLIED LIFE SCIENCES↗